Search tips
Search criteria 


Logo of interfaceThe Royal Society PublishingInterfaceAboutBrowse by SubjectAlertsFree Trial
J R Soc Interface. 2011 May 6; 8(58): 756–760.
Published online 2010 December 1. doi:  10.1098/rsif.2010.0494
PMCID: PMC3061096

Herald waves of cholera in nineteenth century London


Deaths from cholera in London, UK, were recorded weekly from 1824 to 1901. Three features of the time series stand out: (i) cholera deaths were strongly seasonal, with peak mortality almost always in the summer, (ii) the only non-summer outbreaks occurred in the spring of 1832, the autumn of 1848 and the winter of 1853, and (iii) extraordinarily severe summer outbreaks occurred in 1832, 1849, 1854 and 1866 (the four ‘great’ cholera years). The non-summer outbreaks of 1832, 1848 and 1853 appear to have been herald waves of newly invading cholera strains. In addition, a simple mathematical model confirms that a non-summer introduction of a new cholera strain can result in an initial herald wave, followed by a severe outbreak the following summer. Through the analysis of the genomes of nineteenth-century specimens, it may be possible to identify the strains that caused these herald waves and the well-known cholera epidemics that followed.

Keywords: London cholera, herald wave, waterborne disease model, disease seasonality, John Snow

The birth of modern epidemiology is often attributed to John Snow's famous investigation of the 1854 cholera epidemic in London, and his identification of the Broad Street pump as the most important node in the cholera transmission network [1,2]. More than 150 years later, we still do not know what factors triggered the enormous cholera outbreaks in London in the summers of 1832, 1849, 1854 and 1866. In addition to the intrinsic interest of identifying the mechanisms of historical disease invasions, improving our understanding of cholera specifically is important because it remains a serious public health concern in areas where clean water is unavailable [3]. The recent cholera epidemics in Angola in 2006 [4], Zimbabwe in 2008–2009 [5] and Haiti in 2010 [6] are stark examples.

Previous studies of nineteenth century cholera have focused on the ‘great’ cholera years while paying little attention to the years between severe outbreaks [79]. Here, we consider the great cholera years in the context of London's weekly mortality over the course of the nineteenth century as a whole. Cholera deaths were recorded in the Weekly Returns of the Registrar General's Office from 8 January 1842 to 28 December 1901 (the date on which the last cholera death was reported in London). We digitized these and earlier cholera records from the London Bills of Mortality to obtain a contiguous weekly record of cholera deaths for 77 years from 24 August 1824 to 28 December 1901.

Figure 1 displays the London cholera data in two ways. Figure 1a shows the weekly time series, while Figure 1b shows an intensity plot of week-of-year against year to bring out the pattern of seasonality over the years. The three striking features highlighted above are readily apparent: (i) London's cholera epidemics were strongly seasonal and most intense in the summer, (ii) typical outbreaks were far milder than the devastating outbreaks in the summers of 1832, 1849, 1854 and 1866, and (iii) atypically timed outbreaks occurred in the few months preceding the major outbreaks in 1832, 1849 and 1854. The data clearly separate into regular, mild summer outbreaks, together with outliers corresponding to the non-summer outbreaks and the great cholera years. These features persist when cholera deaths are normalized to account for changes in population size and reporting coverage over the century (see figure 2 and the electronic supplementary material).

Figure 1.
Weekly London cholera deaths from 1824 to 1901. (a) Weekly cholera deaths versus time and (b) cholera deaths versus time of year.
Figure 2.
Weekly cholera deaths, normalized by all-cause mortality and plotted against time of year. Open circles correspond to Bills of Mortality data (years prior to 1842), filled circles to Registrar General data (1842 onwards). Filled light blue circles, 1832 ...

What factors might have been responsible for the unusual timing of the non-summer outbreaks in 1832, 1848 and 1853, and the unusual severity of the outbreaks in the following summers? The most obvious hypothesis is that these features of the London cholera time series resulted from the introduction of new cholera strains into the city at random times of year. It is natural to expect that a non-summer introduction would result in an initial outbreak, with severity tempered by the season [10], followed by a severe outbreak in the summer when environmental conditions most strongly promote cholera transmission. The non-summer outbreaks thus appear to have ‘heralded’ the arrival of new cholera strains in London in 1832, 1848 and 1853. The absense of a ‘herald wave’ before the major epidemic in 1866 may simply reflect the invasion of a new strain coincidentally near the start of the normal cholera season that year.

The idea of herald waves has largely been confined to influenza [11], where it has been discussed extensively in relation to the 1918 [1214] and 2009 [15,16] pandemics. To our knowledge, herald waves have never been described for cholera. Our finding of herald waves for London cholera suggests that herald waves may be a common feature of seasonal diseases.

