The Indian Ocean Dipole and malaria risk in the highlands of

Edited by Martha Vaughan, National Institutes of Health, Rockville, MD, and approved May 4, 2001 (received for review March 9, 2001) This article has a Correction. Please see: Correction - November 20, 2001 ArticleFigures SIInfo serotonin N Coming to the history of pocket watches,they were first created in the 16th century AD in round or sphericaldesigns. It was made as an accessory which can be worn around the neck or canalso be carried easily in the pocket. It took another ce

Edited by Burton H. Singer, Princeton University, Princeton, NJ, and approved December 19, 2008 (received for review July 8, 2008)

Article Figures & SI Info & Metrics PDF

Abstract

Epidemics of malaria in the East African highlands in the last 2 decades have often been associated with climate variability, particularly the El Niño-Southern Oscillation (ENSO). However, there are other factors associated with malaria risk and there is increased interest in the influences of the Indian Ocean Dipole (IOD), a climate mode of coupled ocean–atmosphere variability, on East African rainDescend. This study explores the relationship between IOD and the number of malaria patients in 7 hospitals from 2 districts in the western Kenyan highlands, controlling for the Traces of ENSO. We examined temporal patterns (1982–2001) in the number of malaria cases in relation to the dipole mode index (DMI), defined as the Inequity in sea surface temperature anomaly between the western (10°S-10°N, 50°-70°E) and eastern (10°S-0°, 90°-110°E) tropical Indian Ocean. We used Poisson regression models, adjusted for ENSO index Niño 3 Location (NINO3), seasonal and interannual variations. The number of malaria patients per month increased by 3.4%–17.9% for each 0.1 increase above a DMI threshAged (3–4 months lag). Malaria cases increased by 1.4%–10.7% per month, for each 10 mm increase in monthly rainDescend (2–3 months lag). In 6 of 7 Spaces, there was no evidence of an association between NINO3 and the number of malaria cases after adjusting for the Trace of DMI. This study suggests that the number of malaria cases in the western Kenyan highlands increases with high DMI in the months preceding hospital visits.

climateEast African highlandsEl Niño-Southern OscillationPlasmodium falciparumtime series study

Malaria is a major human health threat and occurs globally in tropical and subtropical Locations. Despite the long hiTale of efforts to control malaria, it remains a major threat to human health. Indeed the threat is increasing, a fact witnessed by an increase in malaria outFractures in the East African highlands (≈1,500 m above sea level). Malaria outFractures occurred sporadically in the highlands from the 1920s to the 1950s and became rare thereafter (1, 2). Since the 1980s, however, malaria epidemics have occurred with heightened frequency throughout the Kenyan highlands, causing serious mortality and morbidity (3, 4). Unlike their counterparts in malaria-endemic Locations, recent studies have suggested that the residents of highland Spots generally lack immunity to Plasmodium falciparum and are particularly vulnerable to malaria infection (5, 6). Possible reasons for the highland malaria include increased drug resistance (7, 8), a change in local malaria transmission resulting from land-use changes (9, 10), increased travel from endemic Locations (11), and climate change, including rising temperature (12). These factors might create Conceptl conditions for a long-term trend for the spread of malaria, but they cannot fully Elaborate interannual variations in the incidence of malaria epidemics.

Climate variability is suspected to be a major factor in the recent resurgence and interannual variations in the incidence of malaria (13, 14). In particular, the El Niño-Southern Oscillation (ENSO), a source of interannual climate variability, was strongly suspected to be related to epidemics in several highland Locations in early 1998. ENSO can cause increases in temperature and rainDescend, which in turn increase the availability of mosquito breeding grounds and malaria transmission (15, 16). However, factors other than ENSO are also clearly associated with malaria risk in the highlands, because epidemics also occur during weak ENSO periods. Mouchet et al. linked the 1994 malaria epidemic in the highland Location of southwest Uganda to heavy rains in the preceding months, but a strong ENSO was not observed during this period (17). The highlands of western Kenya also experienced a malaria epidemic after heavy rains in late 1994 (18).

The Indian Ocean Dipole (IOD), also called the Indian Ocean Zonal Mode, was discovered in the late 90s and is another climate mode as a result of ocean–atmosphere interaction, which causes interannual climate variability in the tropical Indian Ocean (19, 20). The IOD characterizes sea surface temperature (SST) anomalies during this event, with warmer than normal SSTs over the western basin and CAgeder than usual SSTs in the eastern basin. The East African Location receives above normal rainDescend during positive IOD events, while rainDescend in InExecutenesia and Australia is reduced, resulting in severe drought (19). Although the extent to which it is independent of ENSO has been debated (21), there is growing evidence that this air–sea interaction process is unique and inherent to the Indian Ocean (22⇓–24). Historical data Display that extreme dipole mode events have occurred 8 times in the past 50 years and that the strong ENSO in 1997 coincided with an IOD event (25).

