Cholera Dynamics and El Niño-Southern Oscillation

See allHide authors and affiliations

Science  08 Sep 2000:
Vol. 289, Issue 5485, pp. 1766-1769
DOI: 10.1126/science.289.5485.1766


Analysis of a monthly 18-year cholera time series from Bangladesh shows that the temporal variability of cholera exhibits an interannual component at the dominant frequency of El Niño–Southern Oscillation (ENSO). Results from nonlinear time series analysis support a role for both ENSO and previous disease levels in the dynamics of cholera. Cholera patterns are linked to the previously described changes in the atmospheric circulation of south Asia and, consistent with these changes, to regional temperature anomalies.

Cholera remains a major public health problem in many areas of the world, including Bangladesh and India. A climate influence on cholera has long been debated (1), and it has been suggested that ENSO, a major source of interannual climate variability, drives the interannual variation of the disease (2, 3). For example, cholera reappeared in Peru with the El Niño event of 1991–92 and seems to fluctuate seasonally in Bangladesh with sea surface temperature (SST) in the Bay of Bengal (2, 4). Recent studies of time series for diarrhoeal diseases in Peruvian children have shown an increase in cases associated with warmer temperatures and the 1997–98 El Niño (5,6). Vibrio cholerae, the bacterium that causes the disease, is now known to inhabit brackish waters and estuarine systems (2) and thus might be sensitive to climate patterns. Here we examine the associations between cholera and ENSO and between cholera and climate at interannual time scales, using an 18-year record from Bangladesh where the disease is endemic. A nonlinear time series approach allows us to consider different hypotheses for the roles of environmental driving variables and the inherent disease dynamics in producing the interannual variability of cholera.

The disease data consist of a monthly time series for cholera incidence between January 1980 and March 1998 in Dhaka, Bangladesh (Fig. 1A). Over the same time span, the monthly SST anomaly in a region of the equatorial Pacific provides an index for ENSO (Fig. 1B). The cholera time series displays the well-known seasonal variation of the disease—typically described as bimodal, with a small peak in the spring and a larger one in the fall or early winter—but also shows a multiyear modulation of the seasonal cycles. The interannual variability of cholera cases has a dominant frequency of 1/3.7 years, as shown by singular spectrum analysis (7, 8) (Fig. 2). The same dominant frequency is found for the ENSO time series, which suggests that climate variability acts as a driver in the dynamics of the disease (Fig. 2). Alternatively, however, this low-frequency variability could arise solely from the seasonal forcing of disease transmission (9). To investigate the role of ENSO in light of this alternative explanation, we consider a nonlinear time series approach that allows us to compare specific alternative hypotheses for the underlying factors in cholera dynamics. Because the null (non-ENSO) hypothesis is a nonlinear interaction between seasonality and cholera dynamics, the use of standard linear time series models would strongly bias the comparison in favor of the ENSO alternative.

Figure 1

(A) Time series of percentage of cholera cases obtained by the ICDDR,B (International Centre for Diarrhoeal Disease Research, Bangladesh, in Dhaka) from a systematic sample of the patients visiting the facility. (B) ENSO time series given by the index Niño 3.4.

Figure 2

Singular spectrum analysis applied to the ENSO index and to the cholera time series with seasonality removed. For both time series, a window length of 60 months was used. The dominant eigenvalues for both data sets are given by a pair, indicating the existence of an underlying oscillatory component. The eigenspectrum of nonlinear signals typically exhibit three distinct regions: first, a group of significant eigenvalues, used here to obtain the reconstructed components (A); then an intermediate slope; and finally the noise floor. Error bars were calculated from the inverse of the lag-one AR coefficient with conservative weights according to (23, 24). Projections of the signal onto the empirical orthogonal functions (EOFs) of the corresponding eigenvalues give the different principal components (PCs). To reconstruct the time series corresponding to several EOFs, the associated PCs are combined. This reconstruction process preserves the phase of the time series. In addition, no information is lost in the reconstruction, because the sum of all the individual reconstructed components gives the original time series. Here the four first components were used, which explain 46 and 76% of the variance for the cholera and ENSO time series, respectively. In (B), spectra were computed with the maximum entropy method (MEM) applied to the reconstructed components in (A), respectively. The dominant peak in both spectra corresponds to a frequency of 0.0225 months (or 1/3.7 years). Results are robust to varying resolutions and MEM orders. Similar results were obtained with the multi-taper method (MTM). The peak corresponding to 0.0225 months was significant at a 99% level and proved robust to varying resolutions and tapers.

