Development of a new version of the Liverpool Malaria Model. II. Calibration and validation for West Africa
© Ermert et al; licensee BioMed Central Ltd. 2011
Received: 17 September 2010
Accepted: 16 March 2011
Published: 16 March 2011
In the first part of this study, an extensive literature survey led to the construction of a new version of the Liverpool Malaria Model (LMM). A new set of parameter settings was provided and a new development of the mathematical formulation of important processes related to the vector population was performed within the LMM. In this part of the study, so far undetermined model parameters are calibrated through the use of data from field studies. The latter are also used to validate the new LMM version, which is furthermore compared against the original LMM version.
For the calibration and validation of the LMM, numerous entomological and parasitological field observations were gathered for West Africa. Continuous and quality-controlled temperature and precipitation time series were constructed using intermittent raw data from 34 weather stations across West Africa. The meteorological time series served as the LMM data input. The skill of LMM simulations was tested for 830 different sets of parameter settings of the undetermined LMM parameters. The model version with the highest skill score in terms of entomological malaria variables was taken as the final setting of the new LMM version.
Validation of the new LMM version in West Africa revealed that the simulations compare well with entomological field observations. The new version reproduces realistic transmission rates and simulated malaria seasons are comparable to field observations. Overall the new model version performs much better than the original model. The new model version enables the detection of the epidemic malaria potential at fringes of endemic areas and, more importantly, it is now applicable to the vast area of malaria endemicity in the humid African tropics.
A review of entomological and parasitological data from West Africa enabled the construction of a new LMM version. This model version represents a significant step forward in the modelling of a weather-driven malaria transmission cycle. The LMM is now more suitable for the use in malaria early warning systems as well as for malaria projections based on climate change scenarios, both in epidemic and endemic malaria areas.
The World Health Organization (WHO) estimated that about two billion people, that is more than 40% of the total world population, are exposed to malaria . Estimates in terms of 2009 revealed that this mosquito-borne disease causes about 225 million cases and 781,000 deaths annually. At least 90% of the worldwide malaria deaths occur in sub-Saharan Africa .
Malaria is a climate-sensitive tropical disease and hence climate exerts a strong impact upon the distribution of malaria transmission in space and time . An assessment of current and future malaria risk is an important topic in the area of research relating climate to disease risk . Reliable forecasts of epidemic malaria outbreaks on seasonal timescales  and assessments of disease vulnerability over decadal timescales are needed . However, this requires the production of a weather-malaria modelling system. World-leading numerical weather forecast centres have already demonstrated useful skill in forecasts far beyond a month lead-time for some tropical regions , and climate projections are becoming more reliable . Advances have also been achieved in weather- or climate-driven malaria modelling. For example, Hoshen and Morse  introduced the Liverpool Malaria Model (LMM), which is a mathematical-biological model of malaria parasite dynamics driven by daily temperature and precipitation data. However, further progress must be obtained in order to enable skilful malaria simulations for epidemic and endemic malaria regions based on meteorological data. For this reason, the present study introduces the LMM version of 2010 (henceforth called LMM2010), which simulates a more realistic spread of malaria in space and time and is hence a useful tool for a weather- or climate-disease modelling system.
LMM parameters and mathematical formulations
humid degree days of the gonotrophic cycle
37.1 degree days
37.1 degree days
dry degree days of the gonotrophic cycle
65.4 degree days
65.4 degree days
humid gonotrophic temperature threshold
dry gonotrophic temperature threshold
10-day accumulated precipitation threshold
rainfall laying multiplier
number of produced eggs per female mosquito
number of oviposited eggs per female mosquito
Eq. Two in 
lower threshold of unsuitable rainfall conditions (fuzzy distribution model)
most suitable rainfall condition (fuzzy distribution model)
upper threshold of unsuitable rainfall conditions (fuzzy distribution model)
cap on the number of fertile mosquitoes
mosquito mature age
rainfall independent immature daily mosquito survival probability
daily immature mosquito survival probability (in %)
Eq. Three in 
Eq. Four in 
daily mosquito survival probability (in %)
Martens I (see )
Martens II (see )
dry season mosquito survival probability shift
degree-days of the sporogonic cycle
111 degree days
111 degree days
sporogonic temperature threshold
human blood index
mosquito-to-human transmission efficiency
adult-child conversion rate
human infectious age
daily human recovery rate
0.0284 day 1
fraction of gametocyte carriers
human-to-mosquito transmission efficiency
trickle of the number of added infectious mosquitoes
Time series of meteorological stations
Regional climate simulations
Driving the LMM with gridded data is needed for the validation of the simulated malaria spread in Africa. Due to the lack of gridded daily temperature and rainfall data from ground observations, data from a so-called regional climate model are used. Such a model is able to more accurately represent spatial variations of the atmosphere than so-called global atmosphere-ocean general circulation models.
Gridded LMM simulations for Africa were based on meteorological data from three ensemble runs between 1960 and 2000 from the REgional MOdel (REMO) . REMO is a hydrostatic, limited area model using primitive equations, which are solved on 20 hybrid atmospheric levels . At the lateral and lower boundaries REMO was nested into ECHAM5/MPI-OM (European Centre HAmburg Model, 5th generation/Max-Planck-Institute-Ocean Model) global coupled climate model simulations that were forced with the observed greenhouse gas increase [17–19]. The space-time resolution of REMO data used in the present study was 0.5° and daily values, respectively.
Expected systematic model errors in annual precipitation and in simulated temperatures were considerably reduced using bias correction, ensuring realistic data input for the malaria simulations . The bias-corrected data are later used to drive gridded malaria runs in order to compare the original and new LMM versions.
Entomological and parasitological observations
In terms of malaria modelling, the entomological and parasitological data are of particular interest since malaria models have to undergo validation procedures. The present study was restricted to West Africa since this area was one of the focus regions of the IMPETUS project (Integrated Approach to the Efficient Management of Scarce Water Resources in West Africa) . Numerous published malaria observations were extracted from the literature.
It was shown that population density is an important factor of malaria transmission [21, 22]. In addition, crop irrigation strongly impacts malaria and its seasonality [23–26]. Also malaria and vector control such as the usage of impregnated bed nets are able to reduce transmission [27–31]. In the present study, field observations were only included for which a rural environment can be assumed from the information given in the publication (see Additional file 2 which includes also data from non-rural sites). Urban and irrigated areas and those when permanent streams influenced transmission were excluded (except for the locations in Cameroon, where this is nearly always the case). The few sites subject to vector control or other local or regional interventions were not used. This procedure ensured that other factors than climate, which influence malaria transmission, were excluded to the extent possible.
Entomological malaria field studies frequently sample biting mosquitoes on humans . Standard measures include the so-called Human Biting Rate (HBR), which is the number of mosquito bites per human per time. However, only female mosquitoes with sporozoites in their salivary glands are able to infect humans. This fraction of the biting females is called CircumSporozoite Protein Rate (CSPR). Multiplying HBR with CSPR results in the so-called Entomological Inoculation Rate (EIR), which is defined as the number of infectious mosquito bites per human per time . Only months revealing infectious mosquito bites (EIR values above zero) are usually used to define the malaria season at a certain location . By contrast, parasitological malaria studies usually measure the asexual parasite ratio (PR) representing the proportion of the survey population, which is positive for the malaria parasite . The literature was, therefore, reviewed in terms of malaria field studies (see Additional file 2).
Wherever applicable, observations were taken into account with regard to eleven different entomological and parasitological variables (these are: the annual Human Biting Rate (HBR a ), the annual Entomological Inoculation Rate (EIR a ), the annual mean CircumSporozoite Protein Rate (CSPR a ), the length, start, and end of the malaria season (Seas, SSeas, and ESeas, respectively), the length of the main malaria season (MSeas; i.e. the number of months in which 75% of EIR a is recorded ), the month of maximum transmission (XSeas; i.e. the month with the largest entomological inoculation rate), the annual mean, maximum, and minimum of the asexual parasite ratio (PR a , PR max; a , and PR min, a , respectively). Included are the literature references as well as some basic features such as the name, geographical position, land-use of the study site, and the time period. Note that all these malaria variables can be computed from the LMM output.
Definition of malaria seasonality
The definition of the malaria season is based on the monthly Entomological Inoculation Rate (EIR m ), which is observed from field studies. The introduction of the malaria parasite into the LMM is assured by a constant influx of new infectious mosquitoes resulting in EIR values that might exceed those of low transmission areas. In order to compare 'truly' modelled EIR values with field observations, the artificial EIR values were removed via two separate LMM runs. The standard run results in a mixture of bites from the added infectious mosquitoes and those which were infected when they bit infectious humans in the simulation. The artificial EIR values were produced in a second run when only the added infectious mosquitoes were biting in the model, achieved by setting the number of produced eggs per female mosquito (#E p ) to zero. Here, the oviposition of Anopheles females is prohibited and the mosquito population is therefore not able to grow. The subtraction of the EIR values of the second run from that of the standard run results finally in the elimination of artificial EIR values. Note that the same procedure was used with regard to the computation of the human biting rate (HBR).
In field studies, the malaria season is mostly determined on a monthly scale and is usually referred to months with EIR m values above zero . For this reason, in the model the malaria season starts (SSeas) by definition in the first month with an EIR m value of at least 0.01 infectious bites per human. According to the formulation of the model, the arbitrary EIR value of 0.01 means that during a 30-day month at least one out of the 100 modelled humans is bitten by an infectious mosquito. In West Africa, the malaria season usually ends when the number of mosquitoes decreases at the end of the rainy season. Therefore, mosquito biting is reduced to such low numbers that the EIR values are reduced to zero indicating that transmission has ceased. Consequently, the last month during the transmission period defines the end of the malaria season (ESeas). Some individual years also reveal year-round or even no malaria transmission for certain locations. The length of the malaria season (Seas) is therefore the number of months with EIR m reaching at least 0.01 infectious bites. For each site or grid box additionally the length of the main transmission season (MSeas) is defined as the number of months in which 75% of EIR a is transmitted, an index which was used by Hay et al. . Where possible the month with the maximum malaria transmission (XSeas) is identified as the month with the highest EIR m value.
The use of the transmission threshold of 0.01 infectious mosquito bites per human per month ensures the attainment of a reasonable transmission level in the model that can be compared with observations from the literature. However, the definition of the malaria season in the model might not be directly comparable to field studies since observations are subject to a certain detection limit. That is due to the fact that field experiments do not continuously measure biting rates (at best two times in each week of the field campaign) and that these field studies do not account for every human of the population of the study site.
Definition of the validation
For the model validation, the LMM runs are driven by weather observations from meteorological stations in West Africa. The malaria model simulations are compared against observations from more than 200 sites using up to eleven different entomological and parasitological variables. The skill of the malaria simulations is measured by a problem-adapted skill score (see below). The performance of every particular set of parameter settings of the model is measured for each of the eleven malaria variables by this skill score (SC(x)). In addition, the performance of a set of parameter settings with regard to all malaria variables is quantified by the sum of the eleven skill scores (SC(all)). This approach enables the production of a ranking of different sets of parameter settings. Finally, the most suitable set of parameter settings in terms of the skill score determines the new LMM version.
The LMM was run with meteorological time series covering 1973 to 2006 discussed above. Therefore, each simulation at a particular station produced 34 annual values for every considered variable (e.g. EIR a ). Since malaria conditions are different between rural and urban areas , only data from rural field sites were used for the LMM validation. On the other hand, entomological data are never measured continuously because of the amount of work required in practice . For this reason, no single long time series of malaria observations exist. However, field observations of a subset of the nine malaria variables were frequently found in the area of the considered weather stations. Unfortunately, such malaria data has not been detected for all the used stations. For example, in the vicinity of Parakou (Benin) no field observations were extracted from the literature.
A further data issue that must be taken into account arises from the fact that the malaria field measurements were not performed at the locations of the weather stations. However, observation sites are available in the vicinity of the weather stations meaning that they are located in same climatic zones (Figure 1). In this context, it must be noted that precipitation can differ greatly between locations just few kilometres apart, but meteorological stations are much more coarsely distributed [39, 40]. Also different environmental conditions (e.g. state of land surface) complicate the comparison between simulations and field observations. Shaman and Day  allude to the mismatch between the scales at which disease vectors respond to hydrologic variability and the scales at which hydrologic variability is actually monitored. For all these reasons, a reproduction of the field values of particular years, for example, of transmission rates is not expected by the LMM. In addition, it implies that paired annual correlations at single stations are not useful. Therefore, the usage of standard skill scores, used to determine the skill in weather forecasting , is not possible with the data sets that are available. In addition, statistical moments of the model and observational malaria data cannot be statistically meaningful tested due to the low sample size of the field observations. In order to overcome these issues, a subjective problem-adapted skill score is defined in the following section.
Definition of a skill score
A special validation procedure must be defined. It must be considered that simulations will not be able to reproduce single field observations. However, it is possible that the LMM simulates about the correct mean and variability of malaria variables, as observed in the field, since the weather station and the assigned field sites encounter similar climate conditions. In the following, location names refer to the name of the weather station, even if malaria observations from a nearby field site are discussed. The assignment of the weather stations to the malaria study sites is given in Additional file 2.
Criteria in terms of the evaluation of LMM simulations
Any observation is included in the simulated range
Every observation is included in the simulated range
The observed median is included in the simulated range
HBR a EIR a CSPR a Seas MSeas PR a PR min,a PR max,a
median quartile enclosure
The observed median is located within the lower and upper quartile of the simulations
HBR a EIR a CSPR a Seas MSeas PR a PR min,a PR max,a
The simulations exceed the one and a half time maximum of all field observations
HBR a EIR a
I: The observations as well as the simulations show the same month maximum occurrence of the monthly entomological inoculation rate (EIR m ).
SSeas ESeas XSeas
II: The majority of the observations and simulations show no or yearround transmission, respectively.
III: The observations and simulations reveal mostly no transmission.
IV: The month showing the most field observations of multiple years of EIR m is identical to the particular month resulting from the simulation.
SSeas/ESeas: Criterion I or II
XSeas: Criterion III or IV
The evaluation of the model runs is performed by means of a skill score (SC), which considers the number of available observations for each of the eleven malaria variables. Stations with at least five observations contribute at maximum three or four points to the calculated skill score. Weather stations with fewer field measurements add one to three points to the skill score (cf. Table 2). No score is added when stations reveal no observations. The validation method also makes allowances for the uncertainty of the year-to-year variability. A proper analysis would only be possible for numerous available observations per station. In fact, an estimate of the frequency distribution of observed malaria parameters is not feasible from the few records available.
The computed skill score is based on the probability that observations, assigned to particular weather stations, fall into certain ranges of simulations (Table 2). It is expected that any observation is included in the range of the model simulations with regard to this station data (criterion 1: overlap). Of course, the model version performs better when every observation is enclosed in the simulated range (criterion 2: enclosure). Where there are at least three available records the observed median is also calculated. The confidence in the run increases when the observed median is contained in the range of the model simulations (criterion 3: median enclosure). The reliability of the median estimate increases for sites where five observations or more field measurements exist. Only for at least five available records the condition has to be met that the observed median falls within the range of the lower (25th percentile) and upper quartile (75th percentile) of the 34 simulated annual values (criterion 4: median quartile enclosure). For model integrations with different parameter settings, every fulfilled criterion at a weather station adds one point to the skill score of a particular malaria variable.
This rating system might favour model versions generating unrealistic high values and a strong interannual variation. This fact is countered by another criterion that eliminates unrealistically high entomological values (criterion 5: penalty). Five penalty points are applied to the skill score of HBR a and EIR a , when any simulated value exceeds one and a half times the maximum of all available field measurements (see Table 2). This threshold seems to be a reasonable measure for the restriction of simulated values. Sets of parameter settings leading to unrealistic high biting rates (HBR a and EIR a values) are rejected.
The third and fourth criteria are not calculable for three of the variables under certain circumstances (cf. Table 2): the start (SSeas) and end (ESeas) of the malaria season (for which the criteria are undefined in the cases of no transmission or year-round transmission), and the month of maximum transmission (XSeas) (for which they are undefined in the case of no transmission). For these variables, the maximum of the frequency distribution is compared instead (criterion 6: frequency). Regarding SSeas and ESeas, a model version is given an additional score of one point when both the observations as well as the simulations show the same month of maximum occurrence. A point is also gained when the majority of the observations and simulations show no or year-round transmission, respectively. For XSeas, a model version can achieve one additional point for a weather station either if both the observations and simulations reveal mostly no transmission or if the month showing the most field observations of multiple years of maximum EIR m values is identical to the particular month resulting from the simulation (e.g. if for a single station the strongest transmission values are both mostly observed and simulated for August).
For each weather station, a skill score is determined for every set of parameter settings and for every one of the eleven entomological and parasitological variables (e.g. for Ouagadougou the skill score in terms of EIR a , i.e. SC(EIR a ), is 4 for the LMM2010; see Figure 3). The skill scores from the weather stations add up to the skill score of a particular variable (SC(x); e.g. SC(EIR a ) is 41 for the LMM2010; see Figure 3). Finally, the sum of all skill scores of a particular set of parameter settings is calculated (i.e. SC(all), which is 279 for the LMM2010).
For each particular station and each variable a certain number of points are achievable (e.g. 54 points can be reached in terms of SC(EIRa)). As stated above this optimum depends on the availability of field observations (see Table 2), which differs for each weather station and each malaria variable. The minimum number of possible points is zero for stations, with no available field observations (e.g. see Parakou in Figures 3, 4, and 5). At maximum four points can be achieved, when at least five observations are available. For example, numerous field observations are available for Bobo-Dioulasso (see Figure 3, 4, and 5). In case of SSeas, ESeas, and XSeas at maximum only three points can be achieved since only three criteria are applied for these variables (see Table 2).
In the following, undefined parameters of the LMM are calibrated by means of West African field observations. A two step approach is used in order to reduce the degrees of freedom of the model within each step. The performance of the model runs is measured by the problem-adapted scoring system. At the end, the set of parameter settings revealing the highest skill score is selected. Finally, the simulations of the LMM2004 and LMM2010 are compared to field observations of the eleven entomological and parasitological malaria variables.
Calibration of the LMM
The majority of the model parameters were determined in the first part of this study . However, some parameters were not allocatable due to a large spread in their published values or due to the lack of data. These undetermined values were two of the three parameters (S and U2) of the fuzzy distribution model. Here, S defines the most suitable 10-day accumulated rainfall conditions for mosquito breeding, meaning rainfall provides open water surfaces, which tend not to be flushed out by heavy rainfall. U2 defines the upper limit of suitable rainfall conditions above which all breeding habitats are flushed out. Note that U1, which is the lower limit of suitable rainfall conditions, was simply set to 0 mm since under prolonged dry conditions ephemeral water bodies will dry up .
In addition, the number of produced eggs per female mosquito (#E p ), the cap on the number of fertile mosquitoes (CAP), as well as the shift off relative to the dry season mosquito survival probability (pd↓) are undetermined. Here, pd↓is used to reduce vector survival during dry atmospheric conditions . In the first part of this study , some parameter values were estimated from the range of published data and few were based on educated guesses, this might have resulted in some false estimation of these predefined parameters. However, the following calibration of the remaining parameters will largely compensate such inaccurate assessments.
As outlined above, the LMM calibration was undertaken for weather conditions at meteorological stations and included two general steps: (i) The initial experiment enabled a rough estimate of realistic parameter values and allowed the setting of two undetermined parameters. (ii) The second set of model runs permitted a final adjustment of the remaining model parameter settings. In order to simplify the calibration procedure, pd↓was initially set to zero. Altogether three different settings were used for U2 (i.e. the upper limit of suitable rainfall conditions in terms of larval breeding of the fuzzy distribution model). Ahumada et al.  defined extreme rainfall as more than 255 mm of cumulative rainfall throughout a period of three days. Their model markedly reduces the mosquito population under excessive rainfall. According to a 10-day period, which is used by the fuzzy distribution model, this suggests a U2 value of about 500-1000 mm. In fact, U2 was either set to 500, 750, or 1000 mm.
For the isolation of particular sets of parameter settings the remaining parameters (S, U2, #E p , and CAP) were simultaneously varied (S: 5, 10, 15, 20, and 30 mm; U2: 500, 750, and 1000 mm; #E p : 50, 75, 100, 125, and 150 eggs; CAP: 250, 500, 750, 1000, and 2000 fertile females). The various LMM simulation runs were compared to data from the eleven entomological and parasitological variables (see Additional file 2), these are: HBR a , EIR a , CSPRa, SSeas, ESeas, Seas, PR a , PR max, a , and PR min, a . The first 375 (originating from the combination of 5·3·5·5 settings) different LMM settings were ranked with regard to the skill score of all variables (SC(all)) and in terms of the two entomological variables HBRa and EIRa (i.e. SC(HBR a , EIR a ) = SC(HBR a )+SC(EIR a )). Here, SC(HBR a , EIR a ) is measuring the skill of the model simulation in terms of both HBR a and EIR a .
Overview of calibration experiments
S ∈ [5, 30]
SC(HBR a , EIR a )
⇒ S : = 10 mm
U 2 ∈ [500, 1000]
⇒ U 2 : = 500 mm
#E p ∈ [50, 150]
SC(HBR a , EIR a )
⇒ #E p ∈ [75, 125]
CAP ∈ [250, 2000]
⇒ CAP ∈ [300, 900]
CAP ∈ [300, 900]
SC(HBR a , EIR a )
⇒ CAP: 400 fertile
#E p ∈ [70, 130]
females, #Ep : = 120 eggs,
pd↓∈ [0, 10]
pd↓: = -10%
The evaluation of the performance of the 375 model versions also enabled a final setting of U2 since highest skill scores were exclusively generated by a value of 500 mm . The lowest skill scores were attained by settings with high values of CAP, #E p , and U2, as well as low values of S, which apparently evoke large mosquito populations (see Additional file 3 and step 1.2 in Table 3).
A closer analysis of the data revealed that malaria transmission rates in the Sahel are fairly sensitive to the setting of #E p (i.e. the number of produced eggs per female mosquito) . The median HBRa value rises, for example, at Linguère, Mopti, and Diourbel from less than 100 to several thousand bites per year, when #E p increases from 50 to 150 eggs. Unfortunately, only nine field observations of HRB a and EIRa are available north of 14°30'N. This fact impeded a proper determination or a further confinement of #E p . However, the LMM underestimates (overestimates) EIR a , when #E p is set to 50 (150) eggs (see Additional file 3 and step 1.3 in Table 3).
After confirming the parameter range of #E p (75-125 eggs) also CAP (i.e. the cap on the number of fertile mosquitoes) was more precisely determined (step 1.4 in Table 3). CAP is only of importance for comparatively large annual rainfall amounts. CAP markedly reduces the number of deposited eggs and hence biting rates under wet conditions. In fact, the reduction is markedly pronounced in the Sudanian zone, along the Guinean coast, and in Cameroon for low CAP values (≤750 fertile females). In contrast, large values of CAP (≥ 1000) cause fairly high numbers of non-infectious and infectious mosquito bites . The ranking relative to SC(all) shows that LMM versions reveal a small skill when CAP is set to 250 fertile females .
The basis of the second iteration of the LMM calibration are the conclusions taken from the first step. For this reason, only the setting of CAP, #E p , and pd↓needed to be varied. For the second set of runs, the #E p values were varied between 70 and 130 eggs (7 steps at increments of 10 eggs) and CAP (i.e. the cap on the number of fertile females) was set between 300 and 900 fertile females (13 steps at 50 fertile females). Five different values for pd↓(i.e. the reduction of dry season mosquito survival) were in addition utilised (0.0, -2.5, -5.0, -7.5, and -10.0%). The second set of model runs included altogether 455 (originating from the combination of 7·13·5 settings) different model sets of parameter settings.
The second set of runs shows that low (high) #E p and high (low) CAP values produce the highest skill scores in terms of the eleven variables (SC(all)). As stated above this is because #E p and CAP tend to compensate each other at the more humid locations (see Additional file 3). Particularly notable is the fact that various model versions exhibited comparatively high skill scores. This fact makes a final objective setting of the remaining parameters difficult. For simplicity, the model version with the highest skill score in terms of HBR a and EIR a (SC(HBR a , EIR a )) was chosen.
The set of parameter settings with #E p = 120 eggs, CAP = 400 fertile females, and pd↓= -10% produced the highest skill score in terms of HBR a and EIR a . A total of 78 from 106 possible points (73.6%) was reached by this model version. With regard to SC(all) altogether 279 from the 440 possible points were gained (Additional file 3).
The #E p value of 120 eggs is in the middle of observations (see Additional file One of Ermert et al. ). The pd↓value of -10% enables the simulation of a realistic end of the malaria season. Malaria transmission stops one to two months earlier in the model when this pd↓value is applied (cf. Figure 4). The correspondence with MARA (mapping Malaria Risk in Afrika project) maps  is improved by the earlier break of malaria transmission.
Validation of the LMM2010
LMM2004 and LMM2010 performance
criteria number (see Table 2)
criteria number (see Table 2)
SC(HBR a )
SC(CSPR a )
SC(EIR a )
SC(PR a )
SC(PR max, a )
SC(PR min, a )
Another result of the calibration is that the simulations capture the variability of malaria transmission. The interannual variability of EIR a , for example, is fairly large. For most stations the number of infectious mosquito bites fluctuates between values of less than 100 and several hundred. Unfortunately, long-term studies are rare and continuous observations from rural sites are only available from Ndiop for four years (Senegal; 13°41'N, 16°23'W). In this Sahelian village, malaria transmission varied in the mid of the 1990 s between seven and 63 infective bites per human per year . For this area, simulated EIR a values range from almost zero to about 158 infectious bites per human per year. Note that the LMM simulation refers to the meteorological data from Kaolack, which is located about 55 km to the northeast of Ndiop.
The simulation of the malaria seasonality by the LMM2010 is consistent with epidemic and endemic field observations, despite the fact that a part of the field observations reveals considerably heterogeneous values. The simulated length of the season agrees well with the observations (SC(Seas) = 31(42)). Only in one out of 16 cases is there no overlap (criterion 1) with observations. In three out of the five cases with at least five field observations, the optimum scoring of four points is achieved (see Additional file 4). In general, the season length shortens with decreasing annual rainfall. Short malaria seasons and no malaria transmission are simulated for the Sahel and year-round transmission is found, for instance, in Cameroon. Also the length of the main transmission season corresponds roughly to the observations (SC(MSeas) = 28(37)). At only two of the 14 weather stations no overlappings (criterion 1) are found. Criterion 3 (median enclosure) is fulfilled for each of the seven cases (Table 4).
The skill score in terms of the maximum transmission month (XSeaS) gains 23 from the 41 achieveable points (SC(XSeas) = 23(41); cf. Additional file 4). The month of maximum transmission overlaps (criterion 1) between simulations and observations in 15 out of 17 cases and all observations are enclosed in the simulation range (criterion 2) in six out of 16 cases. The simulated month with the maximum transmission therefore corresponds frequently with that of the observations, especially in the Sahel. However, in terms of XSeas, the model simulations disagree in Cameroon, at the weather stations of Yaoundé and Koundja for areas with year-round endemic malaria transmission.
LMM2010 versus LMM2004
The validation of the LMM2004 and LMM2010 runs based on weather station data from West Africa clearly demonstrates a significant improvement of the original model version (the P value is 0.0064; see Table 4). Except for ESeas, all skill scores are higher for the LMM2010 than for the LMM2004. Comparison of the LMM2004 simulations with field observations reveal two general features: (i) The LMM2004 fails to simulate malaria transmission in various malarious semi-arid regions and significantly overestimates the transmission in humid areas especially regions with endemic malaria territories. Note here that the LMM2004 was not designed to be used in endemic areas. (ii) Values of the parasite ratio, expressed in the three parasitological variables, are not optimal in the original model version. In terms of ESeas, the improvement of the model, especially in the Sahel, seems not to show up due to the lack of numerous observations of ESeas in this area. Nevertheless, the LMM2010 represents a significant step forward in the modelling of a weather-driven malaria transmission cycle.
The first feature (i) is a result of the simplified linear relationship between rainfall and the deposition of eggs by female mosquitoes in the LMM2004. In predominantly dry areas the model is not able to produce a reasonable size of the mosquito population. The LMM2004 hence fails to simulate malaria transmission in the northern part of the Sahel (Figure 6a & 6c). By contrast, in the LMM2010 runs, malaria transmission is found as far north as about 18°N (Figure 6b & 6d). This seems to be realistic since, for example, a definite malaria season was observed south of Agadez , which is located at about 17°N. Further the MARA project detected epidemic-prone areas up to 20°N .
In the second step, LMM runs were driven by gridded and bias-corrected REMO data for 1960-2000. Based on this input data the LMM2010 seems to reproduce a realistic picture with regard to transmission rates and the malaria seasonality for the whole of Africa including epidemic and endemic malaria areas (cf. Figure 6b & 6d). This is in contrast to the LMM2004 which was originally only constructed for epidemic malaria territories . This former model version simulates a tremendous number of mosquitoes and generates significantly too high EIR a values for humid areas such as the equatorial tropics (Figure 6a). HBR a partly exceeds 1,000,000 bites per annum in these endemic malaria areas . Note that HBR a field values are usually much lower than 50,000 bites (see Additional file 2). However, the original model version also reveals deficiencies in terms of the spread of malaria within epidemic regions. The LMM2004 fails to simulate malarious areas north of about 16°N.
The infectiousness of mosquitoes is underestimated by the LMM2004 since values of the circumsporozoite protein rate (here CSPR a ) are in general lower than 1%  resulting from the comparatively low mosquito survival of the Martens I scheme (see Figure 3 in ). In addition, the start of the malaria season is notably delayed in the LMM2004 simulations. The start of the season occurs about one to two months later under the original than under the new LMM version . The strong growth of the mosquito population causes shorter main transmission seasons (MSeas values) in the LMM2004 runs than in that of the LMM2010. Moreover, a later occurrence of XSeas is found for the LMM2004 simulations.
Deficiencies of the LMM2004 are also found for the three parasitological variables. Almost the whole population clears the malaria parasite during the dry season due to the high recovery rate (r) of the LMM2004 (LMM2004 (LMM2010): r = 0.0284(0.005) day-1; see ). Such a characteristic is, however, not found in parasitological surveys (see Figure 5) and results in a very low skill in terms of PR min, a (see Table 4). The strong recovery from infection, moreover, underestimates the values of PR a and PR max, a . In comparison to LMM2010 runs, parasite ratios are hence much lower in the LMM2004 simulations (see Figure 6e & 6f).
One characteristic of malaria epidemic areas is the sudden and unexpected increase of the parasite ratio in certain years. Epidemic regions therefore should reveal a strong interannual variability of the annual maximum of the asexual parasite ratio (PR max, a ; Figure 6g & 6h) and in general low annual mean asexual parasite ratios (PR a ; Figure 6e & 6f). According to LMM simulations forced with gridded REMO data, epidemic areas are found along a strip within the Sahelian zone as well as for highland territories and fairly dry areas in East Africa, Zambia, and Angola (Figure 6g & 6h). Most notable is the fact that in the LMM2004 runs the strip of high epidemic risk is detected about 2-4° to the south in comparison to that of the LMM2010. Note that the LMM2010 simulations seem to be more realistic in comparison to previous assessments of epidemic territories in the Sahel . Differences between LMM2004 and LMM2010 simulations are also found outside of the Sahelian zone. Overall the LMM2010 seems to result in a stronger interannual variability of PR max, a .
The aim of the present study was the development of an improved weather-driven malaria model, which is able to simulate malaria transmission in both epidemic as well as endemic malaria areas. This section provides a detailed discussion with regard to various aspects of the present study. The present-day climate performance of the new model version are discussed relative to results of former studies and the model calibration is evaluated.
Comparison of LMM2010 runs with those performed by the original formulation of the model reveals a significant improvement (see Table 4) of the model performance both for epidemic as well as endemic malaria areas. In contrast to the LMM2004, the LMM2010 is able to reproduce the low transmission rates in the northern part of the Sahel. This enables, therefore, an improved detection of epidemic areas, in particular, in the Sahel. The LMM2010 validation also shows that the model can now also be applied for endemic malaria areas. The usage of the fuzzy distribution model enables the simulation of realistic sizes of the mosquito population under humid rainfall conditions resulting in reasonable transmission rates. Moreover, the lag in the malaria seasonality has disappeared in the new version of the LMM.
Various published malaria distribution maps [43, 48] correspond well with the simulated spread of malaria by the LMM2010. However, in certain parts the simulated intensity of malaria transmission differs considerably from published malaria maps. The LMM2010 in general seems to predict higher transmission rates than satellite-derived predictions of EIR a from Rogers et al. . The maps of transmission intensity provided by Gemperli et al.  are fairly patchy. In fact, the prediction from Gemperli et al. significantly suffers from the neglected interannual variability of malaria. Based on the few available EIR a observations it is difficult to judge which estimates are closer to reality. However, the validation of the LMM2010 under different climatic conditions provides evidence that the present study generated realistic biting rates and a reasonable interannual variability.
The calibration of the LMM was performed in West Africa for different atmospheric conditions of epidemic and endemic malaria regions. Realistic temperature and precipitation time series were reconstructed from various synoptic weather stations. The comparison with observations from eleven entomological and parasitological variables finally enabled the setting of undetermined model parameters.
The databases (including meteorological, entomological, and parasitological observations) for the LMM calibration are not optimal. There is a mismatch between the scales at which a disease vector responds to hydrologic variability and the scales at which hydrologic variability is actually observed. Systems should be developed that monitor hydrologic variability at scales corresponding to disease system ecologies . In this study, the generation of realistic meteorological station time series enabled the comparison with atmospheric conditions from malaria field studies, which were not conducted directly at the weather stations. These sites therefore in any case exhibit a different temporal variability of rainfall and temperatures. This might be one reason, amongst other factors such as environmental conditions, why year-to-year comparisons between observation and simulation were weakly correlated at single locations. In order to circumvent this problem, the present study refrained from looking at paired annual correlations at single stations but applied a problem-adapted scoring system.
The required historical entomological and parasitological data are rarely available with sufficient coverage. Most locations show only one, two, or even no field measurements. It is therefore likely that a larger set of observations would have an impact on the result of the model calibration. Ideally, model simulations and malaria observations should be compared from year-to-year. However, this would require the simultaneous monitoring of long-term malaria data and meteorological measurements. Such long time series are available for the area of Ndiop/Senegal (S. Louvet, personal communication, 2007), but these data sets are not publicly available.
The close ranking of diverse model runs as well as the lack of sufficient validation data further restricted an objective formal fitting of the model. In fact, various steps of the calibration procedure were subjective. Due to high computational costs it was furthermore not possible to fit all remaining model parameters simultaneously. However, because various settings compensate each other it is likely that the final model formulation conforms as much as possible to reality.
The calibration and validation of the model should also be ideally not only restricted to West Africa. However, such an extension to, for example, East Africa was beyond the scope of this study. The Malaria Atlas Project intends to provide access to various malaria studies . This might provide an efficient access to malaria data beyond that of West Africa. Such an extension would ideally include East African highlands and an estimation of the sporogonic temperature threshold.
This study was naturally not able to account for all processes involved in the spread of malaria. Some factors might be included in a future extension of the LMM. The simulation of the parasitological malaria variables by the LMM2010 is a simplification of real processes. The validation of the LMM2010 by means of parasitological measurements in West Africa revealed shortcomings of the new model version. Lower skill scores were achieved by the three parasitological variables when compared to the results from the eight entomological variables.
In addition to the lack of immunity, the LMM2010 does not account for other malaria factors such as chemoprophylaxis and human activities. However, this could be implemented by means of a variable parameter setting. Observations suggest a greater variability of the parasite ratio. At Bobo-Dioulasso, for example, the ten observed annual mean asexual parasite ratios (PR a ) range between 29.1 and 77.5%. In contrast, the 34 annual values of the LMM2010 only span values between 50 and 70% .
Due to the lack of long-term observations, Kleinschmidt et al.  and Gemperli et al.  were forced to neglect the interannual variability of PR. This fact might again partly be responsible for their projected irregular PR pattern in West Africa. Their maps also show a sharp decrease of PR a north of about 15°N. In contrast to the LMM2010 runs and the Garki model simulations from Gemperli et al. , PR a is frequently lower than 50% south of 15°N. Only few regions exhibit higher PR a values than 70%, which are simulated by the LMM2010.
It should be pointed out here that climate is rarely the only important driver of malaria. Numerous other studies showed [21, 22, 53, 54] that in particular human activities are crucial for the transmission and prevention of malaria across Africa. For example, the modification of the landscape by irrigation , forest clearing , or urbanization  can significantly alter malaria transmission. The present study assessed only the malaria risk for rural areas without the influence of permanent breeding places. The applicability of this analysis is therefore limited when permanent water bodies or urban centres are present. In principle, the calibration of the LMM2010 could also be performed for urban areas. However, such an undertaking seems to be hampered by the lower number of available observations (see Additional file 2).
One of the most comprehensive studies to date in terms of gathering validation data and information from the literature for the development of a new version of an existing malaria model was conducted for West Africa. The first part of this study  provided new parameter settings of the LMM and changed some key processes in the model. The performance of numerous model versions were compared to malaria field observations from rural sites with no malaria measures and irrigation in order to maximise the climate-driven malaria impact. The comparison with observations from eleven entomological and parasitological variables finally enabled the specification of a final set of parameter settings of the model. Validation of the new LMM version in West Africa reveals that the simulations and malaria seasonality compare well with entomological field observations of epidemic and endemic malaria areas. The LMM2010 also demonstrates a fairly realistic simulation of the malaria spread as well as an improved detection of epidemic risk in Africa. Due to model limitations, the performance of the LMM2010 is somewhat weaker with regard to parasitological variables.
It is concluded that the LMM2010 represents a significant step forward in the modelling of a weather-driven malaria transmission cycle. In contrast to the original model, the application of the new model version is not only restricted to epidemic malaria regions but the usability is now extended to endemic malaria areas. Ermert  hence used the LMM2010 for the assessment of malaria risk under the influence of observed and projected climate change using regionalized climate projections for Africa. The LMM2010 simulated transmission rates were passed to the Garki model to form a hybrid malaria model, which is able to consider further aspects of the malaria disease such as age-dependent parasite ratios as well as the immune status of the population. In the near future, this hybrid model will be applied in a health early warning system, which is implemented by the QWeCI (Quantifiying Weather and Climate Impacts on Health in Developing Countries) project from the European Community's Seventh Framework Research Programme. Further work might enable the inclusion of some other processes involved in the spread of malaria into the LMM.
List of abbreviations
DWD: German Weather Service; ECHAM5/MPI-OM: European Centre HAmburg Model, 5th generation/Max-Planck-Institute-Ocean Model; GSOD: Federal climate complex Global Surface Summary of Day version 7; LMM: Liverpool Malaria Model; LMM2004: Liverpool Malaria Model version of 2004; LMM2010: Liverpool Malaria Model version of 2010; MARA: Mapping Malaria Risk in Africa;
List of symbols
#E p : number of produced eggs per female mosquito; CAP: cap on the number of fertile mosquitoes; CSPR: CircumSporozoite Protein Rate; CSPR a : annual mean CircumSporozoite Protein Rate; EIR: Entomological Inoculation Rate; EIR m : monthly Entomological Inoculation Rate; EIR a : annual Entomological Inoculation Rate; ESeas: End month of the malaria Season; HBR: human Biting Rate; HBR a : annual Human Biting Rate; MSeas: length of the Main malaria Season; pd↓: dry season mosquito survival probability shift off; PR: asexual Parasite Ratio; PR a : annual mean asexual Parasite Ratio; PR max, a : annual maximum of the asexual Parasite Ratio; PR min, a : annual minimum of the asexual Parasite Ratio; r: recovery rate; R: rainfall; R a : annual rainfall amount; RΣ10d: 10-day accumulated precipitation; S: most suitable rainfall condition according to the fuzzy distribution model; SC(all): skill score with regard to the used eleven malaria variables; SC(x): skill score with regard to malaria variable x; Seas: length of the malaria Season; SSeas: Start month of the malaria Season; T: temperature; T a : annual mean temperature; U1: lower threshold of unsuitable rainfall conditions with regard to the fuzzy distribution model; U2: upper threshold of unsuitable rainfall conditions regarding the fuzzy distribution model; XSeas: month of maXimum transmission, i.e. the month with the largest EIR value
We are grateful to numerous people who kindly provided useful data, literature references, and special knowledge for this study. These are: Dipl.-Met. C. Baugirdis from DWD/Offenbach (weather station data for West Africa), Dr. S. Louvet from the 'Centre de Recherches de Climatologie (CRC)'/Dijon (malaria literature), as well as S. Stephens and N. Lott from US National Climatic Data Center (interpretation of GSOD precipitation).
This study was part of the IMPETUS West Africa project and was supported by the Federal German Ministry of Education and Research (BMBF) under grant No. 01 LW06001A in the GLOWA programme and by the Ministry of Innovation, Science, Research, and Technology (MIWFT) of the federal state of North Rhine-Westphalia under grant No. 313-21200200. The support by the EU project AMMA (funded by the European Commission's Sixth Framework Research Programme through project No. 004089) and the Natural Environment Research Council (NERC) AMMA is gratefully acknowledged by AEJ and APM. Based on a French initiative, AMMA was built by an international scientific group and was funded by a large number of agencies, especially from France, United Kingdom, United States, and Africa. Detailed information on scientific coordination and funding is available on the AMMA international website http://www.amma-international.org. Finally, we wish to thank one anonymous reviewer whose comments helped to greatly improve the manuscript.
- WHO: World malaria situation in 1994. Wkly Epidemiol Rec. 1997, 72: 269-276.
- Greenwood B, Mutabingwa T: Malaria in 2002. Nature. 2002, 415: 670-672. 10.1038/415670a.View ArticlePubMed
- Rogers DJ, Randolph SE: The global spread of malaria in a future, warmer world. Science. 2000, 289: 1763-1766. 10.1126/science.289.5478.391b. (Corrected in print, 29 Sept. 2000, p. 2283)View ArticlePubMed
- Confalonieri U, Menne B, Akhtar R, Ebi KL, Hauengue M, Kovats RS, Revich B, Woodward A: Human health. Climate change 2007: impacts, adaptation and vulnerability. Edited by: Parry ML, Canziani OF, Palutikof JP, van der Linden PJ, Hanson CE. 2007, Cambridge, United Kingdom: Cambridge University Press, 391-431. [Contribution of working group II to the fourth assessment report of the Intergovernmental Panel on Climate Change]
- 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.View Article
- Patz JA, Campbell-Lendrum D, Holloway T, Foley JA: Impact of regional climate change on human health. Nature. 2005, 438: 310-317. 10.1038/nature04188.View ArticlePubMed
- Tompkins AM, Feudale L: Seasonal ensemble predictions ofWest African monsoon precipitation in the ECMWF system 3 with a focus on the AMMA special observing period in 2006. Weather and Forecasting. 2010, 25: 768-788. 10.1175/2009WAF2222236.1.View Article
- Solomon S, Qin D, Manning M, Chen Z, Marquis M, Averyt KB, Tignor M, Miller HL: Climate change 2007: the physical science basis. 2007, Cambridge, United Kingdom and New York, USA: Cambridge University Press, [Contribution of working group I to the fourth assessment report of the Intergovernmental Panel on Climate Change]
- Hoshen MB, Morse AP: A weather-driven model of malaria transmission. Malar J. 2004, 3: 32-10.1186/1475-2875-3-32.PubMed CentralView ArticlePubMed
- Ermert V, Fink AH, Jones AE, Morse AP: A new version of the Liverpool Malaria Model. I. Review of the parameter setting and model structures. Malar J. 2011, 10 (1): 35-10.1186/1475-2875-10-35.PubMed CentralView ArticlePubMed
- 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, University of Cologne, Cologne, Germany, [http://kups.ub.uni-koeln.de/volltexte/2010/3109/]
- Faye O, Fontenille D, Herve JP, Diack PA, Mouchet SDJ: Le paludisme en zone sahelienne du Sénégal. 1. Données entomologiques sur la transmission. Ann Soc Belg Med Trop. 1993, 73: 21-30.PubMed
- Quakyi IA, Leke RGF, Befidi-Mengue R, Tsafack M, Bomba-Nkolo D, Manga L, Tchinda V, Njeungue E, Kouontchou S, Fogako J, Nyonglema P, Harun LT, Djokam R, Sama G, Eno A, Megnekou R, Metenou S, Ndountse L, Same-Ekobo A, Alake G, Meli J, Ngu J, Tietche F, Lohoue J, Mvondo JL, Wansi E, Leke R, Folefack A, Bigoga J, Bomba-Nkolo C, Titanji V, Walker-Abbey A, Hickey MA, Johnson AH, Taylor DW: The epidemiology of Plasmodium falciparum malaria in two Cameroonian villages: Simbok and Etoa. Am J Trop Med Hyg. 2000, 63: 222-230.PubMed
- Detinova TS: Age-grouping methods in Diptera of medical importance with special reference to some vectors of malaria. No. 47 in Monograph Series, WHO. 1962
- Paeth H, Born K, Girmes R, Podzun R, Jacob D: Regional climate change in Tropical and Northern Africa due to greenhouse forcing and land use changes. J Clim. 2009, 22: 114-132. 10.1175/2008JCLI2390.1.View Article
- Jacob D, van den Hurk BJJM, Andrae U, Elgered G, Fortelius C, Graham LP, Jackson SD, Karstens U, Koepken C, Lindau R, Podzun R, Rockel B, Rubel F, Sass BH, Smith R, Yang X: A comprehensive model intercomparison study investigating the water budget during the PIDCAP period. Meteorology and Atmospheric Physics. 2001, 77: 19-44. 10.1007/s007030170015.View Article
- Roeckner E: IPCC-AR4 MPI-ECHAM5 T63L31 MPI-OM GR1.5L40 20C3 M run no.2: atmosphere 6 HOUR values MPImet/MaD Germany. 2004, World data center for climate. CERA-DB "EH5-T63L31 OM-GR1.5L40 20C 2 6H", [http://cera-www.dkrz.de/WDCC/ui/Compact.jsp?acronym=EH5-T63L31_OM-GR1.5L40_20C_2_6H]
- Roeckner E, Lautenschlager M, Schneider H: IPCC-AR4 MPI-ECHAM5 T63L31 MPI-OM GR1.5L40 20C3 M run no.1: atmosphere 6 HOUR values MPImet/MaD Germany. 2006, World data center for climate
- Roeckner E, Lautenschlager M, Schneider H: IPCC-AR4 MPI-ECHAM5 T63L31 MPI-OM GR1.5L40 20C3 M run no.3: atmosphere 6 HOUR values MPImet/MaD Germany. 2006, World data center for climate
- Speth P, Christoph M, Diekkrüger B, Bollig M, Fink AH, Goldbach H, Heckelei T, Menz G, Reichert B, Rössler M: Impacts of global change on the hydrological cycle in West and Northwest Africa. 2010, Heidelberg, Germany: Springer, [http://www.springer.com/978-3-642-12956-8]View Article
- 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 CentralView ArticlePubMed
- Kelly-Hope LA, 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 CentralView ArticlePubMed
- Carnevale P, Guillet P, Robert V, Fontenille D, Doannio J, Coosemans M, Mouchet J: Diversity of malaria in rice growing areas of the Afrotropical region. Parassitologia. 1999, 41: 273-276.PubMed
- Keiser J, Castro MCD, Maltese MF, Bos R, Tanner M, Singer BH, Utzinger J: Effect of irrigation and large dams on the burden of malaria on a global and regional scale. Am J Trop Med Hyg. 2005, 74: 392-406.
- Konaté L, Diop A, Sy N, Faye MN, Deng Y, Izri A, Faye O, Mouchet J: Comeback of Anopheles funestus in Sahelian Senegal. Lancet. 2001, 358: 336-View ArticlePubMed
- Ijumba JN, Shenton FC, Clarke SE, Mosha FW, Lindsay SW: Irrigated crop production is associated with less malaria than traditional agricultural practices in Tanzania. Trans R Soc Trop Med Hyg. 2002, 96: 476-480. 10.1016/S0035-9203(02)90408-6.View ArticlePubMed
- Lindsay SW, Alonso PL, Schellenberg JRMA, Hemingway J, Thomas PJ, Shenton FC, Greenwood BM: A malaria control trial using insecticide-treated bed nets and targeted chemoprophylaxis in a rural area of The Gambia, West Africa: 3. entomological characteristics of the study area. Trans R Soc Trop Med Hyg. 1993, 87 (Suppl. 2): 19-23. 10.1016/0035-9203(93)90171-L.View ArticlePubMed
- Lindsay SW, Alonso PL, Schellenberg JRMA, Hemingway J, Adiamah JH, Shenton FC, Jawara M, Greenwood BM: A malaria control trial using insecticide-treated bed nets and targeted chemoprophylaxis in a rural area of The Gambia, West Africa: 7. impact of permethrin-impregnated bed nets on malaria vectors. Trans R Soc Trop Med Hyg. 1993, 87 (Suppl. 2): 45-51. 10.1016/0035-9203(93)90175-P.View ArticlePubMed
- Akogbéto M, Nahum A: Impact des moustiquaires imprégnées de deltaméthrine sur la transmission du paludisme dans un milieu côtier lagunaire, Bènin. Bull Soc Pathol Exot. 1996, 89: 291-298.PubMed
- Cuzin-Ouattara N, den Broek AHV, Habluetzel A, Diabate A, Sanogo-Ilboudo E, Diallo DA, Cousens SN, Esposito F: Wide-scale installation of insecticide-treated curtains confers high levels of protection against malaria transmission in a hyperendemic area of Burkina Faso. Trans R Soc Trop Med Hyg. 1999, 93: 473-479. 10.1016/S0035-9203(99)90343-7.View ArticlePubMed
- Le Goff G, Robert V, Fondjo E, Carnevale P: Efficacy of insecticide impregnated bed-nets to control malaria in a rural forested area in southern Cameroon. Mem Inst Oswaldo Cruz. 1992, 87 (Suppl. 3): 355-359.View ArticlePubMed
- Le Goff G, Carnevale P, Fondjo E, Robert V: Comparison of three sampling methods of man-biting anophelines in order to estimate the malaria transmission in a village of south Cameroon. Parasite. 1997, 4: 75-80.View ArticlePubMed
- Awolola TS, Okwa O, Hunt RH, Ogunrinade AF, Coetzee M: Dynamics of the malaria-vector populations in coastal Lagos, south-western Nigeria. Ann Trop Med Parasitol. 2002, 96: 75-82. 10.1179/000349802125000538.View ArticlePubMed
- Fontenille D, Lochouarn L, Diatta M, Sokhna C, Dia I, Diagne N, Lemasson JJ, Ba K, Tall A, Rogier C, Trape JF: Four years' entomological study of the transmission of seasonal malaria in Senegal and the bionomics of Anopheles gambiae and A. arabiensis. Trans R Soc Trop Med Hyg. 1997, 91: 647-652. 10.1016/S0035-9203(97)90506-X.View ArticlePubMed
- Smith DL, Guerra CA, Snow RW, Hay SI: Standardizing estimates of the Plasmodium falciparum parasite rate. Malar J. 2007, 6: 131-10.1186/1475-2875-6-131.PubMed CentralView ArticlePubMed
- Hay SI, Noor AM, Simba M, Busolo M, Guyatt HL, Ochola SA, Snow RW: Clinical epidemiology of malaria in the highlands of Western Kenya. Emerg Infect Dis. 2002, 8: 543-548. 10.3201/eid0804.010529.PubMed CentralView ArticlePubMed
- Hay SI, Rogers DJ, Toomer JF, Snow RW: Annual Plasmodium falciparum entomological inoculation rates (EIR) across Africa: literature survey, internet access and review. Trans R Soc Trop Med Hyg. 2000, 94: 113-127. 10.1016/S0035-9203(00)90246-3.PubMed CentralView ArticlePubMed
- Chalvet-Monfray K, Sabatier P, Bicout DJ: Downscaling modeling of the aggressiveness of mosquitoes vectors of diseases. Ecol Model. 2007, 204: 540-546. 10.1016/j.ecolmodel.2007.01.024.View Article
- Shaman J, Day JF: Achieving operational hydrologic monitoring of mosquitoborne disease. Emerg Infect Dis. 2005, 11: 1343-1350.PubMed CentralView ArticlePubMed
- Balme M, Vischel T, Lebel T, Peugeot C, Galle S: Assessing the water balance in the Sahel: impact of small scale rainfall variability on runoff. Part 1: Rainfall variability analysis. Journal of Hydrology. 2006, 331: 336-348. 10.1016/j.jhydrol.2006.05.020.View Article
- Murphy AH: Skill scores based on the mean square error and their relationships to the correlation coefficient. Monthly Weather Review. 1988, 25: 2417-2424. 10.1175/1520-0493(1988)116<2417:SSBOTM>2.0.CO;2.View Article
- Ahumada JA, Lapointe D, Samuel MD: Modeling the population dynamics of Culex quinquefasciatus (Diptera: Culicidae), along an Elevational Gradient in Hawaii. J Med Entomol. 2004, 41: 1157-70. 10.1603/0022-2585-41.6.1157.View ArticlePubMed
- MARA: Towards an atlas of malaria risk in Africa: first technical report of the MARA/ARMA collaboration. 1998, First technical report, MARA/ARMA, Durban, South Africa
- Le Masson JJ, Fontenille D, Lochouam L, Dia I, Simard F, Ba K, Diop A, Diatta M, Molez JF: Comparison of behavior and vector efficiency of Anopheles gambiae and An. arabiensis ( Diptera: Culicidae) in Barkedji, a Sahelian area of Senegal. J Med Entomol. 1997, 34: 396-403.View Article
- Molez JF, Diop A, Gaye O, Lemasson JJ, Fontenille D: Malaria morbidity in Barkedji, village of Ferlo, in Senegal Sahelian area. Bull Soc Pathol Exot. 2006, 99: 187-190.PubMed
- Stafford Smith DM: Mosquito records from the Republic of Niger, with reference to the construction of the new 'Trans-Sahara Highway'. J Trop Med Hyg. 1981, 84: 95-100.
- Snow RW, Craig M, Deichmann U, Marsh K: Estimating mortality, morbidity and disability due to malaria among Africa's non-pregnant population. Bull World Health Org. 1999, 77: 624-640.PubMed CentralPubMed
- Rogers DJ, Randolph SE, Snow RW, Hay SI: Satellite imagery in the study and forecast of malaria. Nature. 2002, 415: 710-715. 10.1038/415710a.PubMed CentralView ArticlePubMed
- Gemperli A, Sogoba N, Fondjo E, Mabaso M, Bagayoko M, Brit OJT, Anderegg D, Liebe J, Smith T, Vounatsou P: Mapping malaria transmission inWest and Central Africa. Trop Med Int Health. 2006, 11: 1032-1046. 10.1111/j.1365-3156.2006.01640.x.View ArticlePubMed
- Shaman J, Day JF: Reproductive phase locking of mosquito populations in response to rainfall frequency. PLoS One. 2007, 3: e331-10.1371/journal.pone.0000331.View Article
- Guerra CA, Hay SI, Lucioparades LS, Gikandi PW, Tatem AJ, Noor AM, Snow RW: Assembling a global database of malaria parasite prevalence for the Malaria Atlas Project. Malar J. 2007, 6: 17-10.1186/1475-2875-6-17.PubMed CentralView ArticlePubMed
- Kleinschmidt I, Omumbo J, Briët O, van de Giesen N, Sogoba N, Mensah NK, Windmeijer P, Moussa M, Teuscher T: An empirical malaria distribution map for West Africa. Trop Med Int Health. 2001, 6: 779-786. 10.1046/j.1365-3156.2001.00790.x.View ArticlePubMed
- Epstein PR: Global warming and vector-borne disease. Lancet. 1998, 351: 1737-10.1016/S0140-6736(05)77777-1.View ArticlePubMed
- van Lieshout M, Kovats RS, Livermore MTJ, Martens P: Climate change and malaria: analysis of the SRES climate and socio-economic scenarios. Glob Environ Change. 2004, 14: 87-99. 10.1016/j.gloenvcha.2003.10.009.View Article
- Briët OJT: A simple method for calculating mosquito mortality rates, correcting for seasonal variations in recruitment. Med Vet Entomol. 2002, 16: 22-27.View ArticlePubMed
- Munga S, Minakawa N, Zhou G, Mushinzimana E, Barrack OOJ, Githeko AK, Yan G: Association between land cover and habitat productivity of malaria vectors in western Kenyan highlands. Am J Trop Med Hyg. 2006, 74: 69-75.PubMed
- Keiser J, Utzinger J, de Castro MC, Smith TA, Tanner M, Singer BH: Urbanization in sub-Saharan Africa and malaria control. Am J Trop Med Hyg. 2004, 71 (Suppl. 2): 118-127.PubMed
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.