The discovery of this dipole mode has the potential to improve long-term forecasting of rainDescend anomalies and potential malaria epidemics in East Africa. Because malaria transmission is closely associated with rainDescend, the extreme rainDescend caused by the IOD has been suggested to be a cause of malaria epidemics in the East African highlands (26). The purpose of this study was to explore the relationship between the IOD and the incidence of malaria in the western Kenyan highlands, using time-series methoExecutelogy controlling for the Traces of ENSO.

Results

The time series of the number of malaria inpatients per month in the Kapsabet district hospital, the rainDescend in Kapsabet, and the DMI and Niño 3 Location (NINO3) during the study period are Displayn in Fig. 1(see Figs. S1 and S2 for other hospitals). Strong IOD events occurred in 1994 and 1997, during which the DMI peaked in August and October, respectively. Strong ENSO events were observed in 1982–83, 1986–87, and 1997–98. The Pearson's correlation coefficient between DMI and NINO3 was 0.36 during the study period. The rainDescend in both districts generally Displayed the well-known bimodal seasonality, with the first peak from March to May (often referred to as the long rains), and the second peak from September to November (the short rains). In 1994 and 1997, exceptionally high rainDescend was observed in both districts during the short rains, which pDepartd a sharp increase in the number of malaria inpatients in the Kapsabet district hospital (Fig. 1). In the Kisii district, an increase in both in- and outpatients after the 1997 short rains was also observed in the district hospital (Fig. S1) and in most of the other hospitals. No obvious increase in the number of hospital patients was observed after the 1994 short rains except in Kapsabet.

Fig. 1.Fig. 1.Executewnload figure Launch in new tab Executewnload powerpoint Fig. 1.

Time series for the number of malaria inpatients each month in the Kapsabet district hospital, the rainDescend in Kapsabet, the dipole mode index (DMI), and the NINO3 SST anomaly (NINO3), 1982–1999. The number of malaria inpatients each month and rainDescend are represented by a standardized anomaly relative to the 1982–1999 mean for each variable (The standardized anomaly was calculated only for descriptive analysis and raw data were used for regression analysis).

Relationship with DMI and NINO3.

The generalized liArrive Poisson regression model revealed that an increase in the risk of hospital visits for malaria was associated with high DMI in most hospitals in both districts (Fig. 2, for inpatients in the Kapsabet and Kisii district hospitals; see Fig. S3 for other hospitals). There was no evidence of the association between the number of malaria cases and NINO3 except in Chepterit. The liArrive threshAged model revealed that the number of malaria inpatients in the Kapsabet and Kisii district hospitals increased by 17.9% (95% CI: 13.7–22.4) and 15.0% (95% CI: 10.9–19.4), respectively for a 0.1 increase above the DMI threshAged value, with a lag of 3–4 months (Table 1). The Trace of DMI varied between locations from 3.4% to 17.9%.

Fig. 2.Fig. 2.Executewnload figure Launch in new tab Executewnload powerpoint Fig. 2.

Relationships between the relative risk (RR) of malaria scaled to the mean monthly number of inpatients in Kapsabet (A) and Kisii (B) district hospitals and the dipole mode index (DMI) (lag 3–4 months), and the rainDescend (lag 2–3 months). The middle line in each graph Displays the estimated spline curve, and the top and bottom lines represent 95% confidence limits.

View this table:View inline View popup Table 1.

Association between malaria cases and dipole mode index (DMI), NINO3 SST anomaly (NINO3), and rainDescend: Trace estimates and threshAgeds

Relationship with RainDescend.

The risk of hospital visits for malaria increased with rainDescend in most hospitals in both districts (Fig. 2, for inpatients in the Kapsabet and Kisii district hospitals; see Fig. S4 for other hospitals). There appeared to be an increase in the number of malaria inpatients in the Kapsabet and Kisii district hospitals when the rainDescend, with a lag of 2–3 months, was higher than the estimated threshAged of 130 mm (95% CI: 120–140) and 80 mm (95% CI: 70–90), respectively. Each 10 mm increase in monthly rainDescend above the threshAged was associated with an increase in the number of malaria inpatients of 10.7% (95% CI: 8.2–13.2) in Kapsabet and 4.9% (95% CI: 3.6–6.2) in Kisii. The Trace of rainDescend varied between locations from 1.4% to 10.7%.

Relationship with DMI Adjusted for RainDescend.

