A regional-scale, high resolution dynamical malaria model that accounts for population density, climate and surface hydrology
© Tompkins and Ermert; licensee BioMed Central Ltd. 2013
Received: 27 August 2012
Accepted: 17 January 2013
Published: 18 February 2013
The relative roles of climate variability and population related effects in malaria transmission could be better understood if regional-scale dynamical malaria models could account for these factors.
A new dynamical community malaria model is introduced that accounts for the temperature and rainfall influences on the parasite and vector life cycles which are finely resolved in order to correctly represent the delay between the rains and the malaria season. The rainfall drives a simple but physically based representation of the surface hydrology. The model accounts for the population density in the calculation of daily biting rates.
Model simulations of entomological inoculation rate and circumsporozoite protein rate compare well to data from field studies from a wide range of locations in West Africa that encompass both seasonal endemic and epidemic fringe areas. A focus on Bobo-Dioulasso shows the ability of the model to represent the differences in transmission rates between rural and peri-urban areas in addition to the seasonality of malaria. Fine spatial resolution regional integrations for Eastern Africa reproduce the malaria atlas project (MAP) spatial distribution of the parasite ratio, and integrations for West and Eastern Africa show that the model grossly reproduces the reduction in parasite ratio as a function of population density observed in a large number of field surveys, although it underestimates malaria prevalence at high densities probably due to the neglect of population migration.
A new dynamical community malaria model is publicly available that accounts for climate and population density to simulate malaria transmission on a regional scale. The model structure facilitates future development to incorporate migration, immunity and interventions.
Models of malaria transmission are useful tools to aid understanding of the disease dynamic and have long been applied to assess the potential for intervention [2–4]. As weather parameters such as temperature and precipitation are important in determining the disease niche, there is also potential, given accurate forecasts/projections of these parameters, to use malaria models that account for climate in early warning systems in epidemic regions [5, 6] or to assess potential shifts in niche regions under climate change scenarios , although in both cases incorporation of other socio-economic factors is desirable [8–10].
As in other climate-related sectoral fields, predictive systems are often based on statistical models, which relate past case numbers to climate anomalies through regression techniques [5, 11]. The skill of statistical models depends on the data resources available for the region of interest. Such approaches can incorporate interventions, and have increased in complexity to include random errors and relevant socio-economic indicators  as in a recent model developed for dengue fever in Brazil . Poor, sparse or short data records confound the derivation of a reliable statistical model.
Dynamical models instead explicitly model the key equations describing the disease dynamic. A good understanding of the disease is required but if this is available this approach has the advantage of being able to incorporate day-to-day and even sub-daily diurnal variations in climate variables and thus differentiate transmission season risk due to sub-seasonal variability . Such models can potentially be applied to different regions to those in which they were developed, although local calibration may still be required.
Early dynamical disease models used idealized equations to describe the vector and host dynamic as separate discrete compartments [2, 15]. Many models have focused on the disease development within the host, with rate equations determining the passage between uninfected, infectious and recovering states and vectors given a parasite-free or infectious status. These models frequently include categories to account for immune status of the host, either to improve the model per se [16–18], or to improve our understanding of the process of immunity acquisition by determining which model best fits the available data [19–21]. Often the impact of climate was either neglected or incorporated through the force-of-infection related transmission rates.
Some recently developed dynamical models have focused on improving particular climatic aspects of the vector and parasite life-cycles. A dynamical model for the vector life-cycle has been introduced that also included a treatment of water bodies . The model was designed to be run on a local scale with knowledge of the water body in question, and included the temperature sensitivity of larvae while linking water temperature to air temperature through a function of cloud cover. Another recent study also introduced a simple generic treatment of the rainfall influence on breeding sites into a dynamical model applied to single locations . The Liverpool Malaria Model (LMM), which has been applied regionally to forecasting and climate projection issues [24, 25], includes the temperature effect on the sporogonic and gonotrophic cycles, explicitly representing the parasite and egg growth stages, as well as including a treatment of the temperature impact on vector mortality rates. The model was designed to be run over a regional scale but did not include an explicit representation of the surface hydrology; egg laying rates were proportional to the 10-day rain rate. Since flushing of larvae at intense rain rates could be important in reducing malaria [27, 28], this relation was modified to a highly nonlinear function to account for this [29, 30].
While successfully applied to modelling the parasite-vector cycles, one drawback of these models is that they do not explicitly represent the interaction between host and vector. The lower parasite ratio (PR) values in urban areas compared to rural areas is due to a range of socio-economic factors that impact access to health facilities and treatment. Field study and survey based statistical modelling imply a reduction of entomological inoculation rate (EIR) with increasing population density in peri-urban environments [31–33], indicating that surface hydrology and decreasing ratio of vector to host is key. Dynamical models that fail to incorporate host-vector dynamics may not differentiate between peri-urban and rural malaria. Moreover, inclusion of the host dynamic allows for an improved representation of human behaviour in malaria propagation, including interventions, treatments and also economic migration, which can reintroduce the disease into malaria-free zones [34–36].
A dynamical model has recently focused on the interaction between the host and vector, which used climate data from a single location in central Tanzania in a set of idealized simulations to investigate the impact of interventions such as bed net usage, although no comparison to data was conducted for the location . Another malaria model [38, 39] used an agent-based approach to represent vectors and included population density. The study mapped out a rural village at a 10 m spatial scale, modelling growth of breeding sites using a two dimensional hydrological model and followed the evolution of each individual vector within the village domain. Such a modelling approach would be impossible to apply on a regional scale for planning purposes, due to the shear numerical cost and the lack of the detailed input required on a wide scale. While [38, 39] claimed that explicitly modelling each small scale breeding site of the anopheles vector and its interaction with its host was vital to present malaria transmission, we note that in the atmospheric sciences reliable weather predictions are made without explicitly modelling effects occurring below the truncation scale of the model (such as clouds and convection), but instead representing their bulk statistical effects collectively via a so-called parametrization scheme[40, 41]. In the same way, a regional-scale malaria model can incorporate important small-scale effects such as the growth of available breeding sites (ponds and puddles) in the rainy season, representing the net bulk impact in a parametrization scheme.
The aim here is to introduce a new dynamical malaria model that can be used on a relative fine spatial resolution of order 10 km but applied over a regional scale and explicitly incorporates vector-host interaction in addition to a treatment of the aggregate effect of seasonal temporary bodies in a pond model. While idealized, the hydrological component puts a framework into place that will allow future development such as incorporating land-surface related permanent breeding sites and coupling to a more complex surface hydrology. Moreover, the inclusion of the population density will facilitate the future study of such phenomena as host migration, interventions and urbanization. Thus the model presented here provides a tool that takes a first step towards fulfilling two of the recommendations of : to understand the relationships between climate and climate-sensitive diseases and health issues under different environmental conditions by incorporating other data into... health forecast services, for example population, rural vs. urban residence, [and] migration. After introducing the model framework, an evaluation the model is conducted in an idealized framework and also using point and gridded regional integrations comparing the model to field studies of entomological parameters, malaria atlas project (MAP) processed survey data and an existing dynamical model. This evaluation demonstrates the ability of model to predict the seasonality and spatial distribution of malaria transmission, before the future directions are summarized.
Overview of the VECTRI model
The malaria model presented is a grid cell distributed dynamical model and is referred to as VECTRI; the vec tor-borne disease community model of the International Centre for Theoretical Physics, Tri este. This manuscript documents VECTRI release version v1.25 as of June 2012. In as far as possible, the model physics and associated parameters are taken from the literature for the Anopheles gambiae complex and the Plasmodium falciparum malaria parasite. In the present version each location (grid cell) is independent, but the structure of the model will allow communication between grid cells such as vector flight or human population migration to be easily incorporated. The following sections describe the basic structure of the model, with emphasis placed on its novel aspects.
The vector status is also bin resolved, consisting of two properties: the gonotrophic and sporogonic cycles. It is thus represented as a two dimensional array V(N gono ,N sporo ). All vectors in the first gonotrophic bin are in meal-searching mode, and once a meal is obtained, the vectors advance in terms of the egg development state at a rate RV,gono related to ambient temperature until the final bin is reached (using the advection equation similar to eqn. 1). At this point the vector lays a new raft of eggs and is recycled to the first meal-searching bin.
Each timestep, parasite transmission may occur to a proportion of the biting vectors, and the status of these vectors will subsequently additionally progress in the sporogonic dimension, with the rate RV,sporo again determined by temperature. Once vectors reach the final bin they are infective to humans and remain so until death. A third array of bins maps the status of the disease in the human host (H) population (dimension N H ), with the first bin representing the uninfected population. The model does not include age or immunity factors, which is the subject of present model development.
Thus, while the model introduces new relationships regarding the surface hydrology and the explicit interaction between vector and host, the underlying numerical structure is similar to that employed in the LMM . A timestep of one day is used to integrate the model equations, although a shorter timestep could be used if input data (temperature/rainfall) are available on these timescales, and the advection equation is solved using a simple upstream numerical scheme.
The model can be flexibly integrated using a wide range of horizontal resolutions. Since the model presently does not permit vectors to move to neighbouring grid cells, the model resolution is limited to an upper (finest) level on the order of 1-5 km; an indicated range below which mosquito movement can become significant [43–45]. Integrating the model at O(10 km) resolution is desirable, even if corresponding observational data of clinically proven malaria cases is available only on coarser scale health districts. This is because malaria transmission is a highly nonlinear function of its drivers such as climate and land surface and thus it is preferable to use highest possible spatial resolution inputs of population density, land surface and weather to account for this, while the model output is subsequently aggregated to the spatial scale at which health data is available for comparison.
Structure of the VECTRI model
VECTRI default constants
E + I
Total evaporation and infiltration losses
based on 
Larvae growth degree days
Larvae growth degree days
Gonotrophic cycle degree days
Sporogonic cycle degree days
Larvae flushing factor for infinite rain rate
based on 
Pond growth rate factor
set by tuning or local knowledge of basin
K mar 1,0
Constant of Martens I vector survival scheme
K mar 1,1
Constant of Martens I vector survival scheme
K mar 1,2
Constant of Martens I vector survival scheme
K mar 2,0
Constant of Martens II vector survival scheme
K mar 2,1
Constant of Martens II vector survival scheme
K mar 2,2
Constant of Martens II vector survival scheme
Carrying capacity of water bodies
Number of eggs per batch that result in female vectors
P L,surv 0
Larvae base daily survival rate
Probability of transmission from infective host to vector during single bloodmeal
Probability of transmission from infective vector to host during single bloodmeal
Pond water offset from air temperature
Minimum T wat for larvae development
based on 
Maximum T wat for larvae development
based on 
Minimum T2m for egg development
Minimum T2m for sporogonic cycle
Larvae-flushing rainfall e-folding factor
based on 
Population density zoophilic factor
set by tuning
Maximum temporary pond fraction in cell
set by tuning or local knowledge of basin
There is considerable uncertainty in the setting of the rate coefficient K L , however, with K L value of 90.9 degree days, while a linear approximation of the relationship derived by  results in a much slower rate of 200 degree days. A further source of uncertainty is the specification of the water temperature itself, which depends on the shading of the pool and its dimensions in addition to the ambient air temperature, and is described in the hydrology component below. VECTRI also permits the user to avoid this uncertainty by following  in setting a fixed larvae growth rate (i.e. independent of T wat ) to have a cycle of 12 days. The sensitivity to this relationship is investigated later. Irrespective of the scheme used, an upper temperature limit TL,max is specified above which larvae death occurs.
Egg hatching into larvae and the pupae development stage are both typically on the order of one day [48, 49] and thus are poorly resolved by the daily timestep employed by VECTRI and other similar dynamical models. In order to avoid truncation problems the length is fixed in VECTRI to last exactly one day and is temperature independent.
In the first term on the right, M L is the total larvae biomass per unit surface area of a water body, and w is the fraction coverage of a grid cell by potential breeding sites (not open water) and is given by the surface hydrology component described below. If w = 0, the survival rate PL,surv is also zero. The maximum carrying capacity, ML,max, is set to 300 mg m-2 and larvae mass is assumed to increase linearly with a stage 4 larva having a mean mass of 0.45 mg, both following [22, 38] closely. All larvae die above a water temperature of TL,max.
All female vectors are assumed to find a blood meal in the first night of searching, although this fraction can be set as a model parameter. Insecticide treated nets are, for example, able to frustrate host-seeking mosquitoes  and reduce the mosquito population number. It would depend on the human population and availability of animals. In this way both newly emerging and existing adult mosquitoes might be prevented from progressing in the feeding cycle.
At the end of the cycle the female vector lays N egg eggs that will eventually hatch into females; as is usual in such models, the eggs laid that result in the males are neglected. She subsequently cycles to the meal searching box. The number of eggs is highly variable and depends on vector species. The choice of N egg = 120, corresponding to a batch size of 240 assuming equality between the sexes, follows  but [56–58] indicate that this could be an overestimation. Present model development underway will permit stochastic variation of parameters such as N egg in the future to account for their considerable degree of uncertainty in an ensemble modelling framework.
It should be noted that this assumes bites are randomly taken and not influenced by the host’s infective state. Heterogeneous biting can impact the basic reproductive number considerably [59, 60]. Heterogeneity of feeding habits is related to a wide array of factors, including host attractiveness to the vector and their vicinity to breeding sites, and heterogeneity in interventions such as the use of bed nets, in particular the increased use by hosts suffering from clinical malaria. Some of these effects could easily be included in VECTRI if relevant data were available.
After a number of days the sporozoites invade the salivary glands of the mosquito which subsequently becomes infective to humans. The mosquito is assumed to remain in this infective state until its death.
In addition to the larvae, gonotrophic and sporogonic cycles, temperature also plays a role in determining the mortality of vector. High air temperatures increase vector mortality, but the relationship is uncertain, especially at the high and low temperature bounds of transmission. As in the larvae cycle, two schemes are incorporated in the VECTRI model to permit a multi-model approach, with constants given in Table 1:
The results in this paper are obtained using scheme 2, referred to in  as Martens II.
The factor represents the level of vector zoophily. While members of the Anopheles gambiae complex are in general considered anthropophilic to varying degrees , with arabiensis more zoophilic than sensu stricto, vectors take an increasing proportion of blood meals from cattle in lower population density rural areas with high livestock numbers , although the effectiveness of zooprophylaxis is still debated . The exponential factor reflects this, with the e-folding population density for the effect set to τ zoo = 50 km-2. Thus the factor only has a significant impact for rural populations below this number and avoids the model producing excessively high biting rates and EIR for sparsely populated locations. In future, VECTRI will allow vector movement between cells allowing anthropophilic vectors to cluster around population centres.
where is the Poisson distribution for mean . If bed nets are in use, eqn. 11 could be modified to incorporate this, increasing the mean bite rate for a subset of the unprotected population. This involves a number of complications however, since accurate data would be required concerning bed net distribution and use, how this usage correlates to host infective state, and which proportion of bite are taken during the hours of sleeping.
The host population is represented by the vector H(N host ), and each VECTRI timestep a proportion Pv → h of hosts become infected and progress through the array until 20 days later they assume an infective status, an average value for immune and non-immune subjects [76–79]. Non-immune hosts clear infections at an e-folding rate of C ni = 150 days. Even after a century of study of the disease, the paradigm of naturally acquired immunity (NAI) is still hotly debated . Therefore the present version of the model neglects host immunity, and the impact of the various representations of immunity in VECTRI will be the subject of a companion article.
Considering first w perm , converting land use and terrain information into a fractional coverage of breeding areas provided by permanent water bodies is a significant challenge. Wave/ripples action that can drown larvae  and the presence of predators in larger bodies  imply that larvae exist only in a sub-fraction of such water bodies, in pooling that occurs on the edges of lakes and rivers or the shallow edges of ponds. Higher soil moisture in the vicinity of large water bodies can boost available breeding sites by reducing infiltration loss and increasing the lifetime of temporary puddles and ponds. Lastly, rivers and streams can even confound classic relationships between rainfall and vector density by actually providing more breeding sites during drought periods when flow slows or stops altogether . This remains the subject of current research and the default VECTRI model therefore sets w perm to zero in each grid cell, implying that the breeding site availability is dominated by seasonal ephemeral ponds and malaria incidence will be potentially underestimated in the vicinity of larger semi-permanent water bodies. A user operating the model for a local area with knowledge of permanent water bodies can set this value appropriately.
The net aggregated fractional coverage by temporary pools w pond is derived from a simple water balance model. Ponds are replenished by surface runoff Q, while infiltration I (seepage) and evaporation E and pond overflow reduce their water content. An important parameter is the maximum coverage of the temporary ponds w max , which described the extent of depressions that could potentially become water filled at the peak of a wet season, which are often small in scale, with total catchments of the gully systems studies in Niger in the HAPEX-SAHEL experiment measuring 0.2 km-2. Presently it is assumed that the runoff Q that fills the ponds mostly falls within these depressions and thus the runoff is set to Q = w max P, where P is the precipitation rate. Thus sub-surface infiltration occurring within the depressions is also considered a water source for temporary pools . Future developments will introduce an improved runoff treatment accounting for soil texture and slope.
The source from precipitation is balanced by evaporative, infiltration and overflow losses. In high-resolution simulations of surface hydrology in Niger,  found overflow losses to be approximately 20% of total losses, more than three times the losses due to evaporation. It should also be noted that overflow losses in field campaigns are difficult to measure and thus are often incorporated in the infiltration, which is calculated as a residual in the water balance calculation. Losses through pond overflow are assumed to increase linearly with pond fraction in VECTRI, achieved by scaling the runoff by a factor . Once the pond fraction reaches its maximum, all surface runoff overflows and is lost. Infiltration losses vary substantially depending on soil texture and life-scale of the pond in question. Often the infiltration is a highly nonlinear function of water body extent, since silting may significantly reduce infiltration in the lowest part of longer-lived or semi-permanent pools . This results in a fast initial decay after rain events due to high infiltration rates at the pool edges, followed by a slower decay, while temporary, shorter lived water bodies tend to have more uniform infiltration rates. These rates can be very high, exceeding 600 mm day-1. Presently the VECTRI model simply sets a fixed constant infiltration rate per unit pond area.
Finally, it is recalled that the pond dimension limits larvae mortality rates through the availability of breeding sites governed in eqn. 3. This is an approximation of the net affect of crowding which leads to higher mortality rates, longer development times and smaller adults , which in turn have a competitive disadvantage . The biomass is considered to be distributed equally through all available breeding sites and variability between breeding sites in neglected, supported by  who noted that females avoid ponds that are overcrowded with existing larvae.
In addition to pond dimension the other important parameter of water bodies is the temperature of the water near the surface. The Depinay model  developed a complex empirical function for water temperature as a function of ambient relative humidity and water body size. As the VECTRI model is applied regionally, specific information about individual water body size may not be included. The temperature in shallow ponds and puddles is homogeneous to a good approximation and is often one or two degrees warmer than the air temperature [92, 93]. VECTRI therefore assumes that the temperature of pools T wat to have a fixed offset relative to the air temperature. The default value adopted is a positive offset of 2 K, however, in hot locations it is likely that vector will preferentially choose shaded breeding locations and a lower or even negative offset may be more appropriate. If accurate gridded weather information for wind and surface radiation were available, this aspect of the model could be potential improved implementing a single energy balance model along the lines of [94, 95]. While larger permanent water bodies such as lakes and rivers can have complex stratification of the vertical temperature profile, as discussed above, larvae development occurs mostly in the shallow waters and pools that form on the lake/river boundaries and thus the temperature relation for the permanent water fraction is treated in the same way as the temporary ponds.
Evaluation of VECTRI
Evaluation methodology and data
VECTRI is evaluated using a series of three testbed arrangements. The first test determines the sensitivity of the force of infection measured by EIR in relation to climate in an idealized setting of so-called ’equilibrium runs’. The model is driven to equilibrium in an array of experiments using constant values of temperature and rainfall and compared to the reference model described below.
The VECTRI model is next confronted with a wide range of entomological data consisting of EIR and circumsporozoite protein rate (CSPR) measurements collected from field campaigns located throughout West and Central Africa and which were collated in a review by . These campaigns are divided between rural and urban environments. In particular, the availability of a large number of both peri-urban and rural field studies in and nearby Bobo-Dioulasso documented in [97–106] lead to this location being used as a focus site. For each location, the models are driven by daily weather data from the nearest available station, which for points in or nearby population centres is often within 10 km, but can be much further for isolated rural locations. A simple interpolation methodology is used to fill temporal data gaps . While the malaria field campaigns are often only available for a single year, and never more than a few seasons, the malaria models are integrated for over four decades to give an indication of interannual variability, with results from these ’station integrations’ shown from 1973 to 2006.
Lastly, the malaria prevalence is evaluated from ’regional integrations’ conducted with VECTRI, the first in Eastern Africa (domain 27°E-42°E and 12°S-6°N) and the second for West Africa (18°W-9°E and 4.5°N to 19.5°N). Integrations last for 10 years and use a horizontal spatial resolution of 0.1 degrees, equivalent to approximately 11 km. As freely available daily station data is sparse, the model is instead driven with satellite-derived rainfall using the NOAA FEWS RFE product , available from late 2000 to the present on a 11 km resolution which is interpolated to the VECTRI grid using conservative remapping. The gridded 2 metre height temperature data is taken from the ECMWF ERA-Interim reanalysis product . As the resolution of the reanalysis is relatively coarse at approximately 75 km in the tropics, it is statistically downscaled using a constant adiabatic lapse rate of 6.5 K km-1 to correct for the height error between the model topography and the ETOPOv2 topography  bilinearly interpolated to the VECTRI grid. Starting from artificial initial conditions of 5% of the population carrying the parasite, the model is integrated for a spin up period of 3 years cycling the first year of weather data to allow the malaria loading in the population to equilibrate before the integration runs freely. Focus is on the integration for Eastern Africa due to the wide range of topography and the consequential inclusion of malaria free, epidemic and endemic zones, which are compared to the MAP malaria analysis [111–113]. The MAP analysis ingests PR survey data occurring up to 2010 into a statistical model that also incorporates environmental data to provided gridded values of PR. The reliability of the MAP data depends on the local availability of survey data. In Africa, for most areas there is a large uncertainty of the MAP prediction. For example, the prediction is less certain for Uganda and Tanzania, which represent under-surveyed countries. More certain is the MAP analysis for most parts of Kenya, although  marks the coastal and western regions with high uncertainty despite the far greater densities of survey data available there. The two regional runs are thus also compared to direct survey data from 199 field campaigns using the analysis of  in order to assess the impact of population density on EIR.
The VECTRI model is contrasted to the Liverpool Malaria Model (LMM), which has been applied in various forms to regional malaria modelling across Africa for both seasonal timescale forecasting  and to investigate climate change impacts on transmission . The original version of LMM described by  was updated in terms of its surface hydrology by [29, 30], and referred to here as LMM2010. The LMM2010 is a dynamical mathematical-biological malaria model, which is driven by the two variables of daily mean temperature and precipitation. The treatment of surface hydrology in LMM2010 differs significantly to VECTRI. It applies a fuzzy distribution model to relate egg-larvae development to the decadal (10 day accumulated) rainfall amount, with optimal conditions assumed to occur with 1 mm day-1 average rain rate. rain rate. Egg-laying rates are reduced below 1 mm day-1 rain rates to reflect the lack of breeding sites, and also above this threshold since flushing effects are then assumed to dominate. The use of a ten day average rainfall essentially equates to the assumption that rainfall feeds ponds and pools which have a ten day mean lifespan. Water temperatures are not considered by the LMM2010 in the egg-larval development rates.
The LMM2010 uses similar entomological and parasitological numerical components to VECTRI. By contrast, the LMM2010 neglects the population density and is only calibrated to malaria observations from rural locations in West Africa. This means that the LMM2010 is not able to simulate urban malaria conditions. The model is also not valid for irrigated malaria areas, where malaria transmission is not depending solely on rainfall. The vector-host contact is furthermore static in the LMM2010 since a fixed proportion of the mosquito bites is taken up by humans and animals.
Results and discussion
In terms of rainfall sensitivity, the two models differ considerably. The nonlinear function that determines number of egg laid as a function of rainfall in the LMM2010 model is very sensitive to rainfall. As the rain rate deviates from the assessed ideal value of 1 mm/day, the EIR decreases quickly. Thus in the rainy season, when daily rainfall amounts usually exceed 1 mm/day, it is likely that the LMM will produce an anticorrelation between rainfall amount and malaria transmission intensity. Instead, although the VECTRI model includes a flushing parametrization, the effect is much weaker at these moderate rain rates of 0 to 20 mm day-1 and rainfall is therefore positively correlated with transmission transmission intensity. Moreover, above 2 to 3 mm day-1, the sensitivity of EIR to rain rate is much lower in VECTRI than LMM2010, which is likely to result in a lower interannual variability of malaria in endemic regions in the former model.
Intercomparing EIR from the two malaria models, the very high interannual variability of the LMM2010 is obvious, related to its elevated sensitivity to rainfall demonstrated in the earlier constant-input experiments. The LMM2010EIR value ranges by more than an order of magnitude between the lowest and highest year, while the interannual variability is less than a factor of two for the VECTRI model. It is interesting to note that, despite the basic underlying structure being very similar in the two models, the parametrization choices, in particular the implementation of the surface hydrology in the VECTRI model, results in almost a zero correlation between the two models in their representation of interannual variability; indeed, the EIR appear anti-correlated indicating that rainfall variability is determining the interannual variability to a large extent.
The black squares in the figure give EIR values for high population areas in peri-urban Bobo-Dioulasso, which are much smaller than in the rural environments but nevertheless non-zero. The LMM2010 is not designed to simulate these urban cells, however the VECTRI model is seen to reproduce reasonably well the contrast between rural and urban areas. Even though the treatment of the surface hydrology is identical in urban and rural environments - a gross oversimplification - the VECTRI model is able to mimic the drop in EIR which derives merely from the lower ratio of vector to host in urban areas. In contrast with the other West African locations, urban transmission is overestimated in this location. In the second panel, it is seen that while both models are again similar and perform well in reproducing the observed CSPR rates, the VECTRI model does less well in reproducing the distinction between rural and urban environments. Although the CSPR is lower in urban environments, it is still far larger than the observations in VECTRI. An possible implication is that vector lifetimes for urban areas are too long in the model.
The VECTRI model also reproduces the malaria zones in the warmer and more humid coastal regions in Kenya and Tanzania, but with values of PR along the Kenyan coast and northern Tanzania in VECTRI (30-70%) appear to greatly exceed those in the MAP analysis (0-30%). This is a region in which the survey data incorporated into the MAP analysis is dense, although Figure Three of  still identifies this coastal region as relatively uncertain. A large part of the discrepancy between VECTRI and MAP is likely due to the increasing interventions including widespread distribution of insecticide-treated nets (ITNs) that have occurred over the past decade that have greatly decreased parasite ratios and hospital admissions [117, 118]. For example, the survey of 30 villages in Malindi, Kilifi, and Kwale Districts carried out in the late 1990s before ITN distribution started  (included in the MAP analysis) reported a parasite prevalence ranging from 38 to 83%, with a mean slightly exceeding 60% in each district, in close agreement with the VECTRI model. This highlights the importance of the future incorporations of interventions into VECTRI if it is to be applied to the seasonal forecasting task.
It is emphasized that these are single deterministic integrations of the VECTRI model using its default parameter settings, and an important next step will be to document how sensitive this relationship is to the parameter settings in a stochastic integration framework . Nevertheless, Figure 13 includes the sensitivity to the parameter τ zoo , chosen as it is a tuned parameter in the model, with integrations using values of 35, 70 and 100 km-2 also applied. As expected, these altered values impact the biting rate at low population densities, with lower values τ zoo significantly increasing EIR in sparsely populated areas. Above, 200 people per km-2 the parameter setting has no impact, and for all settings the model over-sensitivity relative to  remains. This is because the parameter is designed to account for increasing zoophilicity at very low population densities, and thus are unable to represent the impact of human population movements which impact PR at low and high population densities. Nevertheless, while the model presently neglects differences in surface hydrology in peri-urban versus rural areas , urban micro meteorology and housing differences , the result indicates that a zero order impact of population density is the dilution effect, where increasing population densities reduce the probability of parasite transmission to the vector [31, 122].
A new open source community model (contact the corresponding author) for malaria is presented that specifically represents dynamical equations for the relevant temperature sensitive cycles in the transmission of malaria, namely the larvae growth cycle in addition to the sporogonic and gonotrophic cycles. Each cycle is divided into an array of bins to resolve the respective processes in order to model the delay in the malaria season with respect to the driving cycle of rainfall, while temperature affects development rates. The VECTRI model implements a framework to consider the surface hydrology of temporary water bodies and ponds. It also incorporates the human/vector element by accounting for population density in the calculation of transmission probabilities, enabling the VECTRI model to differentiate between urban, peri-urban and rural environments and permitting the future incorporation of various interventions into the modelling framework in addition to a treatment of urban-rural and transnational migration.
The model has been integrated in an idealized equilibrium environment to compare the sensitivity of weather parameters to an existing dynamical model LMM2010, and the contrasting response to rainfall due to the incorporation of surface hydrology in VECTRI was noted. Two further sets of experiments confronted the model with a range of observational data for infectious biting rates and parasite ratios in the human population taken from several hundred field studies mostly based in West or Eastern Africa. The model reproduces mean rural infectious bite rates for a range of locations throughout West Africa quite well, with a seasonal cycle of PR that appeared reasonable.
The incorporation of the population density in the biting rate calculation enables VECTRI to differentiate peri-urban and rural environments and it was demonstrated that the model was also able to reproduce the much lower EIR and CSPR rates measured in a peri-urban suburb environment of Bobo-Dioulasso less than 20 km distant from a rural site with far higher transmission. This is interesting as it indicates that simply representing the lower ratio of vector to human numbers is sufficient to represent the lower transmission intensities generally found in peri-urban environments without invoking differences in surface hydrology, access to health facilities, housing quality or use of ITNs. Regional integrations for Eastern Africa highlighted this ability of VECTRI to grossly reproduce EIR rates as a function of population density, in addition to its reasonable distribution of endemic and malaria free regions, separated by intermediate altitude, epidemic regions, although relative to the MAP data, the VECTRI model tends to overestimate PR in endemic zones. Whether this is due to missing processes, such as immunity or migration, or poor parameter choices remains unclear and the subject of investigation.
While already representing a useful tool to aid research into the understanding of malaria transmission, the VECTRI is undergoing further development to incorporate further aspects important for malaria transmission. These developments include: a flexible incorporation of various models for host immunity, a stochastic ensemble framework to account for uncertainty in model parameters and observations, improved surface hydrology that accounts for terrain slope, soil type and land use characteristics, a dynamical host model that includes migration, population demographics and urbanisation, explicit incorporation of commonly employed interventions, and lastly multiple vector types and vector dispersion.
This work contributes towards two EU seventh framework projects: Quantifying Weather and Climate impacts on health in developing countries (QWeCI, grant 243964) and HEALTHY FUTURES (grant 266327). AMT received support from both projects, while VE was fully funded by QWeCI. Dr. Catherine Linard is thanked for providing an Africa-wide version of the Afripop dataset for use in the model, while the MAP team at Oxford University kindly provided original files of the MAP malaria analysis. Comments on the manuscript by Prof. David Taylor were appreciated, in addition to discussion with Profs. Mercedes Pascual and Andy Dobson at the ICTP workshop on infectious diseases in Arusha, Tanzania, 2012. The reviewers also gave suggestions that greatly improved the manuscript.
- McKenzie F: Why model malaria?. Parasitol Today. 2000, 16: 511-516. 10.1016/S0169-4758(00)01789-0.PubMedGoogle Scholar
- MacDonald G, Cuellar CB, Foll CV: The dynamics of malaria. Bull World Health Organ. 1968, 38: 743-755.PubMed CentralPubMedGoogle Scholar
- Dietz K, Molineaux L, Thomas A: A malaria model tested in the African savannah. Bull World Health Organ. 1974, 50: 347-357.PubMed CentralPubMedGoogle Scholar
- Koella JC: On the use of mathematical models of malaria transmission. Acta Trop. 1991, 49: 1-25. 10.1016/0001-706X(91)90026-G.PubMedGoogle Scholar
- Thomson MC, Doblas-Reyes FJ, Mason SJ, Hagedorn R, Connor SJ, Phindela T, Morse AP, Palmer TN: Malaria early warnings based on seasonal climate forecasts from multi-model ensembles. Nature. 2006, 439: 576-579. 10.1038/nature04503.PubMedGoogle Scholar
- Ceccato P, Ghebremeskel T, Jaiteh M, Graves P, Levy M, Ghebreselassie S, Ogbamariam A, Barnston A, Bell M, Del Corral J, Connor SJ, Fesseha I, Brantly EP, Thomson MC: Malaria stratification, climate, and epidemic early warning in Eritrea. Am J Trop Med Hyg. 2007, 77: 61-68.PubMedGoogle Scholar
- Lafferty KD: The ecology of climate change and infectious diseases. Ecol. 2009, 90: 888-900. 10.1890/08-0079.1.Google Scholar
- Ghebreyesus TA, Haile M, Witten KH, Getachew A, Yohannes M, Lindsay SW, Byass P: Household risk factors for malaria among children in the Ethiopian highlands. Trans R Soc Trop Med Hyg. 2000, 94: 17-21. 10.1016/S0035-9203(00)90424-3.PubMedGoogle Scholar
- Ijumba JN, Mosha FW, Lindsay SW: Malaria transmission risk variations derived from different agricultural practices in an irrigated area of northern Tanzania. Med Vet Entomol. 2002, 16: 28-38. 10.1046/j.0269-283x.2002.00337.x.PubMedGoogle Scholar
- Gething PW, Smith DL, Patil AP, Tatem AJ, Snow RW, Hay SI: Climate change and the global malaria recession. Nature. 2010, 465: 342-345. 10.1038/nature09098.PubMed CentralPubMedGoogle Scholar
- Githeko AK, Ndegwa W: Predicting malaria epidemics in the Kenyan highlands using climate data: a tool for decision makers. Glob Chan Hum Health. 2001, 2: 54-63. 10.1023/A:1011943131643.Google Scholar
- Ngom R, Siegmund A: Urban malaria in Africa: an environmental and socio-economic modelling approach for Yaoundé, Cameroon. Nat Hazards. 2010, 55: 599-619. 10.1007/s11069-009-9485-x.Google Scholar
- Lowe R, Bailey TC, Stephenson DB, Graham RJ, Coelho CAS, Sá Carvalho M, Barcellos C: Spatio-temporal modelling of climate-sensitive disease risk: Towards an early warning system for dengue in Brazil. Comput Geosci. 2011, 37: 371-381. 10.1016/j.cageo.2010.01.008.Google Scholar
- Paaijmans KP, Blanford S, Bell AS, Blanford JI, Read AF, Thomas MB: Influence of climate on malaria transmission depends on daily temperature variation. Proc Nat Acad Sci. 2010, 107: 15135-15139. 10.1073/pnas.1006422107.PubMed CentralPubMedGoogle Scholar
- Hethcote HW: Qualitative analyses of communicable disease models. Math Biosci. 1976, 28: 335-356. 10.1016/0025-5564(76)90132-2.Google Scholar
- Ngwa GA, Shu WS: A mathematical model for endemic malaria with variable human and mosquito populations. Math Comput Model. 2000, 32: 747-763. 10.1016/S0895-7177(00)00169-2.Google Scholar
- Pascual M, Cazelles B, Bouma MJ, Chaves LF, Koelle K: Shifting patterns: malaria dynamics and rainfall variability in an African highland. Proc R Soc Lond B Biol Sci. 2008, 275: 123-132. 10.1098/rspb.2007.1068.Google Scholar
- Chitnis N, Hyman JM, Cushing JM: Determining important parameters in the spread of malaria through the sensitivity analysis of a mathematical model. Bull Math Biol. 2008, 70: 1272-1296. 10.1007/s11538-008-9299-0.PubMedGoogle Scholar
- Aron JL: Mathematical modelling of immunity to malaria. Math Biosci. 1988, 90: 385-396. 10.1016/0025-5564(88)90076-4.Google Scholar
- Filipe JAN, Riley EM, Drakeley CJ, Sutherland CJ, Ghani AC: Determination of the processes driving the acquisition of immunity to malaria using a mathematical transmission model. PLoS Comput Biol. 2007, 3: e255-10.1371/journal.pcbi.0030255.PubMed CentralPubMedGoogle Scholar
- Doolan DL, Dobano C, Baird JK: Acquired immunity to malaria. Clin Microbiol Rev. January 2009, 22 (1): 13-36. 10.1128/CMR.00025-08.Google Scholar
- Depinay JMO, Mbogo CM, Killeen G, Knols B, Beier J, Carlson J, Dushoff J, Billingsley P, Mwambi H, Githure J, Toure AM, McKenzie FE: A simulation model of African Anopheles ecology and population dynamics for the analysis of malaria transmission. Malar J. 2004, 3: 10.1186/1475-2875-3-29.Google Scholar
- Alonso D, Bouma M, Pascual M: Epidemic malaria and warmer temperatures in recent decades in an East African highland. Proc R Soc Lond B Biol Sci. 2011, 278: 1661-1669. 10.1098/rspb.2010.2020.Google Scholar
- Jones AE, Morse AP: Application and validation of a seasonal ensemble prediction system using a dynamic malaria model. J Clim. 2010, 23: 4202-4215. 10.1175/2010JCLI3208.1.Google Scholar
- Jones AE, Morse AP: Skill of ENSEMBLES seasonal re-forecasts for malaria prediction in West Africa. Geophys Res Lett. 2012, 39: L23707-10.1029/2012GL054040.Google Scholar
- Hoshen MB, Morse AP: A weather-driven model of malaria transmission. Malar J. 2004, 3: 32-10.1029/2012GL054040.PubMed CentralPubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- Dieng H, Rahman GMS, Abu Hassan A, Che Salmah MR, Satho T, Miake F, Boots M, Sazaly AB: The effects of simulated rainfall on immature population dynamics of Aedes albopictus and female oviposition. Int J Biometeorol. 2012, 56: 113-120. 10.1007/s00484-011-0402-0.PubMedGoogle Scholar
- 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. Malar J. 2011, 10: 10.1186/1475-2875-10-35.Google Scholar
- 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. Malar J. 2011, 10: 62-10.1186/1475-2875-10-62.PubMed CentralPubMedGoogle Scholar
- Robert V, Macintyre K, Keating J, McWilson W, Trappe JP, Duchemin JB, Beier JC: Malaria transmission in urban sub-Saharan Africa. Am J Trop Med Hyg. 2003, 68: 169-176.PubMedGoogle Scholar
- Hay SI, Guerra CA, Tatem AJ, Atkinson PM, Snow RW: Urbanization, malaria transmission and disease burden in Africa. Nat Rev Microbiol. 2005, 3: 81-90. 10.1038/nrmicro1069.PubMed CentralPubMedGoogle Scholar
- Kelly-Hope L, McKenzie FE: The multiplicity of malaria transmission: a review of entomological inoculation rate measurements and methods across sub-Saharan Africa. Malar J. 2009, 8: 19-10.1186/1475-2875-8-19.PubMed CentralPubMedGoogle Scholar
- Martens P, Hall L: Malaria on the move: human population movement and malaria transmission. Emerg Infect Dis. 2000, 6: 103-109. 10.3201/eid0602.000202.PubMed CentralPubMedGoogle Scholar
- Stoddard ST, Morrison AC, Vazquez-Prokopec GM, Soldan VP, Kochel T, Kitron U, Elder JP, Scott TW: The role of human movement in the transmission of vector-borne pathogens. PLoS Negl Trop Dis. 2009, 3: e481-10.1371/journal.pntd.0000481.PubMed CentralPubMedGoogle Scholar
- Tatem AJ, Smith DL: International population movements and regional Plasmodium falciparum malaria elimination strategies. Proc Nat Acad Sci. 2010, 107: 12222-12227. 10.1073/pnas.1002971107.PubMed CentralPubMedGoogle Scholar
- Eckhoff PA: A malaria transmission-directed model of mosquito life cycle and ecology. Malar J. 2011, 10: 10-10.1186/1475-2875-10-303.Google Scholar
- Bomblies A, Duchemin JB, Eltahir EAB: Hydrology of malaria: Model development and application to a Sahelian village. Water Resour Res. 2008, 44: W12445-10.1029/2008WR006917.Google Scholar
- Bomblies A, Duchemin JB, Eltahir EAB: A mechanistic approach for accurate simulation of village scale malaria transmission. Malar J. 2009, 8: 223-10.1186/1475-2875-8-223.PubMed CentralPubMedGoogle Scholar
- Tiedtke M: A comprehensive mass flux scheme for cumulus parameterization in large-scale models. Mon Wea Rev. 1989, 117: 1779-1800. 10.1175/1520-0493(1989)117<1779:ACMFSF>2.0.CO;2.Google Scholar
- Tompkins AM: A prognostic parameterization for the subgrid-scale variability of water vapor and clouds in large-scale models and its use to diagnose cloud cover. J Atmos Sci. 2002, 59: 1917-1942. 10.1175/1520-0469(2002)059<1917:APPFTS>2.0.CO;2.Google Scholar
- Omumbo J, Platzer B, Girma A, Connor S J: Climate and health in, Africa: 10 years on - Workshop report. Tech. rep., International Research Institute for climate and society, IRI, 2011 available at [http://iri.columbia.edu/publications/download=1090]
- Gillies MT: Studies on the dispersion and survival of Anopheles gambiae Giles in East Africa, by means of marking and release experiments. Bull Entomol Res. 1961, 52: 99-127. 10.1017/S0007485300055309.Google Scholar
- Rowley WA, Graham CL: The effect of age on the flight performance of female Aedes aegypti mosquitoes. J Insect Physiol. 1968, 14: 719-728. 10.1016/0022-1910(68)90230-8.PubMedGoogle Scholar
- Thomson MC, Connor SJ, Quinones ML, Jawara M, Todd J, Greenwood BM: Movement of Anopheles gambiae sl malaria vectors between villages in The Gambia. Med Vet Entomol. 1995, 9: 413-419. 10.1111/j.1365-2915.1995.tb00015.x.PubMedGoogle Scholar
- Detinova TS: Age-grouping Methods in Diptera of Medical Importance with Special Reference to Some Vectors of Malaria. Monogr Ser World Health Organ. 1962, 47: 13-191.PubMedGoogle Scholar
- Jepson WF, Moutia A, Courtois C: The malaria problem in Mauritius: The binomics of Mauritian anophelines. Bull Entomol Res. 1947, 38: 177-208. 10.1017/S0007485300030273.PubMedGoogle Scholar
- 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: 375-381.PubMedGoogle Scholar
- 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. Entomol Exp Appl. 1992, 63: 265-271. 10.1111/j.1570-7458.1992.tb01583.x.Google Scholar
- 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 Trop. 2009, 109: 118-123. 10.1016/j.actatropica.2008.09.025.PubMedGoogle Scholar
- Armstrong JA, Bransby-Williams WR: The maintenance of a colony of Anopheles gambiae, with observations on the effects of changes in temperature. Bull World Health Organ. 1961, 24: 427-435.PubMed CentralPubMedGoogle Scholar
- Martens WJM, Niessen LW, Rotmans J, Jetten TH, McMichael AJ: Potential impact of global climate change on malaria risk. Environ Health Perspect. 1995, 103: 458-464. 10.1289/ehp.95103458.PubMed CentralPubMedGoogle Scholar
- Thomson MC, Mason SJ, Phindela T, Connor SJ: Use of rainfall and sea surface temperature monitoring for malaria early warning in Botswana. Am J Trop Med Hyg. 2005, 73: 214-221.PubMedGoogle Scholar
- Muirhead-Thomson RC: The ecology of vector snail habitats and mosquito breeding-places: The experimental approach to basic problems. Bull World Health Organ. 1958, 19: 637-659.PubMed CentralPubMedGoogle Scholar
- 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. Malar J. 2007, 6: 10-10.1186/1475-2875-6-10.PubMed CentralPubMedGoogle Scholar
- 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: 328-332. 10.1111/j.1365-2915.1993.tb00700.x.PubMedGoogle Scholar
- Hogg JC, Thompson MC, Hurd H: Comparative fecundity and associated factors for two sibling species of the Anopheles gambiae complex occuring sympatrically in The Gambia. Med Vet Entomol. 1996, 10: 385-391. 10.1111/j.1365-2915.1996.tb00761.x.PubMedGoogle Scholar
- 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: 639-645.PubMedGoogle Scholar
- Smith DL, Dushoff J, Snow RW, Hay SI: The entomological inoculation rate and Plasmodium falciparum infection in African children. Nature. 2005, 438: 492-495. 10.1038/nature04024.PubMed CentralPubMedGoogle Scholar
- Smith DL, McKenzie FE, Snow RW, Hay SI: Revisiting the basic reproductive number for malaria and its implications for malaria control. PLoS Biol. 2007, 5: e42-10.1371/journal.pbio.0050042.PubMed CentralPubMedGoogle Scholar
- Martens WJM, Jetten TH, Rottmans J, Niessen LW: Climate change and vector-borne diseases: A global modelling perspective. Glob Environ Change. 1995, 5: 195-209. 10.1016/0959-3780(95)00051-O.Google Scholar
- Martens WJM, Jetten TH, Focks DA: Sensitivity of malaria, schistomiasis and dengue to global warming. Clim Change. 1997, 35: 145-156. 10.1023/A:1005365413932.Google Scholar
- Craig MH, Snow RW, le Sueur D: A climate-based distribution model of malaria transmission in sub-Saharan Africa. Parasitol Today. 1999, 15: 105-111. 10.1016/S0169-4758(99)01396-4.PubMedGoogle Scholar
- Linard C, Gilbert M, Snow RW, Noor AM, Tatem AJ: Population distribution, settlement patterns and accessibility across Africa in 2010. PLoS ONE. 2012, 7: e31743-10.1371/journal.pone.0031743.PubMed CentralPubMedGoogle Scholar
- Balk DL, Deichmann U, Yetman G, Pozzi F, Hay SI, Nelson A: Determining global population distribution: methods, applications and data. Adv Parasitol. 2006, 62: 119-156.PubMed CentralPubMedGoogle Scholar
- Dekker T, Steib B, Carde RT, Geier M: L-lactic acid: a human-signifying host cue for the anthropophilic mosquito Anopheles gambiae. Med Vet Entomol. 2002, 16: 91-98. 10.1046/j.0269-283x.2002.00345.x.PubMedGoogle Scholar
- Mahande A, Mosha F, Mahande J, Kweka E: Feeding and resting behaviour of malaria vector, Anopheles arabiensis with reference to zooprophylaxis. Malar J. 2007, 6: 100-10.1186/1475-2875-6-100.PubMed CentralPubMedGoogle Scholar
- Killeen GF, McKenzie FE, Foy BD, Bøgh C, Beier JC: The availability of potential hosts as a determinant of feeding behaviours and malaria transmission by African mosquito populations. Trans R Soc Trop Med Hyg. 2001, 95: 469-476. 10.1016/S0035-9203(01)90005-7.PubMed CentralPubMedGoogle Scholar
- Bøgh C, Clarke SE, Pinder M, Sanyang F, Lindsay SW: Effect of passive zooprophylaxis on malaria transmission in The Gambia. J Med Entomol. 2001, 38: 822-828. 10.1603/0022-2585-38.6.822.PubMedGoogle Scholar
- Koella JC, SÖrensen FL, Anderson RA: The malaria parasite, Plasmodium falciparum, increases the frequency of multiple feeding of its mosquito vector, Anopheles gambiae. Proc R Soc Lond B Biol Sci. 1998, 265: 763-768. 10.1098/rspb.1998.0358.Google Scholar
- Dye C, Hasibeder G: Population dynamics of mosquito-borne disease: effects of flies which bite some people more frequently than others. Trans R Soc Trop Med Hyg. 1986, 80: 69-77. 10.1016/0035-9203(86)90199-9.PubMedGoogle Scholar
- Knols BGJ, de Jong R, Takken W: Differential attractiveness of isolated humans to mosquitoes in Tanzania. Trans R Soc Trop Med Hyg. 1995, 89: 604-606. 10.1016/0035-9203(95)90406-9.PubMedGoogle Scholar
- Mukabana WR, Takken W, Coe R, Knols BGJ: Host-specific cues cause differential attractiveness of Kenyan men to the African malaria vector Anopheles gambiae. Malar J. 2002, 1: 17-10.1186/1475-2875-1-17.PubMed CentralPubMedGoogle Scholar
- Rickman LS, Jones TR, Long GW, Paparello S, Schneider I, Paul CF, Beaudoin RL, Hoffman SL: Plasmodium falciparum-infected Anopheles stephensi inconsistently transmit malaria to humans. Am J Trop Med Hyg. 1990, 43: 441-445.PubMedGoogle Scholar
- Beier JC, Oster CN, Onyango FK, Bales JD, Sherwood JA, Perkins PV, Chumo DK, Koech DV, Whitmire RE, Roberts CR, Diggs CL, Hoffman SL: Plasmodium falciparum incidence relative to entomologic inoculation rates at a site proposed for testing malaria vaccines in western Kenya. Am J Trop Med Hyg. 1994, 50: 529-536.PubMedGoogle Scholar
- Shute PG, Maryon M: A study of gametocytes in a West African strain of Plasmodium falciparum. Trans R Soc Trop Med Hyg. 1951, 44: 421-438. 10.1016/S0035-9203(51)80020-8.PubMedGoogle Scholar
- Miller MJ: Observations on the natural history of malaria in the semi-resistant West African. Trans R Soc Trop Med Hyg. 1958, 52: 152-168. 10.1016/0035-9203(58)90036-1.PubMedGoogle Scholar
- Hawking F, Wilson ME, Gammage K: Evidence for cyclic development and short-lived maturity in the gametocytes of Plasmodium falciparum. Trans R Soc Trop Med Hyg. 1971, 65: 549-555. 10.1016/0035-9203(71)90036-8.PubMedGoogle Scholar
- Day KP, Hayward RE, Dyer M: The biology of Plasmodium falciparum transmission stages. Parasitol. 1998, 116 (Suppl): S95—S109-Google Scholar
- Le Prince JA: Mosquito control in relation to impounded water supply. J Am Water Works Assoc. 1927, 17: 31-36.Google Scholar
- 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 Journal. 2009, 8: 62-10.1186/1475-2875-8-62.PubMed CentralPubMedGoogle Scholar
- Haque U, Hashizume M, Glass GE, Dewan AM, Overgaard HJ, Yamamoto T: The role of climate variability in the spread of malaria in Bangladeshi highlands. PLoS ONE. 2010, 5: e14341-10.1371/journal.pone.0014341.PubMed CentralPubMedGoogle Scholar
- Desconnets JC, Taupin J, Lebel T, Leduc C: Hydrology of the HAPEX-Sahel Central Super-Site: surface water drainage and aquifer recharge through the pool systems. J Hydrol. 1997, 188: 155-178.Google Scholar
- Gianotti RL, Bomblies A, Eltahir EAB: Hydrologic modeling to screen potential environmental management methods for malaria vector control in Niger. Water Resour Res. 2009, 45: W08438-10.1029/2008WR007567.Google Scholar
- Martin-Rosales W, Leduc C: Dynamiques de vidange d’une mare temporaire au Sahel: l’exemple de Banizoumbou (Sud-Ouest du Niger). Comptes Rendus Geoscience. 2003, 335 (5): 461-468. 10.1016/S1631-0713(03)00059-2.Google Scholar
- Hayashi M, Kamp Van der G: Simple equations to represent the volume-area-depth relations of shallow wetlands in small topographic depressions. J Hydrol. 2000, 237: 74-85. 10.1016/S0022-1694(00)00300-0.Google Scholar
- Brooks RT, Hayashi M: Depth-area-volume and hydroperiod relationships of ephemeral (vernal) forest pools in southern New England. Wetlands. 2002, 22 (2): 247-255. 10.1672/0277-5212(2002)022[0247:DAVAHR]2.0.CO;2.Google Scholar
- Soti V, Puech C, Lo Seen D, Bertran A, Vignolles C, Mondet B, Dessay N, Tran A: The potential for remote sensing and hydrologic modelling to assess the spatio-temporal dynamics of ponds in the Ferlo Region (Senegal). Hydrol Earth Sys Sci Discuss. 2010, 14: 1449-1464. 10.5194/hess-14-1449-2010.Google Scholar
- Wanjala CL, Githeko AK, Waitumbi JN: Assessing the impact of topography on malaria exposure and malaria epidemic sensitivity in the Western Kenya highlands. Malar J. 2010, 9 (Suppl 2): P59-10.1186/1475-2875-9-S2-P59.PubMed CentralGoogle Scholar
- 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: 162-172. 10.1603/0022-2585-39.1.162.PubMedGoogle Scholar
- Munga S, Minakawa N, Zhou G, Barrack OOJ, Githeko AK, Yan G: Effects of larval competitors and predators on oviposition site selection of Anopheles gambiae sensu stricto. J Med Entomol. 2006, 43: 221-224. 10.1603/0022-2585(2006)043[0221:EOLCAP]2.0.CO;2.PubMedGoogle Scholar
- Paaijmans KP, Jacobs AFG, Takken W, Heusinkveld BG, Githeko AK, Dicke M, Holtslag AAM: Observations and model estimates of diurnal water temperature dynamics in mosquito breeding sites in western Kenya. Hydrol Process. 2008, 22: 4789-4801. 10.1002/hyp.7099.Google Scholar
- Paaijmans KP, Takken W, Githeko AK, Jacobs AFG: The effect of water turbidity on the near-surface water temperature of larval habitats of the malaria mosquito Anopheles gambiae. Int J Biometeorol. 2008, 52: 747-753. 10.1007/s00484-008-0167-2.PubMedGoogle Scholar
- Jacobs AFG, Heusinkveld BG, Kraai A, Paaijmans KP: Diurnal temperature fluctuations in an artificial small shallow water body. Int J Biometeorol. 2008, 52: 271-280. 10.1007/s00484-007-0121-8.PubMed CentralPubMedGoogle Scholar
- Paaijmans KP, Heusinkveld BG, Jacobs AFG: A simplified model to predict diurnal water temperature dynamics in a shallow tropical water pool. Int J Biometeorol. 2008, 52: 797-803. 10.1007/s00484-008-0173-4.PubMedGoogle Scholar
- Ermert V: Risk assessment with regard to the occurrence of malaria in Africa under the influence of observed and projected climate change. PhD thesis. 2010, Cologne, Germany: University of Cologne, Available online at [http://kups.ub.uni-koeln.de/3109/]Google Scholar
- Carnevale P, Robert V: Introduction of irrigation in Burkina Faso and its effect on malaria transmission. Effects of Agricultural Development on Vector-borne Diseases. Edited by: FAO. 1987, Rome: FAO, 57-67.Google Scholar
- Lochouarn L, Gazin I: La transmission du paludisme dans la ville de Bobo-Dioulasso (Burkina Faso). Ann Soc Belg Med Trop. 1993, 73: 287-296.PubMedGoogle Scholar
- Robert V, Gazin P, Boudin C, Molez JF, Oueddraogo V, Carnevale P: La transmission du paludisme en zone de savane arboreé et en zone rizicole des Bobo-Dioulasso (Burkina Faso). Ann Soc Belg Med Trop. 1985, 65: 201-214.PubMedGoogle Scholar
- Robert V, Gazin P, Benasseni R, Carnevale P: Le paludisme urbain à Bobo-Dioulasso (Burkina Faso). Urbanisation et santé dans le Tiers Monde: transition épidémiologique, changement social et soins de santé primaires. Edited by: Salem G, Emile J. 1989, Paris: ORSTROM, 181-185.Google Scholar
- Gazin P, Goncalves K, Koné B, Lochouam L: Incidence des accès palustres dans un quartier de la ville de Bobo-Dioulasso (Burkina Faso). Bull Soc Pathol Exot. 1996, 89: 200-203.PubMedGoogle Scholar
- Robert V, Gazin P, Boudin C, Molez JF, Oueddraogo V, Carnevale P: La transmission du paludisme en zone de savane arboreé et en zone rizicole des Bobo-Dioulasso (Burkina Faso). Ann Soc Belg Med Trop. 1985, 65: 201-214.PubMedGoogle Scholar
- Boudin C, Robert V, Carnevale P, Ambroise TP: Epidemiology of Plasmodium falciparum in a rice field and a savanna area in Burkina Faso: seasonal fluctuations of gametocytaemia and malaria infectivity. Ann Trop Med Parasitol. 1991, 85: 377-385.PubMedGoogle Scholar
- Boudin C, Robert V, Carnevale P, Ambroise-Thomas P: Epidemiology of Plasmodium falciparum in a rice field and a savanna area in Burkina Faso. Comparative study on the acquired immunoprotection in native populations. Acta Trop. 1992, 51: 103-111. 10.1016/0001-706X(92)90052-Y.PubMedGoogle Scholar
- Boudin C, Robert V, Verhave JP, Carnevale P, Ambroise-Thomas P: Plasmodium falciparum and P. malariae epidemiology in a West African village. Bull World Health Organ. 1991, 69: 199-205.PubMed CentralPubMedGoogle Scholar
- Robert V, Carnevale P: Influence of deltamethrin treatment of bed nets on malaria transmission in the Kou valley, Burkina Faso. Bull World Health Organ. 1991, 69: 735-740.PubMed CentralPubMedGoogle Scholar
- Love T: The climate prediction center rainfall algorithm version 2. Tech. rep., NOAA Climate Prediction Center, 2002. available at [http://www.cpc.noaa.gov/products/fews/]
- Dee DP, Uppala SM, Simmons AJ, Berrisford P, Poli P, Kobayashi S, Andrae U, Balmaseda MA, Balsamo G, Bauer P coauthors: The ERA-Interim reanalysis: configuration and performance of the data assimilation system. Q J R Meteorol Soc. 2011, 137 (656): 553-597. 10.1002/qj.828.Google Scholar
- Giorgi F, Francisco R, Pal J, Effects of a subgrid-scale topography and land use scheme on the simulation of surface climate and hydrology. Part I: Effects of temperature and water vapor disaggregation: Effects of a subgrid-scale topography and land use scheme on the simulation of surface climate and hydrology. Part I: Effects of temperature and water vapor disaggregation. J Hydrometeor. 2003, 4: 317-333. 10.1175/1525-7541(2003)4<317:EOASTA>2.0.CO;2.Google Scholar
- Smith WHF, Sandwell DT: Global sea floor topography from satellite altimetry and ship depth soundings. Science. 1997, 277 (5334): 1956-1962. 10.1126/science.277.5334.1956.Google Scholar
- Hay SI, Snow RW: The malaria atlas project: developing global maps of malaria risk. PLoS Med. 2006, 3: e473-10.1371/journal.pmed.0030473.PubMed CentralPubMedGoogle Scholar
- Hay SI, Guerra CA, Gething PW, Patil AP, Tatem AJ, Noor AM, Kabaria CW, Manh BH, Elyazar IRF, Brooker SJ, Smith DL, Moyeed RA, Snow RW: A world malaria map Plasmodium falciparum endemicity in 2007. PLoS Med. 2009, 6: e1000048-10.1371/journal.pmed.1000048.PubMed CentralPubMedGoogle Scholar
- Gething PW, Patil AP, Smith DL, Guerra CA, Elyazar IRF, Johnston GL, Tatem AJ, Hay SI: A new world malaria map: Plasmodium falciparum endemicity in 2010. Malar J. 2011, 10: 378-10.1186/1475-2875-10-378.PubMed CentralPubMedGoogle Scholar
- Patil AP, Gething PW, Piel FB, Hay SI: Bayesian geostatistics in health cartography: the perspective of malaria. Trends Parasitol. 2011, 27: 246-253. 10.1016/j.pt.2011.01.003.PubMed CentralPubMedGoogle Scholar
- Ermert V, Fink AH, Morse AP, Paeth H: Environ Health Perspect. 2012, 120: 77-84.Google Scholar
- Thomson RC: The reactions of mosquitoes to temperature and humidity. Bull Entomol Res. 1938, 29: 125-140. 10.1017/S0007485300026158.Google Scholar
- Okiro EA, Hay SI, Gikandi PW, Sharif SK, Noor AM, Peshu N, Marsh K, Snow RW: The decline in paediatric malaria admissions on the coast of Kenya. Malar J. 2007, 6: 151-10.1186/1475-2875-11-2.PubMed CentralPubMedGoogle Scholar
- O’Meara WP, Bejon P, Mwangi TW, Okiro EA, Peshu N, Snow RW, Newton CRJC, Marsh K: Effect of a fall in malaria transmission on morbidity and mortality in Kilifi, Kenya. J Lancet. 2008, 372: 1555-1562. 10.1016/S0140-6736(08)61655-4.Google Scholar
- Mbogo CM, Mwangangi JM, Nzovu J, Gu W, Yan G, Gunter JT, Swalm C, Keating J, Regens JL, Shililu JI, Githure JI, Beier JC: Spatial and temporal heterogeneity of Anopheles mosquitoes and Plasmodium falciparum transmission along the Kenyan coast. Am J Trop Med Hyg. 2003, 68: 734-742.PubMedGoogle Scholar
- Dongus S, Nyika D, Kannady K, Mtasiwa D, Mshinda H, Gosoniu L, Drescher A, Fillinger U, Tanner M, Killeen G, Castro MC: Urban agriculture and Anopheles habitats in Dar es Salaam, Tanzania. Geospat Health. 2009, 3: 189-210.PubMedGoogle Scholar
- Atieli H, Menya D, Githeko A, Scott T: House design modifications reduce indoor resting malaria vector densities in rice irrigation scheme area in western Kenya. Malar J. 2009, 8: 108-10.1186/1475-2875-8-108.PubMed CentralPubMedGoogle Scholar
- Manga L, Fondjo E, Carnevale P, Robert V: Importance of low dispersion of Anopheles gambiae (Diptera: Culicidae) on malaria transmission in hilly towns in south Cameroon. J Med Entomol. 1993, 30: 936-938.PubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.