- Open Access
Effect of climatic variability on malaria trends in Baringo County, Kenya
© The Author(s) 2017
- Received: 26 October 2016
- Accepted: 2 May 2017
- Published: 25 May 2017
Malaria transmission in arid and semi-arid regions of Kenya such as Baringo County, is seasonal and often influenced by climatic factors. Unravelling the relationship between climate variables and malaria transmission dynamics is therefore instrumental in developing effective malaria control strategies. The main aim of this study was to describe the effects of variability of rainfall, maximum temperature and vegetation indices on seasonal trends of malaria in selected health facilities within Baringo County, Kenya.
Climate variables sourced from the International Research Institute (IRI)/Lamont-Doherty Earth Observatory (LDEO) climate database and malaria cases reported in 10 health facilities spread across four ecological zones (riverine, lowland, mid-altitude and highland) between 2004 and 2014 were subjected to a time series analysis. A negative binomial regression model with lagged climate variables was used to model long-term monthly malaria cases. The seasonal Mann–Kendall trend test was then used to detect overall monotonic trends in malaria cases.
Malaria cases increased significantly in the highland and midland zones over the study period. Changes in malaria prevalence corresponded to variations in rainfall and maximum temperature. Rainfall at a time lag of 2 months resulted in an increase in malaria transmission across the four zones while an increase in temperature at time lags of 0 and 1 month resulted in an increase in malaria cases in the riverine and highland zones, respectively.
Given the existence of a time lag between climatic variables more so rainfall and peak malaria transmission, appropriate control measures can be initiated at the onset of short and after long rains seasons.
- Baringo County
- Malaria transmission
- Seasonal trends
Malaria is a global health problem that causes an estimated 438,000 deaths annually; 88% of which occur in the sub-Saharan Africa . Seventy-five percent of the malaria clinical episodes worldwide occur in Africa with a corresponding high public health burden . Up to 35.4 million disability adjusted life years (DALYs) are lost in the sub-Saharan Africa region alone due to malaria mortality and morbidity .
In Kenya, malaria is among the leading causes of morbidity and mortality and is responsible for almost half of all outpatient attendance and 20% of all admissions to health facilities . Pregnant women and children under five years old are most vulnerable to malaria infections  with an estimated 170 million working days being lost to malaria in Kenya each year .
The high burden of malaria in Kenya and the larger sub-Saharan Africa region may be associated with a number of factors among them climatic and environmental . Given that malaria is vector-transmitted, with a complex life cycle in both the mosquito and human, transmission and patterns of malaria infection are dependent on both environmental and climatic factors . A study by Githeko et al.  showed that inter-annual and inter-decadal climate variability influences the epidemiology of vector-borne diseases directly, while temperature and rainfall have long been known to influence seasonal and inter-annual variability of malaria .
The effect of temperature on the life history traits of mosquitoes and malaria transmission has been reported. Temperature can affect the development time of mosquito larvae, the probability of mosquito survival and the development time of malaria parasite (Plasmodium falciparum) in infected mosquitoes either positively or negatively . A rise in temperature to a certain threshold can accelerate the metabolic rate of vectors, increase egg production and increase frequency of blood meals, while temperatures below or above these thresholds can be detrimental to mosquitoes and parasite development . Several mechanistic models concur that effect of temperature on malaria transmission is non-linear, limited to temperatures between 16 and 34 °C with a peak at 25 °C [13–15]. The non-linear temperature sensitivities throughout the mosquito life cycle have a large impact on the adult population dynamics and, therefore, on the mosquitoes’ ability to act effectively as malaria vectors.
Rainfall influences vector longevity indirectly by creating wet conditions that favour vector breeding. This in turn influences the geographical range and seasonal variability of disease vectors . The relationship between malaria incidence and rainfall is non-linear, implying that an increase in precipitation would not necessarily increase malaria cases . Moderate rainfall has a positive effect on mosquito abundance, while intense precipitation can wash away mosquito breeding sites, and therefore reduce malaria transmission shortly following heavy rains . The Normalized Difference Vegetation Index (NDVI) is a spectral measure of amount, relative greenness, phenological characteristics and productivity of vegetation . It is defined as the difference between the visible (RED) and near-infrared (NIR) bands over their sum, (NIR−RED)/(NIR + RED). It is a robust indicator of vegetation condition which allows valid comparisons of seasonal and inter-annual variations in vegetation growth and activity . In the study area, the seasonal NDVI variations are linked to rainfall. NDVI values range between −1 to +1 An NDVI value of zero means no green vegetation and close to +1 (0.8–0.9) indicates the highest possible density of green leaves. NDVI can be used as a surrogate for precipitation based on their close correlation . The capability of NDVI time-series to monitor and predict vector-borne diseases depends on the correlation between disease incidence, vegetation greenness and precipitation .
The nature of vector biological processes and the degree to which the vectors depend on environmental and climatic factors makes malaria transmission somewhat region specific . In Kenya, there are four epidemiological zones whose diversity in malaria transmission and risk is determined by altitude, rainfall patterns and temperature. The zones include: the endemic Lake Victoria and coastal regions, epidemic-prone Western highlands, seasonal transmission arid and semi-arid areas, and low risk central highlands . Parts of Baringo County being semi-arid experience seasonal malaria transmission , while the presence of numerous seasonal and permanent water bodies provides suitable breeding microhabitats for malaria vectors at certain times of the year.
This study was thus necessitated by the lack of information on the interactions of climatic and environmental factors, and their role in driving the transmission and prevalence of malaria across the different ecological zones of Baringo County, which has been hampering the planning of intervention strategies against malaria. The study modelled the effect of climatic variations on the prevalence and long-term trend of malaria so as to identify the seasonal climatic drivers of malaria transmission in the different ecological zones of Baringo County.
Selection of health facilities
Ten health facilities were selected across the four zones based on availability of health records and the catchment population served; with those serving larger populations being selected. The facilities included: Kabarnet County Hospital, Kituro Health Centre and Kabartonjo sub-county Hospital in the highlands; Sabor Dispensary and Kipcherere Health Centre in the mid-altitudes; Marigat sub-county Hospital and Kampi ya Samaki Health Centre in the lowlands; and Barwessa, Salawa and Keturwo Health Centres in the riverine zone. The available datasets varied per region. The highland and mid-altitude zones had data spanning for a period of 10 years (2004–2013) while lowland (2005–2013) and riverine (2006–2014) zone had data spanning 9 years.
Retrospective data extraction
Retrospective health records of clinically-diagnosed and treated malaria cases between 2004 and 2014 were extracted from the registers in the 10 health facilities. Daily counts of malaria cases among all age groups and gender were entered in MS Excel and computed into monthly data sets for analyses. Only two out of the 10 health facilities that were studied had some missing data, i.e. Marigat sub-county Hospital which had 13% (14/108 months) of the data missing and Keturwo Health Centre which had 8.3% (9/108 months) of the data missing. The missing data were imputed on aggregated monthly data sets using predictive mean matching method using mice package in R 3.0.3 statistical software .
Imputation of missing data
The unavailable monthly data was assumed to be missing at random (MAR) before multiple imputation was conducted. Predictive mean matching (PMM) method was used to enable imputation of the missing values based on the observed values. Through this method the missing values were imputed by means of nearest neighbour values with distances based on the expected values of the missing variables . Unlike other imputation methods, the PPM method produces acceptable estimates and preserves the underlying distribution of the observed data especially for quantitative variables that are not normally distributed. Compared with other methods, PMM produces imputed values that are much more like real values. The major pitfall with PMM however is that there is no mathematical theory to justify it and only a handful studies have evaluated its performance, meaning that it is still not clear how well it compares with alternative methods .
Climate and environmental data
Monthly average rainfall, inferred maximum air temperatures and enhanced vegetation indices over the study period were sourced from the International Research Institutes of Climate and Society’s database . The monthly average rainfall data used was obtained from University of California Santa Barbara (UCSB) Climate Hazards Group InfraRed Precipitation with Station Data (CHIRPS) v2p0 . Annual precipitation was averaged from monthly estimates for each year. Inferred maximum air temperature (Tmax) and minimum land surface temperature (Tmin) used as proxy for minimum and maximum air temperatures were obtained from the United States Geological Survey (USGS) LandDAAC MODIS 1 km 8 day version 005 datasets [34, 35]. Sixteen-day MODIS Enhanced Vegetation Index (EVI) composites from MODIS-Terra MOD13Q1 at 250 m spatial resolution  were used to derive 15 years temporal profiles.
For all data sets, the spatial averages limited to the spatial extent of the ecological zones were downloaded in expert mode. The 8- and 16-day composites were averaged into monthly means for each year.
Exploratory data analysis was used to visualize seasonal patterns of climate variables in relation to malaria cases across the four zones. A smoothing line was added to the patterns to get good visual information. Locally weighted regression was used to smoothen the data points with a smoother span of 0.67, tricube as weighted function, number of iterations for robust fitting 3 and an order of the polynomial of 1.
The Seasonal Mann–Kendall trend test was used to detect malaria monthly trends over the study period. Pair wise comparison of monthly means of malaria cases was carried out in order to draw simultaneous inference about the dominant malaria seasons using Tukey multiple comparison procedure. Additive decomposition of malaria cases and climatic variables were conducted in order to estimate the trend component. The trend was determined using moving averages as the smoothing method. The sample cross correlation function was used to identify lags of climate variables that were useful predictors of malaria cases. Dominant cross correlations between malaria cases and climate variables were selected and included in the regression model. A negative binomial regression model with lagged climate variables was used to model the monthly malaria cases. Model fit was assessed by checking the autocorrelation function and partial autocorrelation function of the model residuals. Variance inflation factors (VIF) were computed to check multi-collinearity. The analyses were done using dyn and trend packages in R [36, 37].
Ethical approval to access hospital registers for data extraction was obtained from the Kenyatta National Hospital/University of Nairobi Ethical Review Committee (P70/02/2013) and also from the Ministry of Health (Ref.: CNTY/GEN/Vol.1/83) and Department of Medical Services (Ref: BCG/CDH/GEN/VOL.II/2015), Baringo County.
Long-term mean monthly malaria cases (2004–2014)
Long-term malaria cases (2004–2014)
Highest cases (month, year)
Lowest cases (month, year)
1799 (Oct 2012)
285 (May 2009)
197 (Oct 2009)
20 (Apr 2006)
1318 (Aug 2006)
126 (Jan 2005)
979 (Nov 2009)
110 (May 2014)
Yearly cumulative malaria cases (2004–2014) and trends
Malaria seasonality over 10 year period
Long-term trends of malaria cases against climatic variables
Identification of possible climatic and environmental predictors of malaria cases
Rainfall, temperature and EVI lags in relation to malaria cases in the four zones
Rainfall (lag 1)
Rainfall (lag 2)
Temperature (lag 1)
EVI (lag 0)
EVI (lag 1)
EVI (lag 2)
Rainfall (lag 1)
Rainfall (lag 2)
Temperature (lag 1)
EVI (lag 0)
Rainfall (lag 0)
Rainfall (lag 1)
Rainfall (lag 2)
Temperature (lag 0)
EVI (lag 0)
Rainfall (lag 1)
Rainfall (lag 2)
Temperature (lag 0)
Temperature (lag 1)
EVI (lag 1)
EVI (lag 2)
Climatic factors are considered important in the spatial and temporal distribution of vector borne diseases as they determine vector distribution, and influence inter-annual variability, epidemics and long-term trends . There is a strong discernible link between malaria outbreaks, temperature  and rainfall . In the current study, malaria cases generally increased in highland and mid-altitude zones but decreased in the riverine and lowland zones from the year 2011 onwards during the study period (Fig. 2). Two malaria peak seasons were identified in the lowland zone while three malaria peak seasons were identified in the other zones, largely following climatic seasons in the study area. Statistically significant differences in monthly malaria peaks was recorded in the highlands and mid-altitude zones suggesting seasonal malaria transmission. However there was no statistical significance in malaria peaks in the lowland and riverine zones, suggesting that malaria transmission in these two zones is perennial rather than seasonal as previously thought .
The mean monthly rainfall for this study showed positive significant correlation with malaria cases at 2 months lag across all zones while the mean maximum temperature showed positive significant correlation with malaria cases in two zones, the highlands and the riverine zones. Previous studies examining the link between climate and malaria established lagged associations between climate variables (temperature and rainfall) and malaria cases over time periods ranging from weeks to months [10, 18, 39–42]. These studies attributed the lags to the creation of mosquito breeding habitats, the time required by mosquitoes to develop to adulthood, acquire and transmit malarial infection, and for symptoms to arise in the human host as the most probable cause.
According to Confalonieri et al. , periods of unusually high rainfall, altered humidity or warmer temperatures can result in modified distribution and duration of malaria, as well as increased transmission; even in areas where control is strong. Consistent with the current findings, Small et al.  cited precipitation and temperature as key drivers of malaria case variations across Africa, while acknowledging the complexities of some climatic factors.
The difference in environmental relationship to malaria cases across the zones is attributed to variations in environmental factors between the zones. The mid-altitude zone has no rivers or water bodies and rainfall is therefore the only source of surface water that serves as breeding points for malaria vectors. Although the other zones have permanent water bodies in the form of lakes, rivers, swamps, dams and water pans, rainfall still contributes to malaria cases through creation of additional seasonal breeding sites for malaria vectors.
Consistent with our study findings, Chaves et al.  cited increased microhabitats resulting from relative humidity caused by moderate rainfall. These conditions increase the longevity of adult mosquitoes by prolonging vector life span. Paaijmans et al.  highlighted the complex interrelationship between precipitation and vectors, noting that drought may eliminate mosquito habitats, while floods could create isolated pools suitable for vector breeding. The relatively low annual rainfall in the mid-altitude zone and the general absence of permanent water bodies contributed to the observed low but varying numbers of recorded malaria cases; possibly due to the varying climatic conditions. All in all, there is a general consensus that rainfall can influence malaria transmission either positively by creating suitable habitats or negatively by flushing breeding sites depending on its intensity [46, 47].
Temperature plays a key role in malaria transmission by influencing vector and parasite life cycles. Studies have highlighted the biological amplification nature of temperature on mosquitoes . This study showed that the mean maximum temperatures within the four zones varied. While the mean maximum temperature significantly influenced malaria cases at lag 0 in the riverine zone and lag 1 in the highlands, it was non-significant in the mid-altitude and lowland zones. The difference in the contribution of maximum temperature to malaria cases between zones is attributed to the differences in prevailing temperatures in the four zones. Being colder, temperature was probably the limiting factor in malaria vector development in the highland and riverine zones; hence a rise in the maximum temperature increased vector and parasite development rates . Since temperature influences the development and survival rates of both vectors and parasites, malaria transmission rates tend to increase with increasing temperature but up to a given threshold .
Craig et al.  put the optimal temperatures for malaria transmission at between 22 and 32 °C, while Bi et al.  reported temperatures of between 20 and 30 °C as being optimal for Anopheles survival and that temperatures below 16 °C and above 30 °C have a negative impact on mosquitoes survival. Chikodzi  noted that temperatures above 32 °C can cause high vector population turnover, with thermal death for mosquitoes expected to occur around 41–42 °C.
Vegetation index often acts as a surrogate for precipitation and surface temperatures and has been correlated to vector borne diseases . In this study, vegetation cover followed a positive trend with the amount of precipitation received. In this study, EVI did not play any significant role in malaria transmission across the four zones.
This study established seasonality in malaria transmission over the study period (2004–2014) in the highland and mid-altitude zone. Malaria transmission in the lowland and riverine zone was shown to be perennial. Peak malaria cases followed increased rainfall with a time lag of 2 months across the study area and increased maximum temperatures with a time lag of 0 and 1 months in the riverine and highland zones respectively. The observed time lags between peak malaria cases and climatic variables are particularly important in forecasting malaria outbreak using local weather data. Therefore, monitoring rainfall and temperature trends and early recognition of anomalies in weather patterns can provide a fairly accurate forecast of transmission risk within Baringo County, and hence inform timely action including vector control measures.
BE, IK were study PI and Co-PI respectively and conceptualized and supervised the work. EK, AO, DA, MM, DO, EM developed the tools and collected the data. EK, AO, DA and MM analysed the data and drafted the manuscript. BE and IK reviewed the manuscript and authorized submission for publication. All authors read and approved the final manuscript.
This research was carried out within a larger 3-year study on early warning systems for improved human health and resilience to climate sensitive vector borne diseases in Kenya through a Grant by WHO/TDR/IDRC Project ID B20278. The authors acknowledge the contributions of the entire project team, the Baringo County Ministry of Health, the County Director of Health, the Public Administration and community members without whose support this research would not have been conducted.
The authors declare that they have no competing interests.
This research received both national and the World Health Organization (WHO) ethical clearance reference P70/02/2013 and Protocol ID B20278, respectively.
This study received financial assistance from the WHO’s Special Programme for Research and Training in Tropical Diseases (TDR) through a grant agreement with the International Development Research Centre of Canada (106905–00).
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.
- WHO. World malaria report 2013. Geneva: World health Organization; 2013.Google Scholar
- Brooker S, Akhwale W, Pullan R, Estambale B, Clarke SE, Snow RW, et al. Epidemiology of plasmodium-helminth co-infection in Africa: populations at risk, potential impact on anemia, and prospects for combining control. Am J Trop Med Hyg. 2007;77:88–98.PubMedPubMed CentralGoogle Scholar
- Fernando SD, Rodrigo C, Rajapakse S. The ‘hidden’ burden of malaria: cognitive impairment following infection. Malar J. 2010;9:366.View ArticlePubMedPubMed CentralGoogle Scholar
- Ministry of Health. Towards a malaria free Kenya: Kenya revised malaria strategy 2009–2018. National Malaria Control Programme; 2014.Google Scholar
- Ministry of Health. National guidelines for the diagnosis, treatment and prevention of malaria for health workers in Kenya. Kenya: Div Malar Control, MoH; 2006.Google Scholar
- Ministry of Health. National malaria strategy 2041–2010. Kenya: Div Malar Control, MoH; 2001.Google Scholar
- Tonui WK, Otor SCJ, Kabiru EW, Kiplagat WK. Patterns and trends of malaria morbidity in western highlands of Kenya. Int J Educ Res. 2013;1:8.View ArticleGoogle Scholar
- Guthmann JP, Llanos-Cuentas A, Palacios A, Hall AJ. Environmental factors as determinants of malaria risk. A descriptive study on the northern coast of Peru. Trop Med Int Health. 2002;7:518–25.View ArticlePubMedGoogle Scholar
- Githeko AK, Lindsay SW, Confalonieri UE, Patz JA. Climate change and vector-borne diseases: a regional analysis. Bull World Health Organ. 2000;78:1136–47.PubMedPubMed CentralGoogle Scholar
- Zhou G, Minakawa N, Githeko AK, Yan G. Association between climate variability and malaria epidemics in the East African highlands. Proc Natl Acad Sci USA. 2004;101:2375–80.View ArticlePubMedPubMed CentralGoogle Scholar
- Bayoh MN, Lindsay SW. Effect of temperature on the development of the aquatic stages of Anopheles gambiae sensu stricto (Diptera: Culicidae). Bull Entomol Res. 2003;93:375–81.View ArticlePubMedGoogle Scholar
- Mellor PS, Leake CJ. Climatic and geographic influences on arboviral infections and vectors. Rev Sci Tech. 2000;19:41–54.View ArticlePubMedGoogle Scholar
- Mordecai EA, Paaijmans KP, Johnson LR, Balzer C, Ben-Horin T, de Moor v, et al. Optimal temperature for malaria transmission is dramatically lower than previously predicted. Ecol Lett. 2013;16:22–30.View ArticlePubMedGoogle Scholar
- Craig MH, Snow RW, le Sueur D. A climate-based distribution model of malaria transmission in sub-Saharan Africa. Parasitol Today. 1999;15:105–11.View ArticlePubMedGoogle Scholar
- Parham PE, Michael E. Modeling the effects of weather and climate change on malaria transmission. Environ Health Perspect. 2010;118:620–6.View ArticlePubMedGoogle Scholar
- Patz JA, Strzepek K, Lele S, Hedden M, Greene S, Noden B, et al. Predicting key malaria transmission factors, biting and entomological inoculation rates, using modelled soil moisture in Kenya. Trop Med Int Health. 1998;3:818–27.View ArticlePubMedGoogle Scholar
- Huang F, Zhou S, Zhang S, Wang H, Tang L. Temporal correlation analysis between malaria and meteorological factors in Motuo County. Tibet Malar J. 2011;10:54.View ArticlePubMedGoogle Scholar
- Briet OJ, Vounatsou P, Amerasinghe PH. Malaria seasonality and rainfall seasonality in Sri Lanka are correlated in space. Geospat Health. 2008;2:183–90.View ArticlePubMedGoogle Scholar
- Evrendilek F, Gulbeyaz O. Deriving vegetation dynamics of natural terrestrial ecosystems from MODIS NDVI/EVI data over Turkey. Sensors. 2008;8:5270–302.View ArticlePubMedPubMed CentralGoogle Scholar
- Pettorelli N, Vik JO, Mysterud A, Gaillard JM, Tucker CJ, Stenseth NC. Using the satellite-derived NDVI to assess ecological responses to environmental change. Trends Ecol Evol. 2005;20:503–10.View ArticlePubMedGoogle Scholar
- Schultz P, Halpert M. Global correlation of temperature, NDVI and precipitation. Adv Space Res. 1993;13:277–80.View ArticleGoogle Scholar
- Wayant NM, Maldonado D, Rojas de Arias A, Cousino B, Goodin DG. Correlation between normalized difference vegetation index and malaria in a subtropical rain forest undergoing rapid anthropogenic alteration. Geospat Health. 2010;4:179–90.View ArticlePubMedGoogle Scholar
- Woube M. Geographical distribution and dramatic increases in incidences of malaria: consequences of the resettlement scheme in Gambela, SW Ethiopia. Indian J Malariol. 1997;34:140–63.PubMedGoogle Scholar
- Division of malaria control. 2010 Kenya malaria indicator survey. Nairobi: Ministry of Public Health and Sanitation; 2011.Google Scholar
- Mohajan HK. Improvement of health sector in Kenya. Am J Public Health Res. 2014;2:159–69.View ArticleGoogle Scholar
- KNBS. The 2009 Kenya population and housing census, vol. IC. In: Population distribution by age, sex and administrative units. 2010.Google Scholar
- Odada EO, Onyando JO, Obudho PA. Lake Baringo: addressing threatened biodiversity and livelihoods. Lakes Reserv Res Manag. 2006;11:287–99.View ArticleGoogle Scholar
- Ochieng AO, Nanyingi M, Kipruto E, Ondiba IM, Amimo FA, Oludhe C, et al. Ecological niche modelling of Rift Valley fever virus vectors in Baringo, Kenya. Infect Ecol Epidemiol. 2016;6:32322.View ArticlePubMedGoogle Scholar
- Buuren S, Groothuis-Oudshoorn K. Mice: multivariate imputation by chained equations in R. J Stat Softw. 2011;45:1–67.View ArticleGoogle Scholar
- Vink G, Frank LE, Pannekoek J, Buuren S. Predictive mean matching imputation of semicontinuous variables. Stat Neerl. 2014;68:61–90.View ArticleGoogle Scholar
- Molenberghs G, Fitzmaurice G, Kenward MG, Tsiatis A, Verbeke G. Handbook of missing data methodology. Boca Raton: CRC Press; 2014.Google Scholar
- International Research Institute for Climate and Society. SOURCES: the IRI/LDEO collection of climate data. http://iridl.ldeo.columbia.edu/SOURCES/. Accessed 25 Nove 2014.
- Funk CC, Peterson PJ, Landsfeld MF, Pedreros DH, Verdin JP, Rowland JD, Romero BE, Husak GJ, Michaelsen JC, Verdin AP. A quasi-global precipitation time series for drought monitoring: US Geological Survey Data Series 832, p. 4. 2014. http://dx.doi.org/10.3133/ds832. Accessed 1 May 2016.
- Vancutsem C, Ceccato P, Dinku T, Connor SJ. Evaluation of MODIS land surface temperature data to estimate air temperature in different ecosystems over Africa. Remote Sens Environ. 2010;114:449–65.View ArticleGoogle Scholar
- Huete A, Didan K, Miura T, Rodriguez EP, Gao X, Ferreira LG. Overview of the radiometric and biophysical performance of the MODIS vegetation indices. Remote Sens Environ. 2002;83:195–213.View ArticleGoogle Scholar
- Grothendieck G. dyn: Time series regression. R package version 0.2-8.1. 2017. https://cran.r-project.org/src/contrib/Archive/dyn/dyn_0.2-8.1.tar.gz. Accessed 1 May 2016.
- Pohlert T. Trend: non-parametric trend tests and change-point detection. R package version 0.2.0. 2015. https://cran.r-project.org/web/packages/trend/index.html. Accessed 20 May 2016.
- 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
- Midekisa A, Senay G, Henebry GM, Semuniguse P, Wimberly MC. Remote sensing-based time series models for malaria early warning in the highlands of Ethiopia. Malar J. 2012;11:165.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Pascual M, Cazelles B, Bouma MJ, Chaves LF, Koelle K. Shifting patterns: malaria dynamics and rainfall variability in an African highland. Proc Biol Sci. 2008;275:123–32.View ArticlePubMedGoogle Scholar
- Wyse APP, Bevilacqua L, Rafikov M. Simulating malaria model for different treatment intensities in a variable environment. Ecol Model. 2007;206:322–30.View ArticleGoogle Scholar
- Confalonieri U, Menne B, Akhtar R, Ebi KL, Hauengue M, Kovats RS, et al. Human health. Climate change 2007: impacts, adaptation and vulnerability. Contribution of working group II to the fourth assessment report of the intergovernmental panel on climate change; 2007.Google Scholar
- Small J, Goetz SJ, Hay SI. Climatic suitability for malaria transmission in Africa, 1911–1995. Proc Natl Acad Sci USA. 2003;100:15341–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Paaijmans KP, Imbahale SS, Thomas MB, Takken W. Relevant microclimate for determining the development rate of malaria mosquitoes and possible implications of climate change. Malar J. 2010;9:196.View ArticlePubMedPubMed CentralGoogle Scholar
- Smith M, Macklin M, Thomas C. Hydrological and geomorphological controls of malaria transmission. Earth Sci Rev. 2013;116:109–27.View ArticleGoogle Scholar
- Hardy AJ, Gamarra JG, Cross DE, Macklin MG, Smith MW, Kihonda J, et al. Habitat hydrology and geomorphology control the distribution of malaria vector larvae in rural Africa. PLoS ONE. 2013;8:e81931.View ArticlePubMedPubMed CentralGoogle Scholar
- Patz JA, Olson SH. Malaria risk and temperature: influences from global climate change and local land use practices. Proc Natl Acad Sci USA. 2006;103:5635–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Beck-Johnson LM, Nelson WA, Paaijmans KP, Read AF, Thomas MB, Bjornstad ON. The effect of temperature on Anopheles mosquito population dynamics and the potential for malaria transmission. PLoS ONE. 2013;8:e79276.View ArticlePubMedPubMed CentralGoogle Scholar
- Bi P, Tong S, Donald K, Parton KA, Ni J. Climatic variables and transmission of malaria: a 12-year data analysis in Shuchen County, China. Public Health Rep. 2003;118:65–71.View ArticlePubMedPubMed CentralGoogle Scholar
- Chikodzi D. Spatial modelling of malaria risk zones using environmental, anthropogenic variables and geographical information systems techniques. J Geosci Geomat. 2013;1:8–14.Google Scholar