Previous Article |
Table of Contents
| Next Article
APPLIED BIOLOGICAL SCIENCES
Climatic suitability for malaria transmission in Africa, 19111995


¶
*Department of Geography, University of Maryland, College Park, MD 20742-8225;
Woods Hole Research Center, P.O. Box 296, Woods Hole, MA 02543-0296;
Trypanosomiasis and Land-Use in Africa Research Group, Department of Zoology, University of Oxford, South Parks Road, Oxford OX1 3PS, United Kingdom; and ¶Kenya Medical Research Institute/Wellcome Trust Collaborative Programme, P.O. Box 43640, 00100 Nairobi GPO, Kenya
Communicated by George M. Woodwell, Woods Hole Research Center, Woods Hole, MA, October 27, 2003 (received for review May 6, 2003)
| Abstract |
|---|
|
|
|---|
Here we examine trends in a climate-driven biological model (19) of malaria transmission for the entire African continent between 1911 and 1995. We specifically address the role of climate change in the African malaria resurgence (12) by using a spatially and temporally extensive gridded climate data-set (ref. 20; Appendix: Methods and Data) to identify locations where the malaria transmission climate suitability index (MTCSI; Appendix: Methods and Data) has changed significantly. In those areas of change, we examine more closely the underlying climate forcing of transmission suitability.
Recent investigations of the influence of climate on malaria transmission in Africa have reached different conclusions. For example, in the East African highlands several current studies (16, 21, 22) found no significant climate trends, whereas others associated the proportion of total admissions caused by malaria at a site with regional temperature anomalies (23) and claim to have shown evidence of regional warming (24). Disagreement remains (24, 25) over the appropriate use of coarse-resolution climate data in drawing inferences at the facility level for which malaria incidence data are available. Further, it has been suggested that even statistically insignificant changes in climate might result in significant changes in malaria transmission potential (24). Concerns have also been raised over the effect of seasonal noise on long-term trend analyses (25). These analyses address these issues. First, scaling concerns are eliminated by performing a comprehensive continental analysis at the original 0.5° latitude and longitude resolution of the climate data. Second, overlooking potentially significant changes in transmission is avoided by examining the MTCSI as opposed to analyzing the climate data directly. Last, there are no issues of seasonal variability confounding the results because the MTCSI is calculated on an annual basis.
| Methods and Data |
|---|
|
|
|---|
The MTCSIt for t = 19111995 was calculated for each of the 10,246 half-degree grid cells covering Africa. Excluding areas that were never suitable (maximum MTCSIt
0.1) or perennially suitable (minimum MTCSIt
0.9) for transmission reduced further analyses to 45% of the continent (4,603 of 10,246 grid cells) for which trends in MTCSI were examined. Autoregressive models were fitted with ordinary least-squares regression to each TS. The presence of trends was assessed with augmented DickeyFuller (refs. 26 and 27; Appendix: Methods and Data) tests, and each series was classified according to behavior (stationary or nonstationary) and the statistical significance of the trend (zero or nonzero). Grid cells classified as either "stationary with trend" or "random walk with drift" contained significant deterministic and stochastic trends in climatic suitability, respectively, whereas those classified as "stationary with no trend" and "random walk with no drift" had no significant trend (Appendix: Methods and Data).
| Results |
|---|
|
|
|---|
|
|
To determine the underlying climate forcings on MTCSI in these specific areas with pronounced trends (Fig. 1D), we further examined the mean annual MTCSI in relation to mean air temperature T and rainfall R during the month most limiting to transmission (Fig. 2). In southern Mozambique (Fig. 2 A), transmission was limited primarily by rainfall. Air temperature, which fluctuated around the 22°C optimum, was the limiting factor in just 14 of the 85 years, mostly before 1940. Precipitation was usually lower than the 80-mm optimum. Rainfall levels were stable from 1945 to 1965 but variable in the earlier and later time periods. The positive MTCSI trend was thus forced by the prevalence of wet, warm years since the mid-1960s.
In the southeastern Ethiopian highlands, MTCSI was always forced by rainfall (Fig. 2B). Mean monthly precipitation was persistently lower than the optimum, with a consistent decline over the time period, resulting in reduced malaria suitability in the region. In contrast, mean MTCSI, temperature, and precipitation for southern Ivory Coast and Ghana (Fig. 2C) suggest transmission in this area was primarily temperature-limited (68 of the 85 years). There was, however, a shift toward precipitation limitation in later years. From 1911 to 1955 there were just 6 rainfall-limited years, but after 1955 there were 11, with 4 of these occurring from 1990 to 1995.
The TS for northwestern Tanzania (Fig. 2D) shows that temperature and rainfall fluctuated around the optimal values, with more temperature-limited, wetter years from 1930 to 1950 and more rainfall-limited, warmer years between 1910 and 1930 and 1950 and 1995. The MTSCI showed a high and increasing interannual variability indicative of a stochastic trend. This region has unstable malaria transmission dominated by variability in the timing and amount of rainfall (28).
The southwestern Congo basin TS (Fig. 2E) showed marked fluctuation in MTCSI in the middle of the time period (19301975) but little variation in the earlier (19101930) or later (19751995) years. This same variability occurred with rainfall, whereas the temperature increased steadily over the time period. The years limited by temperature occurred mostly in the earlier and later periods, when rainfall was high and relatively stable. In contrast, the years with low MTCSI in the middle of the time period were primarily rainfall-limited. This change in variability was likely an artifact of the interpolation method used to generate the climate data (Appendix: Materials and Data), i.e., low variability occurred when rainfall measurements were unavailable within 450 km of the region.
In the Sahel (Fig. 2F), MTCSI was clearly limited by rainfall, with just 3 of the 85 years forced by temperature. The significant negative trend in MTCSI was a result of decreased precipitation over the time period, with a concentration of dry years since the late 1950s.
| Discussion |
|---|
|
|
|---|
| Appendix: Methods and Data |
|---|
|
|
|---|
Additionally, diurnal temperature range records were sparse before 1941 for the entire continent and before 1961 for all but extreme southern Africa. For time periods where no meteorological station data existed within a specified distance (450 km for rainfall), the grid cell value is set to the long-term mean from the TS. Therefore, the rainfall TS in the Congo basin shows no variability for the time periods where observations within 450 km were not available. The lack of diurnal temperature range data is less problematic as we would expect the Tmin over a half-degree grid cell to have little affect on limiting malaria transmission for most of the African continent.
MTCSI Model. The MTCSI is defined by a series of curves y = cos2{[(x U)/(S U)]·(
/2)}, where x is a climate parameter, U is the value of x when conditions are unsuitable, and S is the value of x when conditions are suitable. When S > U, the suitability (1 y) increases with x; when S < U, the suitability y decreases as x increases. The model defines monthly increasing (S = 22°C, U = 18°C) and decreasing (S = 22°C, U = 32°C) curves for T, a monthly increasing (S = 80 mm, U = 0 mm) curve for R, and a single increasing (S = 6°C, U = 4°C) curve for annual minimum temperature. For each month m = 1, 2,..., 12, we calculated the suitabilities yTm and yRm resulting from temperature and rainfall constraints, respectively. Monthly suitability ym was then computed as ym = min(yTm, yRm). For each year, the suitability yTmin because of annual minimum temperature was estimated by using Tmin = min(Tm 0.5DTRm). The suitability index for year t is defined as MTCSIt = min(ymax, yTmin), where ymax = max(ym) persisting for 3 months poleward of 8° north latitude and 5 months elsewhere (19).
Augmented DickeyFuller Test. Each TS yt was initially fit by ordinary least-squares to a pth-order autoregressive model
, where
yn represents the differenced series at a lag of n years;
,
,
, and the
i are constants; yt1 is the series lagged 1 year; and the
t are random shocks. Lag order was increased stepwise from p = 0 to p = 4; p for the accepted model was that which maximized the adjusted r2 statistic for the regression. The standard errors SE(
), SE(
), and SE(
) were estimated for parameters
,
,
. The residuals
t were tested for autocorrelation by using the Q statistic
![]() |
where acfk is the error autocorrelation function at lag k and P represents the maximum number of lags to consider (30). We used P = 24 for each fitted model, based on recommendations of P > 20 lags (31). Autocorrelation was sufficiently represented by the inclusion of four or fewer lags; for the 4,603 grid cells analyzed, only 115 had a Q statistic less than the 5% critical value, and these cells did not exhibit any obvious spatial pattern.
Each series was examined to determine its tendency toward being stationary by using the augmented DickeyFuller (26, 27) test statistic
t =
/SE(
) with null model (
= 1,
= 0) and alternative (
1). If
t exceeded its order-adjusted 5% critical value (32) then the series was stationary, i.e., the effect of a random shock
i diminished over time and the series tended to revert to its mean value,
+
t. In this case we tested for a deterministic trend by comparing the test statistic 
=
/SE(
) to the 5% critical value for the student's t distribution. If 
exceeded the critical value, the series was classified as stationary with trend; otherwise the series was classified as stationary with no trend. For a stationary with trend series, the expected change in MTCSI over t years is
t.
If the null model (
= 1,
= 0) could not be rejected, then the model was refit by ordinary least-squares but with the linear trend term
t omitted. The test statistic in this case was 
=
/SE(
), with null model (
= 0,
= 0,
= 1) and alternative (
0,
= 0,
1). If 
exceeded the 5% critical value (33), then the series was classified as a random walk with drift; otherwise the series was classified as a random walk with no drift. For a random walk with drift
, the expected value of the time series at time t is
t.
| Acknowledgements |
|---|
| Footnotes |
|---|
To whom correspondence should be addressed at: Department of Geography, University of Maryland, 2181 Lefrak Hall, College Park, MD 20742-8225. E-mail: jsmall{at}geog.umd.edu.
| References |
|---|
|
|
|---|
This article has been cited by other articles in HighWire Press-hosted journals:
![]() |
Ying Zhang, Peng Bi, and J. E. Hiller Climate Change and the Transmission of Vector-Borne Diseases: A Review Asia Pac J Public Health, January 1, 2008; 20(1): 64 - 76. [Abstract] [PDF] |
||||
![]() |
S. J. Goetz, A. G. Bunn, G. J. Fiske, and R. A. Houghton Satellite-observed photosynthetic trends across boreal North America associated with climate and fire disturbance PNAS, September 20, 2005; 102(38): 13521 - 13525. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||