- Open Access
Spatio-temporal analysis of the role of climate in inter-annual variation of malaria incidence in Zimbabwe
International Journal of Health Geographicsvolume 5, Article number: 20 (2006)
On the fringes of endemic zones climate is a major determinant of inter-annual variation in malaria incidence. Quantitative description of the space-time effect of this association has practical implications for the development of operational malaria early warning system (MEWS) and malaria control. We used Bayesian negative binomial models for spatio-temporal analysis of the relationship between annual malaria incidence and selected climatic covariates at a district level in Zimbabwe from 1988–1999.
Considerable inter-annual variations were observed in the timing and intensity of malaria incidence. Annual mean values of average temperature, rainfall and vapour pressure were strong positive predictors of increased annual incidence whereas maximum and minimum temperature had the opposite effects. Our modelling approach adjusted for unmeasured space-time varying risk factors and showed that while year to year variation in malaria incidence is driven mainly by climate, the resultant spatial risk pattern may to large extent be influenced by other risk factors except during high and low risk years following the occurrence of extremely wet and dry conditions, respectively.
Our model revealed a spatially varying risk pattern that is not attributable only to climate. We postulate that only years characterized by extreme climatic conditions may be important for developing climate based MEWS and for delineating areas prone to climate driven epidemics. However, the predictive value of climatic risk factors identified in this study still needs to be evaluated.
The risk of malaria infection varies widely according geographic region, season and year [1, 2]. On the fringes of endemic zones, particularly at the southernmost latitudes in southern Africa, across arid regions of northern Africa and among the highlands of east and central horn of Africa climate is a major determinant of seasonal and inter-annual (year to year) variation in malaria transmission .
In Southern Africa annual variation in climatic conditions and associated changes in malaria infection affect the timing and intensity of malaria incidence. This has an impact on the effectiveness of interventions . As a result there is a need for the development of climate-based malaria early warning systems (MEWS) capable of predicting seasonal to inter-annual variations with a lead time that allows health authorities to respond in a timely manner with preparatory/preventative measures [5, 6]. However, despite the fact that climate data are often used to account for spatial, seasonal and inter-annual variation in malaria risk in Africa, there is often little or no consensus about the relative importance and predictive value of different factors involved . The disagreements seem to stem from differences in perspective and methods used .
In Southern Africa, few studies or models of the relationship between malaria and climatic factors have been published. In a recent meeting of the Southern African Inter-Country Programme on Malaria Control (SAMC) a number of countries acknowledged having poor empirical basis on which to develop and test climate-based early warning and detection indicators . Existing models include a non-spatial model for Zimbabwe which identified temperature as the main determinant of increased malaria risk years . Spatial and temporal models which used both temperature and rainfall for analysis and mapping of malaria risk in KwaZulu-Natal, South Africa [10, 11]. An exploratory analysis of 30 years worth of data in the epidemic prone area of KwaZulu-Natal, South Africa, which showed that certain aspects of climate appear to drive inter-annual variation of malaria incidence but not its overall level . In a recent study Thomson et al  showed that rainfall and sea surface temperature (SST) have a potential for application in the development of seasonal forecasts . We also recently developed a space-time seasonality model based on the relationship between monthly clinical malaria case data and environmental factors in Zimbabwe .
Year to year predictability of malaria incidence still remains a challenge, and more work is required before malaria climate based forecasting models can realize their full potential in the region. In this study we use Bayesian spatio-temporal analysis to describe year to year variation of malaria incidence data from Zimbabwe in relation to variation in climatic risk factors to enhance our ability of developing an operational MEWS and determine areas prone to climate-driven epidemics.
In common with Angola, Namibia, Botswana, Zambia, Mozambique and South Africa Zimbabwe lies at the southern limits of malaria distribution in Africa. Malaria remains a major cause of mortality and morbidity despite more than four decades of sustained national control programme . Moreover, as a result of reduced level of transmission there is propensity for malaria epidemics unless adequately controlled or prevented. Overall, about 45–50% of the 12.5 million people of Zimbabwe are at risk of malaria. In 1998, it was estimated that approximately 8% of all deaths and 12% of all outpatient cases were due to malaria http://www.malaria.org/zw/countries/zimbabwe.htm. Recently, substantial socio-economic changes have further compromised the malaria control programme .
Climate is another major factor that determines the extent of malaria transmission in Zimbabwe, and its variability may work with or against efforts to bring malaria under control . The most important factors governing malaria epidemiology in the country are season, altitude, and associated rainfall and temperature changes [14, 16–18]. Malaria is found mainly within the low and mid altitude zones and rarely at higher altitude. However, this can vary tremendously from one year to another.
We used normalized difference vegetation index (NDVI) available between 1988 and 1999 from Advanced Very High Resolution Radiometer (AVHRR) sensor onboard the National Oceanic and Atmospheric Administration (NOAA) satellite http://daac.gsfc.nasa.gov/data/dataset/AVHRR/index.html. We also used mean annual values of rainfall, vapour pressure, minimum, maximum and mean temperature obtained for each district and year for the 12 year period. These were obtained from the climate research unit (CRU) climate surfaces derived from interpolated weather station data as a function of latitude, longitude, and elevation using thin-plate splines .
In Zimbabwe malaria is a notifiable disease and records from hospitals and clinics are compiled at a district, provincial and national levels to describe the malaria situation and trends. We used annual clinical malaria case data for children under the age of five reported in 58 districts covering the whole country between 1988 and 1999 . This is the highest risk group, with relatively little protective immunity and is therefore expected to be more sensitive to changes in malaria transmission. The data included both microscopically confirmed and unconfirmed clinically diagnosed cases. District population projections based on the 1982 and 1992 census were used to calculate incidence rate per 1000 person years .
We used the annual proportion of monthly malaria cases and Markham's seasonality index [22, 23] to display between-year variation in the data. The seasonality index has been described in detail elsewhere . Briefly, this method calculates the seasonal concentration of the malaria case load and the peak month in a given year.
A preliminary negative binomial regression analysis was carried out in STATA 9.0 (Stata Corp., College Station, TX, USA) to assess the relationship between annual malaria incidence and annual values of each climatic covariate. Thereafter, Bayesian negative binomial models were fitted in WinBUGS  to examine the association between inter-annual variation in malaria incidence and a combination of climatic covariates selected from the preliminary analysis (see the appendix for more details). Basically, spatial random effects were used at a district level to take into account spatial correlation present in the data. Temporal random effects were also used at yearly intervals to account for temporal correlation. Spatial correlation was incorporated by assuming a conditional autoregressive (CAR) process in the random effects. A first order autoregressive process was applied for temporal random effects .
Markov Chain Monte Carlo simulation (MCMC) was applied to estimate model parameters . After the initial burn-in of 5000 the number of iterations thereafter depended on convergence which was assessed using ergodic averages. After convergence a final sample of 5000 was collected to obtain summaries of the posterior distribution of the parameters. The Deviance Information Criterion (DIC)  was used for the comparison of model fit. Small values of DIC indicate superior model fit. Model estimates were exponentiated to represent incidence rate ratios (IRR), that is, per unit change in incidence for each covariate.
Figure 1A–D shows that malaria transmission in Zimbabwe is characterized by considerable between and within year variations. The highest malaria incidence recorded during the 12 year period was in 1988 with 15.5 cases per 1000 person years and the lowest was in 1992 with 5.2 cases per 1000 person years (Figure 1A). In addition since 1996, there is a rising trend in annual incidence with reported malaria cases remaining at high levels. The intensity and timing of the seasonal peak also varies from year to year (Figure 1B &1D). From 1988 to 1999 the peak month fluctuated between March and April with the exception of 1992 and 1995 which were characterized by peaks in January and May respectively.
High annual malaria incidence coincide with high rainfall and relatively warm conditions while low incidence years coincide only with low rainfall (Figure 2). Vapour pressure and NDVI follow the rainfall pattern. Temperature derived covariates seem to be important only in the presence of sufficient rainfall. The intensity and timing of the seasonal peak in each year appears to follow variability in rainfall.
In the bivariate analysis all selected covariates showed a significant relationship (P < 0.001) with malaria incidence (Table 1). All models (Table 2) indicated that mean annual temperature, rainfall, vapour pressure and NDVI were strong positive predictors of increased annual incidence rate in contrast maximum and minimum temperature had a reducing effect. However, in the spatial model (model 2) rainfall had no significant effect and in the spatio-temporal model (Model 3) only NDVI was not significant. Model comparison showed that the spatial-temporal model had a small DIC value and therefore was the best fitting model. This model had large spatially correlated random effects.
Figure 3 show differences in the spatial pattern of modelled malaria incidence during the 12 year study period. In 1988, 1993 and 1996 onwards high to moderate incidence rates were more widespread with the highest incidence rates in the north western and eastern part of the country. In 1991 to 1992 and 1994 to 1995 incidence rates were predominantly moderate to very low levels across the country with pockets of high incidence rates in the districts situated along of the Zambezi river system in the north western part and on the border with Mozambique in the eastern part as well as along the Limpopo river system in the south eastern part.
Our observations confirm that malaria transmission in Zimbabwe is characterized by substantial inter-annual variation. Incidence rates show a rising, albeit fluctuating, trend with reported cases remaining relatively high after 1996 (Figure 1). The intensity and timing of seasonal transmission also varies from year to year. The highest and lowest risk years were recorded during one of the most wet and severe drought periods in Zimbabwe, respectively . From a regional perspective 1988, 1993, 1996 and 1997 were the most serious epidemic years .
Besides the limitation imposed by the effect of misdiagnosis as a consequence of unconfirmed clinically diagnosed cases, which would probably be to smooth over differences since this is not influenced by covariates. There are multiple explanations for the observed trend and it is difficult to identify true causes . In general, there is more support for non-climatic explanations of recent trends especially in the fringe areas of endemic zones in Africa. These include deterioration of malaria control efforts, development of drug and insecticide resistance and a rise in co-infection with HIV/AIDS [8, 29–31]. In Zimbabwe, none of these have been adequately quantified for inclusion into the modelling framework. Hence our modelling approach adjusts for unmeasured spatially and temporally structured sources of variation while investigating the association between annual malaria incidence and climatic covariates. However, the large residual spatial variation observed in the spatial-temporal model suggests that there are other important covariates not accounted for in the analysis that could explain most of the spatial variation in malaria incidence. Most of the temporal variation in the data appears to be explained by the selected covariates.
Our results showed that annual mean values of temperature, rainfall and vapour pressure are strong positive predictors of increased annual malaria incidence. Annual minimum and maximum temperature had the opposite effects probably due to fact that in an average annual cycle these climatic covariates are associated with the cold dry period (June-August) and hot dry period (August to October), respectively, both of which have low transmission . The association between vapour pressure (a measure of humidity) and annual incidence reflect the importance of the interaction between rainfall and temperature which modulates the ambient air humidity which in turn affect the survival and activity of Anopheles mosquitoes [32, 33].
Furthermore, the spatio-temporal model (Table 2) showed that while NDVI, a surrogate for the response of vegetation to rainfall, appears to have a positive effect, it had no statistically significant association with annual malaria incidence. This may be due to the fact that vegetation greenness is to a large extent dependant on the amount of rainfall available in a given year [34, 35]. Moreover, total annual rainfall in Zimbabwe is characterized by strong variability . Hence, rainfall is a stronger predictor of malaria incidence at an inter-annual time scale.
The modelled spatial pattern showed that areas in the low lying north-western and south-eastern part of the country have high but less variable incidence, whereas areas of lower incidence in the middle to highland areas showed greater year to year variability. Eco-epidemiological conditions in these areas determine the stability of between year variation in malaria transmission [16, 34, 36]. Our analysis therefore corroborates the fact that year to year variation in malaria incidence in Zimbabwe is driven mainly by climatic covariates, and further demonstrates that resultant incidence in a given area may also be a function of other unmeasured risk factors. The spatial effect of climatic conditions during the study period appears to be more evident only during the much reported drought years from 1990–1992 and following heavy rainfall in 1998 .
In our previous analysis optimum ranges of climatic risk factors identified in the present study were also significantly associated with malaria transmission at a seasonal time scale . It is likely therefore that increase in mean annual temperature, rainfall and humidity is linked to within year changes in average climatic conditions and in turn lead to changes in the full annual cycle of malaria transmission. The magnitude of which probably varies from one year another. For example, Freeman and Bradley  found that higher than average annual mean monthly temperatures during the critical period of malaria transmission (April and September) in the previous year were associated with an increase in the severity of malaria in the following year in some areas. This becomes more pronounced and widespread if both extremes of warm and wet conditions coincide and vice versa as observed during high and low incidence years (Figure 2 and 3). We deduce that extreme climatic events associated with covariates identified in this study may be useful for developing climate based malaria forecasting models operational at both seasonal and inter-annual time scales.
In conclusion, our modelling approach adjusted for unobserved spatial and temporal varying risk factors, and showed that while inter-annual variation in malaria incidence is driven mainly by climatic conditions, the resultant spatially varying risk pattern may also be influenced by other risk factors. Nevertheless, high and low incidence years following the occurrence of extreme climatic conditions may be useful for developing climate based MEWS and for delineating areas prone to climate driven epidemics. However, the predictive value of climatic risk factors identified in this study still needs to be evaluated.
Appendix: statistical model
We assumed that the observed counts of malaria cases Y it in district i (i = 1,..., 58) and year t (1988–1999) follow a negative binomial distribution with parameters p it and r, that is, Y it ~ NB (pit, r), where pit relates to average number of cases via the formula (μ it ) = p it /r and r is the overdispersion parameter. We modelled average number of new cases (μ it ) as a function of potential risk factors as follows:
(1) log(μ it ) = log N it + α + β (non-spatial model)
(2) log(μ it ) = log N it + α + β+ ø i (spatial model)
(3) log(μ it ) = log N it + α + β+ ø i + ω it (spatio-temporal model)
N it denotes population at risk in district i and year t, α is the incidence rate when all covariates have zero value, X it is a vector of climatic covariate effect in district i and year t, βa vector of the regression coefficients, ø i is the spatial random effect for district i and ω it is the temporal random effect for year t and district i. District specific random effects were modelled via a conditional autoregressive (CAR) process, which implies that each φ i conditional on the neighbour φ j follows a normal distribution with mean equal to the average of neighbouring φ j and variance inversely proportional to the number of neighbours n i , that is, ø i |ø j , j neighbouring of i ~ Normal () where γ is a parameter that quantifies the amount of spatial correlation present in the data and measures the spatial variance. The temporal effects were modelled by first order autoregressive process with temporal variance , which allows correlation between consecutive time periods for each district and year. We assumed inverse gamma hyper-prior distributions for the variance parameters of the spatial and temporal random effects, non-informative Uniform prior distributions U(-∞, ∞) for the regression coefficients β, and a Uniform distribution U(a, b) for γ, with limits a and b specified as described in Gelfand and Vounatsou .
Advanced Very High Resolution Radiometer
Climate research unit
Deviance Information Criterion
Malaria early warning system
Markov Chain Monte Carlo simulation
National Oceanic and Atmospheric Administration
Normalized Difference Vegetation Index
Sea surface temperature
Southern African Inter-Country Programme on Malaria Control
Baird JK, Agyel SO, Utz GC, Koram K, Barcus MJ, Jones TR, Fryauff DJ, Binka FN, Hoffman SL, Nkrumah FN: Seasonal malaria attack rates in infants and young children in Northern Ghana. Am J Trop Med Hyg. 2002, 66: 280-286.
Thomson M, Indeje M, Connor S, Dilley M, Ward N: Malaria early warning and seasonal climate forecasts. Lancet. 2003, 362: 280-10.1016/S0140-6736(03)14135-9.
WHO: Climate change and human health, risks and responses. 2003, Geneva: World Health Organization
Da Silva J, Garanganga B, Teveredzi V, Marx SM, Mason SJ, Connor SJ: Improving epidemic malaria planning, preparedness and response in Southern Africa. Malar J. 2004, 3: 37-10.1186/1475-2875-3-37.
le Sueur D, Sharp BL: Malaria Forecasting Project: Report on Workshop on Reducing Climate Related Vulnerability in Southern Africa. 1996, Victoria Falls, Zimbabwe. Washington, D.C, USA: Office of Global Programs, National Oceanic and Atmospheric Administration, 65-73.
le Sueur D, Sharp BL, Gouws E, Ngxongo S: Malaria in South Africa. S Afr Med J. 1996, 8: 936-939.
Teklehaimanot HD, Lipsitch M, Teklehaimanot A, Schwartz J: Weather-based prediction of Plasmodium falciparum malaria in epidemic-prone regions of Ethiopia I. Patterns of lagged weather effects reflect biological mechanisms. Malar J. 2004, 3: 41-10.1186/1475-2875-3-41.
Bouma MJ: Climate change and tropical disease: methodological problems and amendments to demonstrate effects of temperature on the epidemiology of malaria, a new perspective on the highland epidemics in Madagascar, 1972–89. Trans R Soc Trop Med Hyg. 2003, 97: 133-139. 10.1016/S0035-9203(03)90099-X.
Freeman T, Bradley M: Temperature is predictive of severe malaria years in Zimbabwe. Trans R Soc Trop Med Hyg. 1996, 90: 232-10.1016/S0035-9203(96)90224-2.
Kleinschmidt I, Sharp BL, Clarke GPY, Curtis B, Fraser C: Use of generalized linear mixed models in the spatial analysis of small area malaria incidence rates in KwaZulu Natal, South Africa. Am J Epidemiol. 2001, 153: 1213-12121. 10.1093/aje/153.12.1213.
Kleinschmidt I, Sharp BL, Mueller I, Vounatsou P: Rise in malaria incidence rates in South Africa: a small area spatial analysis of variation in time trends. Am J Epidemiol. 2002, 155: 257-267. 10.1093/aje/155.3.257.
Craig MH, Kleinschmidt I, le Sueur D, Sharp BL: Exploring 30 years of malaria case data in KwaZulu-Natal, South Africa: part I. The impact of climatic factors. Trop Med Int Health. 2004, 9: 1247-1257. 10.1111/j.1365-3156.2004.01340.x.
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.
Mabaso MLH, Craig M, Vounatsou P, Smith T: Towards empirical description of malaria seasonality in southern Africa: the example of Zimbabwe. Trop Med Int Health. 2005, 10: 909-918. 10.1111/j.1365-3156.2005.01462.x.
Makono R, Sibanda S: Review of prevalence of malaria in Zimbabwe with specific references to parasite drug resistance (1984–96). Trans R Soc Trop Med Hyg. 1999, 93: 449-452. 10.1016/S0035-9203(99)90331-0.
Hartman J, Ebi K, McConnell K, Chan N, Weyant J: Climate suitability for stable malaria transmission in Zimbabwe under different climate change scenarios. Glob Change Hum Health. 2002, 3: 42-54. 10.1023/A:1019621231796.
Taylor P: The malaria problem in Zimbabwe epidemiology. Centr Afr J Med. 1985, 31: 63-165.
Taylor P, Mutambu SL: A review of the malaria situation in Zimbabwe with a special reference to the period 1972–1248. Trans R Soc Trop Med Hyg. 1986, 80: 12-19. 10.1016/0035-9203(86)90185-9.
Mitchell TD, Jones PD: An improved method of constructing a database of monthly climate observations and associated high-resolution grids. Int J Climatol. 2005, 25: 693-712. 10.1002/joc.1181.
Ministry of Health and Child Welfare: National Malaria Control Programme. Disease Prevention and control. 2000, Harare, Zimbabwe: Ministry of Health and Child Walfare
Central Statistics Office: District population projections in Zimbabwe. 2002, Causeway, Harare, Zimbabwe: Central Statistics Office
Markham CG: Seasonality of precipitation in the United States. Ann Assoc Am Geogr. 1970, 60: 593-597. 10.1111/j.1467-8306.1970.tb00743.x.
McGee OS: The determination of rainfall seasons in South Africa using Makharm's technique. S Afr Geogr. 1977, 5: 390-395.
WinBUGS: Software version 1.4. 2000, MRC Biostatistics Unit, Institute of public Health, UK: Robinson Way, Cambridge CB2 2SR
Bernardinelli L, Clayton D, Montomali C: Bayesian estimates of disease maps: How important are priors?. Stat Med. 1995, 14: 2411-2431.
Gelfand AE, Smith AMF: Sampling-based approach to calculating marginal densities. J Am Stat Assoc. 1990, 85: 398-409. 10.2307/2289776.
Spiegelhalter DJ, Best NG, Carlin BP, van der Linde A: Bayesian measures of model complexity and fit. J R Stat Soc. 2002, 64 (B): 583-639. 10.1111/1467-9868.00353.
Kovats S: El Nino and human health. Bull World Health Organ. 2000, 78: 1127-1135.
Hay SI, Cox J, Rogers DJ, Randolph SE, Stern DI, Shanks GD, Myers MF, Snow RW: Climate change and resurgence of malaria in East African Highlands. Nature. 2002, 415: 905-909. 10.1038/415905a.
Hay SI, Rogers DJ, Randolph SE, Stern DI, Cox J, Shanks GD, Snow RW: Climate change and resurgence of malaria in East African Highlands. Trends Parasitol. 2002, 18: 530-534. 10.1016/S1471-4922(02)02374-7.
Craig MH, Kleinschmidt I, Nawn JB, le Sueur D, Sharp BL: Exploring 30 years of malaria case data in KwaZulu-Natal, South Africa: part II. The impact of non climatic factors. Trop Med Int Health. 2004, 9: 1258-1266. 10.1111/j.1365-3156.2004.01341.x.
Bruce-Chwatt LJ: Essential Malariology. 1980, London: Heineman medical books Ltd
Molineaux L, Eds: The epidemiology of human malaria as an explanation of its distribution, including some implications for its control. Malaria: Principles and practices of malariology. Edited by: Wernsdorfer WH, McGregor I. 1988, London: Churchill Livingstone, 913-998.
Reid ET, Woods RW: Anopheline mosquitoes of southern Rhodesia, a general survey. Proc Trans Rhod Sci Assoc. 1957, XLV: 47-72.
Ramsdale CD: The effect of residual hut spraying with HCH on mixed populations of Anopheles gambiae complex in Rhodesia. WHO/Mal/508.65. 1965
Craig HM, 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.
Gelfand AE, Vounatsou P: Proper multivariate conditional autoregressive models for spatial data analysis. Biostat. 2003, 4: 11-25. 10.1093/biostatistics/4.1.11.
This work is part of the Mapping Malaria Risk in Africa (MARA) collaboration between the Malaria Research Lead Programme at the Medical Research Council in South Africa and the Swiss Tropical Institute in Basel Switzerland. We are indebted to the National Malaria Control Programme in Zimbabwe and the World Health Organization Southern Africa Inter-Country Programme for Malaria Control for assistance during data collection. This work was financed partially by the Rudolf Geigy Stiftung zu Gunsten des Schweizerischen Tropeninstituts and partially by the Swiss National Science Foundation Project 3252B0-102136/1.
The authors of this paper do not have any commercial interest or other associations which pose a conflict of interest with this paper.
MLHM conceptualized, analyzed and drafted the manuscript. PV and TM participated in the conception, analysis and drafting of the manuscript. SM and JDS critically reviewed the manuscript.