Lacking information that could be used to specify a valid mechanistic model for the ENSO effect, we use time series models that are both nonlinear and nonparametric and are effective at modeling high-dimensional relationships. The dynamics of a variable of interest,Nt , a measure of cholera levels, are modeled with a nonlinear equation of the formEmbedded Image Embedded Image(1)where T p is a prediction time,f is a nonlinear function, andEt is the environmental forcing under consideration (10, 11). The sin and cos functions implement a seasonal clock and et represents the IID random noise variables. The parameters τ, τf, and d denote, respectively, two different time lags and the number of time delay variables. Time delay coordinates are used in the model as surrogates for unobserved variables influencing the endogenous dynamics of the disease, such as the fraction of susceptible individuals in the population (12, 13). The functional form off is not specified in a rigid form. Instead, the shape off is determined by the data, using an objective model selection criterion: generalized cross-validation (GCV) (14). We used the GCV criterion to compare models with and without seasonality and with and without the environmental covariate Et (Table 1). The selected model is low-dimensional and incorporates both seasonality and ENSO as external forcings (Fig. 3). The model also accounts for the largest fraction of the variance in the data.

Figure 3

The (square-root–transformed) cholera data (black line) and the 2-months-ahead prediction of the fitted model incorporating both seasonality and ENSO at a lag τf = 11 (red line).

Table 1

Comparison of models using the GCV model selection criterion V c. The model (defined by Eqs. 1and 2) is fitted to the cholera time series after the data are square-root–transformed, which normalizes the residuals and stabilizes the variance. A time lag t of 2 months is used, based on the well-known rule of thumb of choosing the lag for which the autocorrelation function first crosses 0.5. Thus, Eq. 1 takes the general formEmbedded Image Embedded Image(4)where Yt =Embedded Image. When seasonality is incorporated but not the ENSO index, low-dimensional models (d = 1 and d = 2) are selected (that is, they have the smallest values of the cross-validation criterionV c). The importance of seasonality is demonstrated by comparing these models to their autonomous counterparts. Models with an equivalent or larger number of independent variables (for example, d = 4 and d = 6) but no seasonality have larger values ofV c. The importance of ENSO is examined by incorporating the ENSO index into the simplest seasonal model (d = 1) at different time lags (τf between 0 and 12 months). The smallestV c value is obtained for the model with τf = 11. This model is also selected over the seasonal one with an equal number of independent variables (d = 2). It has a higherr 2 value than the seasonal and autonomous models, accounting for a larger fraction of the variance. (For all models, the number of neurons k used in fitting the function f of Eq. 2 varied between 1 and 5. Only the model with the smallest V c is reported here.) A value of k larger than one in the selected ENSO model indicates that a linear time series model is not adequate to fit the causal relationships and would give a less reliable conclusion about the role of ENSO.

View this table:

Could the better fit (the higherr 2) of the model incorporating ENSO simply result from the larger number of variables and parameters in that model? To address this question, we conducted parametric bootstrap tests for the significance of the improvement in fit (15). The result of these comparisions (Table 2) shows that the addition of the ENSO index as a predictor variable is highly significant. But so is the addition of previous disease levels. These results support a role for both extrinsic factors (ENSO and seasonality) and intrinsic ones (previous disease levels) in the dynamics of cholera.

Table 2

Parametric bootstrap test for the significance of additional predictor variables. Comparisons are made only across the best models previously selected with the GVC criterion. The models with previous disease levels and seasonality as the predictor variables are labeled as “dynamic seasonal,” and the model that in addition incorporates the ENSO index is labeled as “seasonal ENSO” (τf = 11). The mean, SD, and maximum (Max.) value of Δi r 2 are computed from the bootstrap data with n = 1000. We assess first the significance of adding the ENSO index as a predictor variable by considering the dynamic seasonal models as the null hypothesis. In all cases, the addition of ENSO as a predictor variable is highly significant (α = 0.001). We then assess the significance of considering previous disease levels as a predictor variable. We compare the seasonal ENSO model, which does incorporate previous disease levels, to the corresponding model built only from the seasonal clock and the ENSO index, which we label “environmental” (τf = 11, k = 3) . The environmental model performs poorly, withr 2 = 44% (even for an equal number of neurons k). The addition of previous disease levels as a predictor variable is highly significant (α = 0.001). For all comparisons, we also recorded the number of times that the cross-validation criterion failed to select the correct model (the null model). This frequency is given as the probability of Δi V c being negative in the last column. Results show that this frequency is extremely small for all comparisons.

View this table:

A positive influence of previous disease levels indicates a density-dependent effect in cholera transmission. In the fall, when the largest values generally occur, this effect is of the compensatory type: Cholera incidence increases monotonically with previous levels but with an ever decreasing slope. Compensatory density dependence can result from a decrease in the number of susceptibles in the population after a large number of cases and/or from the complexities of transmission through water in the environment. Overcompensatory density dependence, a negative effect of previous levels on current levels, has been observed in other infectious diseases (16) but is not apparent here for a time lag of 2 months. Although we cannot rule out overcompensation for disease levels higher than the ones observed, this difference can result from a less effective immunity in the case of cholera.

