Effect of meteorological factors on clinical malaria risk among children: an assessment using village-based meteorological stations and community-based parasitological survey

Background Temperature, rainfall and humidity have been widely associated with the dynamics of malaria vector population and, therefore, with spread of the disease. However, at the local scale, there is a lack of a systematic quantification of the effect of these factors on malaria transmission. Further, most attempts to quantify this effect are based on proxy meteorological data acquired from satellites or interpolated from a different scale. This has led to controversies about the contribution of climate change to malaria transmission risk among others. Our study addresses the original question of relating meteorological factors measured at the local scale with malaria infection, using data collected at the same time and scale. Methods 676 children (6–59 months) were selected randomly from three ecologically different sites (urban and rural). During weekly home visits between December 1, 2003, and November 30, 2004, fieldworkers tested children with fever for clinical malaria. They also collected data on possible confounders monthly. Digital meteorological stations measured ambient temperature, humidity, and rainfall in each site. Logistic regression was used to estimate the risk of clinical malaria given the previous month's meteorological conditions. Results The overall incidence of clinical malaria over the study period was 1.07 episodes per child. Meteorological factors were associated with clinical malaria with mean temperature having the largest effect. Conclusion Temperature was the best predictor for clinical malaria among children under five. A systematic measurement of local temperature through ground stations and integration of such data in the routine health information system could support assessment of malaria transmission risk at the district level for well-targeted control efforts.


Background
Meteorological factors are important drivers of malaria transmission. Temperature, rainfall and humidity have been associated with the dynamics of malaria vector population and, therefore; with spread of the disease. Ambient temperature plays a major role in the life cycle of the malaria vector. The development of the parasite within the mosquito (sporogonic cycle) is dependent on temperature. The sporogonic cycle takes about 9 to 10 days at temperatures of 28°C, but stops at temperatures below 16°C [1][2][3]. The daily survival of the vector is dependent on temperature as well. At temperatures between 16°C and 36°C, the daily survival is about 90%. This survival drops rapidly at temperature above 36°C. The highest proportion of vectors surviving the incubation period is observed at temperatures between 28° and 32°C [4]. The gonotrophic cycle, which is the time between two blood meals of the vector, is short at higher temperatures because the digestion speed increases [5]. Therefore, higher temperatures result in more frequent vector-host contact.
Several field studies have reported the impact of ambient temperature on malaria outcomes. [6][7][8][9][10]. In South Africa, Craig and colleagues [8] identified a significant correlation between temperature and the number of malaria cases. Mean maximum daily temperatures from January to October of the preceding season were positively associated with the incidence of clinical cases of malaria. In Ethiopia, minimum temperature was associated with malaria in a cold district (minimum temperature below 12°C); while in a hot district (minimum temperature above 12°C) the effect was not significant [9].
Rainfall provides breeding sites for mosquitoes to lay their eggs, and ensures a suitable relative humidity of at least 50 to 60% to prolong mosquito survival. Relative humidity below 60% shortens the life span of the mosquitoes. The onset of the rainy season was associated with an increase in vector abundance [11,12]. The 1958 malaria epidemic in Ethiopia was associated with unusual high amounts of rain. Similarly, in Nairobi, outbreaks of malaria occurred in 1940 after heavy rains [13]. In the Ugandan highlands, rainfall anomalies (differences from the mean) because of El Nino were positively correlated with vector density one month later, and this may have started the resulting epidemic [14].