We note that the identity and origins of new cholera strains in the nineteenth century are not known. Current evolution of new cholera strains appears to be facilitated by horizontal gene transfer among strains in different serogroups ([17]; see the electronic supplementary material). This mechanism might also account for the invasion of antigenically novel strains in the nineteenth century. Cyclical replacement of the predominant cholera serotype has also been observed [18], and may be relevant for London cholera.

We used a simple mathematical model to verify theoretically that a non-summer introduction of a new cholera strain can result in a herald wave followed by a severe outbreak in the summer. The model extends the classical ‘susceptible–infectious–recovered’ (SIR) framework to include a water compartment (W), with transmission occurring through both person–person and person–water–person pathways ([19]; see appendix A for model equations).

As an example, figure 3 compares the London cholera mortality time series from 1845 to 1852 (figure 3a) with an SIWR simulation (figure 3b). Assuming the transmission rate in the water varies seasonally (peaking in the summer, as in figure 3c), simulated cholera time series resemble the typical pattern in London, with mild annual summer outbreaks (as in the period 1845–1848 in figure 3b). The arrow in figure 3b indicates the time at which we introduced a new cholera strain into the SIWR model (the end of September 1848). The population is completely susceptible to the new strain, resulting in an initial outbreak when the strain is introduced. However, owing to the season [10], disease transmissibility is waning at the time of introduction. Eventually transmission decreases past the point where an outbreak can be sustained (and the initial herald wave terminates), but the new cholera strain persists at low levels in the water (figure 3d). When transmissibility from the water rises again the following summer, it triggers an unusually severe epidemic owing to the large number of susceptible individuals (simulation details given in appendix A).

Figure 3.
Model simulations of London cholera, compared with observed cholera deaths. (a) Observed weekly deaths from 1845 to 1852. (b) Simulation results corresponding to endemic cholera, together with the introduction of a new cholera strain in late September ...

While the strain invasion hypothesis is simple—and appealing from the point of view of parsimony—many other explanations are possible, including misdiagnosis involving other diarrhoeal diseases. A direct test of the strain invasion hypothesis would require genetic information from the circulating cholera strains. Relevant tissue samples from patients who died of cholera do exist in museum collections, and recent advances in the recovery and sequencing of DNA as well as the reconstruction of complete genomes from fossil materials [20,21] make sequencing substantial portions of the genomes of nineteenth century cholera strains a realistic goal (see the electronic supplementary material).

Given the fame and historical importance of the four major London cholera epidemics in the nineteenth century, it is surprising that the herald waves we have identified here have not been highlighted previously. Unravelling the mechanisms behind these herald waves will deepen our understanding of the evolutionary and ecological history of this important disease, and in turn help us understand the factors underlying severe cholera outbreaks in modern times. Our study of London cholera also suggests that herald waves may occur for more diseases than has been previously realized, and emphasizes the need for further work examining the relationship between the timing and magnitude of seasonal outbreaks [10]. The systematic digitization of lengthy historical records of disease incidence and mortality will be invaluable for this endeavour. In particular, evidence that herald waves have preceded major epidemics of other diseases may be hidden in untapped historical sources.


Many people deserve thanks for contributing to the acquisition and digitization of London's weekly mortality records, especially Kelly Hancock, Susan Marsh-Rollo, James McDonald and David Richardson. We thank Alison Devault for valuable discussions regarding cholera evolution and archival cholera samples. This project was funded primarily by a grant to D.J.D.E. from the J. S. McDonnell Foundation. D.J.D.E. and H.N.P. also thank NSERC and CIHR.

Appendix A. mathematical model

A.1. Model equations

Our SIWR model for waterborne disease modelling is expressed as a simple system of ordinary differential equations [19],

equation image
A 1

where the host population is divided into susceptible (S), infectious (I) and recovered (R) compartments. The variable W tracks pathogen concentration in a water compartment (e.g. the River Thames and natural wells), and D is the number of individuals killed by the disease. Recovered individuals are immune to further infection. The total host population size is N = S + I + R. The parameter ν is the birth rate, μ is the natural death rate, 1/γ is the mean infectious period and 1/ξ is the mean pathogen lifetime in the water compartment. The parameter α is the pathogen shedding rate into the water, and f is the case fatality proportion. Disease transmission can occur either through person–person contact, with rate parameter bI, or through the water, with rate parameter bW(t). Seasonality in waterborne transmission is modelled using sinusoidal forcing,

equation image
A 2

A.2. Simulation details