A positive effect of ENSO is observed on the 2-month-ahead predictions of cholera incidence in the fall (November, for example). This is the time of the year during which the second and largest outbreak of cholera typically occurs in the bimodal seasonal cycle. Cholera incidence 2 months ahead increases with both SST anomaly and present incidence. The sensitivity to changes in the ENSO index is largest for increasing negative values of the SST anomaly; that is, for changes that anticipate an El Niño, and for increasing and large positive values—changes that occur during an El Niño. Later in the spring (April, for example), the levels of cholera incidence are typically lower. Disease levels can either increase or decrease with the ENSO index. Incidence goes up with increasing negative values of the SST anomaly, as for changes that would anticipate an El Niño. It decreases, however, with increasing positive values of the SST anomaly. One possible interpretation of this pattern is that large peaks in the fall are typically followed by small increases of incidence in the spring (from February into April). This inverse relation is indeed observed in the data. These results indicate that the dynamics of cholera in Bangladesh are consistent with a remote forcing by ENSO.

ENSO is thought to affect the atmospheric circulation in the Indian Ocean and south Asia: Changes in cloud cover and evaporation associated with a weakening of the local Hadley cell increase the heat flux entering this region a few months after warming of the Pacific during an El Niño (17). This provides a possible mechanistic connection between El Niño and regional climate variables potentially having an impact on cholera in Bangladesh. We considered three interrelated climate variables: upper-tropospheric humidity, cloud cover, and top-of-atmosphere absorbed solar radiation. From global satellite–retrieved fields for these variables (18–20), we computed global correlation maps between the cholera time series in Dhaka and the climate time series at the different points of the global grid, with temporal lags ranging from ±12 months. The strongest associations are found for negative lags; that is, for the climate variables leading cholera by 4 to 6 months (Fig. 4). These maps also show coherence with those obtained from correlations of the climate anomaly fields with an ENSO index, instead of cholera, and for a positive lag of 6 months [(21) and figure 3 in (17)]. Thus, the same changes in the atmospheric circulation in south Asia that trail the warming of the Pacific appear to anticipate changes in the interannual variation of cholera in Bangladesh. The time lags in these associations are consistent with the 11-month delay found between cholera and ENSO in our time series analysis.

Figure 4

Maps of correlation coefficients between cholera in Dhaka and anomalies in two satellite-retrieved fields 4 months earlier: (A) upper tropospheric humidity and (B) top-of-atmosphere absorbed solar radiation. The colored regions indicate areas of positive and negative correlations, and the black boundaries indicate the areas where correlations are significant at the 0.002 level. (For each particular point in the map, the significance was calculated with a Monte Carlo test by randomly rearranging the elements of each time series at all points in 999 permutations.)

The above results suggest that an increase in local temperature ultimately mediates the influence of ENSO on cholera. Higher ambient temperatures would correspond to higher water temperatures in shallow bodies of water, such as ponds and rivers in the large estuary of Bangladesh and shallow coastal waters of the Bay of Bengal. Satellite data show that a positive association between cholera and temperature is first observed to the north of Bangladesh over the Himalayas, where temperature leads cholera increases by 6 months (Fig. 5). The pattern then moves south, though it weakens, as the lag to cholera decreases. Ambient temperatures have also been implicated in the dynamics of diarrhoeal diseases and ofV. cholerae in the environment in Peru (5,22), and SSTs have been shown to display a bimodal seasonal cycle similar to that of cholera cases in Bangladesh (2,4).

Figure 5

Maps of correlation coefficients between cholera in Dhaka and temperatures on land and at sea for different lags, with the environmental variable anticipating disease. Monthly temperature data for the period from 1980 to 1995 were extracted from the Global Ocean Surface Temperature Atlas (GOSTAplus) for a grid of 5° latitude and longitude. Data were provided by the NASA Physical Oceanography Distributed Active Archive Center at the Jet Propulsion Laboratory/California Institute of Technology. Black boundaries indicate regions where the correlations are significant at the 0.002 level (see the legend to Fig. 4 for details).

Another mediating factor in the ENSO-cholera relation might be the melting of the snowpack in the Himalayas, through its effect on the monsoons, precipitation, and river discharge. This scenario, which remains to be investigated, is suggested by the strong but reduced pattern appearing to the north of Bangladesh (Fig. 5, first and second panels). Floods and droughts can affect not only human interactions with water resources and therefore exposure to the pathogen, but also sanitary conditions and susceptibility to disease.

  • * To whom correspondence should be addressed. E-mail: mercedes{at}


View Abstract

Navigate This Article