Research | Open | Published:
A geographic analysis of population density thresholds in the influenza pandemic of 1918–19
International Journal of Health Geographicsvolume 12, Article number: 9 (2013)
Geographic variables play an important role in the study of epidemics. The role of one such variable, population density, in the spread of influenza is controversial. Prior studies have tested for such a role using arbitrary thresholds for population density above or below which places are hypothesized to have higher or lower mortality. The results of such studies are mixed. The objective of this study is to estimate, rather than assume, a threshold level of population density that separates low-density regions from high-density regions on the basis of population loss during an influenza pandemic. We study the case of the influenza pandemic of 1918–19 in India, where over 15 million people died in the short span of less than one year.
Using data from six censuses for 199 districts of India (n=1194), the country with the largest number of deaths from the influenza of 1918–19, we use a sample-splitting method embedded within a population growth model that explicitly quantifies population loss from the pandemic to estimate a threshold level of population density that separates low-density districts from high-density districts.
The results demonstrate a threshold level of population density of 175 people per square mile. A concurrent finding is that districts on the low side of the threshold experienced rates of population loss (3.72%) that were lower than districts on the high side of the threshold (4.69%).
This paper introduces a useful analytic tool to the health geographic literature. It illustrates an application of the tool to demonstrate that it can be useful for pandemic awareness and preparedness efforts. Specifically, it estimates a level of population density above which policies to socially distance, redistribute or quarantine populations are likely to be more effective than they are for areas with population densities that lie below the threshold.
Studying influenza pandemics of the past may offer valuable lessons for preparedness for the next great pandemic [1, 2]. In this paper, we analyze the 1918 influenza pandemic, during which up to 100 million people worldwide and 670,000 people in the U.S.A. are estimated to have died [3–6]. India, where over 15 million perished in the short span of one year, was the single worst-affected country in terms of total mortality [7, 8]. Factors affecting mortality from the pandemic included a lack of immunity to the virus, which resulted in the infection of higher-than-normal numbers of people, its highly contagious nature , and the prior presence of other infections such as tuberculosis or subsequent development of pneumonia as the result of infection [10, 11]. Because influenza viruses spread through human contact, geography and population density in particular are potential factors for transmission and, indirectly, human mortality. The aim of this paper is to analyze the role of population density in the influenza pandemic of 1918. Because the statistics on influenza mortality for India are deeply flawed, following Davis , we estimate population growth trajectories allowing for a break between 1918 and 1919 to capture population loss from the disease, and use a threshold estimation method to test whether low population density districts in India experienced rates of population loss that were different from high density districts.
Epidemiology of the influenza pandemic of 1918–19 in India
The 1918–19 influenza pandemic was one of the worst epidemics in history with an estimated global mortality between 20 and 100 million [9, 12]. The pandemic occurred in two or three waves [13, 14]. The first was a mild wave in the spring of 1918, followed by a second more severe wave in the following autumn that was responsible for the majority of deaths. The third wave was sporadic . The influenza pandemic of 1918–19 in India shared many characteristics with the pandemic in other parts of Asia and the world. In this account, we draw heavily on the Report of the Sanitary Commissioner of India for 1918  and Chandra . The virus is believed to have entered India in the early months of 1918 through the port of Bombay on the west coast. This first wave was relatively mild. The same early and mild first wave was observed in other countries including Indonesia , England, Scotland, and Wales , Portugal and Spain , Mexico , and Peru , and cities including New York City  and Copenhagen . It subsided by August, only to be followed by a second and far more virulent wave that peaked between September and November 1918 in various parts of India. Indonesia , England, Scotland and Wales , Portugal and Spain  and Mexico  experienced a similar pattern of timing and relative severity. Populations that were exposed to the virus early have been shown to have benefited from acquired immunity against the deadly second wave in Denmark , Norway , and Britain [26, 27]. A distinctive characteristic of this epidemic was its disproportionate impact on victims aged between 15 and 35 years [6, 14, 15, 28–30]. This phenomenon was also observed in India .
Epidemics and urban responses in history
It is no coincidence that the public health community worries about new and virulent infectious diseases . As recently as 2009, an episode of pandemic influenza is estimated to have claimed between 151,700 and 575,500 lives worldwide . While vaccinations are frequently the first line of defense against influenza viruses, development of vaccinations for a new strain of virus may take months. In such a situation, only short-term measures including social distancing, and in extreme cases evacuations and quarantines, can protect citizens from a severe epidemic outbreak. The early introduction of social distancing measures, such as school and church closures and banning of mass gatherings, significantly reduced excess mortality during the 1918–19 influenza pandemic . With similar contagious diseases, such as SARS, plague and cholera, quarantines and evacuations have been used to counter severe outbreaks. In the 15th century, for example, the government of Venice combated outbreaks of plague by establishing the “Lazzaretto Vecchio” on a small island off the coast of the Piazza San Marco . More recently, in the late 19th and early 20th centuries, New York City quarantined travelers on Swinburne and Hoffman islands to prevent the spread of cholera . At about the same time, the city established the Metropolitan Board of Health to develop zoning codes to prevent overcrowding in the city and to establish standards for sanitary conditions . With reference to the influenza pandemic of 1918-19 in India, the Sanitary Commissioner of India wrote: “As the striking distance of the influenza virus is probably short the obvious ideal is free ventilation and isolation of sufferers with a view to increase the air space between infected and uninfected” , p.66].
Over the past two decades, sudden virus outbreaks that could have led to widespread human pandemics, including H1N1 , SARS , and H5N1  prompted a series of studies on non-pharmaceutical interventions [40–48], including measures to increase social distance, such as the creation of spatial barriers through quarantine [49, 50], relocating populations to ‘safe’ areas [51, 52], or imposing travel restrictions . Evacuations, a last resort among social distancing measures, are still used; in the aftermath of the Haiti hurricane of November, 2009, and the subsequent earthquake of January 12, 2010, the government ordered the evacuation of the capital, Port-au-Prince, to prevent the spread of epidemic cholera . In sum, it is widely believed that public health interventions, including social distancing measures and the controlled movement of people to either sequester those infected or as a means to lower population density below some critical threshold can significantly decrease the likelihood of a contagious disease spreading. Yet, to date, there is little if any guidance as to what such a population density threshold might be for any disease. Using the influenza pandemic of 1918 as a case, this paper presents an approach to identifying such a threshold value as a guideline for public health policy.
Population density as a factor in influenza population loss
Studies examining the potential relationship between population density and mortality during the 1918 influenza pandemic have produced mixed results. Garrett  found a positive relationship between mortality rates and population densities measured on a state-wide scale in the USA. Once cities were introduced into the equation, normalizing the mortality rate of the cities with those of the states also showed a positive relationship with population density. For Nigeria, crowding contributed to comparatively higher mortality than less-crowded areas, and “there is enough evidence to support the view that large towns suffered more than small and remote villages” . During the 2009–10 influenza pandemic, it was observed that H1N1 infections were sustained over longer periods of time in Taiwanese areas with higher population densities . Theoretical models for influenza and other transmissible respiratory diseases consisting of agent, host, and environment interactions usually require a high host density .
Chowell et al.  also identified a link between population density and mortality for the 1918 pandemic, though it is the opposite of Garrett’s  findings; in Wales and England, low population density in rural areas was positively associated with mortality. On a larger county scale, however, they found no connection between population density or residential crowding and mortality or transmissibility. Supporting this stream of research, Mills et al.  as well as Nishiura and Chowell  could not identify an association between mortality and population size or density, measured as household size. While a review of the literature suggests that the evidence linking higher population densities with higher mortality rates is mixed, intrinsically this relationship makes sense, because influenza viruses spread via human interactions . With rapidly rising population densities around the world, the creation of mega cities, and growing international connectivity, there is, therefore, a dire need for more research on this phenomenon. Given its high population density and rapid urbanization, India is of particular interest for the study of the emergence and spread of viruses posing significant pandemic threats [61, 62].
Despite the interest in the link between population density and influenza morbidity and mortality, little is known about critical turning points or population density thresholds above which the demographic cost of the pandemic may have exceeded that of low-population density areas. Previous studies on the 1918 pandemic have emphasized the size of populations, be they in cities, towns, or rural areas, by selecting arbitrary thresholds on the basis of jurisdictional sizes or reporting boundaries. In some instances the choice of threshold appears to have been made on the basis of convenience (i.e., using categories that may have been originally created by the producers of the data using some criterion other than epidemiology (see Table 1)). Therefore, there is a hitherto unfulfilled need for a deliberate exercise to determine thresholds that focus on the epidemiologic phenomenon at hand. The aim of this paper is to apply a threshold estimation method to identify a population density threshold separating high-density districts in British India from low-density districts. In estimating this threshold, we simultaneously test the hypothesis that the low-density districts so identified differed from high-density districts in terms of population loss. As the world becomes increasingly urbanized, knowledge of how to estimate critical levels of population density above which populations may be at graver risk of contracting or succumbing to influenza than populations in lower density areas can play an important role in fostering pandemic preparedness. The findings of our study contribute directly to the field of spatial epidemiology, which is concerned with “the study of spatial variation in disease risk or incidence” to assist public health decision making , p.328; .
Results and discussion
Estimates from the initial set of models, described below, indicated multiple possible thresholds for population density. The point estimate of the threshold in model (1), described in the methods section below, was 19,067 PPSM (people per square mile), for which only one district, Calcutta, lay above the threshold (see the table in the Appendix). The results from this model suggest that Calcutta, the most densely populated district in British India, with a population density of 35,025 PPSM, is an outlier. Therefore, Calcutta was removed from the dataset and the subsequent analyses were conducted using a dataset of 198 districts and 1,188 observations. The results for datasets containing the Calcutta outlier in the Appendix are broadly consistent with the results presented in the paper.
Table 2 contains the parameter estimates for the models without the Calcutta outlier. These models indicated the presence of two discrete intervals of possible threshold values. Therefore, we present two models (columns 1 and 2) corresponding to the threshold value of population density that minimized the threshold test statistic in each interval (Figure 1). Figure 1 and Column 1 in Table 2 shows the results of the threshold estimation procedure from model (1). The point estimate of the threshold value of population density was 175 PPSM. At the 5% level of significance, a discontinuous set of threshold values that could not be rejected was obtained, corresponding to the intervals 148–209 and 381–464 PPSM (see Figure 1). These ranges are analogous to a 95% confidence interval. The alternate threshold estimates, presented in Column 2 in Table 2, represent values of the threshold that (a) cannot be rejected at the 5% level of significance and (b) yield the minimum threshold test statistic over the alternate range of contiguous possible (in the sense that they cannot be rejected) threshold values within which they occur. In this case, the point estimate is 435 PPSM, which produced the lowest test statistic for the 381–464 PPSM interval (see Figure 1).
An important characteristic of the threshold estimation procedure is the ability to simultaneously test for differences in population loss between below-threshold and above-threshold districts. Interestingly, the difference between low-density and high-density districts is large. For the 175 PPSM threshold, below-threshold districts experienced a population growth rate of −3.72% between 1918 and 1919, while above-threshold districts experienced a growth rate of −4.69%, for a net difference of 0.97%. This difference is, moreover, statistically significant, as denoted by the significance of the parameter estimate for β30. For the 435 PPSM threshold, the corresponding figures are −3.51%, −5.85%, and 2.31% respectively.
In the above models, district-specific intercepts and coefficients on the time trend were also estimated. The intercepts correspond to the logarithm of district-specific population in 1891 and the coefficients on the time trend correspond to the annual rate of population growth. Because of the large number of estimates (198 each), these are not reported in Table 2.
Figures 2 and 3 are a spatial illustration of the relationship between population density and influenza-attributable population loss. Figure 2 is a map of India that contrasts the districts that are above and below the estimated population density threshold. Figure 3 is a map of the districts arranged by quintile of population loss using estimates computed in Chandra et al. . In general, the coastal areas and Gangetic plain of India show coinciding areas of above-threshold population density and high rates of population loss.
The results of this study suggest the presence of population density thresholds that can be used to separate low population loss districts from high population loss districts in British India during the influenza pandemic of 1918. Using Hansen’s  method of threshold estimation, we identified a threshold of 175 PPSM. Below this threshold, districts experienced a decline in population of 3.72%, compared with a 4.69% decline for districts above the threshold, indicating a role for population density in understanding population loss from the epidemic. This evidence is significant in light of the often mixed findings of scholars on the relationship between population density and influenza mortality. The results from the other models presented in this paper suggest even greater differences in population loss between low- and high-density districts, so the chosen point estimates and corresponding difference in this paper should be considered conservative estimates.
While this study advances the literature in a number of ways, it has a number of limitations. The data do not contain information about cases of influenza or mortality that resulted directly from influenza. Therefore, the types of available data preclude analysis of case-mortality, transmission (with or without mortality), and certain other phenomena of epidemiologic significance. It is also not possible to ascribe the entire population loss from the influenza pandemic to mortality alone because of likely depressing effects of the disease on fertility in the immediate aftermath of the pandemic. Even though many studies on the 1918–1919 influenza in New Zealand, the USA and European countries have attempted to factor fertility into either post-pandemic growth rates or depression of conception during the pandemic [68–72], in the context of India, Davis  estimated underreporting of births  by as much as 50% of the true figures. While likely minimal, migration effects also cannot be accounted for due to the lack of appropriate data. In light of these limitations, we have followed the approach of Davis , using changes in population size between 1918 and 1919 that are not explained by the normal population growth trajectory, obtained from relatively accurate population census data, to create a picture of severity of the pandemic in the different districts of India. This paper takes the position that the estimated population loss is an indicator of the severity of the pandemic as a whole, be that loss the direct result of influenza or of other conditions resulting from influenza, including pneumonia, decreased fertility as a consequence of influenza or starvation due to pandemic-influenced famine. The data also do not contain information about socioeconomic status (i.e., poverty, social class or income, or ethnicity) [5, 25, 55, 74, 75] or remoteness , precluding the inclusion of other potential threshold variables that have been shown to be associated with mortality outcomes in other contexts.
The last data limitation described above also raises the issue of the constraints that the methodology places on the analysis. Thresholds may have temporal and spatial variability depending on the above variables and population structure, viral virulence, or transmissibility, though there appears to be little by way of theory or evidence on this subject. In addition, the threshold estimation method used, which has been developed relatively recently, is designed to detect a single threshold based on a single threshold variable (in this case, population density). A more versatile version of the method which allows for multiple thresholds estimated by simultaneously using multiple threshold variables would have been ideal. Finally, to the extent that the conditions in British India were different from conditions encountered in many developed and developing countries today, the results need to be interpreted in their specific context.
The above limitations do not detract conceptually from the utility of the sample-splitting methodology for the identification of demographic or other thresholds for health and potentially other phenomena, and for its value in establishing the critical role for population density in separating high-loss districts from districts that did not suffer as much during the single-worst epidemic of the 20th century in the single-worst-affected country. In this light, this study may be viewed as a prototype on which health geographers can build using the more sophisticated data that are available in modern contexts and as the sample-splitting methodology evolves. In addition, the data used in this study cover a larger area and population than any other single-country study of the pandemic, and thus have the strength of a large sample size.
More generally, the findings presented above have a number of implications for researchers and policy makers in demography, epidemiology, planning, and public health. Most importantly, they introduce a new analytic method, threshold estimation, to the study of epidemics and their effects on populations and population growth. Subject to possible caveats about data and methodology, they also demonstrate that low population density districts in British India may not have suffered as much as high population density districts from the influenza pandemic of 1918–19. At a broader level, the mixed nature of results of studies of this question and their possible connection to contextual factors is an interesting one, and merits further study. The results of this study suggest that plans for pandemic preparedness and adaptation can be informed by the results of studies using this method, especially where high-quality data are available. Threshold estimates can be used to inform the public about location-based risk in times of such epidemics where such risk is found to be present. In addition to introducing a new analytic tool to the study of the geography of health, therefore, it is hoped that this study will be used as a template to inform guidelines for pandemic preparedness issued by public health agencies with a view to minimizing the impact of such events in the future.
Following earlier studies that used census data [7, 17, 77], the data used in this study were obtained from six decennial censuses held in India, for the years 1891, 1901, 1911, 1921, 1931, and 1941 . We focus on these censuses, and not the two censuses preceding 1891, those of 1872 and 1881, on the basis of Davis’  diagnosis that the earlier censuses undercounted the population by over 1%. In addition, these two earlier censuses were conducted using methods that had significantly changed by 1891.
Within these data, we focus on population figures for the districts that were directly ruled by the British India government. We do not use parallel data from the princely states of India, which were nominally under the control of local princes and kings. Data collection by administrative authorities was in some cases significantly different from the British India administration in capacity and function. The coverage of the dataset is broad, encompassing 199 districts for each of the censuses for a total sample size of 1,194 observations. In addition, to ensure comparability over time, in the 1941 census, the population statistics for each census were reported after having been adjusted to conform to the district boundaries as of the 1941 census , providing a convenient dataset on population that is comparable across all the censuses. The district areas used in the computation of population density are, therefore, based on the 1941 boundaries.
In order to estimate population density thresholds that separate low- and high-density districts on the basis of population loss from the pandemic, we use the threshold estimation technique of Hansen  embedded in a population growth model as follows. First, we compute the mean population density of each district across the six censuses as the population of the district divided by the area of the district in square miles to yield persons per square mile (henceforth PPSM). This mean population density provides an estimate of population density at the onset of the pandemic, and is used as the ordering variable for the threshold estimation procedure. It is also highly correlated (r = 0.999) with an alternate estimate of population density computed using the estimate of population in 1918 computed from the standard exponential population growth model allowing for a break in 1918-19 . After ascertaining that the correlation between this mean population density and the dependent variable in the population growth model, namely the log of population, is not significantly different from zero (r = 0.036, P = 0.219), a condition for the threshold estimation procedure, we estimate the population density threshold. The use of the mean density rather than census-specific density ensures that, during the implementation of the sample-splitting algorithm (see below), all six observations for each district lie on the same side of the threshold.
The model used for this exercise was a fixed effects model allowing for each district to have different intercept and time trend terms to absorb heterogeneity in (the log of) population size and population growth, and allowing for a drop in the population between 1918 and 1919, the year of the pandemic. This broad approach was also employed in Davis’  classic study and developed elsewhere [7, 17, 77]. The general model, which follows the approach of the latter three studies, can be expressed as
where LPOP it is the log of population in district i in year t, T t is a time trend, FLU t is a year-specific indicator variable defined as
ε it is a random error term, β 0i , β1i, and β3i, are vectors of district-specific parameters, and β2 is a (fixed) parameter. To this model, for each of the 199 possible values of the threshold level of population density (corresponding to the 199 different districts in the sample), an indicator variable was added to the data such that the variable took on the value 1 if the observation was drawn from a district with a population density lower than the threshold and 0 otherwise. This indicator variable was interacted with the variable of interest, namely the term capturing the drop in population from the influenza pandemic to produce model (1):
Here, I D is the indicator variable corresponding to the threshold population density D. For each of the 199 possible values of D, the above equation was estimated using the 1,194 available observations. The point estimate of the threshold value of D was the one for which the sum of squared errors for the above model was minimized. As a robustness check of the above model, we also estimated models without district-level heterogeneity in the coefficient estimates corresponding to population growth. Thus model (2) was
Appendix: Threshold models for influenza population loss with Calcutta outlier
National Institute of Allergy and Infectious Diseases: "Reconstructed 1918 Influenza Virus Has Yielded Key Insights, Scientists Say.". 2012, Retrieved 6 October, 2012, from http://www.niaid.nih.gov/news/newsreleases/2012/Pages/Reconstructed1918Flu.aspx
Hobday RAP, Cason JWP: The Open-Air Treatment of Pandemic Influenza. Am J Public Health. 2009, 99 (S2): S236-242. 10.2105/AJPH.2008.134627.
Lezzoni L: Influenza 1918: The Worst Epidemic in American History. 1918, New York, NY: TV Books, 1999-
U.S. Department of Health & Human Services: "Pandemic Flu History.". 2012, Retrieved 2012, 7 October, from http://www.flu.gov/pandemic/history/index.html
Mamelund S-E: A socially neutral disease? Individual social class, household wealth and mortality from Spanish influenza in two socially contrasting parishes in Kristiania 1918–1919. Social Science & Medicine. 2006, 62 (4): 923-940. 10.1016/j.socscimed.2005.06.051.
Johnson NPAS, Mueller J: Updating the accounts: global mortality of the 1918-1920" Spanish" influenza pandemic. Bulletin of the History of Medicine. 2002, 76 (1): 105-115. 10.1353/bhm.2002.0022.
Chandra S, Kuljanin G, Wray J: Mortality from the Influenza Pandemic of 1918–1919: The Case of India. Demography. 2012, 49 (3): 857-865. 10.1007/s13524-012-0116-x.
Davis K: The Population of India and Pakistan. 1951, Princeton, NJ: Princeton University Press
Santibañez S, Fiore A, Merlin T, Redd S: A Primer on Strategies for Prevention and Control of Seasonal and Pandemic Influenza. Am J Public Health. 2009, 99 (S2): S216-224. 10.2105/AJPH.2009.164848.
Brundage JF, Shanks GD: Deaths from bacterial pneumonia during 1918–19 influenza pandemic. Emerg Infect Dis. 2008, 14 (8): 1193-1199. 10.3201/eid1408.071313.
Brundage JF, Shanks GD: What really happened during the 1918 influenza pandemic? The importance of bacterial secondary infections. J Infect Dis. 2007, 196 (11): 1717-1718. 10.1086/522355.
Curson P, McCracken K: An Australian perspective of the 1918–1919 influenza pandemic. New South Wales Public Health Bulletin. 2006, 17 (8): 103-107.
Chowell G, Luís MAB, Niall J, Wladimir JA, Cécile V: The 1918–1919 influenza pandemic in England and Wales: spatial patterns in transmissibility and mortality impact. Proceedings of the Royal Society B. 2008, 275 (1634): 501-509. 10.1098/rspb.2007.1477.
Langford C: The age pattern of mortality in the 1918–19 influenza pandemic: an attempted explanation based on data for England and Wales. Medical history. 2002, 46 (1): 1-14.
Patterson KD, Pyle GF: The geography and mortality of the 1918 influenza pandemic. Bulletin of the History of Medicine. 1991, 65 (1): 4-21.
Sanitary Commissioner: Annual report of the Sanitary Commissioner with the Government of India. 1920, Calcutta: Superintendent of Government Printing
Chandra S: Mortality from the influenza pandemic of 1918–19 in Indonesia. Popul Stud (Camb). 2013, [Epub ahead of print]
Burgerlijk Geneeskundigen Dienst: Reports of the Netherlands-Indian Medical Civil Service. Mededeelingen van den Burgerlijken Geneeskundigen dienst in Nederlandsch-Indië . 1912, Batavia: Landsdrukkerij, 1912-1925.
Johnson N: Britain and the 1918–1919 Influenza Pandemic: A Dark Epilogue. 2006, Oxford: Routledge
Erkoreka A: The Spanish influenza pandemic in occidental Europe (1918–1920) and victim age. Influenza and other respiratory viruses. 2010, 4 (2): 81-89. 10.1111/j.1750-2659.2009.00125.x.
Chowell G, Viboud C, Simonsen L, Miller MA, Acuna-Soto R: 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. 2010, 202 (4): 567-575. 10.1086/654897.
Chowell G, Viboud C, Simonsen L, Miller MA, Hurtado J, Soto G, Vargas R, Guzman M, Ulloa M, Munayco C: The 1918–1920 influenza pandemic in Peru. Vaccine. 2011, 29: B21-B26.
Olson DR, Simonsen L, Edelson PJ, Morse SS: Epidemiological evidence of an early wave of the 1918 influenza pandemic in New York City. Proc Natl Acad Sci U S A. 2005, 102 (31): 11059-11063. 10.1073/pnas.0408290102.
Andreasen V, Viboud C, Simonsen L: Epidemiologic characterization of the 1918 influenza pandemic summer wave in Copenhagen: implications for pandemic control strategies. J Infect Dis. 2008, 197 (2): 270-278. 10.1086/524065.
Mamelund SE: Spanish influenza mortality of ethnic minorities in Norway 1918–1919. European Journal of Population/Revue européenne de Démographie. 2003, 19 (1): 83-102. 10.1023/A:1022179025843.
Barry JM, Viboud C, Simonsen L: Cross-protection between successive waves of the 1918–1919 influenza pandemic: epidemiological evidence from US Army camps and from Britain. J Infect Dis. 2008, 198 (10): 1427-1434. 10.1086/592454.
Mathews JD, McBryde ES, McVernon J, Pallaghy PK, McCaw JM: Prior immunity helps to explain wave-like behaviour of pandemic influenza in 1918–9. BMC Infect Dis. 2010, 10 (128): 1-9.
Crosby AW: America's forgotten pandemic: the influenza of 1918. 2003, Cambridge: Cambridge University Press
Pyle GF: The diffusion of influenza: patterns and paradigms. 1986, Totowa, NJ: Rowman & Littlefield Pub Incorporated
Barry J: The Great Influenza: The Epic Story of the Deadliest Plague in History. 2004, New York: Penguin Group
Epstein P: Climate Change and Infectious Disease: Stormy Weather Ahead?. Epidemiology. 2002, 13 (4): 373-375. 10.1097/00001648-200207000-00001.
Dawood F-S, Iuliano A-D, Reed C, Meltzer M-I, Shay D-K, Cheng P-Y, Bandaranayake D, Breiman R-F, Brooks W-A, Buchy P: Estimated global mortality associated with the first 12 months of 2009 pandemic influenza A H1N1 virus circulation: a modelling study.Lancet Infect Dis. 2012, 12 (9): 687-95. 10.1016/S1473-3099(12)70121-4.
Bootsma MCJ, Ferguson NM: The effect of public health measures on the 1918 influenza pandemic in US cities. Proc Natl Acad Sci. 2007, 104 (18): 7588-7593. 10.1073/pnas.0611071104.
Valsecchi M: Mass Plague Graves Found On Venice 'Quarantine’ Island. 2007, Washington, DC: National Geographic News
Seitz S, Miller S: The Other Islands of New York City. 2003, Woodstock, VT: Countryman Press
Sussna S: Bulk control and zoning :the New York city experience. Land economics. 1967, 43 (2): 158-171. 10.2307/3145239.
Alexandrov BB, Hardman J, Janies DA, Studer J, Treseder TW, Valson C, Voronkin IO: Selection for resistance to oseltamivir in seasonal and pandemic H1N1 influenza and widespread co-circulation of the lineages. Int J Heal Geogr. 2010, 9: 13-10.1186/1476-072X-9-13.
Lai PC, Wong CM, Hedley AJ, Lo SV, Leung PY, Kong J, Leung GM: Understanding the Spatial Clustering of Severe Acute Respiratory Syndrome (SARS) in Hong Kong. Environ Heal Perspect. 2004, 112 (15): 1550-1556. 10.1289/ehp.7117.
Peterson AT, Williams RAJ: Ecology and geography of avian influenza (HPAI H5N1) transmission in the Middle East and northeastern Africa. Int J Heal Geogr. 2009, 8: 47-10.1186/1476-072X-8-47.
Anderson RM, May RM: Infectious Diseases of Humans: Dynamics and Control. 1992, Oxford: Oxford Univ. Press
Ferguson NM, Cummings DAT, Cauchemez S, Fraser C: Strategies for containing an emerging influenza pandemic in Southeast Asia. Nature. 2005, 437 (7056): 209-214. 10.1038/nature04017.
Ferguson NM, Cummings DAT, Fraser C, Cajka JC, Cooley PC, Burke DS: Strategies for mitigating an influenza pandemic. Nature. 2006, 442 (7101): 448-452. 10.1038/nature04795.
Eubank S, Guclu H, Anil Kumar VS, Marathe MV: Modelling disease outbreaks in realistic urban social networks. Nature. 2004, 429 (6988): 180-184. 10.1038/nature02541.
Germann TC, Kadau K, Longini IM, Macken CA: Mitigation strategies for pandemic influenza in the United States. Proc Natl Acad Sci U S A. 2006, 103 (15): 5935-5940. 10.1073/pnas.0601266103.
Glass RJ, Glass LM, Beyeler WE, Min HJ: Targeted social distancing design for pandemic influenza. Emerg Infect Dis. 2006, 12: 1671-1681. 10.3201/eid1211.060255.
Halloran ME, Ferguson NM, Eubank S, Longini IM, Cummings DAT, Lewis B, Xu S, Fraser C, Vullikanti A, Germann TC: Modeling targeted layered containment of an influenza pandemic in the United States. Proc Natl Acad Sci U S A. 2008, 105 (12): 4639-10.1073/pnas.0706849105.
Atkinson M, Wein L: Quantifying the Routes of Transmission for Pandemic Influenza. Bull Math Biol. 2008, 70: 820-867. 10.1007/s11538-007-9281-2.
Wein LM, Atkinson MP: Assessing Infection Control Measures for Pandemic Influenza. Risk Analysis. 2009, 29 (7): 949-962. 10.1111/j.1539-6924.2009.01232.x.
Cliff AD, Haggett P: Spatial aspects of epidemic control. Prog Hum Geogr. 1989, 13 (3): 315-347.
Bensimon CM, Upshur REG: Evidence and Effectiveness in Decisionmaking for Quarantine. Am J Public Health. 2007, 97 (1): 44-48. 10.2105/AJPH.2005.077305.
Frenkel S, Western J: Pretext or prophylaxis? Racial segregation and malarial mosquitos in a British tropical colony: Sierra Leone. Ann Assoc Am Geogr. 1988, 78 (2): 211-228. 10.1111/j.1467-8306.1988.tb00203.x.
Curtin PD: Medical knowledge and urban planning in tropical Africa. Am Hist Rev. 1985, 90 (3): 594-613. 10.2307/1860958.
McLafferty S: Placing pandemics: geographical dimensions of vulnerability and spread. Eurasian Geography and Economics. 2010, 51 (2): 143-161. 10.2747/1539-722.214.171.124.
Margesson R: Haiti Earthquake: Crisis and Response. Congressional Research Service. 2010, 7–5700 (R41023): 1-58.
Garrett TA: Economic Effects of the 1918 Influenza Pandemic: Implications for a Modern-day Pandemic. Working paper CA0721. 2007, : Federal Reserve Bank of St. Louis Web site, 2010
Ohadike DC: Diffusion and Physiological Responses to the Influenza Pandemic of 1918–19 in Nigeria. Social Science & Medicine. 1991, 32 (12): 1393-1399. 10.1016/0277-9536(91)90200-V.
Kao C-L, Chan T-C, Tsai C-H, Chu K-Y, Chuang S-F, Lee C-C, Li Z-RT, Wu K-W, Chang L-Y, Shen Y-H: Emerged HA and NA Mutants of the Pandemic Influenza H1N1 Viruses with Increasing Epidemiological Significance in Taipei and Kaohsiung, Taiwan, 2009–10. PLoS One. 2012, 7: e31162-10.1371/journal.pone.0031162.
Spicknall IH, Koopman JS, Nicas M, Pujol JM, Li S, Eisenberg JNS: Informing Optimal Environmental Influenza Interventions: How the Host, Agent, and Environment Alter Dominant Routes of Transmission. PLoS Computational Biology. 2010, 6 (10): 1-11.
Mills CE, Robins JM, Lipsitch M: Transmissibility of 1918 pandemic influenza. Nature. 2004, 432 (7019): 904-906. 10.1038/nature03063.
Nishiura H, Chowell G: Rurality and pandemic influenza: geographic heterogeneity in the risks of infection and death in Kanagawa, Japan (1918–1919). New Zealand Medical Journal. 2008, 121 (1284): 6-10.
Gani SR, Ali ST, Kadi AS: he transmission dynamics of pandemic influenza A/H1N1 2009–2010 in India. Current Science (Bangalore). 2011, 101 (8): 1065-1072.
Kumar S, Quinn SC: Existing health inequalities in India: informing preparedness planning for an influenza pandemic. Health Policy and Planning. 2012, 27 (6): 516-10.1093/heapol/czr075.
Ostfeld RS, Glass GE, Keesing F: Spatial epidemiology: an emerging (or re-emerging) discipline. Trends Ecol Evol. 2005, 20: 328-336. 10.1016/j.tree.2005.03.009.
Linard C, Tatem AJ: Large-scale spatial population databases in infectious disease research. Int J Heal Geogr. 2012, 11: 7-10.1186/1476-072X-11-7.
McSweeny K, Colman A, Fancourt N, Parnell M, Stantiall S, Rice G, Baker M, Wilson N: Was rurality protective in the 1918 influenza pandemic in New Zealand?. New Zealand Medical Journal. 2007, 120 (1256): 2579-
Kolte V, Skinho P, Keiding N, Lynge E: The Spanish flu in Denmark. Scand J Infect Dis. 2008, 40: 538-645. 10.1080/00365540701870903.
Hansen BE: Sample Splitting and Threshold Estimation. Econometrica. 2000, 68 (3): 575-603. 10.1111/1468-0262.00124.
Mamelund SE: Can the Spanish Influenza Pandemic of 1918 Explain the Baby Boom of 1920 in Neutral Norway?. Popul Engl Ed. 2004, 59 (2): 229-260.
Pool D: The effects of the 1918 pandemic of influenza on the Maori population of New Zealand. Bulletin of the History of Medicine. 1973, 47 (3): 273-
Bloom-Feshbach K, Simonsen L, Viboud C, Mølbak K, Miller MA, Gottfredsson M, Andreasen V: Natality decline and miscarriages associated with the 1918 influenza pandemic: the scandinavian and United States experiences. J Infect Dis. 2011, 204 (8): 1157-1164. 10.1093/infdis/jir510.
Bloom-Feshbach K, Simonsen L, Viboud C, Mølbak K, Miller MA, Gottfredsson M, Andreasen V: Reply to Mamelund. J Infect Dis. 2012, 206 (1): 141-143. 10.1093/infdis/jis316.
Mamelund SE: Fertility fluctuations in times of war and pandemic influenza. J Infect Dis. 2012, 206 (1): 140-141. 10.1093/infdis/jis315.
Mills ID: The 1918–1919 Influenza Pandemic – The Indian Experience. The Indian Economic and Social History Review. 1986, XXIII (1): 1-40.
Murray CJL, Lopez AD, Chin B, Feehan D, Hill KH: Estimation of potential global pandemic influenza mortality on the basis of vital registry data from the 1918–20 pandemic: a quantitative analysis. Lancet. 2006, 368 (9554): 2211-2218. 10.1016/S0140-6736(06)69895-4.
Mamelund SE: An egalitarian disease? Socioeconomic status and individual survival of the Spanish Influenza pandemic of 1918–19 in the Norwegian capital of Kristiania. 2004, Oslo: Memorandum, Department of Economics, University of Oslo
Mamelund SE: Geography may explain adult mortality from the 1918–20 influenza pandemic. Epidemics. 2011, 3 (1): 46-60. 10.1016/j.epidem.2011.02.001.
Chandra S: Deaths Associated with Influenza Pandemic of 1918–19, Japan. Emerg Infect Dis. 2013, 19 (4): In press. http://wwwnc.cdc.gov/eid/article/19/4/12–0103_article.htm
Census of India: Census of India. Calcutta: Superintendent: Government Printing, Various years
The data used in this research were collected with support from the National Institute on Drug Abuse (NIDA, grant number 1R21DA025917-01A1 (PI: Chandra)), National Institutes of Health. Its contents are solely the responsibility of the authors and do not necessarily represent the official views of NIDA.
None of the authors has any competing interests.
SC: Project leadership, design and execution of quantitative analysis; writing of entire manuscript. EK-N: Project leadership, writing of background, results, conclusion, and reference sections. GK: Design and execution of quantitative analysis; writing of analytic section. JV: Preparation of maps. All authors read and approved the final manuscript.