The changing burden of malaria and association with vector control interventions in Zambia using district-level surveillance data, 2006–2011
© Kamuliwo et al.; licensee BioMed Central Ltd. 2013
Received: 3 July 2013
Accepted: 25 November 2013
Published: 1 December 2013
Malaria control was strengthened in Zambia over the past decade. The two primary interventions for vector control are indoor residual spraying (IRS) and long-lasting insecticide-treated nets (LLINs). Using passive malaria surveillance data collected from 2006 to 2011 through the Zambian District Health Information System, the associations between increased coverage with LLINs and IRS and the burden of malaria in Zambia were evaluated.
National passive malaria surveillance data from 2006 to 2011 were analysed. A district-level, random-effects model with Poisson regression was used to explore the association between malaria cases and coverage with LLINs and IRS. Malaria cases and LLINs and IRS coverage were mapped to visualize spatiotemporal variation in malaria for each year.
From 2006–2011, 24.6 million LLINs were distributed and 6.4 million houses were sprayed with insecticide. Coverage with LLINs was not uniformly distributed over the study period and IRS was targeted to central and southern districts where malaria transmission was low. LLIN coverage was associated with a reduction in malaria cases, although an increase in the number of malaria cases was reported in some districts over the study period. A high burden of malaria persisted in north-eastern Zambia, whereas a reduction in the number of reported malaria cases was observed in western and southern Zambia.
Enhanced and targeted interventions in north-eastern Zambia where the burden of malaria remains high, as well as efforts to sustain low malaria transmission in the south-west, will be necessary for Zambia to achieve the national goal of being malaria free by 2030.
An estimated 219 million episodes of malaria and 660,000 malaria deaths occurred worldwide in 2010 . Approximately 80% of malaria episodes and 90% of the deaths were reported from the African continent . International funding for malaria control rose to a peak of USD 1.84 billion in 2012 . These financial commitments to malaria control made increased coverage of four key interventions possible: distribution of long-lasting insecticide-treated nets (LLINs), scale-up of indoor residual spraying (IRS), early diagnosis and treatment with artemisinin-combination therapy, and intermittent preventive treatment in pregnancy (IPTp) [1, 2]. Coverage with LLINs and IRS increased rapidly in some countries of sub-Saharan Africa, with household LLIN ownership across the region reaching 53% by 2012 and IRS protecting 11% of the population at risk .
Zambia has a history of highly endemic malaria transmission with varying levels of malaria control efforts over the past 50 years . In the past decade, malaria control was strengthened, leading to substantial increases in coverage following the availability of external funding . Zambia was successfully awarded the Global Fund Rounds 1, 4 and 7 for malaria control and was successful in obtaining financial support from the US President’s Malaria Initiative, the UK’s Department for International Development and the World Bank to scale-up LLIN and IRS coverage. However, through 2011, malaria remained endemic in all of the country’s 72 administrative districts [3, 5].
The two primary interventions for vector control in Zambia are IRS and LLINs. Using passive malaria surveillance data collected from 2006 to 2011 through the Zambian District Health Information System (DHIS), the district-level associations between IRS and LLIN coverage and the number of reported malaria cases were evaluated.
Zambia is a land-locked country in southern Africa, bordering eight malaria-endemic countries: Angola, Botswana, Democratic Republic of Congo, Malawi, Mozambique, Namibia, Tanzania, and Zimbabwe. The population of Zambia is approximately 13 million, with two-thirds of the population residing in rural areas . The country has a land area of 752,618 sq km, with an overall population density of 17 persons per sq km. The average life expectancy in Zambia is 49 years for men and 50 years for women .
Zambia is stratified into three malaria epidemiological zones: a low-transmission region in south-eastern Zambia with parasite prevalence <1%; a low stable-transmission zone in north-western/south-central Zambia with a parasite prevalence of 10%; and a high-transmission zone in northern and eastern Zambia with a parasite prevalence of >20% . Zambia experiences three distinct seasons, with a rainy season extending from November to May during which malaria transmission peaks, a cool dry season from late May to August, and a hot dry season from September to November. Persistent year-round malaria transmission occurs in focal locations across the country .
Malaria in Zambia is almost entirely caused by Plasmodium falciparum infection (~98%), with a low frequency of infections due to Plasmodium malariae and Plasmodium ovale, and little or no transmission of Plasmodium vivax. The main malaria vectors in Zambia are Anopheles gambiae s.s., Anopheles arabiensis and Anopheles funestus s.s[7, 8]. The epidemiology of malaria in Zambia is the result of two major forces: geographic, climatic and social features that establish ecological conditions conducive or restrictive to malaria transmission, and the introduction and scaling-up of malaria control interventions. Both forces alter malaria transmission and infection, as well as associated morbidity and mortality.
Malaria morbidity and mortality by district
Annual reported, district-level, aggregated malaria case data were obtained from the Zambian DHIS, the equivalent of the national health management information system. DHIS data are collected at heath facilities (including district hospitals and health centres) in paper form and are sent to the District Health Office (DHO) for data capture and validation in the DHIS. The DHO transmits the district-level data to the provincial office for further processing and aggregation for the province. The consolidated provincial data are then transmitted to the National Malaria Control Centre (NMCC), which in turn provides feedback and technical support to the provinces, districts and health facilities.
Severe malaria, a set of clinical and laboratory parameters associated with an increased risk of death in the presence of P. falciparum parasitaemia, was reported from hospital admission records [9, 10]. Malaria episodes, including deaths, were reported by age group: <five and ≥ five years. From 2006 to 2008, only the number of malaria episodes was reported. Clinical and confirmed malaria by rapid diagnostic test (RDT) or microscopy were reported separately from 2009 to 2011.
Spatial and population data
District polygon boundaries (shape files) were obtained from the Zambian Ministry of Environment and Statistics in Lusaka. For 2010, demographic and housing data were obtained from the Zambian Bureau of Statistics (ZBS) . District-level demographic data were projected for 2006, 2007, 2008, 2009, and 2011 using an exponential population growth model based on ZBS reports for 2000 and 2010. The number of houses in each district was projected based on 2000 and 2010 estimates, assuming linear growth. These estimates served as denominators for the calculation of rates [11, 12]. Morbidity and mortality rates attributable to malaria were calculated by age group (<five and ≥ five years) per 1,000 population based on data from the DHIS and population estimates. This study was deemed not be human subjects research by the Johns Hopkins Bloomberg School of Public Health Institutional Review Board.
IRS and LLIN coverage by district
Data on IRS coverage were captured using daily spray forms that were annually consolidated at the NMCC. Data were also available on the number of LLINs distributed annually through various distribution channels, including antenatal and under-five clinics and mass vaccination campaigns. Both IRS and LLIN coverage data were aggregated at the district level. LLIN rates were calculated as the number of nets per two household residents and an assumed three-year life span of LLINs [13–15]. IRS rates were calculated per 1,000 houses. IRS and LLIN coverage rates also were estimated at the provincial level.
Poisson regression was used to explore the association between the burden of malaria, as measured by the total number of malaria cases either diagnosed clinically or confirmed by microscopy or RDT, and LLIN and IRS coverage at the district level using 2006 as the baseline year. The random effect of the district variable was used to model variation between districts and by year. All statistical analyses were performed using STATA 11 (STATA Corp. 2003, College Station, TX, USA).
Annual malaria cases and annual IRS and LLIN coverage data were linked with district shape files. Spatial autocorrelation was tested for total yearly malaria cases using incremental spatial autocorrelation. The local Getis-OrdGi* statistic  was used to generate hotspot maps for malaria prevalence and IRS and LLIN coverage from 2006 to 2011. The local Getis-OrdGi* statistic compares the local mean rate (the rates for a district and its nearest neighbouring districts) to the global mean rate (the rates for all districts), producing a z-score and p-value for each district, reflecting whether the differences between the local and global means are statistically significant . A statistically significant positive z-score indicates a hotspot for high rates. Similarly, a statistically significant negative z-score for a district indicates local spatial clustering of low rates [16, 18, 19]. To identify persistent malaria hotspots, the number of malaria cases in each year from 2007 to 2011 was subtracted from the number of cases in 2006. This change in the number of malaria cases for each year also was used in a hot spot analysis. These analyses were conducted using ArcGIS 10 (ESRI, Redlands, CA, USA).
Changes in the burden of malaria
Reported national annual cases for malaria indicators, 2006 to 2011
Malaria prevalence/1,000 persons
Malaria prevalence <5 years/1,000 persons
Malaria prevalence ≥5 years/1,000 persons
Severe malaria/1,000 persons
Severe malaria <5 years/1,000 persons
Severe malaria ≥5 years/1,000 persons
Malaria related mortality/1,000 persons
Malaria related mortality <5 years/1,000 persons
Malaria related mortality ≥5 years/1000 persons
Clinical malaria cases (%)
Confirmed* malaria cases (%)
LLIN distribution/per two household residents
IRS coverage/1,000 houses
Reporting of both clinical diagnoses and confirmed cases by RDT or microscopy started in Zambia from 2009. Only half (50.2%) of the reported cases were confirmed by either microscopy or RDT in 2011 (Figure 1).
LLIN and IRS coverage
Associations between vector control activities and the burden of malaria
Impact of IRS and LLIN on the prevalence of malaria in Zambia, 2006 to 2011
Deaths due to malaria
PRR* (95% CI)
PRR (95% CI)
PRR (95% CI)
IRS (per 1,000 household coverage)
LLIN (per two household residents coverage)
0 .54 (0.51–0.57)
Impact of IRS and LLIN on the prevalence of malaria by age group in Zambia, 2006 to 2011
PRR* (95% CI)
PRR (95% CI)
PRR (95% CI)
PRR (95% CI)
PRR (95% CI)
PRR (95% CI)
IRS (per 1,000 household coverage)
LLIN (per two household residents coverage)
Zambia made much progress in reducing the burden of malaria between 2006 and 2011, in part through widespread coverage with vector control interventions. This reduction in the burden of malaria was associated with LLIN coverage at the district level. However, malaria remains a major cause of morbidity and mortality in parts of the country, and particularly among children younger than five years of age. After six years of malaria control interventions, both resurgent and persistently stable malaria hotspots were observed in the north-eastern districts of Zambia. Additional malaria control strategies will be needed to reduce the burden of malaria in these regions.
Overall, Zambia achieved distribution of one LLIN for every two persons. However, LLINs may lose effectiveness with the emergence of insecticide resistance [5, 8]. In 2010, significant resistance to deltamethrin, permethrin and DDT was detected in both An. gambiae s.s. and An. funestus s.s. in several districts in Zambia . Vector and insecticide resistance mapping in malaria hotspots could help guide vector control strategies. Larval control has never been widely implemented in Zambia. However, with the emergence of insecticide resistance, larval control through environmental management and larviciding could provide alternatives for vector control in some areas. The selection of larvicides, insecticides for IRS and deployment of LLINs must be guided by detailed knowledge of vector species, their biting and oviposition behaviours and patterns of resistance.
High coverage with IRS targeted a subset of districts in central Zambia, with little extension into the southwest of the country, and was conducted in mostly urban areas. Further progress will require additional strategies, with a focus on north-eastern districts that consistently reported significant and stable malaria hotspots .
All countries sharing a border with Zambia have endemic malaria. Cross-border collaborations with neighbouring countries should be strengthened, such as the recent Zambia Zimbabwe Cross-Border Malaria Initiative [3, 22]. Persistently high malaria transmission along the Mozambique and Malawi borders indicates the need for additional cross-border collaborations, with monitoring of population movement and malaria transmission. The low level of malaria transmission and decline in the number of malaria cases in southern Zambia, which shares a border with Namibia, Botswana and Zimbabwe, demonstrates the unprecedented achievement of the National Malaria Control Programme.
This study has several limitations, most notably the fact that these analyses were based on routinely collected data within the DHIS with the potential for both over and under reporting of malaria cases. Secular changes in reporting could have biased prevalence estimates and errors could have been introduced as data were aggregated at higher levels of the health information system. However, data accuracy and completeness were not systematically assessed. Many cases reported in the DHIS were not confirmed by RDT and microscopy but were based on clinical signs and symptoms, with the potential for misclassification. Both RDT and microscopy have limited sensitivity and specificity, and are particularly likely to misclassify individuals with low levels of parasitaemia. The observed associations between LLIN and IRS coverage and the burden of malaria were ecologic and not at the level of individuals.
Sustained efforts are critical to maintain the gains that have already been achieved and further progress may bring malaria transmission to near-zero in all four districts that reported fewer than ten cases per 1,000 population in 2011, thus creating malaria-free zones in Zambia. Enhanced and targeted interventions in north-eastern Zambia where the burden of malaria remains high, as well as sustaining low malaria transmission in the southwest, will be necessary for Zambia to achieve the national goal of being malaria free by 2030.
UH is supported by an A Ralph and Sylvia E Barr Fellowship from the W Harry Feinstone Department of Molecular Microbiology and Immunology, Johns Hopkins Bloomberg School of Public Health.
- WHO: World malaria report. 2012, Geneva: World Health Organization, Available from: http://www.who.int/malaria/publications/world_malaria_report_2012/en/index.html [Accessed 15 June 2013]Google Scholar
- Sutcliffe CG, Kobayashi T, Hamapumbu H, Shields T, Mharakurwa S, Thuma PE, Louis TA, Glass G, Moss WJ: Reduced risk of malaria parasitemia following household screening and treatment: a cross-sectional and longitudinal cohort study. PLoS One. 2012, 7: e31396-10.1371/journal.pone.0031396.PubMed CentralView ArticlePubMedGoogle Scholar
- NMCC: Zambia national malaria control program performance review. 2011, Lusaka, Zambia: National Malaria Control Centre, Ministry of HealthGoogle Scholar
- Chizema-Kawesha E, Miller JM, Steketee RW, Mukonka VM, Mukuka C, Mohamed AD, Miti SK, Campbell CC: Scaling up malaria control in Zambia: progress and impact 2005–2008. Am J Trop Med Hyg. 2010, 83: 480-488. 10.4269/ajtmh.2010.10-0035.PubMed CentralView ArticlePubMedGoogle Scholar
- NMCC: Achievements in malaria control, the Zambian story 2000–2010. 2011, Lusaka, Zambia: National Malaria Control Centre, Ministry of HealthGoogle Scholar
- United Nations: World mortality report. 2011, Available from: http://www.un.org/en/development/desa/population/publications/pdf/mortality/worldMortalityReport2011.pdf [Accessed 15 June, 2013]Google Scholar
- Chanda E, Hemingway J, Kleinschmidt I, Rehman AM, Ramdeen V, Phiri FN, Coetzer S, Mthembu D, Shinondo CJ, Chizema-Kawesha E, Kamuliwo M, Mukonka V, Baboo KS, Coleman M: Insecticide resistance and the future of malaria control in Zambia. PLoS One. 2011, 6: e24336-10.1371/journal.pone.0024336.PubMed CentralView ArticlePubMedGoogle Scholar
- Norris LC, Norris DE: Efficacy of long-lasting insecticidal nets in use in Macha, Zambia, against the local Anopheles arabiensis population. Malar J. 2011, 10: 254-10.1186/1475-2875-10-254.PubMed CentralView ArticlePubMedGoogle Scholar
- Marsh K, Forster D, Waruiru C, Mwangi I, Winstanley M, Marsh V, Newton C, Winstanley P, Warn P, Peshu N: Indicators of life-threatening malaria in African children. N Engl J Med. 1995, 332: 1399-1404. 10.1056/NEJM199505253322102.View ArticlePubMedGoogle Scholar
- WHO: Severe falciparum malaria. Trans R Soc Trop Med Hyg. 2000, 94 (Suppl 1): S1-S90.Google Scholar
- Central statistical office: Census of population and housing. 2011, Lusaka, Zambia: Ministry of Finance and National PlanningGoogle Scholar
- Central statistical office: Census of population and housing preliminary report 2001. 2001, Lusaka, Zambia: Ministry of Finance and National PlanningGoogle Scholar
- WHO: Insecticide-Treated Mosquito Nets: a WHO Position Statement. 2007, http://www.who.int/malaria/publications/atoz/itnspospaperfinal.pdf [Accessed 10.30.2013]Google Scholar
- Pulkki-Brannstrom AM, Wolff C, Brannstrom N, Skordis-Worrall J: Cost and costeffectiveness of long-lasting insecticide-treated bed nets - a model-based analysis. Cost Eff Resour Alloc. 2012, 10: 5-10.1186/1478-7547-10-5.PubMed CentralView ArticlePubMedGoogle Scholar
- Kilian A, Byamukama W, Pigeon O, Gimnig J, Atieli F, Koekemoer L, Protopopoff N: Evidence for a useful life of more than three years for a polyester-based long-lasting insecticidal mosquito net in Western Uganda. Malar J. 2011, 10: 299-10.1186/1475-2875-10-299.PubMed CentralView ArticlePubMedGoogle Scholar
- Ord JK, Getis A: Local spatial autocorrelation statistics: distributional issues and an application. Geogr Anal. 1995, 27: 286-306.View ArticleGoogle Scholar
- Haque U, Scott LM, Hashizume M, Fisher E, Haque R, Yamamoto T, Glass GE: Modelling malaria treatment practices in Bangladesh using spatial statistics. Malar J. 2012, 11: 63-10.1186/1475-2875-11-63.PubMed CentralView ArticlePubMedGoogle Scholar
- Getis A, Ord J: The analysis of spatial association by use of distance statistics. Geogr Anal. 1995, 24: 189-206.View ArticleGoogle Scholar
- ESRI: http://www.esriurl.com/spatialstats [Accessed 15 June, 2013]
- Coleman PG, Goodman CA, Mills A: Rebound mortality and the cost-effectiveness of malaria control: potential impact of increased mortality in late childhood following the introduction of insecticide treated nets. Trop Med Int Health. 1999, 4: 175-186. 10.1046/j.1365-3156.1999.43382.x.View ArticlePubMedGoogle Scholar
- Bousema T, Griffin JT, Sauerwein RW, Smith DL, Churcher TS, Takken W, Ghani A, Drakeley C, Gosling R: Hitting hotspots: spatial targeting of malaria for control and elimination. PLoS Med. 2012, 9: e1001165-10.1371/journal.pmed.1001165.PubMed CentralView ArticlePubMedGoogle Scholar
- Macepa M: A cross-border partnership in Zambia and Zimbabwe: A step toward elimination on World Malaria Day. Available from: http://www.makingmalariahistory.org/a-cross-border-partnership-in-zambia-and-zimbabwe-a-step-toward-elimination-on-world-malaria-day/ [Assessed 25 June 2013]
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.