Model variables and parameters for system (A 1) are summarized in tables 1 and and2,2, together with initial conditions and parameter values for simulating endemic and introduced cholera (figure 3). The birth and natural death rates were chosen to match London's population growth between 1801 and 1901 ( An expected infectious period of 3 days was used in the model (the typical infectious period is 1–5 days for cholera patients [22]). The ability of Vibrio cholerae to persist outside of human hosts depends upon environmental factors such as salinity [23] and temperature [24]. Under appropriate conditions, V. cholerae can persist for extended periods of time in environmental water sources [25]. Here, we model the expected pathogen lifetime in the water to be of the order of one to two weeks. Case fatality rates for cholera in modern times range from a few per cent to as high as 50 per cent [26], and was fixed at 10 per cent for our model. We set α, the rate at which infected individuals shed pathogen into the water compartment, to 10 cells ml−1 d−1, a value that has been used in previous modelling efforts [27,28]. The transmission parameters were tuned to give reasonable fits to the cholera mortality data. This tuning was accomplished by first locating a periodic orbit for the model when ν = μ, such that this orbit matched the general seasonal pattern of London cholera in ‘typical’ years. Initial conditions for endemic cholera were taken from this periodic orbit. An initial population size of 100 000 was used in the simulations, rather than the population size of the entire city of London, since only a portion of the city's population was at risk for cholera (e.g. John Snow's finding that Vauxhall and Southwark Waterworks customers were at greatly elevated risk of infection [29]).

Table 1.

Variables for system (A 1), together with initial conditions used for simulations of endemic and introduced London cholera (figure 3).

initial conditions
Ssusceptible individuals19 006100 000
Iinfected individuals30
Rrecovered individuals80 9910
Wpathogen concentration in water reservoir (cells ml −1)61014 000
Dindividuals who have died from the disease00
Ntotal population size100 000100 000

Table 2.

Parameters for system (A 1), including values used for simulations of endemic and introduced London cholera (figure 3).

demographic parameters
νbirth rateyr−10.0440.044
μnatural death rateyr−10.0330.033
pathogen parameters
γ−1mean infectious perioddays33
ξ−1mean pathogen lifetime in water reservoirweeks21
fcase fatality proportion0.10.1
contact parameters
bIperson–person contact rateindividuals−1 yr−19.12 × 10−43.65 × 10−4
αrate of pathogen shedding into reservoircells ml−1 yr−1 individuals−136503650
bW(t)reservoir–person contact rateml cells−1 yr−1
 Baverage value of bW(t)ml cells−1 yr−13.9 × 10−52.14 × 10−5
 Aamplitude of seasonality of bW(t)0.50.88
 t1time of maximum seasonal transmissibilityyear0.41 (May 31)0.47 (June 20)
 Tperiod of seasonal forcingyears11


1. Snow J. 1854. The cholera near Golden-square, and at Deptford. Med. Times Gazette 9, 321–322
2. Johnson S. 2006. The ghost map. New York, NY: Riverhead Books
3. Prüss-Üstün A., Bos R., Gore F., Bartram J. 2008. Safer water, better health: costs, benefits and sustainability of interventions to protect and promote health. Geneva, Switzerland: World Health Organization
4. WHO 2007. Cholera annual report 2006. Wkly Epidemiol. Rec. 82, 273–284 [PubMed]
5. WHO 2009. Cholera: global surveillance summary, 2008. Wkly Epidemiol. Rec. 84, 309–324 [PubMed]
6. Enserink M. 2010. Haiti's outbreak is latest in cholera's new global assault. Science 330, 738–739 (doi:10.1126/science.330.6005.738)10.1126/science.330.6005.738 [PubMed] [Cross Ref]
7. Brody H., Rip M., Vinten-Johansen P., Paneth N., Rachman S. 2000. Map-making and myth-making in broad street: the London cholera epidemic, 1854. Lancet 356, 64–68 (doi:10.1016/S0140-6736(00)02442-9)10.1016/S0140-6736(00)02442-9 [PubMed] [Cross Ref]
8. Paneth N., Vinten-Johansen P., Brody H., Rip M. 1998. A rivalry of foulness: official and unoffical investigations of the London cholera epidemic of 1854. Am. J. Public Health 88, 1545–1553 (doi:10.2105/AJPH.88.10.1545)10.2105/AJPH.88.10.1545 [PubMed] [Cross Ref]
9. Rosenberg C. E. 1987. The cholera years: the United States in 1832, 1849, and 1866. Chicago, IL: University of Chicago Press
10. Stone L., Olinsky R., Huppert A. 2007. Seasonal dynamics of recurrent epidemics. Nature 446, 533–536 (doi:10.1038/nature05638)10.1038/nature05638 [PubMed] [Cross Ref]
11. Glezen W. P., Couch R. B., Six H. R. 1982. The influenza herald wave. Am. J. Epidemiol. 116, 589–598 [PubMed]
12. Andreasen V., Viboud C., Simonsen L. 2008. Epidemiologic characterization of the 1918 influenza pandemic summer wave in Copenhagen: implications for pandemic control strategies. J. Infect. Dis. 197, 270–278 (doi:10.1086/524065)10.1086/524065 [PMC free article] [PubMed] [Cross Ref]
13. Chowell G., Viboud C., Simonsen L., Miller M. A., Acuna-Soto R. 2010. Mortality patterns associated with the 1918 influenza pandemic in Mexico: evidence for a spring herald wave and lack of preexisting immunity in older populations. J. Infect. Dis. 202, 567–575 (doi:10.1086/654897)10.1086/654897 [PMC free article] [PubMed] [Cross Ref]
14. Olson D. R., Simonsen L., Edelson P. J., Morse S. S. 2005. Epidemiologic evidence of an early wave of the 1918 influenza pandemic in New York City. Proc. Natl Acad. Sci. USA 102, 11 059–11 063 (doi:10.1073/pnas.0408290102)10.1073/pnas.0408290102 [PubMed] [Cross Ref]
15. CDC 2009. Outbreak of swine-origin influenza A (H1N1) virus infection—Mexico, March–April 2009. Morb. Mortal. Wkly Rep. 58, 467–470 [PubMed]
16. Morens D. M., Taubenberger J. K. 2009. Understanding influenza backward. J. Am. Med. Assoc. 302, 679–680 [PMC free article] [PubMed]
17. Faruque S. M., Sack D. A., Sack R. B., Colwell R. R., Takeda Y., Nair G. B. 2003. Emergence and evolution of Vibrio cholerae O139. Proc. Natl Acad. Sci. USA 100, 1304–1309 (doi:10.1073/pnas.0337468100)10.1073/pnas.0337468100 [PubMed] [Cross Ref]
18. Koelle K., Pascual M., Yunus M. 2006. Serotype cycles in cholera dynamics. Proc. R. Soc. B 273, 2879–2886 (doi:10.1098/rspb.2006.3668)10.1098/rspb.2006.3668 [PMC free article] [PubMed] [Cross Ref]
19. Tien J. H., Earn D. J. D. 2010. Multiple transmission pathways and disease dynamics in a waterborne pathogen model. Bull. Math. Biol. 72, 1506–1533 (doi:10.1007/s11538-010-9507-6)10.1007/s11538-010-9507-6 [PubMed] [Cross Ref]
20. Poinar H. N., et al. 2006. Metagenomics to paleogenomics: large-scale sequencing of mammoth DNA. Science 311, 392–394 (doi:10.1126/science.1123360)10.1126/science.1123360 [PubMed] [Cross Ref]
21. Briggs A. W., et al. 2009. Targeted retrieval and analysis of five Neandertal mtDNA genomes. Science 325, 318–321 (doi:10.1126/science.1174462)10.1126/science.1174462 [PubMed] [Cross Ref]
22. Kaper J. B., Morris J. G., Jr, Levine M. M. 1995. Cholera. Clin. Microbiol. Rev. 8, 48–86 [PMC free article] [PubMed]
23. Miller C. J., Drasar B. S., Feachem R. G. 1982. Cholera and estuarine salinity in Calcutta and London. Lancet 319, 1216–1218 (doi:10.1016/S0140-6736(82)92340-6)10.1016/S0140-6736(82)92340-6 [PubMed] [Cross Ref]
24. Lacey S. W. 1995. Cholera: calamitous past, ominous future. Clin. Infect. Dis. 20, 1409–1419 [PubMed]
25. Rollins D. M., Colwell R. R. 1986. Viable but nonculturable stage of Campylobacter jejuni and its role in survival in the natural aquatic environment. Appl. Environ. Microbiol. 52, 531–538 [PMC free article] [PubMed]
26. Siddique A. K., Salam A., Islam M. S., Akram K., Majumdar R. N., Zaman K., Fronczak N., Laston S. 1995. Why treatment centers failed to prevent cholera deaths among Rwandan refugees in Goma, Zaire. Lancet 345, 359–361 (doi:10.1016/S0140-6736(95)90344-5)10.1016/S0140-6736(95)90344-5 [PubMed] [Cross Ref]
27. Codeco C. T. 2001. Endemic and epidemic dynamics of cholera: the role of the aquatic reservoir. BMC Infect. Dis. 1, 1. (doi:10.1186/1471-2334-1-1)10.1186/1471-2334-1-1 [PMC free article] [PubMed] [Cross Ref]
28. Hartley D. M., Morris J. G., Smith D. L. 2006. Hyperinfectivity: a critical element in the ability of V. cholerae to cause epidemics? PLoS Med. 3, 63–69 (doi:10.1371/journal.pmed.0030007)10.1371/journal.pmed.0030007 [PMC free article] [PubMed] [Cross Ref]
29. Snow J. 1936. Snow on cholera: being a reprint of two papers. New York, NY: The Commonwealth Fund

Articles from Journal of the Royal Society Interface are provided here courtesy of The Royal Society