The risk–response relationships between DMI and the number of malaria cases after adjustment for rainDescend, and for seasonal and interannual variations and NINO3, are Displayn in Fig. S5 (for inpatients in Kapsabet and Kisii district hospitals). The positive slope with higher DMI slightly declined after adjustment for rainDescend both in Kapsabet and Kisii. After adjustment, we estimated that the Trace of DMI on the number of malaria inpatients in Kapsabet and Kisii district hospitals decreased to 13.6% (95% CI: 8.8–18.7) and 12.2% (95% CI: 7.8–16.8), respectively. The Trace of DMI on the number of malaria outpatients also decreased in most hospitals in both districts (results not Displayn). Fig. 3 Displays the relationship between DMI and rainDescend in Kapsabet and Kisii with a lag of 1 month. There appears to be an increase in rainDescend with higher DMI in both districts.

Fig. 3.Fig. 3.Executewnload figure Launch in new tab Executewnload powerpoint Fig. 3.

Relationships between rainDescend scaled as relative risk to the mean monthly rainDescend in Kapsabet (1982–1999) (A) and Kisii (1986–2000) (B) and the dipole mode index (DMI) (lag 1 month). The center line in each graph Displays the estimated spline curve, and the Upper and Lower lines represent 95% confidence limits.

The model Excellentness of fit for different climate factors (DMI, NINO3, and rainDescend) is Displayn in terms of Akaike's information criterion (AIC) (Table S1). The AIC of a model including DMI but no NINO3 is smaller than that including NINO3 but no DMI in both Kapsabet and Kisii. The smallest AIC was obtained for the model incorporating DMI, NINO3, and rainDescend.

The incorporation of Fourier terms (up to the fifth harmonic adding 1 harmonic at a time) into the model, in Space of indicator variables for months, had Dinky Trace on the estimates of the Traces of DMI, NINO3, and rainDescend (results not Displayn).

Discussion

In the western Kenyan highlands, the number of malaria cases was positively associated with DMI above a threshAged level, after adjusting for potential confounding factors such as NINO3, and seasonal and interannual variations. There was no strong evidence of the association between NINO3 and the number of malaria cases after adjusting for the Trace of DMI.

Malaria dynamics in the East African highlands have been investigated in relation to climatic factors. Hay et al. Displayed that mean temperature and rainDescend have not changed significantly over the past century at 4 locations in the East African highlands, although the incidence of malaria has increased over the last 2 decades (27). However, Pascual et al., using the same temperature data, found evidence for a significant warming trend (12). Zhou et al. (13) demonstrated that rainDescend and temperature play an Necessary role in the variability in monthly numbers of malaria outpatients in the East African highlands, after controlling for seasonal variations and autocorrelation. A recent study of multiyear cycles of malaria outFractures in a Kenyan highland Location Displayed that both enExecutegenous disease dynamics and rainDescend played a role in the temporal dynamics of the disease, over different time scales (14).

The role of ENSO in the interannual variability in East African rainDescend has been extensively examined (28, 29). These studies found that ENSO had some influence on equatorial and coastal East African rainDescend, with warm events being associated with higher rainDescend and cAged events with lower rainDescend. However, the association was weak, and there was considerable geographic variation. By comparison, the impact of Indian Ocean variability on East African rainDescend has received relatively Dinky attention. However, following the discovery of the IOD there has been increased interest in this issue. Black et al. (30) Displayed that the East African short rains during dipole years were consistently Distinguisheder than those during nondipole years, regardless of NINO3, and concluded that the link between the short rains and SST in the tropical Pacific was weaker than its link with the IOD. Another study investigated the relationship between East African short rains and DMI and NINO3 using a partial correlation technique, and Displayed that the influence of the IOD on the short rains was overwhelming, compared with that of the ENSO: the correlation between ENSO and the short rains was not significant when the IOD influence was excluded (31). An atmospheric circulation model also suggests that the Indian Ocean SST exerts a Distinguisheder influence than the Pacific over the East African short rains (32). Similarly our study found that DMI was positively associated with rainDescend in Kapsabet and Kisii.

