Conceived and designed the experiments: AA KJL. Performed the experiments: AA KJL JLS KMC EWP. Analyzed the data: AA KJL JLS KMC EWP. Contributed reagents/materials/analysis tools: AA JLS KMC KJL CJT JRE JEP KLR. Wrote the paper: AA KJL JLS KMC CJT SCB JRE JEP KLR. Contributed to data processing and graphic design: EWP JLS KMC AA.
Recent clusters of outbreaks of mosquito-borne diseases (Rift Valley fever and chikungunya) in Africa and parts of the Indian Ocean islands illustrate how interannual climate variability influences the changing risk patterns of disease outbreaks. Although Rift Valley fever outbreaks have been known to follow periods of above-normal rainfall, the timing of the outbreak events has largely been unknown. Similarly, there is inadequate knowledge on climate drivers of chikungunya outbreaks. We analyze a variety of climate and satellite-derived vegetation measurements to explain the coupling between patterns of climate variability and disease outbreaks of Rift Valley fever and chikungunya.
We derived a teleconnections map by correlating long-term monthly global precipitation data with the NINO3.4 sea surface temperature (SST) anomaly index. This map identifies regional hot-spots where rainfall variability may have an influence on the ecology of vector borne disease. Among the regions are Eastern and Southern Africa where outbreaks of chikungunya and Rift Valley fever occurred 2004–2009. Chikungunya and Rift Valley fever case locations were mapped to corresponding climate data anomalies to understand associations between specific anomaly patterns in ecological and climate variables and disease outbreak patterns through space and time. From these maps we explored associations among Rift Valley fever disease occurrence locations and cumulative rainfall and vegetation index anomalies. We illustrated the time lag between the driving climate conditions and the timing of the first case of Rift Valley fever. Results showed that reported outbreaks of Rift Valley fever occurred after ∼3–4 months of sustained above-normal rainfall and associated green-up in vegetation, conditions ideal for Rift Valley fever mosquito vectors. For chikungunya we explored associations among surface air temperature, precipitation anomalies, and chikungunya outbreak locations. We found that chikungunya outbreaks occurred under conditions of anomalously high temperatures and drought over Eastern Africa. However, in Southeast Asia, chikungunya outbreaks were negatively correlated (
Extremes in climate conditions forced by the
Interannual climate variability associated with the
Climate fluctuations leading to extreme temperatures, storm surges, flooding, and droughts produce conditions that precipitate mosquito-borne disease epidemics directly affecting global public health. Abnormally high temperatures affect populations of mosquito disease vectors by influencing: mosquito survival; susceptibility of mosquitoes to viruses; mosquito population growth rate, distribution, and seasonality; replication and extrinsic incubation period of a virus in the mosquito; and virus transmission patterns and seasonality
The
Correlation of sea surface temperatures and rainfall anomalies illustrate ENSO teleconnection patterns. There is a tendency for above (below) normal rainfall during
Periodic anomalous warming and cooling in the tropical central to eastern Pacific Ocean region (ENSO) can trigger a tropospheric bridge effect that propagates globally. This propagation is what triggers teleconnections and is postulated to have a lagged response of between 3–5 months. The effect of these teleconnections is that they can trigger anomalous convective activity, or lack of, at huge distances from the original site of warming
At a gross global scale, the
Symbols indicate distribution of recent outbreaks of chikungunya (2004–2006) shown by yellow dots and Rift Valley fever (2006–2009) shown by red, blue and green dots over eastern and southern Africa and the Indian Ocean islands.
To determine the ecological and climatic conditions leading to and associated with Rift Valley fever and chikungunya mosquito-borne disease outbreaks, we analyzed relationships between locations of outbreaks and patterns of change in vegetation, rainfall, and temperature.
The baseline disease case location data used in this study were based on human and livestock epidemiological surveys by different institutions in various countries. These data are limited to outbreaks between 2006–2009 for Rift Valley fever in East Africa, Sudan, and Southern Africa. For chikungunya data, we compiled data from two sources. The first source was data from the recent epidemic in Eastern Africa and Western Indian Ocean islands, covering 2004–2006. The second source was a historical 1952–2010 dataset compiled from various sources, including literature in the online archives of the U S Centers for Disease Control and Prevention (CDC), the WHO, and other relevant sources to create a statistically robust sample for this study. We were primarily concerned with geographic locations of the disease cases in order to compare them with environmental data. At each historical outbreak location an approximate geographic latitude/longitude location was determined using place names. Historical outbreaks were located in East Africa, Central Africa, South Asia (primarily India and Bangladesh), and Southeast Asia, providing a broader historical context in which to analyze chikungunya-climate relationships. Most comprehensive geo-referenced records of chikungunya were limited to the recent epidemic period (2004–2010). Full details of chikungunya and Rift Valley fever case data are provided in sections 3 and 5 of
We utilized a number of environmental and climate data sets described in detail in
| Data | Source | Coverage | Climatology period | Purpose |
| NINO 3.4 SST | NOAA/CPC | 5°N–5°S, 170°W–120°W, monthly | 1971–2000 | teleconnections |
| WIO | NOAA/CPC | 10°N–10°S, 40°–64°E, monthly | 1971–2000 | teleconnections |
| GPCP Rainfall | NOAA/CPC | global, monthly, 1° | 1979–2009 | teleconnections, chikungunya |
| NCEP Air Temperature | NOAA/CPC | global, monthly, 2.5° | 1968–1996 | chikungunya |
| ARC Rainfall | NOAA/CPC | Africa, monthly, 10 km spatial resolution | 1995–2006 | Rift Valley fever, chikungunya |
| SPOT Vegetation AVHRR NDVI | VITONASA/GIMMS | Africa, monthly, 1 km, 8 km spatial resolution | May 1998–April 2008 | Rift Valley fever |
| Disease Data (Rift Valley fever, chikungunya) | CDC-K, WHO, FAO, OIE and various national governments | episodic (Rift Valley fever: 2006–2009, chikungunya: 1952–2010) | N/A | climate-ecology-disease teleconnections and relationships |
All anomaly indices were computed as monthly departures from their respective climatological values (long-term means) defined by the periods shown above. NINO 3.4 SST index was computed by the National Oceanic and Atmospheric Administration Climate Prediction Center (NOAA/CPC) as part of operational ENSO monitoring activities. We computed the WIO index directly from the global SST data based on previous research by Linthicum et al. (1999). SPOT Vegetation data were processed by Vlaamse Instelling voor Technologisch Onderzoek (VITO) in Belgium into 10-day composite data. Monthly composites, long-term means, and anomalies from these data were processed by the NASA/Global Inventory Modeling and Mapping Studies (GIMMS) group.
We mapped disease location data on corresponding NDVI and climate data anomalies in order to understand associations between specific anomaly patterns in ecological and climate variables and disease outbreak patterns through space and time. We further explored the associations by plotting and comparing disease data against cumulative rainfall and vegetation index anomalies to illustrate the lag time between the driving climate conditions and the timing of first case disease occurrence for Rift Valley fever. For chikungunya we further investigated the relationships among surface air temperature, precipitation anomalies, and chikungunya outbreaks through correlation analysis. For Rift Valley fever we further investigated relationships between precipitation and Rift Valley fever outbreaks through logistic regression. Results were interpreted in terms of vector biology and population dynamics. We suggested caveats for non-outbreak years when climate and ecological conditions would indicate an imminent outbreak.
In general, for each climate/environmental variable, we calculated anomalies as follows:
Negative rainfall anomalies correspond to the large-scale regional drought in Eastern Africa during October–December, 2005. Anomalies were calculated with reference to the 1995–2000 long term mean. Epicenters of chikungunya outbreaks during this period are shown by the four open black dots.
Patterns of rainfall anomalies preceding outbreaks of Rift Valley fever in (A) East Africa: September–December, 2006, (C) Sudan: June–September, 2007, and (E) Southern Africa: October, 2007–January, 2008. Each outbreak was preceded by persistent and above-normal rain on the order of +200 mm for a period of ∼2–4 months (
To illustrate the concept of teleconnections globally we calculated monthly rainfall anomalies for the GPCP data set based on 1979–2008 long term means. The rainfall anomalies were then correlated with the NINO3.4 sea surface temperature (SST) anomaly index by calculating Pearson's correlation coefficient over the monthly time series to produce the map shown in
To illustrate the ecoclimatic teleconnection connection patterns in relation to Rift Valley fever outbreaks we plotted a Hovmöller diagram of NDVI anomalies for each outbreak region (East Africa, Sudan, and South Africa) against the NINO3.4 SST anomalies. This diagram shows the spatial and temporal dynamics of NDVI anomalies in relation to temporal dynamics of ENSO (represented by the NINO 3.4 index), Western equatorial Indian Ocean (WIO) SST anomalies, and the Rift Valley fever outbreak patterns. Details of how the Hovmöller diagram was derived are given in section 6 of
Spatial and temporal anomaly patterns in normalized difference vegetation index for selected areas of South Africa (A: 29°E and 32.5°E, averaged from 23°S to 27°S), Sudan (B: 32.5°E and 34°E, averaged from 11°N to 15°N), Tanzania (C: 34°E and 37°E, averaged from 4.5°S to 8.5°S) and Kenya (D: 37°E and 42.5°E, averaged from 2°S to 2°N). Regions were plotted by geographic position west to east and represent areas with dense concentrations of Rift Valley fever cases. NDVI anomalies are depicted as percent departures from the 2002–2008 long-term mean, and show the response of vegetation to variations in rainfall. Periods shaded in green to purple indicate above-normal vegetation conditions associated with above-normal rainfall. Periods of persistent drought or below normal rainfall are shown in shades of yellow to red. Each Rift Valley fever outbreak was preceded by above-normal vegetation conditions resulting from persistent above-normal rainfall in the Horn of Africa and Sudan in 2006–2007. Chikungunya epidemics occurred over East Africa and Indian Ocean islands during the 2005–2006 drought period shown by negative NDVI anomalies from 2005–2006 [D: red boxed area]. Clusters of epidemics/epizootics of Rift Valley fever in East Africa (2006–2007) and Sudan (2007) occurred during the
Since Rift Valley fever outbreaks are known to follow periods of extended above-normal rainfall, we calculated a cumulative rainfall anomaly index based on the ARC data set as follows:
In order to quantify the relationship between rainfall anomaly and the occurrence of Rift Valley fever, we used logistic regression. For each region (East Africa, Sudan, Southern Africa, Madagascar), we calculated the cumulative rainfall anomaly for the four months immediately prior to and including the onset of Rift Valley fever activity. These regional reference periods were as follows: (1) East Africa, September–December, 2006, for the December, 2006 outbreak; (2) Sudan, June–October, 2007, for the October, 2007 outbreak; (3) South Africa, October, 2007–January, 2008, for the January, 2008 outbreak; and (4) Madagascar, December, 2007–March, 2008, for the March, 2008 outbreak. For each outbreak site given in
| Region | coefficients | std error | z value | p(>|z|) | confidence | |
| East Africa (n = 383) | −3.1202 | 0.3050 | −10.231 | <2×10−16 | >99.9% | |
| 2.8096 | 0.3264 | 8.608 | <2×10−16 | >99.9% | ||
| Sudan (n = 257) | −4.6153 | 0.8952 | −5.156 | 2.53×10−7 | >99.9% | |
| 26.5603 | 5.1397 | 5.168 | 2.53×10−7 | >99.9% | ||
| South Africa (n = 185) | −2.022 | 0.260 | −7.777 | 7.40×10−15 | >99.9% | |
| 4.575 | 1.106 | 4.135 | 3.55×10−5 | >99.9% | ||
| Madagascar (n = 65) | −2.0897 | 0.5849 | −3.573 | 0.000353 | >99.9% | |
| −10.6091 | 3.8901 | −2.727 | 0.006387 | 99.9% |
Logistic regression of Rift Valley fever presence/absence on cumulative rainfall anomalies over a 4 month period. For each region, the top row presents results for the intercept term (
For chikungunya, we tested the hypothesis that disease outbreaks were correlated with elevated temperature and/or drought by plotting occurrences against surface air temperature anomalies and precipitation anomalies under two scenarios. In the first scenario, meant to simulate high temperatures and moderate drought, cumulative anomalies were made for a 3-month period prior to and then including the actual month of the case, for a total of 4 months. An example for April would be: January, February, March, and April anomalies aggregated to create the 4-month cumulative anomaly. Using the same method, a second scenario was constructed for a prolonged period of high temperatures and severe drought, except this time for a 7-month total. A frequency of occurrence was obtained by extracting values from each of the two scenarios for each of the chikungunya locations, for each region. For the temperature anomalies, the occurrences were classified as higher-than-normal if the anomalies were >0, or cooler-than-normal if <0. The precipitation anomalies were classified similarly as wetter-than-normal or drought if >0 or <0, respectively. Because each geographic region had an occurrence sample size of ≥30, the binomial test of significance was used. This test assumes a normal distribution, random sampling, and mutually exclusive data. A confidence level of 95% was assumed for the 2-tailed test. If the region's occurrences passed the test of significance, the nature of the relationship was tested by calculating the Pearson's correlation coefficient. Additional details of chikungunya correlation analyses are provided in sections 2 and 5 of
Frequency distributions of chikungunya outbreak events and 4-month cumulative temperature anomalies for East Africa (A), Central Africa (B), South Asia (C), and Southeast Asia (D). The 4-month anomaly threshold was used to represent periods of either cool temperatures or drought and extreme high temperatures The dashed line at zero depicts the 1979–2009 long-term mean temperature, with warmer-than-normal temperatures shown to the right (red) and cooler-than-normal temperatures shown to the left (blue) of the line. Cases shown to the right of the dashed line occurred during periods of elevated temperature with a persistence of 4 months.
Frequency distributions of chikungunya outbreak events and 4-month cumulative precipitation anomalies for East Africa (A), Central Africa (B), South Asia (C), and Southeast Asia (D). The 4-month anomaly threshold was used to represent periods of either persistent above-normal rainfall/wetness or persistent drought conditions. The dashed line at zero depicts the (1979–2009) long term rainfall, with greater-than-normal precipitation shown to the right (blue) and lower-than-normal precipitation shown to the left (red) of the line. Cases shown to the left of the dashed line occurred during periods of drought with a persistence of 4 months.
The map in
The 2004–2009 period of analysis is an illustration of the above teleconnection patterns and demonstrates how different climate and ecological anomaly extremes resulting from these teleconnections influence vector-borne disease outbreaks through variations in temperature, rainfall, and ecology. The clusters of outbreaks and epidemics/epizootics of chikungunya and Rift Valley fever across Africa (
Historically, chikungunya is known to be enzootic in west and central Africa in a sylvatic cycle involving wild non-human primates and forest species of
In time the chikungunya outbreak spread and impacted the western Indian Ocean islands including Seychelles, Comoros, Mayotte, Mauritius, and La Reunion, all in 2005, infecting between 30–75% of the populations in affected areas
Cessation of the drought in eastern Africa was marked by the development of contrasting patterns of SST anomalies in the equatorial Indian Ocean, with positive SST anomalies in the equatorial western sector and negative SST anomalies in the eastern sector, and a warm ENSO event in the eastern Pacific Ocean. These patterns of SST anomalies (
The enhancement of
The shift from
An analysis of the January, 1979–February, 2010 relationships between chikungunya activity and surface air temperature anomalies and precipitation anomalies is shown in
The relationship between precipitation and Rift Valley fever was determined through a logistic regression of Rift Valley fever presence/absence on cumulative precipitation anomalies for the 4 months immediately preceding each Rift Valley fever outbreak
In contrast, for Madagascar a negative relationship was found (
We have shown that inter-annual climate variability, as expressed by the ENSO phenomenon in association with regional climatic circulation mechanisms in the equatorial Indian Ocean, had broad influence on two mosquito-borne disease outbreaks over the greater Eastern Africa region, Southern Africa, and western Indian Ocean islands through opposite spatial shifts in precipitation and vegetation anomaly patterns (
Historically, large scale outbreaks of chikungunya have been in large highly populated urban settings of tropical Asia, transmitted by
In South Asia the occurrence of the majority of chikungunya cases during elevated temperatures during both rainy and drought periods strongly suggests that temperature alone was a major driving factor, involving both urban and rural transmission by
We have also illustrated that there was a spatial shift in the area at risk of Rift Valley fever activity from Eastern to Southern Africa in tandem with a phase shift from
Our analyses of the 2006–2009 Rift Valley fever outbreaks confirmed that there was a very close correlation between outbreaks and persistent (i.e., 3–4 months) above-normal rainfall
Our knowledge of teleconnection events and the quasi-cyclical nature of climate variability may allow parts of Africa, the Indian Ocean basin islands, and elsewhere within the greater tropics to have more than a year warning prior to Rift Valley fever outbreaks
Outbreaks of mosquito-borne diseases on epidemic scales, such as those experienced during 2005–2009 in Africa and the western Indian Ocean islands, place a huge burden on public healthcare systems and the economy. Outbreaks of chikungunya are also an impediment to tourism, a major contributor to the gross national product of countries and island nation states in the region. The costs to the economies of East Africa in lost trade in livestock due to Rift Valley fever outbreaks were estimated to be $65 million
It is apparent from our analyses that in changing and variable climate, arboviruses and their mosquito vectors are going to adapt to the existing climatic and ecological conditions in a new region, and the resultant disease transmission will vary accordingly and may not be the same manifestation as observed in the original endemic regions. Combining satellite-derived measurements and analyses of climate and ecology with an understanding of mosquito vector biology and human and animal population immunity status can contribute substantially towards reducing the global burden of vector-borne diseases.
All data analyzed were anonymized. We only used GPS latitude-longitude coordinates for cases, and we did not handle or deal with any human or animal specimens.
(DOC)
Click here for additional data file.
(TIFF)
Click here for additional data file.
(TIFF)
Click here for additional data file.
(TIFF)
Click here for additional data file.
(TIF)
Click here for additional data file.
We acknowledge the contributions of various personnel from the U.S. Army Medical Research Unit – Kenya, U.S. Centers for Disease and Control – Kenya, Kenya Medical Research Foundation, the World Health Organization's Department of Epidemic and Pandemic Alert and Response, the Food and Agricultural Organization of the United Nations, World Organization for Animal Health (OIE), Ministries of Medical Services, Public Health in Kenya, Tanzania and Sudan, Special Pathogens Unit, National Institute for Communicable Diseases and Ministry of Agriculture-Livestock Department, Republic of South Africa, and the Institute Pasteur, Madagascar, for collecting and providing access to georeferenced human and livestock case data used in this study.
The authors have declared that no competing interests exist.
This research is supported in part by funding from the Department of Defense - Armed Forces Health Surveillance Center, Division of GEIS Operations, the United States Department of Agriculture - Agricultural Research Service, Google and Betty-Moore Foundations and the National Aeronautics and Space Administration Grant # NNH08CD31C/ROSES 2007. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.