Temporal dynamic of malaria in a suburban area along the Niger River
© The Author(s) 2017
Received: 7 April 2017
Accepted: 16 October 2017
Published: 23 October 2017
Even if rainfall and temperature are factors classically associated to malaria, little is known about other meteorological factors, their variability and combinations related to malaria, in association with river height variations. Furthermore, in suburban area, urbanization and growing population density should be assessed in relation to these environmental factors. The aim of this study was to assess the impact of combined environmental, meteorological and hydrological factors on malaria incidence through time in the context of urbanization.
Population observational data were prospectively collected. Clinical malaria was defined as the presence of parasites in addition to clinical symptoms. Meteorological and hydrological factors were measured daily. For each factors variation indices were estimated. Urbanization was yearly estimated assessing satellite imaging and field investigations. Principal component analysis was used for dimension reduction and factors combination. Lags between malaria incidences and the main components were assessed by cross-correlation functions. Generalized additive model was used to assess relative impact of different environmental components, taking into account lags, and modelling non-linear relationships. Change-point analysis was used to determine transmission periods within years.
Malaria incidences were dominated by annual periodicity and varied through time without modification of the dynamic, with no impact of the urbanization. The main meteorological factor associated with malaria was a combination of evaporation, humidity and rainfall, with a lag of 3 months. The relationship between combined temperature factors showed a linear impact until reaching high temperatures limiting malaria incidence, with a lag 3.25 months. Height and variation of the river were related to malaria incidence (respectively 6 week lag and no lag).
The study emphasizes no decreasing trend of malaria incidence despite accurate access to care and control strategies in accordance to international recommendations. Furthermore, no decreasing trend was showed despite the urbanization of the area. Malaria transmission remain increase 3 months after the beginning of the dry season. Addition to evaporation versus humidity/rainfall, nonlinear relationship for temperature and river height and variations have to be taken into account when implementing malaria control programmes.
In Mali, malaria remains the leading cause of mortality and morbidity representing 42% of consultations in health centres . The average national prevalence of malaria parasitaemia was 38% in 2010 and 52% in 2012–2013 in children under-5 with 10% in Bamako . In urban and suburban areas where malaria is hypo-endemic with prevalence of less than 10% in children aged from 0 to 9, the burden of clinical malaria is classically low considering socio-economic, development and environmental factors . The proportion of children under 5 years with fever in Bamako was 14.6% [4 ] 2 years before and 8.2% in 2012 . The use of insecticide-treated mosquito nets, the treatment of malaria cases by artemisinin-based combination therapy, and the chemoprevention of malaria in pregnant women are the main malaria control strategies in the area. The biological diagnosis of malaria is free at Sotuba Malaria Research Centre. Malaria cases treatment in children under 5 years and chemoprevention in pregnant women are also free. Insecticide-treated mosquito nets for children under 5 and pregnant women is distributed freely by the National Malaria Control Programme (PNLP). In general, Sotuba has not been affected by large movements of refugees during the study period from 2008 to 2012. Only the last year of the study coincided with the problem in the north of Mali, which did not significantly affect the study area. Urbanization has been reported as a factor which contributes to the decline in risk of malaria .
In Sahelian countries, such as Mali, rainfall is globally weak (200–1300 mm of rain) with variable period of drought. Meteorological factors are known to play an important role in the intra and inter-annual evolution of malaria vectors and incidences [6–9]. Among the numerous available meteorological factors only few of them are classically studied (cumulated rainfall, temperature, humidity). However, other specific factors, such as number of rainfall events, evaporation, temperature at different height, sunlight, wind speed, their variations and combinations are not assessed, nor the association with hydrologic factors, in the context of urbanization.
Malaria transmission is markedly heterogenic in endemo-epidemic areas according to the season and the year . This variability has an impact on control activity planning against malaria. In the context of elimination of malaria, knowledge of malaria risk periods (high transmission periods) is essential for proper use of resources .
Observation data on malaria incidence and burden of clinical malaria are available in many health centres in Mali. The hypothesis was that the burden of clinical malaria cases may change according to environmental and urbanization risk factors. The aim of the study was to assess the impact of urbanization combined with meteorological and hydrological factors on the malaria dynamic in a growing suburban area along the Niger River.
Population observational data were prospectively collected using the outpatient consultation and laboratory registers in the SMRC. From May 2008 to December 2012, all patients residing in Sotuba who consulted at the SMRC for clinical care with malaria compatible signs or symptoms were included in the study. SMRC was the only health centre in the village during the study period.
During outpatient consultations, clinical malaria was defined as the presence of any parasitaemia determined by microscopy and rapid diagnosis test (RDT) (without any changes in the diagnosis strategy during the study period) in addition to clinical variables as fever, headache, vomiting, diarrhoea, abdominal pain, nausea, body-aches, anorexia, convulsion, pallor, cough, myalgia and chills. Fever was defined as axillary temperature over 37.5 °C using an electronic thermometer.
In order to highlight the most important factors related to malaria incidence variations, the following environmental variables were measured daily in the village by the Rural Economic Institute (Institut d’Economie Rurale), located at Sotuba and the National Direction of Hydraulic: maximum inside temperature (t1M); minimum inside temperature (t1m); maximum ground temperature (t2M); minimum ground temperature (t2m); temperature at 10 cm (t10); temperature at 20 cm (t20); outside extreme humidity (hx); relative humidity (hr); evaporation using an experimental plug (evg); evaporation using an experimental pan (evn); cumulated rainfall (r); number of rainfall events (nre); sunlight (sl); tension value (tv); wind speed (ws); river height (rvh). The population density variation was estimated by counting buildings using satellite images from Google earth (© 2015 Google Inc. Menlo Park, California) from 2009 to 2012 based on the population census in 2008 (census done by MRTC team). Estimates were further completed by field investigations and interviews of village leaders to assess the number of inhabitants in the new buildings of the area.
The weekly cumulative malaria incidences, cumulated rainfall and rainfall events, and weekly mean, maximum, minimum, standard deviation (SD), range, and coefficient of variation (CV) of all the other meteorological and hydrological factors were estimated.
Meteorological and hydrological factors (19) and their variations were analysed through two principal component analysis (PCA). This method was used in order to combine factors, taking into account collinearity and for dimension reduction. The main components, determined by eigenvalues, were then characterized, for meteorological factors and hydrological factors.
Lags between each component and malaria incidence were assessed by cross-correlation function estimations, after modelling time processes for each time series, following the Box–Jenkins approach . Each component showing a significant cross-correlation was then selected and shifted (according to the estimated lag) for the multivariate analysis.
- IV.For the multivariate analysis, a regressive approach was developed, following a modelling strategy determined in order to take into account the different situations that may bias classical regression models:
Taking into account over-dispersion in the number of malaria cases, a quasi-Poisson distribution was used, replacing a classical Poisson distribution for count data;
In order to interpret the parameters (after exponential transformations) as incidence ratios, the population was used as an offset of the regression model after log transformation;
Assessing non-linear relationships between malaria cases and each cofactors (the shifted components selected at the first step), spline functions were used;
Taking into account collinearities between the different cofactors, we did use the main components of the previous PCAs (first step), which were, by definition, independent.
A stepwise approach was also used based on the generalized Cross-validation criterion (GCV).
These modelling strategy led us to the use of a general additive model (GAM), following Simon Wood recommendations .
Population socio demographic characteristics
The yearly rainfall varied from 1 year to another (Fig. 2). The highest yearly cumulative rainfall was seen in 2010 (1186.4 mm) and the lowest in 2011 (867 mm). The maximum (Mt1M) and minimum (mt1m) inside temperatures were similar throughout years, but varied intra-annually. The mean river height (rvh) changed annually (annual peak) and the dynamic remain similar from year to year. The years 2009 and 2010 showed the highest malaria incidence peaks, together with the highest population growth rates and the highest rainfall peaks.
PCA of malaria risk factors
River profiles (Fig. 3 B1, B2): For river factors the first three primary components explained 99.84% of the inertia. The first component (Riv1) was constituted by the height of river and the second and the third components were constituted by the variations of river height (Riv2 and Riv3).
Determination of PCA axes elapsing times by ARIMA of series according Box–Jenkins approach
According to the hydrological factors, the river height component (Riv1) was significantly positively cross-correlated with malaria incidences (0.30) with a lag of 6 weeks (Fig. 4). The second component (first river variation component, Riv2) was significantly negatively cross-correlated with no lag (− 0.12), and the third component (second river variation component, Riv3) was significantly positively cross-correlated with a lag of 7 weeks (0.17) (Fig. 4).
Modelling malaria incidences (GAM)
General additive model results for univariate and multivariate analysis
Change point analysis determining transmission periods
The change point analysis clearly showed the classical pattern of seasonal malaria, with low and high transmission periods (Fig. 1). Low transmission periods were not totally associated with dry seasons, and high transmission periods were not totally associated with rainy seasons, but less shifted with river height (Fig. 1).
Clinical malaria incidence and potential environmental risk factors, including temporal variations were investigated in suburban area of Bamako, Mali over 5 years from 2008 to 2012. Despite good access to care, diagnosis and treatment procedure in accordance to international recommendation, and the presence of a well-trained team, the temporal pattern of clinical malaria incidence did not change over years (Fig. 1). In particular, no reduction in the incidence of malaria was observed despite increasing urbanization compared to a phenomenon commonly reported in the literature that urbanization decrease the risk of malaria [5, 18]. The intra annual [12, 19–21] and inter annual [12, 20] variations of malaria are usual, with a pattern already described in many transmission settings. A lag between low transmission periods of malaria and dry seasons was observed, as well as between high transmission periods and rainy seasons. Conversely, a better concordance and a shorter lag were found between malaria incidence and river height patterns (6 weeks lag, Fig. 1).
The main meteorological factor associated with malaria is a combination of evaporation, humidity and rainfall (Fig. 5A1) described by several studies before. This linear correlation occurred with a lag of about 3 months. Even if in other area the lag between rainfall and malaria cases could be different, rainfall is not the only factor related to malaria transmission. Indeed, humidity, temperature, but also land cover/land use are important factors. As malaria transmission is known to be heterogeneous, according to geographical/environmental contexts, the lag difference between different areas is not surprising. Furthermore, Niakhar is a small city, and the data assessed by Ndiaye et al. aggregated informations from the whole rural area, showing millet and peanuts fields . But Sotuba is a sub-urban area, very close to Bamako, the capital city. Note that, the SMRC team did not report any refugees camp near the Sotuba area during the study period.
The relationship between river high and rainfall at the same geographical area is not direct. Indeed, the Niger river high, in Bamako, is mainly due to rainfall in Conakry Guinea, at the source of the Niger river, and less related to the rainfall at Bamako. Furthermore, it is not the river high that is related to malaria incidence but also river variations, which could create, or destroy, breeding sites. This has also been observed in Bandiagara, Mali .
Most of the studies in the literature found that rainfall and humidity increase the risk of malaria by developing suitable breeding sites and increasing density of mosquitoes [9, 10, 20–28]. However, the positive association between humidity/rainfall and malaria incidence is inconsistent with other authors . Furthermore, the length of the lag between rainfall and malaria incidence should be taken into account by prevention policies, including traveler recommendations on prophylaxis, malaria risk is still important, 3 months after rain stops.
The relationship between combined temperature factors was not linear (Fig. 5A2). This association was positively correlated and when temperature attempted a threshold where malaria incidence was constant, and then slightly decreased with temperature . The lag time was 3.25 months between combined temperature factors and the incidence of malaria. This nonlinear association is concordant with malaria vector and parasite development cycles. This type of relationship is not taken into account in the major part of the literature. It can be explained by the existence of a temperature threshold fatal to the development of malaria vector (malaria transmission is optimal at 25 °C and decrease drastically at > 28 °C) [23, 30–32]. The main vector species are Anopheles gambiae s.l. and Anopheles funestus, with predominance of An. gambiae, which accounted for 99.7% of transmission . Other studies found a positive linear association between malaria incidence and maximum temperature [9, 24, 27, 28, 33]. However, some studies observed negative association . This apparent contradiction may be explained by the difference between the transmission periods, statistical methods (not exploring nonlinear relationships) and the diversity of the risk factors combined for the analysis.
Height of the Niger River was positively linearly correlated to malaria incidence (Fig. 5B1 and B2) with a lag of about 1.5 months. The variation of the river height plays an important role in the creation of anopheles breeding sites. When the river is full, the breeding sites on the banks of the river are completely submerged with water. However, the periods of flood recession are generally a source of breeding sites along the river. Humidity of the area could be driven also by river variation, and have an effect on mosquito survival . Publications focused much more on the presence or absence of vectors or parasites carriers in the areas near the rivers [35, 36].
In conclusion, the present study emphasizes temporal variation of malaria incidence without change in the dynamics through the time. The lag between rainfall and malaria is classically acknowledged. These has to be considered for the implementation of malaria control strategies such as seasonal malaria chemoprevention, mass distribution of insecticide impregnated bed nets and mass sensitization campaigns to cover 3 months after rain stops. Nonlinear relationships for combined meteorological factors and other source of breeding sites (such as river) have to be taken into account when developing environmental based control strategies.
Furthermore, the persistence of clinical malaria during the low transmission period should be assessed, determining if different symptoms, clinical profiles and parasitaemia are associated with the different transmission periods.
MSS participated in the design of the study, coordinate data collection and data management, undertook the statistical analysis, interpretation of results and drafted the manuscript. KS, BK, YS and SG have collected data and assisted in interpretation of results. AD, ND and RP assisted in interpretation of results and was involved in the critical revision of the manuscript. OKD participated in the design of the study, assisted in interpretation of results and critical revision of the manuscript. JG participated in the design of the study, undertook the statistical analysis, interpretation of results and drafted the manuscript. MY collected hydrologic data and assisted in interpretation of results. All authors read and approved the final manuscript.
Many thanks to the population of Sotuba. We would like to thank also the Institut d’Economie Rurale (IER), Bamako, Mali for helping us to extract the meteorological data. We thanks the National Direction of Hydraulic in Bamako, Mali [Direction Nationale de l’Hydraulique Bamako (MALI)] for collecting River data. Thanks also to UMI 3189 and REDOH/grant challenge. Furthermore, we would like to thanks the MRTC parasitology, clinical and data management teams in Sotuba.
The authors declare that they have no competing interests.
This work was supported by the AMMA Consortium (African Monsoon Multidisciplinary Analysis) and by the Prospective & Cooperation NGO. The Sotuba Malaria Research Centre, was built by NIAID/NIH Grants and US-Embassy Bamako. The centre has been recently strengthen by the MSP3 Grants and WANECAM Grant. The field activities are supported by the UMI3189 Grant, the IRBA/MRTC Grant and FMERIEUX Grants. The funding partners did not interfere in the design, analysis and interpretation of this work.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Ministère de la Santé. Annuaire Système Local d’Information Sanitaire (SLIS). Mali. 2012.Google Scholar
- Cellule de Planification et de Statistique (CPS/SSDSPF), Institut National de la Statistique (INSTAT/MPATP), INFO-STAT et ICF International. Enquête Démographique et de Santé au Mali 2012–2013. Rockville: CPS, INSTAT, INFO-STAT et ICF International; 2014.Google Scholar
- Ceesay SJ, Bojang KA, Nwakanma D, Conway DJ, Koita OA, Doumbia SO, et al. Sahel, savana, riverine and urban malaria in West Africa: similar control policies with different outcomes. Acta Trop. 2012;121:166–74.View ArticlePubMedGoogle Scholar
- Cellule de Planification et de Statistique du Ministère de la Santé (CPS/MS), Direction Nationale de la Statistique et de l’Informatique du Ministère de l’Economie, de l’Industrie et du Commerce (DNSI/MEIC) et Macro International Inc. Enquête Démographique et de Santé du Mali 2006. Calverton: CPS/DNSI et Macro International Inc; 2007.Google Scholar
- Tatem AJ, Gething PW, Smith DL, Hay SI. Urbanization and the global malaria recession. Malar J. 2013;12:133.View ArticlePubMedPubMed CentralGoogle Scholar
- Breman JG, Alilio MS, Mills A. Conquering the intolerable burden of malaria: what’s new, what’s needed: a summary. Am J Trop Med Hyg. 2004;71(Suppl 2):1–15.PubMedGoogle Scholar
- Al-Mansoob MA, Al-Mazzah MM. The role of climate on malaria incidence rate in four governorates of Yemen. Med J Malays. 2005;60:349–57.Google Scholar
- Chaves LF, Koenraadt CJ. Climate change and highland malaria: fresh air for a hot debate. Q Rev Biol. 2010;85:27–55.View ArticlePubMedGoogle Scholar
- Reid HL, Haque U, Roy S, Islam N, Clements AC. Characterizing the spatial and temporal variation of malaria incidence in Bangladesh, 2007. Malar J. 2012;11:170.View ArticlePubMedPubMed CentralGoogle Scholar
- Coulibaly D, Rebaudet S, Travassos M, Tolo Y, Laurens M, Kone AK, et al. Spatio-temporal analysis of malaria within a transmission season in Bandiagara, Mali. Malar J. 2013;12:82.View ArticlePubMedPubMed CentralGoogle Scholar
- Hay SI, Snow RW. The malaria atlas project: developing global maps of malaria risk. PLoS Med. 2006;3:e473.View ArticlePubMedPubMed CentralGoogle Scholar
- Dicko A, Sagara I, Diemert D, Sogoba M, Niambele MB, Dao A, et al. Year-to-year variation in the age-specific incidence of clinical malaria in two potential vaccine testing sites in Mali with different levels of malaria transmission intensity. Am J Trop Med Hyg. 2007;77:1028–33.PubMedGoogle Scholar
- Shumway RH, Stoffer DS. Time series analysis and its applications. 3rd ed. Springer-Verlag: New York; 2011: XI, 596. p. 202.Google Scholar
- Wood SN. Generalized additive models: an introduction with R. London: Chapman and Hall/CRC; 2006.Google Scholar
- Killick R, Fearnhead P, Eckley IA. Optimal detection of changepoints with a linear computational cost. JASA. 2012;107:1590–8.View ArticleGoogle Scholar
- Chen J, Gupta AK. Parametric statistical change point analysis. 2nd ed. Birkhauser Basel; 2012: XIII, 273.Google Scholar
- Lê S, Josse J, Husson F. FactoMineR: an R package for multivariate analysis. J Stat Softw. 2008;25:1–18.View ArticleGoogle Scholar
- Trape JF, Zoulani A. Malaria and urbanization in Central Africa: the example of Brazzaville. Part III: relationships between urbanization and the intensity of malaria transmission. Trans R Soc Trop Med. 1987;81(Suppl 2):19–25.View ArticleGoogle Scholar
- Yé Y, Hoshen M, Kyobutungi C, Louis VR, Sauerborn R. Local scale prediction of Plasmodium falciparum malaria transmission in an endemic region using temperature and rainfall. Glob Health Action. 2009. doi:10.3402/gha.v2i0.1923.Google Scholar
- Stefani A, Hanf M, Nacher M, Girod R, Carme B. Environmental, entomological, socioeconomic and behavioural risk factors for malaria attacks in Amerindian children of Camopi, French Guiana. Malar J. 2011;10:246.View ArticlePubMedPubMed CentralGoogle Scholar
- Donovan C, Siadat B, Frimpong J. Seasonal and socio-economic variations in clinical and self-reported malaria in Accra, Ghana: evidence from facility data and a community survey. Ghana Med J. 2012;46(2):85–94.PubMedPubMed CentralGoogle Scholar
- Ndiaye O, Le Hesran JY, Etard JF, Diallo A, Simondon F, Neil Ward M, et al. Variations climatiques et mortalité attribuée au paludisme dans la zone de Niakhar, Senegal, de 1984 à 1996. Cahiers d’études et de recherches francophones/Santé. 2001;11:25–33.Google Scholar
- Mordecai EA, Paaijmans KP, Johnson LR, Balzer C, Ben-Horin T, de Moor E, et al. Optimal temperature for malaria transmission is dramatically lower than previously predicted. Ecol Lett. 2013;16:22–30.View ArticlePubMedGoogle Scholar
- Orlando PZ, Mikael A. Spatial and temporal patterns of malaria incidence in Mozambique. Malar J. 2011;10:189.View ArticleGoogle Scholar
- Gillioli G, Mariani L. Sensitivity of Anopheles gambiae population dynamics to meteo-hydrological variability: a mechanistic approach. Malar J. 2011;10:294.View ArticleGoogle Scholar
- Chaves LF, Satake A, Hashizume M, Minakawa N. Indian Ocean Dipole and rainfall drive a Moran effect in East Africa malaria transmission. J Infect Dis. 2012;205:1885–91.View ArticlePubMedGoogle Scholar
- Zhang Y, Liu QY, Luan RS, Liu XB, Zhou GC, Jiang JY, et al. Spatial-temporal analysis of malaria and the effect of environmental factors on its incidence in Yongcheng, China, 2006–2010. BMC Public Health. 2012;12:544.View ArticlePubMedPubMed CentralGoogle Scholar
- Fang H, Shuisen Z, Shaosen Z, Hongju W, Linhua T. Temporal correlation analysis between malaria and meteorological factors in Motuo County, Tibet. Malar J. 2011;10:54.View ArticleGoogle Scholar
- Ali A, Monica CJ, Cesar K. Modelling the effects of weather and climate on malaria distributions in West Africa. Malar J. 2014;13:126.View ArticleGoogle Scholar
- Pascual M, Ahumada JA, Chaves LF, Rodo X, Bouma M. Malaria resurgence in the East African highlands: temperature trends revisited. Proc Natl Acad Sci USA. 2006;103:5829–34.View ArticlePubMedPubMed CentralGoogle Scholar
- Gosoniu L, Veta AM, Vounatsou P. Bayesian geostatistical modeling of malaria indicator survey data in Angola. PLoS ONE. 2010;5:e9322.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhou SS, Huang F, Wang JJ, Zhang SS, Su YP, Tang LH. Geographical, meteorological and vectorial factors related to malaria re-emergence in Huang-Huai River of central China. Malar J. 2010;9:337.View ArticlePubMedPubMed CentralGoogle Scholar
- Yé Y, Louis VR, Simboro S, Sauerborn R. Effect of meteorological factors on clinical malaria risk among children: an assessment using village-based meteorological stations and community-based parasitological survey. BMC Public Health. 2007;7:101.View ArticlePubMedPubMed CentralGoogle Scholar
- Yamana TK, Eltahir EAB. Incorporating the effects of humidity in a mechanistic model of Anopheles gambiae mosquito population dynamics in the Sahel region of Africa. Parasite Vectors. 2013;6:235.View ArticleGoogle Scholar
- Labbo R, Fandeur T, Jeanne I, Czeher C, Williams E, Arzika I, et al. Ecology of urban malaria vectors in Niamey, Republic of Niger. Malar J. 2016;15:314.View ArticlePubMedPubMed CentralGoogle Scholar
- Sissoko MS, van den Hoogen LL, Samake Y, Tapily A, Diarra AZ, Coulibaly M, et al. Spatial patterns of Plasmodium falciparum clinical incidence, asymptomatic parasite carriage and Anopheles density in two villages in Mali. Am J Trop Med Hyg. 2015;93:790–7.View ArticlePubMedPubMed CentralGoogle Scholar