Study site location
Despite consistently reported evidence that meteorological factors affect malaria transmission, conflicting opinions persist about their relative importance. Shanks and colleagues [15] could not find a link between ambient temperature and malaria admissions over a 30-year period in Kenya. Similarly, Hay and colleagues [16] attributed the increase in malaria morbidity at four sites in East Africa to drug resistance, rather than changes in temperature. Patz argued that this discrepancy might arise from methodological issues. In these studies, climate data were obtained by interpolating broad-scale grid regional climate data, based on sparse historical meteorological station data, which are unsuitable for individual village sites [17]. This was described as an absence of evidence, rather than evidence of absence of a climate effect. To assess climate factors and infectious disease outcomes, Kovats and colleagues recommended that enough measurements on each aspect should be performed for a specific study region [18].
The purpose of the current study was to estimate the effects of temperature, rainfall, and relative humidity on the incidence of clinical malaria, by examining data on these meteorological factors and laboratory-confirmed cases of clinical malaria. A binary-response logistic regression model with fractional polynomial (FP) transformations was used to assess these effects. The FP approach allowed us to explore better the nature of the relationships between these meteorological parameters and the risk of clinical malaria. This is because presumed linearity of these relationships is not justified [1][2][3].

Study site
The study was conducted in the town of Nouna and the villages of Cissé and Goni. These three sites are part of the Nouna Demographic Surveillance Systems (DSS) area [19], which is located in Kossi province in the north-western part of Burkina Faso ( Figure 1).

Study population
Children were selected at each site by cluster sampling of households, using a sampling frame produced from the DSS database. The sample size was set to detect inter-site differences in clinical malaria incidence of at least 10%, with 80% power and 95% confidence, and accounting for 15% loss to follow up. Children were involved after getting informed consent from their parents. In total, 676 children (Cissé: 171, Goni: 240 and Nouna: 265), aged 6 to 59 months, took part in the study. The Nouna ethical committee approved the study.

Clinical malaria detection
Three trained interviewers, one based at each site, visited each child at home every week. At each visit, they measured the child's axillary temperature and collected a blood film (by finger prick) from any child who was febrile (axillary temperature > = 37.5 C). In addition, interviewers collected information on bed net use and housing conditions. Blood films were read in the Laboratory of the Nouna Health Research Centre. Details on the clinical malaria detection procedure are provided elsewhere [20]. The outcome measure was clinical malaria episode, defined as being febrile with parasiteamia.
For purposes of this study, incidence was defined as the number of clinical malaria episodes detected per 1000 units of person time bearing in mind that some illness episodes which occurred between visits may have been missed.
For ethical reasons, children with fever were treated presumptively with cholorquine (CQ) according to the then national treatment guidelines for malaria. When fever persisted for two days or other symptoms surfaced, the interviewers immediately referred the child to the nearest health centre. The project covered all related costs.

Meteorological data
Rainfall, temperature, and humidity were measured on the ground using meteorological units installed at each of the three sites. Each meteorological unit consisted of a Digital Datalogger (THIES Datalogger, MeteoLOG TDL 14), to which three sensors (temperature, humidity, and Meteorological station (a) and datalogger (b) in one of the study sites rainfall) were connected ( Figure 2). Following standard meteorological conventions, the temperature and humidity sensors were set at a height of 2.5 metres and the rainfall sensor 1.5 metres above the ground. The Dataloggers were set for 10-second measurement cycles and 10minute recording cycles. Every month, a meteorological supervisor visited each site, downloaded the data from the Datalogger into a memory card, and transferred them into a meteorological database. Daily minimum, mean, and maximum temperatures, relative humidity, and total rainfall were then calculated.