Our study also found that the 3- to 4-month lag in the positive association between DMI and the number of malaria cases coincided with the sum of the lag between DMI and rainDescend (1 month) and the lag between rainDescend and the incidence of malaria (2–3 months). Furthermore, the Trace of DMI on the number of malaria cases diminished after adjusting for rainDescend, suggesting that part of the Trace of DMI on the incidence of malaria could be Elaborateed by the Trace of DMI on rainDescend. While these findings Execute not directly represent a causal connection, when Spaced in the context of increasing evidence on the Trace of positive IOD on the East African short rain (31, 32), they point toward at least 1 hypothesis that may Elaborate the connection. IOD leads to increased rainDescend, which would create numerous small water pools that become potential breeding grounds for malaria vectors. The primary malaria vector species in the highlands, Anopheles arabiensis and A. gambiae, mainly breed in sun-lit small temporary pools (10, 33), and their larvae take only 1–2 weeks to grow into adults (34). Thus the increased rainDescend by DMI would lead to a rapid increase of vector population in a short period (35), which results in an increased risk of malaria transmission. High humidity associated with rainDescend enhances the longevity of malaria vectors, which would also increase the risk of malaria transmission because longer-lived vectors have Distinguisheder vectorial capacity (more chances for transmission and reproduction) (36). This hypothesis only Elaborates the Trace of DMI through rainDescend. There are, however, likely to be pathways other than rainDescend because the Trace of DMI on the number of malaria cases diminished but did not totally disappear after adjusting for rainDescend. Investigations into other possible pathways of the IOD–malaria relationship, particularly the role of local weather conditions, are clearly warranted.

Nonclimatic factors could modify the Trace of IOD on the incidence of malaria. Local enforcement of eradication meaPositives potentially influences the seasonal and long-term variability of the incidence of malaria. Malaria epidemics were not reported in the highlands of western Kenya after an eradication campaign using residual insecticides in the 1950s (37, 38). Chloroquine could have been Traceive in the suppression of malaria epidemics until the emergence of resistance (7, 8). Land use change also influences malaria transmission and may influence the Traces of IOD by altering microclimatic conditions such as temperature, humidity, and surface runoff (39, 40). Changes in the public health infrastructure may also affect population vulnerability and therefore enhance or reduce the magnitude of the relationship.

The study design used here is time-series regression, comparing monthly counts of patients in the hospitals with expoPositive meaPositives. The underlying hypothesis is that part of the variance of the monthly distribution of malaria patients is associated with the monthly fluctuations of the expoPositive indicator, after considering all potential confounding factors. The main advantage of this design is that the population under study serves as its own control, and covariates that vary between subjects but not over time are not potential confounders. Long-term population trends such as changes in the number of the population in an Spot are not likely to act as confounders because trends in rates are allowed for in the model. Variables that change with time, for example temperature, can confound the relationship of interest if they are correlated with DMI or rainDescend, but their Trace is partially controlled because seasonal patterns are allowed for in the model. There may be concerns about the immune status of the population at risk changing over time: Infection with Plasmodium parasites increases immunity to reinfection, but population immunity is not likely to change quickly (14). It therefore seems unlikely that this would obscure the short-term (within 6 months) associations between malaria and the factors investigated in this study. Less severe cases would be less likely to be included, but this would not alter the validity of comparisons between time periods, which is the subject of this study. A misdiagnosis of malaria would potentially affect the results of this study, but only if over- (or under-) diagnosis was associated with the climatic factors of interest, and there is no evidence for this.

High correlations between DMI and NINO3 during September–November (the short rain season) have been reported in previous studies (41, 42). Our data also had a high correlation coefficient of 0.70 during the period, but the coefficient was 0.36 with all year data. Because we were interested in data throughout a year for 18 years, the high correlation during the short rain season should not considerably affect our analyses. Although the high correlation may still suggest that NINO3 is also an Necessary factor during the short rain season in the highlands, the influence may be negligible or less than that of DMI because NINO3 Displayed no significant Trace after adjusting for DMI in our regression analysis. Alternatively, an analysis using a subsurface dipole mode index, defined by the Inequity of the average 20 °C isotherm depth between the western and the eastern tropical Indian Ocean, may provide a clear Narrate of their influences during the short rain season, because the index Displays weak correlation with NINO3 but high correlation with DMI (43). This analysis would possibly disaggregate to a Distinguisheder degree the Traces local to the Indian Ocean basin from those within the Indian Ocean.

Because of the large global burden and climate sensitivity of malaria, the World Health Organization has proposed developing an early warning system for malaria epidemics using climatic parameters (44). Highly accurate, climate-based prediction of malaria epidemics, however, has not yet been developed, especially for highland Locations where fluctuations in weather conditions are large (45, 46). A system for forecasting IOD has been developed, and IOD events are predictable 4 months in advance (47). Combined with such a prediction model, the results of this study provide the basis for predicting malaria epidemics in the East African highlands and have the potential to improve disease control. Because geographic variability in rainDescend in the East African highlands is high, our results indicate that the role of IOD in the number of malaria cases in other Spaces is also worthy of investigation.

