Effect of climatic variability on malaria trends in Baringo County, Kenya
Malaria Journalvolume 16, Article number: 220 (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.
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,14,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.
This study focused on Baringo County, Kenya, located between 35.602°E, 36.277°E and 0.541°N, 0.723°N, at altitudes ranging between 870 and 2499 m above sea level (asl). The county covers 11,075 km2 with an estimated population of 555,561 as at 2009 . The average annual rainfall ranges between 300 and 1500 mm while air temperature varies between 16 °C in the highland areas and 42 °C in the lake ecosystem . Rainfall pattern is tri-modal, with the long rains received during March–May (MAM) and two short rains seasons that are experienced between June–August (JJA) and October–December (OND). This study focused on the central part of Baringo County (Fig. 1) encompassing three sub-counties (Baringo North, Baringo Central and Baringo South). The study area was sub-divided into four ecological zones on the basis of hydrology, altitude, vegetation cover, soil types, and precipitation. The four zones were the lowland zone lying at an elevation of 1000 m asl and surrounding the permanent water bodies (L. Baringo, L. Bogoria and L. 94), the mid-altitude zone lying between 1000 and 1500 m asl, the highland zone, lying between 1,500 and 2,300 m asl and the riverine zone (bordering the Kerio River) at altitudes ranging between 1100 and 1200 m asl .
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)
The long-term (2004–2014) mean monthly malaria cases varied across zones. On average, the observed monthly malaria cases were highest in the highlands and lowest in mid-altitudes. Low variability was observed in midland zone while high variability was observed in the other three zones. Highest and lowest malaria cases per zone were however recorded in different months and years over the study period. Table 1 summarizes the data for each zone over the study period.
Yearly cumulative malaria cases (2004–2014) and trends
The highest yearly cumulative malaria cases were recorded in 2012 (11,249 cases) in the highland, 2008 (971 cases) in the mid-altitude, 2009 (9275 cases—over 9 year period) in the lowland and 2010 (7119 cases—over 9 year period) in the riverine zone. Health records for the year 2004 in the lowland and the year 2004 and 2005 in the riverine zone were however not available. The Seasonal Mann–Kendall trend test showed a significant increasing trend in malaria cases over the study period in the highland (z = 2.5, p = 0.0142) and mid-altitude (z = 2.5, p = 0.0141) zones while no change was observed in the lowland (z = −0.5, p = 0.609) and riverine zones (z = −0.6, p = 0.567) (Fig. 2). However, further analysis showed a significant decrease in malaria cases in the riverine zone between 2011 and 2014 (z = −3.2, p = 0.0012).
Malaria seasonality over 10 year period
Three malaria peak seasons were observed in the riverine, highland and mid-altitude zones, while two seasons were observed in the lowland zone (Fig. 3). Pair wise mean comparison tests for the monthly malaria cases in each zone yielded different results. In the highland zone, the month of October experienced significantly higher malaria cases compared to December (estimate = 315.8, CI 9.95–621.65). Mean differences for the other months were statistically insignificant. In the mid-altitude zone, there were significant differences in mean monthly malaria cases for the month of July and April (estimate = 42.3, CI 3.06–81.54), July and December (estimate = 40.7, CI 1.46–79.94), October and February (estimate = 43.4, CI 4.16–82.64), October and March (estimate = 41.7, CI 2.46–80.94) and, October and August (estimate = 50.2, CI 10.96–89.44). The months of July and October recorded significantly higher malaria cases in the highland and mid-altitude zones, a period corresponding to the short rainy seasons of June–August and October–November. There were however no significant differences in mean monthly malaria cases in the lowland and riverine zones.
Relating malaria cases to climatic variables
Long-term trends of malaria cases against climatic variables
In the highland zone, peaks in malaria cases followed low peaks of rainfall while declines were associated with high peaks of rainfall and decrease in maximum temperatures. Decline in malaria cases between September and December 2006 corresponded to heavy precipitation and reduced maximum temperature recorded the previous months (Fig. 4). In the mid-altitude zone, malaria cases and rainfall followed a similar upward trend between 2006 and 2007, before declining between 2008 and 2009 following extremely low rainfall amounts and increased maximum temperatures recorded over the same period. Again in 2012, malaria cases rose following increased rainfall and slight increase in maximum temperatures recorded that year (Fig. 5). In the lowland zone, malaria cases increased following an increase in rainfall and a decrease in maximum temperature between 2006 and 2007. A sharp decline in malaria cases was however noted between 2007 and 2008 following high rainfall and low temperatures (Fig. 6). In the riverine zone, malaria cases increased between 2006 and 2007 and again between 2009 and 2010 following increased rainfall, while a decrease in malaria cases was reported between 2011 and 2013, corresponding to high rainfall amounts and low temperatures received over the same period (Fig. 7).
Identification of possible climatic and environmental predictors of malaria cases
Using the sample cross correlation function, it was established that rainfall at lags 1 and 2 months across all zones and also at lag 0 in the lowland zone had the most dominant cross correlation with malaria cases. Temperature at lag 0 in the lowland and riverine zones, and lag 1 in the highland, mid-altitude and riverine zones had the most dominant cross correlation with malaria cases. Enhanced Vegetation Index at lag 0 in the highland, mid-altitude and lowland zones and lag 1 and 2 in the highland and riverine zones had the most dominant cross correlations with malaria cases. Based on the CCF, a negative binomial regression model was built, and it was observed that rainfall lagging 2 months in all the zones were statistically significant at 5% level of significant. However, rainfall in midland was highly significant as compared to the other zones. Further it was observed that maximum temperature at lag 0 in the riverine zone and lag 1 in the highland zone were significant. Joint significance of rainfall and temperatures in highland and riverine can be attributed to high variations in malaria cases in these zones (Table 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,40,41,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.
WHO. World malaria report 2013. Geneva: World health Organization; 2013.
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.
Fernando SD, Rodrigo C, Rajapakse S. The ‘hidden’ burden of malaria: cognitive impairment following infection. Malar J. 2010;9:366.
Ministry of Health. Towards a malaria free Kenya: Kenya revised malaria strategy 2009–2018. National Malaria Control Programme; 2014.
Ministry of Health. National guidelines for the diagnosis, treatment and prevention of malaria for health workers in Kenya. Kenya: Div Malar Control, MoH; 2006.
Ministry of Health. National malaria strategy 2041–2010. Kenya: Div Malar Control, MoH; 2001.
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.
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.
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.
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.
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.
Mellor PS, Leake CJ. Climatic and geographic influences on arboviral infections and vectors. Rev Sci Tech. 2000;19:41–54.
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.
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.
Parham PE, Michael E. Modeling the effects of weather and climate change on malaria transmission. Environ Health Perspect. 2010;118:620–6.
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.
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.
Briet OJ, Vounatsou P, Amerasinghe PH. Malaria seasonality and rainfall seasonality in Sri Lanka are correlated in space. Geospat Health. 2008;2:183–90.
Evrendilek F, Gulbeyaz O. Deriving vegetation dynamics of natural terrestrial ecosystems from MODIS NDVI/EVI data over Turkey. Sensors. 2008;8:5270–302.
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.
Schultz P, Halpert M. Global correlation of temperature, NDVI and precipitation. Adv Space Res. 1993;13:277–80.
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.
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.
Division of malaria control. 2010 Kenya malaria indicator survey. Nairobi: Ministry of Public Health and Sanitation; 2011.
Mohajan HK. Improvement of health sector in Kenya. Am J Public Health Res. 2014;2:159–69.
KNBS. The 2009 Kenya population and housing census, vol. IC. In: Population distribution by age, sex and administrative units. 2010.
Odada EO, Onyando JO, Obudho PA. Lake Baringo: addressing threatened biodiversity and livelihoods. Lakes Reserv Res Manag. 2006;11:287–99.
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.
Buuren S, Groothuis-Oudshoorn K. Mice: multivariate imputation by chained equations in R. J Stat Softw. 2011;45:1–67.
Vink G, Frank LE, Pannekoek J, Buuren S. Predictive mean matching imputation of semicontinuous variables. Stat Neerl. 2014;68:61–90.
Molenberghs G, Fitzmaurice G, Kenward MG, Tsiatis A, Verbeke G. Handbook of missing data methodology. Boca Raton: CRC Press; 2014.
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.
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.
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.
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.
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.
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.
Wyse APP, Bevilacqua L, Rafikov M. Simulating malaria model for different treatment intensities in a variable environment. Ecol Model. 2007;206:322–30.
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.
Small J, Goetz SJ, Hay SI. Climatic suitability for malaria transmission in Africa, 1911–1995. Proc Natl Acad Sci USA. 2003;100:15341–5.
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.
Smith M, Macklin M, Thomas C. Hydrological and geomorphological controls of malaria transmission. Earth Sci Rev. 2013;116:109–27.
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.
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.
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.
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.
Chikodzi D. Spatial modelling of malaria risk zones using environmental, anthropogenic variables and geographical information systems techniques. J Geosci Geomat. 2013;1:8–14.
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.