Statistical Analysis
A conventional binary outcome logistic regression model was used to assess the effects of mean, minimum, and maximum temperatures (T°C), relative humidity (RH), and amount of rainfall (Pmm) of the previous month on clinical malaria rates among study participants. The interaction terms (T°C*RH', 'T°C* Pmm' and 'RH*Pmm') were included in the model to control for the high correlation of the three individual variables. Other covariates included in the model were site of the study; sex; age; use of bed net; type of house; presence of a well; presence of a farm (presence of any farming; which included cereal and vegetable farms); presence of animals; and presence of a mosquito breeding site within a 30-metre radius of a participant's house. As the temperature, rainfall and relative humidity are continuous variables and their relationships with clinical malaria might not be linear, multivariate FP procedures [21,22] were used to determine the best-fitting relationship they had with clinical malaria. The fitting procedure involves transforming a continuous variable, using a class of eight possible functions to identify the one that provides the best fit. These functions are H 1 (X) = X p , where p takes eight possible values: -2, -1, -0.5, 0, 0.5, 1, 2, 3. The linear model is represented by X 1 , and X 0 repre-sents the logarithm of X. The transformation can be either first or second-degree [21,22].
For the simplest model, the one with one continuous variable (e.g. temperature) and one binary variable (e.g. bed net use, Yes or No), the logistic regression model using the first-degree FP is ... Where π ib is the predicted probability of a child testing positive for clinical malaria at a temperature i, and for bed net use b (yes or no); H 1 (X) is the functional form to which the co-variable X is transformed; β 1 its coefficient; and γ is the coefficient of the co-variable; W b , the use of a bed net. The rate ratio for the first-degree PF transformation is given by the formula: Second-degree transformation uses a combination of two powers from the list of eight. In total, 36 combinations are possible. This is calculated as follows: , where n is the number of possible powers (8), and k the number of powers in each combination (2); then . In addition, the eight combinations with the same power are added. The eight powers and the 36 combinations are tested consecutively. For the first-degree FP, the differences in deviance of each model from the linear one are calculated and compared with the chi-square distribution, with one degree of freedom at α = 0.05. For the second-degree FP, the differences in deviance of each model from the best fitting first-degree FP are calculated and compared with the chi-square distribution with two degrees of freedom. The model with the largest significant deviance difference compared to the best firstdegree FP is selected as the best second-degree FP. The second-degree FP is therefore implicitly a better fit than either the linear model or the first-degree FP. Using the same example as in the first-degree FP; the second-degree FP model is defined mathematically as follows: .where, H 1 (X) and H 2 (X) are the respective functions to which co-variable X is transformed.
The rate ratio for the second-degree FP is given by the formula: All the models were run and fit using STATA ® software [23]. The output of each model (coefficients and transformed variables) was used to calculate rate ratios (RR), only including variables with a significant effect on the risk of clinical malaria. These are mean temperature, rainfall, and relative humidity. Other variables, including minimum and maximum temperature, were removed from the final model by an in-built backward elimination procedure in the FP algorithm. For each variable, reference points (x 0 ) were defined. These are mean temperature = 27°C, the value at which clinical malaria risk peaked; rainfall = 164 mm, the highest value observed; and rela-tive humidity = 60%, the minimum value required for malaria vector survival.

Study population
Characteristics Table 1 provides a summary of the study population characteristics at the beginning of the observation period. In total, 676 children distributed in 350 households took part in the study. The average number of children per household was 1.9. Overall, females were more (52.8%) than males, but the sex distribution was not significantly Monthly rainfall, mean temperature and relative humidity for Nouna Figure 5 Monthly rainfall, mean temperature and relative humidity for Nouna. Monthly rainfall, mean temperature and relative humidity for Cissé Figure 3 Monthly rainfall, mean temperature and relative humidity for Cissé. Monthly rainfall, mean temperature and relative humidity for Goni Figure 4 Monthly rainfall, mean temperature and relative humidity for Goni.

Rainfall
Relative humidity Mean Temperature different across the sites (p = 0.420). Similarly, the age distribution did not differ between sites (p = 0.938).

Follow up status
During the one-year follow up, 20 children (3.0% of 676) left the cohort, either because of death (11) or migration out of the study sites (9). Although 52 home visits were planned for each child, on average, each child was seen for 45.4 times, because children were not always present at each visit. In total 594.9 person-years (PYs) were observed.

Meteorological conditions
The three sites presented a similar pattern of meteorological conditions (Figure 3, 4, 5). Rainfall was observed in all sites from May to November; the highest amounts were observed in August, September, and October. The total amount of rainfall was higher in Nouna (508. 3 mm in 68 days) than in Cissé (334 mm in 49 days) or Goni (408.5 mm in 54 days). The pattern of relative humidity followed the one of rainfall. At all sites, relative humidity was low during the first six months of the year, and then increased with the onset of the rains.    Table 3 shows the transformation functions for each variable. The relationships between meteorological parame-ters and clinical malaria incidence were not linear since none of the functions had the power of +1. Minimum and maximum temperatures, which were included in the initial model, were removed by backward elimination. Only the interaction terms of 'temperature and humidity', and 'temperature and rainfall' had a linear relationship with first-degree transformation. The best fit for mean temperature, rainfall and relative humidity alone were obtained with second-degree transformations. For the interaction term 'rainfall and humidity', the first-degree transformation with a power of -1 yielded the best fit.

Model estimates of coefficients and confidence limits
Although other variables were included in the model for adjustment purposes, their estimates are not shown in Table 4, because the focus was on meteorological parameters. For each power of transformation, an estimate is given. All variables, individually or combined, had a significant effect on clinical malaria rates. The highest effect was observed with mean temperature. The combined effects of rainfall-humidity, temperature-humidity, and temperature-rainfall, although significant, were small.   Mean temperature ranged from 21°C to 34°C. The relationship between mean temperature and clinical malaria incidence displayed a bell-shape curve ( Figure 6). The risk of clinical malaria increased with an increase in mean temperature up to 27°C. At 23°C, the risk for clinical malaria was 53% of the risk at 27°C. Temperatures above 27°C led to a significant decrease in clinical malaria risk. The risk was therefore least at the lower and higher extremes of the temperature range.
The effect of rainfall on clinical malaria risk was observed only for values above 100 mm (Figure 7). Below that level, the risk reduction compared to 164 mm was almost 100%. There was no difference in risk between the rainfall levels of 60 mm and 90-mm. Above 100 mm, the risk increased significantly for each increase of 10 mm of rain. For example, the RR at 150 mm was 0.3 compared to 0.7 at 160 mm.
Relative humidity observed in the field ranged between 15% and 80%. The clinical malaria risk increased with an increase in relative humidity, but not linearly (Figure 8). Below 60% humidity, there was low risk. For example at 55% relative humidity, the risk was 25% lower than at 60% relative humidity. The risk of clinical malaria increased exponentially at relative humidity exceeding 60%.

Discussion
The individual effects of temperature, rainfall, and relative humidity of the preceding month on clinical malaria risk among children under age five were explored. The results confirmed that all meteorological factors are predictors for malaria risk. However, mean temperature was the best predictor.
Mean temperature alone was a strong predictor of clinical malaria rates among the children in our sample. Its rela-Effect of relative humidity on clinical malaria risk (RR = rate ratio) among study children  Rainfall mm (lag of one month) Rate ratio of clinical malaria tionship with clinical malaria was bell-shaped such that the risk was lowest at low and high temperatures. The infection risk was highest at temperatures around 27°C. There is no direct biological link between human infection and ambient temperature. However, the temperature-dependence of the Anopheles gambiae (the primary malaria vector in this region) life cycle could explain these results. Indeed, the sporogonic cycle (the development of the parasite within the vector) is shorter (9-10 days) at 28°C and longer (about 100 days) at temperatures below 20°C [1,24]. Therefore, a short incubation period allows the vector to live long enough to become infectious. The lifespan of the Anopheles gambiae is about 21 days. At high temperatures, vector survival decreases and temperatures of 40°C will result in zero vector survival. Craig and colleagues [4] reported that the proportion of vectors surviving the incubation period is high at temperatures between 28°C and 32°C. This explains the high-risk of clinical malaria observed in this range of temperature. In addition, temperatures above 27°C increase the feeding frequency of the vector to every two days, as the blood digestion rate increases. This results in more frequent vector-human contact. These findings agree with those from several other studies [25,26,7,9].
Cumulative rainfall from the previous month had a positive effect on clinical malaria rates as well. Several other studies reported similar results [26,6,9,10]. The relationship was 'J-shaped', suggesting that the relationship is non-linear, contrary to the results reported by Teklehiemanot and colleagues [9]. In addition, a minimum amount of rainfall is necessary for transmission to take place. In our study, a minimum of 100 mm was necessary before any case of clinical malaria could be observed. This could be explained by the hot and dry climate of the study area leading to high evaporation and water infiltration into the soil. Small amounts of rainfall will evaporate or infiltrate over a shorter time. Conversely, with higher amounts of rainfall, some water remains long enough for the vector to complete its development cycle. At higher amounts of rainfall, further increases may have little effect on clinical malaria risk, suggesting that a saturation level exists, as reported by Teklehiemanot and colleagues [9]. However, because of the relatively low amount of maximum rainfall (160 mm), the model could not capture this saturation phenomenon. As with temperature, the effect of rainfall is not directly related to human infection; but rainfall influences vector population abundance by providing open surface water for breeding. A positive effect of rainfall on malaria transmission has been reported widely in the literature, even though some studies have shown a converse or negligible effect [27][28][29][30]7]. This might be explained by the overwhelming effect of temperature. For rainfall to have a positive effect on malaria incidence, the temperature must be warm enough to support mosquito and parasite development [31]. This explains the low rate of malaria transmission in highlands, where rainfall is abundant, but temperatures are relatively low (below 20°C).
The individual effect of relative humidity on risk of clinical malaria, compared to temperature and rainfall, has been reported less often in the literature. Mean relative humidity of the previous month had a positive effect on malaria incidence. The effect was greater than that of rainfall. This is in line with the findings of Bi and colleagues [32]. To account for the strong correlation between rainfall and relative humidity an interaction term was included in the model. Thus, the estimate represents the individual effect of relative humidity on clinical malaria risk. The relative humidity observed ranged from 15% to 80%, with an average of 50%, mainly driven by rainfall and temperature. At a relative humidity below 60%, there was a decline in the risk of clinical malaria. This could be explained by a decrease in vector lifespan under these conditions, therefore, reducing the proportion of vector survival [33]. In contrast, at relative humidity above 60%, the infection rate increases substantially as humidity increases. The risk at 80% humidity is twice the risk at 60%. Again, this can be explained by improved vector survival.
The study was performed only for one year; thus, it cannot claim to be representative of all years. In fact, rainfall and humidity were both below the average of the last 10 years (data from the national meteorological station), but average temperatures were similar. Clinical data from several years (at least three years) of observation would have been ideal; but this was not feasible, given the costs involved in collecting the malaria infection data. In addition, the reported malaria incidence was not significantly different from that reported in a study conducted in the preceding year. Without a catastrophic event (like flooding or drought), the effect of meteorological variables on malaria infection should not change much from one year to another.
In such a semi-immune population, malaria may be transient. Therefore, weekly visits may not be sensitive enough to capture all the episodes. Infection acquired between visits could have been treated and cured and therefore missed. Indeed, Snow and colleagues reported a sensitivity of 74% of weekly visits in capturing clinical malaria episodes compared to daily visits [34]. In this setting, we expect a slightly higher sensitivity than the 74% because the most commonly used drug at the time was chloroquine which was reported to have a treatment failure of 36% [35]. Also, unless the child was severely ill, parents preferred to wait for the next visit of the interviewer because treatment of study participants was free. It should also be noted that Snow and colleagues did not find statistically significant differences in the risk of detecting malaria episodes when comparing daily visits to weekly ones [34]. Considering the cost and logistics involved in daily visits, the more cost-effective alternative was the weekly visits. The probable underestimation of the clinical malaria episodes notwithstanding, we do not expect the nature of the relationships with meteorological parameters to differ greatly from our findings.

Conclusion
The findings of this study confirm that meteorological conditions have a significant influence on clinical malaria rates among children under-five years old. Although several individual meteorological parameters have an impact on clinical malaria incidence, mean temperature is the best predictor and the main driving force, at least in the region. This is also observed in other regions. For example in the highlands region the preventing factor for malaria transmission is the low temperature because of the high altitude. Changes in ecological settings because of climate change are to be expected in this sub-Saharan region, especially the rain patterns. These changes may modify local mosquito microhabitats and affect transmission widely. A health information system including systematic monitoring of temperature and rainfall could yield an early warning system to support malaria control efforts at district level.