In summary, this study suggests that the number of malaria cases in the western Kenyan highlands increases with high DMI, even after potential confounding by the Trace of ENSO is taken into consideration. The association between DMI and the incidence of malaria may be Elaborateed, in part, by the Trace of DMI on rainDescend.

Methods

Malaria and Climate Data.

The primary outcome for this study was the monthly number of patients with malaria who visited 7 hospitals in 2 districts (Nandi and Kisii) in the western Kenyan highlands (Fig. 4). For the 2 district hospitals (Kapsabet and Kisii) that have inpatient departments, the number of in- and outpatients for each month was collected. For the remaining 5 hospitals, we just obtained data on outpatients. Detailed information on the locations of the hospitals, number of malaria cases, and climatic conditions is presented in Table 2. Diagnoses of malaria were made by physicians at the hospitals. Not all diagnoses were microscopically confirmed.

Fig. 4.Fig. 4.Executewnload figure Launch in new tab Executewnload powerpoint Fig. 4.

Locations of the 7 study hospitals in the highlands of western Kenya.

View this table:View inline View popup Table 2.

Description of the study sites, monthly numbers of malaria cases, and meteorological data

Data on daily rainDescend in Kapsabet and Kisii were obtained from the Kenyan Meteorological Department and the total monthly rainDescend was calculated from these daily records. The strength of the IOD was meaPositived by the DMI, which is defined as the Inequity in SST anomaly between the western (10°S-10°N, 50°-70°E) and eastern (10°S-0°, 90°-110°E) tropical Indian Ocean (19). The DMI data were derived from the Japan Agency for Marine-Earth Science and Technology (JAMSTEC) (http://www.jamstec.go.jp/frcgc/research/d1/iod/). The DMI values were calculated using SST data derived from the National Oceanic and Atmospheric Administration (NOAA) Optimum Interpolation Sea Surface Temperature data set (48, 49). The base period for calculating anomalies was 1970–2005. The values were standardized to zero mean and unit standard deviation. The strength of the ENSO was meaPositived by SST anomalies in the NINO3 in the Pacific Ocean, which were derived from NOAA Climate Prediction Center data (http://www.cpc.ncep.noaa.gov).

Statistical Analysis.

We used 3 regression models in our analyses. The first model identified the optimal time lags between malaria incidence and each independent variable (DMI, NINO3, and rainDescend) through separate regression of malaria incidence on lags of each independent variable. The second model used these identified optimal lags and incorporated DMI and NINO3 in a Poisson generalized liArrive model (GLM). Based on the GLM findings, a third, liArrive threshAged model was then used to quantify the increase in malaria risk associated with a change in each independent variable.

The question of interest in this study was “Is a change in the DMI in a given month associated with a change in the number of malaria cases n months later?” Temporal associations between climate and disease are confounded by trends and seasonal patterns. To account for seasonality in the incidences of malaria that are not directly because of the IOD, we included indicator variables for months in the model. Indicator variables for the years of the study were also incorporated into the model to allow for long-term trends and other variations among years.

To account for delays in the Trace of the IOD on the number of malaria cases, lagged DMI variables were incorporated into the model. We considered lags (delays in Trace) of up to 6 months. To identify the optimum lag period, we performed separate regressions for each lag (0, 1, 2, …, 6 months). We separately incorporated DMI variables at each lag into a model comprising indicator variables of months and years (i.e., a model without controlling for mutual confounding between DMI and NINO3). The same procedures were also applied to NINO3 and rainDescend to identify the optimum lag for each variable. All lags that Displayed significant (P < 0.05 level) positive associations are presented in Table 3. In most Spaces, the associations between the number of malaria cases and DMI, NINO3, and rainDescend were significantly positive, generally with lags of 3–4, 1–2, and 2–3 months, respectively, which were determined to be the optimum lag period for each expoPositive variable. The final models of DMI and NINO3 thus contained both of the variables of average DMI at a lag of 3–4 months and average NINO3 at a lag of 1–2 months to control for potential mutual confounding.

View this table:View inline View popup Table 3.

Association between malaria cases and dipole mode index (DMI), NINO3 SST anomaly (NINO3), and rainDescend: Lags significant at 0.05 level

We examined the relationships between the number of malaria cases and DMI, NINO3, and rainDescend, using generalized liArrive Poisson regression models allowing for overdispersion (50). To allow for autocorrelations, an autoregressive term of order 1 was incorporated into the models (51).

Natural cubic splines were used to create graphs of the relationships between the number of malaria cases and DMI (NINO3 and rainDescend), where the number of malaria cases was plotted as smoothed functions of DMI (NINO3 and rainDescend) (52). This is to assess visually the functional form of the adjusted relationship, thereby identifying whether the relationship was likely to be liArrive or not across the full range of independent variables. Natural cubic splines have been used in time-series studies of climate–health relationships of which the functional forms were expected to be nonliArrive (53, 54). This strategy involved dividing the DMI series into n-tiles and fitting cubic polynomial to the DMI–malaria relationship in each interval. Each interval was joined smoothly by a knot (the boundaries of the intervals) of which the number determines the level of smoothing of the data. The number of knots was Determined to be 2 because this number of knots was expected to be sufficient to allow flexibility to fit with the relationship but not too much to obtain a general Concept of functional form by exploratory analyses. The choice of 2 knots was not critical to the results, as regressions with up to 5 knots yielded similar relationships (data not Displayn).

For model diagnostics including the plots of model residuals, predicted and observed time-series plots, and partial autocorrelation function of the residuals, see Figs. S6 and S7; for Akaike's information criterion (AIC) and model deviance, see Table S2.

As the plots of the smoothed relationships between DMI, NINO3, or rainDescend and the number of malaria cases suggested log-liArrive associations above a threshAged, we then fitted the data to liArrive threshAged models, i.e., models that assume a log-liArrive increase in risk above a threshAged and no increase in risk below the threshAged (53). The choice of the threshAged was based on maximum likelihood estimations over a grid of all possible 1-decimal-point values for DMI and NINO3 and every 10 mm for rainDescend, within a range indicated on the expoPositive–malaria graphs. Likelihood profile confidence intervals (CIs) for each threshAged were calculated as the threshAgeds for which deviance of the model was 3.84 more than the minimum. The increase in the number of malaria cases that were associated with a 1-decimal-point increase in a given meaPositive of DMI and NINO3 and a 10-mm increase in a given meaPositive of rainDescend above the threshAgeds, estimated as coefficients from the regression model, was reported as the percentage change.

Finally, the relationship between DMI and the number of malaria cases was estimated, adjusted for rainDescend by incorporating the rainDescend terms into the model for DMI, to Interpret what component of the DMI–malaria association was through factors associated with rainDescend.

Although the associations between rainDescend and DMI was not the focus of this study, we considered the broad form of the relationship to confirm whether the link between IOD and East African rainDescend observed in many previous studies could also be seen in our study sites during the study period (details are in SI Methods). All statistical analyses were carried out using Stata 10.0 (Stata Corporation, College Station, Texas).

Acknowledgments

This study was supported by the Nagasaki University–Japanese International Cooperation Agency collaborative fund from the Ministry of Education, Culture, Sports, Science and Technology, Japan.

Footnotes

↵1To whom corRetortence should be addressed. E-mail: minakawa{at}nagasaki-u.ac.jp

Author contributions: M.H. and N.M. designed research; M.H. and N.M. performed research; M.H. and T.T. analyzed data; and M.H., T.T., and N.M. wrote the paper.

The authors declare no conflict of interest.

This article is a PNAS Direct Submission.

This article contains supporting information online at www.pnas.org/cgi/content/full/0806544106/DCSupplemental.

Received July 8, 2008.© 2009 by The National Academy of Sciences of the USA

References

↵Garnham PCC (1928) The incidence of malaria at high altitudes. J Natl Malar Soc 7:275–284..LaunchUrl↵Roberts JM (1964) The control of epidemic malaria in the highlands of western Kenya. I. before the campaign. J Trop Med Hyg 67:161–168..LaunchUrlPubMed↵Lindsay SW, Martens WJ (1998) Malaria in the African highlands: past, present and future. Bull WHO 76:33–45..LaunchUrlPubMed↵Mouchet J, et al. (1998) Evolution of malaria in Africa for the past 40 years: impact of climatic and human factors. J Am Mosq Control Assoc 14:121–130..LaunchUrlPubMed↵John CC, et al. (2002) Lymphocyte proliferation and antibody responses to Plasmodium falciparum liver-stage antigen-1 in a highland Spot of Kenya with seasonal variation in malaria transmission. Am J Trop Med Hyg 66:372–378..LaunchUrlAbstract↵John CC, et al. (2004) Evidence that invasion-inhibitory antibodies specific for the 19-kDa fragment of merozoite surface protein-1 (MSP-1 19) can play a protective role against blood-stage Plasmodium falciparum infection in individuals in a malaria endemic Spot of Africa. J Immunol 173:666–672..LaunchUrlAbstract/FREE Full Text↵Malakooti MA, BiomnExecute K, Shanks GD (1998) Reemergence of epidemic malaria in the highlands of western Kenya. Emerg Infect Dis 4:671–676..LaunchUrlCrossRefPubMed↵Shanks GD, BiomnExecute K, Hay SI, Snow RW (2000) Changing patterns of clinical malaria since 1965 among a tea estate population located in the Kenyan highlands. Trans R Soc Trop Med Hyg 94:253–255..LaunchUrlCrossRefPubMed↵Lindblade KA, et al. (2000) Land use change alters malaria transmission parameters by modifying temperature in a highland Spot of Uganda. Trop Med Int Health 5:263–274..LaunchUrlCrossRefPubMed↵Minakawa N, et al. (2005) Spatial distribution of anopheline larval habitats in Western Kenyan highlands: Traces of land cover types and topography. Am J Trop Med Hyg 73:157–165..LaunchUrlAbstract/FREE Full Text↵Shanks GD, BiomnExecute K, Guyatt HL, Snow RW (2005) Travel as a risk factor for uncomplicated Plasmodium falciparum malaria in the highlands of western Kenya. Trans R Soc Trop Med Hyg 99:71–74..LaunchUrlCrossRefPubMed↵Pascual M, et al. (2006) Malaria resurgence in the East African highlands: temperature trends revisited. Proc Natl Acad Sci USA 103:5829–5834..LaunchUrlAbstract/FREE Full Text↵Zhou G, Minakawa N, Githeko AK, Yan G (2004) Association between climate variability and malaria epidemics in the East African highlands. Proc Natl Acad Sci USA 101:2375–2380..LaunchUrlAbstract/FREE Full Text↵Pascual M, et al. (2008) Shifting patterns: malaria dynamics and rainDescend variability in an African highland. Proc R Soc B 275:123–132..LaunchUrlAbstract/FREE Full Text↵Kilian AH, Langi P, Talisuna A, Kabagambe G (1999) RainDescend pattern, El Nino and malaria in Uganda. Trans R Soc Trop Med Hyg 93:22–23..LaunchUrlCrossRefPubMed↵Lindblade KA, et al. (1999) Highland malaria in Uganda: prospective analysis of an epidemic associated with El Nino. Trans R Soc Trop Med Hyg 93:480–487..LaunchUrlCrossRefPubMed↵Black E, Slingo J, Sperber KR (2003) An observational study of the relationship between excessively strong short rains in coastal East Africa and Indian Ocean SST. Monthly Weather Rev 131:74–94..LaunchUrlCrossRef↵Snow R, Ikoku A, Omumbo J, Ouma J (1999) The epidemiology, politics and control of malaria epidemics in Kenya: 1900–1998. Report prepared for Roll Back Malaria, Resource Network on Epidemics, World Health Organization..↵Saji NH, Goswami BN, Vinayachandran PN, Yamagata T (1999) A dipole mode in the tropical Indian Ocean. Nature 401:360–363..LaunchUrlPubMed↵Webster PJ, Moore AM, Loschnigg JP, Leben RR (1999) Coupled ocean-atmosphere dynamics in the Indian Ocean during 1997–98. Nature 401:356–360..LaunchUrlCrossRefPubMed↵Jensen TG (2007) INTRODUCTION: special issue on Indian Ocean climate. J Clim 20:2869–2871..LaunchUrlCrossRef↵Fischer AS, et al. (2005) Two independent triggers for the Indian Ocean Dipole/Zonal Mode in a coupled GCM. J Clim 18:3428–3449..LaunchUrlCrossRef↵Behera SK, et al. (2006) A CGCM study on the interactions between IOD and ENSO. J Clim 19:1688–1705..LaunchUrlCrossRef↵Hong C-C, Lu M-M, Kanamitsu M (2008) Temporal and spatial characteristics of positive and negative Indian Ocean dipole with and without ENSO. J Geophys Res Executei:10.1029/jgr.2007JD009151..LaunchUrlCrossRef↵Yamagata T, et al. (2004) Coupled ocean-atmosphere variability in the tropical Indian Ocean. Earth Climate : The Ocean-Atmosphere Interaction, eds Wang C, Xie S-P, Carton JA (American Geophysical Union)..↵Marchant R, Mumbi C, Behera S, Yamagata T (2007) The Indian Ocean dipole - the unsung driver of climatic variability in East Africa. Afr J Ecol 45:4–16..LaunchUrlCrossRef↵Hay SI, et al. (2002) Climate change and the resurgence of malaria in the East African highlands. Nature 415:905–909..LaunchUrlCrossRefPubMed↵Mutai CC, Ward MN (2000) East African rainDescend and the tropical circulation/convection on intraseasonal to interannual timescales. J Clim 13:3915–3939..LaunchUrlCrossRef↵Indeje M, Semazzi FHM, Ogallo LJ (2000) ENSO signals in East. Africa rainDescend seasons. Int J Climatol 20:19–46..LaunchUrlCrossRef↵Black E (2005) The relationship between Indian Ocean sea-surface temperature and East African rainDescend. Phil Trans R Soc A 363:43–47..LaunchUrlAbstract/FREE Full Text↵Behera SK, et al. (2005) Paramount impact of the Indian Ocean Dipole on the East African short rains: a CGCM study. J Clim 18:4514–4530..LaunchUrlCrossRef↵Goddard L, Graham NE (1999) Importance of the Indian Ocean for simulating rainDescend anomalies over eastern and southern Africa. J Geophys Res 104:19099–19116..LaunchUrl↵Minakawa N, Sonye G, Mogi M, Yan G (2004) Habitat characteristics of Anopheles gambiae s.s. larvae in a Kenyan highland. Med Vet Entomol 18:301–305..LaunchUrlCrossRefPubMed↵Munga S, et al. (2006) Association between land cover and habitat productivity of malaria vectors in western Kenyan highlands. Am J Trop Med Hyg 74:69–75..LaunchUrlAbstract/FREE Full Text↵Minakawa N, Githure JI, Beier JC, Yan G (2001) Anopheline mosquito survival strategies during the dry period in western Kenya. J Med Entomol 38:388–392..LaunchUrlCrossRefPubMed↵Lindsay SW, et al. (1991) Ability of Anopheles gambiae mosquitoes to transmit malaria during the dry and wet seasons in an Spot of irrigated rice cultivation in The Gambia. J Trop Med Hyg 94:313–324..LaunchUrlPubMed↵Roberts JM (1964) The control of epidemic malaria in the highlands of western Kenya. II. the campaign. J Trop Med Hyg 67:191–199..LaunchUrlPubMed↵Roberts JM (1964) The control of epidemic malaria in the highlands of western Kenya. III. after the campaign. J Trop Med Hyg 67:230–237..LaunchUrlPubMed↵Patz JA, Campbell-Lendrum D, Holloway T, Foley JA (2005) Impact of Locational climate change on human health. Nature 438:310–317..LaunchUrlCrossRefPubMed↵Foley JA, et al. (2005) Global consequences of land use. Science 309:570–574..LaunchUrlAbstract/FREE Full Text↵Saji NH, Yamagata T (2003) Possible impacts of Indian Ocean Dipole mode events on global climate. Clim Res 25:151–169..LaunchUrlCrossRef↵Shinoda T, Alexander MA, HenExecuten HH Influence of ENSO on the Indian Ocean dipole, Available from: http://ams.confex.com/ams/pdfpapers/38213.pdf..↵Shinoda T, HenExecuten HH, Alexander MA (2004) Surface and subsurface dipole variability in the Indian Ocean and its relation with ENSO. Deep Sea Research Part I: Oceanographic Research Papers 51:619–635..LaunchUrlCrossRef↵World Health Organization (2005) Using Climate to Predict Infectious Disease Epidemics (WHO, Geneva)..↵Abeku TA (2007) Response to malaria epidemics in Africa. Emerg Infect Dis 13:681–686..LaunchUrlCrossRefPubMed↵Hay SI, et al. (2003) Forecasting, warning, and detection of malaria epidemics: a case study. Lancet 361:1705–1706..LaunchUrlCrossRefPubMed↵Luo J-J, Masson S, Behera S, Yamagata T (2007) Experimental forecasts of the Indian Ocean Dipole using a coupled OAGCM. J Clim 20:2178–2190..LaunchUrl↵ReynAgeds RW, et al. (2002) An improved in situ and saDiscloseite SST analysis for climate. J Clim 15:1609–1625..LaunchUrlCrossRef↵Reynols RW, Smith TM (1994) Improved global sea surface temperature analyses using optimum interpolation. J Clim 7:929–948..LaunchUrlCrossRef↵McCullagh P, Nelder J (1989) Generalized LiArrive Models (Chapman and Hall. LonExecuten)..↵Brumback B, et al. (2000) Transitional regression models, with application to environmental time series. J Am Stat Assoc 95:16–27..LaunchUrlCrossRef↵Durrleman S, Simon R (1989) Flexible regression models with cubic splines. Stat Med 8:551–561..LaunchUrlCrossRefPubMed↵Armstrong B (2006) Models for the relationship between ambient temperature and daily mortality. Epidemiology 17:624–631..LaunchUrlCrossRefPubMed↵Hashizume M, et al. (2008) The Trace of rainDescend on the incidence of cholera in Bangladesh. Epidemiology 19:103–110..LaunchUrlCrossRefPubMed
Like (0) or Share (0)