The equitable impact of sugary drink taxation structures on sugary drink consumption among Canadians: a modelling study using the 2015 Canadian Community Health Survey-Nutrition.

OBJECTIVE: Estimate the impact of 20 % flat-rate and tiered sugary drink tax structures on the consumption of sugary drinks, sugar-sweetened beverages and 100 % juice by age, sex and socio-economic position. DESIGN: We modelled the impact of price changes - for each tax structure - on the demand for sugary drinks by applying own- and cross-price elasticities to self-report sugary drink consumption measured using single-day 24-h dietary recalls from the cross-sectional, nationally representative 2015 Canadian Community Health Survey-Nutrition. For both 20 % flat-rate and tiered sugary drink tax scenarios, we used linear regression to estimate differences in mean energy intake and proportion of energy intake from sugary drinks by age, sex, education, food security and income. SETTING: Canada. PARTICIPANTS: 19 742 respondents aged 2 and over. RESULTS: In the 20 % flat-rate scenario, we estimated mean energy intake and proportion of daily energy intake from sugary drinks on a given day would be reduced by 29 kcal/d (95 % UI: 18, 41) and 1·3 % (95 % UI: 0·8, 1·8), respectively. Similarly, in the tiered tax scenario, additional small, but meaningful reductions were estimated in mean energy intake (40 kcal/d, 95 % UI: 24, 55) and proportion of daily energy intake (1·8 %, 95 % UI: 1·1, 2·5). Both tax structures reduced, but did not eliminate, inequities in mean energy intake from sugary drinks despite larger consumption reductions in children/adolescents, males and individuals with lower education, food security and income. CONCLUSIONS: Sugary drink taxation, including the additional benefit of taxing 100 % juice, could reduce overall and inequities in mean energy intake from sugary drinks in Canada.

The Impacts of Selling Alcohol in Grocery Stores in Ontario, Canada: A Before-After Study.

OBJECTIVE: From 2015 to 2019, the Government of Ontario expanded privatized sales of alcohol, licensing 450 grocery stores to sell beer, cider, and wine. The impacts of a nearby grocery store gaining an alcohol license on adults' alcohol use in Ontario are examined, including whether impacts differed by gender. METHOD: Data from 2015-2019 Canadian Community Health Survey participants in Ontario (age ≥ 20 years), living within 1,000 m and 1,500 m of grocery stores that gained a license to sell alcohol and propensity-matched controls were included (1,000 m n = 14,052, 1,500 m n = 30,486). Alcohol use outcomes included past-7-day number of standard drinks consumed, near-daily drinking (≥4 days/week), and heavy drinking (5+ drinks in men/4+ in women, at least once/month). Gender-specific difference-in-differences (DiD) analyses compared changes in alcohol use before and after intervention in intervention and control populations. RESULTS: Decreases in past-7-day drinks, near-daily drinking, and heavy drinking were observed after intervention in both intervention and control populations. At the 1,000 m level, adjusted DiD analyses showed past-7-day drinking in women (risk ratio = 1.21, 95% CI [0.88, 1.60]) and heavy drinking in men (odds ratio = 1.38, 95% CI [0.92, 2.08]) had effect sizes above 1, a relative increase over controls, although confidence intervals crossed 1. Findings did not indicate significant differences in alcohol use in intervention relative to controls for other alcohol use measures and at 1,500 m. CONCLUSIONS: Findings suggest no association between a partial alcohol deregulation initiative in Ontario and alcohol use from 2015 to 2019. It is important to monitor the impacts on alcohol use over time as further alcohol deregulation plans in Ontario and other jurisdictions are considered.

PCR Test Positivity and Viral Loads during Three SARS-CoV-2 Viral Waves in Mumbai, India.

SARS-CoV-2 polymerase chain reaction (PCR) tests generally report only binary (positive or negative) outcomes. Quantitative PCR tests can provide epidemiological information on viral transmission patterns in populations. SARS-CoV-2 transmission patterns during India's SARS-CoV-2 viral waves remain largely undocumented. We analyzed 2.7 million real-time PCR testing records collected in Mumbai, a bellwether for other Indian cities. We used the inverse of cycle threshold (Ct) values to determine the community-level viral load. We quantified wave-specific differences by age, sex, and slum population density. Overall, PCR positivity was 3.4% during non-outbreak periods, rising to 23.2% and 42.8% during the original (June-November 2020) and Omicron waves (January 2022), respectively, but was a surprisingly low 9.9% during the Delta wave (March-June 2021; which had the largest increase in COVID deaths). The community-level median Ct values fell and rose ~7-14 days prior to PCR positivity rates. Viral loads were four-fold higher during the Delta and Omicron waves than during non-outbreak months. The Delta wave had high viral loads at older ages, in women, and in areas of higher slum density. During the Omicron wave, differences in viral load by sex and slum density had disappeared, but older adults continued to show a higher viral load. Mumbai's viral waves had markedly high viral loads representing an early signal of the pandemic trajectory. Ct values are practicable monitoring tools.

Mortality Associated with Ambient PM2.5 Exposure in India: Results from the Million Death Study.

BACKGROUND: Studies on the extent to which long-term exposure to ambient particulate matter (PM) with aerodynamic diameter ≤2.5µm (PM2.5) contributes to adult mortality in India are few, despite over 99% of Indians being exposed to levels that the World Health Organization (WHO) considers unsafe. OBJECTIVE: We conducted a retrospective cohort study within the Million Death Study (MDS) to provide the first-ever quantification of national mortality from exposure to PM2.5 in India from 1999 to 2014. METHODS: We calculated relative risks (RRs) by linking a total of ten 3-y intervals of satellite-based estimated PM2.5 exposure to deaths 3 to 5 y later in over 7,400 small villages or urban blocks covering a total population of 6.8 million. We applied using a model-based geostatistical model, adjusted for individual age, sex, and year of death; smoking prevalence, rural/urban residency, area-level female illiteracy, languages, and spatial clustering and unit-level variation. RESULTS: PM2.5 exposure levels increased from 1999 to 2014, particularly in central and eastern India. Among 212,573 deaths at ages 15-69 y, after spatial adjustment, we found a significant RR of 1.09 [95% credible interval (CI): 1.04, 1.14] for stroke deaths per 10-µg/m3 increase in PM2.5 exposure, but no significant excess for deaths from chronic respiratory disease and ischemic heart disease (IHD), all nonaccidental causes, and total mortality (after excluding stroke). Spatial adjustment attenuated the RRs for chronic respiratory disease and IHD but raised those for stroke. The RRs were consistent in various sensitivity analyses with spatial adjustment, including stratifying by levels of solid fuel exposure, by sex, and by age group, addition of climatic variables, and in supplementary case-control analyses using injury deaths as controls. DISCUSSION: Direct epidemiological measurements, despite inherent limitations, yielded associations between mortality and long-term PM2.5 inconsistent with those reported in earlier models used by the WHO to derive estimates of PM2.5 mortality in India. The modest RRs in our study are consistent with near or null mortality effects. They suggest suitable caution in estimating deaths from PM2.5 exposure based on MDS results and even more caution in extrapolating model-based associations of risk derived mostly from high-income countries to India. https://doi.org/10.1289/EHP9538.

Spatio-temporal modelling of malaria mortality in India from 2004 to 2013 from the Million Death Study.

BACKGROUND: India has a substantial burden of malaria, concentrated in specific areas and population groups. Spatio-temporal modelling of deaths due to malaria in India is a critical tool for identifying high-risk groups for effective resource allocation and disease control policy-making, and subsequently for the country's progress towards United Nations 2030 Sustainable Development Goals. METHODS: In this study, a spatio-temporal model with the objective of understanding the spatial distribution of malaria mortality rates and the rate of temporal decline, across the country, has been constructed. A spatio-temporal "random slope" model was used, with malaria risk depending on a spatial relative risk surface and a linear time effect with a spatially-varying coefficient. The models were adjusted for urban/rural status (residence of the deceased) and Normalized Difference Vegetation Index (NDVI), using 2004-13 data from the Million Death Study (MDS) (the most recent data available), with nationwide geographic coverage. Previous studies based on MDS had focused only on aggregated analyses. RESULTS: The rural population had twice the risk of death due to malaria compared to the urban population. Malaria mortality in some of the highest-risk regions, namely the states of Odisha and Jharkhand, are declining faster than other areas; however, the rate of decline was not uniformly correlated with the level of risk. The overall decline was faster after 2010. CONCLUSION: The results suggest a need for increased attention in high-risk rural populations, which already face challenges like inadequate infrastructure, inaccessibility to health care facilities, awareness, and education around malaria mortality and prevalence. It also points to the urgent need to restart the MDS to document changes since 2013, to develop appropriate malaria control measures.

Assessment of SARS-CoV-2 Seropositivity During the First and Second Viral Waves in 2020 and 2021 Among Canadian Adults.

Importance: The incidence of infection during SARS-CoV-2 viral waves, the factors associated with infection, and the durability of antibody responses to infection among Canadian adults remain undocumented. Objective: To assess the cumulative incidence of SARS-CoV-2 infection during the first 2 viral waves in Canada by measuring seropositivity among adults. Design, Setting, and Participants: The Action to Beat Coronavirus study conducted 2 rounds of an online survey about COVID-19 experience and analyzed immunoglobulin G levels based on participant-collected dried blood spots (DBS) to assess the cumulative incidence of SARS-CoV-2 infection during the first and second viral waves in Canada. A sample of 19 994 Canadian adults (aged ≥18 years) was recruited from established members of the Angus Reid Forum, a public polling organization. The study comprised 2 phases (phase 1 from May 1 to September 30, 2020, and phase 2 from December 1, 2020, to March 31, 2021) that generally corresponded to the first (April 1 to July 31, 2020) and second (October 1, 2020, to March 1, 2021) viral waves. Main Outcomes and Measures: SARS-CoV-2 immunoglobulin G seropositivity (using a chemiluminescence assay) by major geographic and demographic variables and correlation with COVID-19 symptom reporting. Results: Among 19 994 adults who completed the online questionnaire in phase 1, the mean (SD) age was 50.9 (15.4) years, and 10 522 participants (51.9%) were female; 2948 participants (14.5%) had self-identified racial and ethnic minority group status, and 1578 participants (8.2%) were self-identified Indigenous Canadians. Among participants in phase 1, 8967 had DBS testing. In phase 2, 14 621 adults completed online questionnaires, and 7102 of those had DBS testing. Of 19 994 adults who completed the online survey in phase 1, fewer had an educational level of some college or less (4747 individuals [33.1%]) compared with the general population in Canada (45.0%). Survey respondents were otherwise representative of the general population, including in prevalence of known risk factors associated with SARS-CoV-2 infection. The cumulative incidence of SARS-CoV-2 infection among unvaccinated adults increased from 1.9% in phase 1 to 6.5% in phase 2. The seropositivity pattern was demographically and geographically heterogeneous during phase 1 but more homogeneous by phase 2 (with a cumulative incidence ranging from 6.4% to 7.0% in most regions). The exception was the Atlantic region, in which cumulative incidence reached only 3.3% (odds ratio [OR] vs Ontario, 0.46; 95% CI, 0.21-1.02). A total of 47 of 188 adults (25.3%) reporting COVID-19 symptoms during phase 2 were seropositive, and the OR of seropositivity for COVID-19 symptoms was 6.15 (95% CI, 2.02-18.69). In phase 2, 94 of 444 seropositive adults (22.2%) reported having no symptoms. Of 134 seropositive adults in phase 1 who were retested in phase 2, 111 individuals (81.8%) remained seropositive. Participants who had a history of diabetes (OR, 0.58; 95% CI, 0.38-0.90) had lower odds of having detectable antibodies in phase 2. Conclusions and Relevance: The Action to Beat Coronavirus study found that the incidence of SARS-CoV-2 infection in Canada was modest until March 2021, and this incidence was lower than the levels of population immunity required to substantially reduce transmission of the virus. Ongoing vaccination efforts remain central to reducing viral transmission and mortality. Assessment of future infection-induced and vaccine-induced immunity is practicable through the use of serial online surveys and participant-collected DBS.

COVID mortality in India: National survey data and health facility deaths.

India's national COVID death totals remain undetermined. Using an independent nationally representative survey of 0.14 million (M) adults, we compared COVID mortality during the 2020 and 2021 viral waves to expected all-cause mortality. COVID constituted 29% (95%CI 28-31%) of deaths from June 2020-July 2021, corresponding to 3.2M (3.1-3.4) deaths, of which 2.7M (2.6-2.9) occurred in April-July 2021 (when COVID doubled all-cause mortality). A sub-survey of 57,000 adults showed similar temporal increases in mortality with COVID and non-COVID deaths peaking similarly. Two government data sources found that, when compared to pre-pandemic periods, all-cause mortality was 27% (23-32%) higher in 0.2M health facilities and 26% (21-31%) higher in civil registration deaths in ten states; both increases occurred mostly in 2021. The analyses find that India's cumulative COVID deaths by September 2021 were 6-7 times higher than reported officially.

Trends in snakebite deaths in India from 2000 to 2019 in a nationally representative mortality study.

The World Health Organization call to halve global snakebite deaths by 2030 will require substantial progress in India. We analyzed 2833 snakebite deaths from 611,483 verbal autopsies in the nationally representative Indian Million Death Study from 2001 to 2014, and conducted a systematic literature review from 2000 to 2019 covering 87,590 snakebites. We estimate that India had 1.2 million snakebite deaths (average 58,000/year) from 2000 to 2019. Nearly half occurred at ages 30-69 years and over a quarter in children < 15 years. Most occurred at home in the rural areas. About 70% occurred in eight higher burden states and half during the rainy season and at low altitude. The risk of an Indian dying from snakebite before age 70 is about 1 in 250, but notably higher in some areas. More crudely, we estimate 1.11-1.77 million bites in 2015, of which 70% showed symptoms of envenomation. Prevention and treatment strategies might substantially reduce snakebite mortality in India.

Seasonal variation and etiologic inferences of childhood pneumonia and diarrhea mortality in India.

Control of pneumonia and diarrhea mortality in India requires understanding of their etiologies. We combined time series analysis of seasonality, climate region, and clinical syndromes from 243,000 verbal autopsies in the nationally representative Million Death Study. Pneumonia mortality at 1 month-14 years was greatest in January (Rate ratio (RR) 1.66, 99% CI 1.51-1.82; versus the April minimum). Higher RRs at 1-11 months suggested respiratory syncytial virus (RSV) etiology. India's humid subtropical region experienced a unique summer pneumonia mortality. Diarrhea mortality peaked in July (RR 1.66, 1.48-1.85) and January (RR 1.37, 1.23-1.48), while deaths with fever and bloody diarrhea (indicating enteroinvasive bacterial etiology) showed little seasonality. Combining mortality at ages 1-59 months with prevalence surveys, we estimate 40,600 pneumonia deaths from Streptococcus pneumoniae, 20,700 from RSV, 12,600 from influenza, and 7200 from Haemophilus influenzae type b and 24,700 diarrheal deaths from rotavirus occurred in 2015. Careful mortality studies can elucidate etiologies and inform vaccine introduction.

Mortality attributable to hot and cold ambient temperatures in India: a nationally representative case-crossover study.

BACKGROUND: Most of the epidemiological studies that have examined the detrimental effects of ambient hot and cold temperatures on human health have been conducted in high-income countries. In India, the limited evidence on temperature and health risks has focused mostly on the effects of heat waves and has mostly been from small scale studies. Here, we quantify heat and cold effects on mortality in India using a nationally representative study of the causes of death and daily temperature data for 2001-2013. METHODS AND FINDINGS: We applied distributed-lag nonlinear models with case-crossover models to assess the effects of heat and cold on all medical causes of death for all ages from birth (n = 411,613) as well as on stroke (n = 19,753), ischaemic heart disease (IHD) (n = 40,003), and respiratory diseases (n = 23,595) among adults aged 30-69. We calculated the attributable risk fractions by mortality cause for extremely cold (0.4 to 13.8°C), moderately cold (13.8°C to cause-specific minimum mortality temperatures), moderately hot (cause-specific minimum mortality temperatures to 34.2°C), and extremely hot temperatures (34.2 to 39.7°C). We further calculated the temperature-attributable deaths using the United Nations' death estimates for India in 2015. Mortality from all medical causes, stroke, and respiratory diseases showed excess risks at moderately cold temperature and hot temperature. For all examined causes, moderately cold temperature was estimated to have higher attributable risks (6.3% [95% empirical confidence interval (eCI) 1.1 to 11.1] for all medical deaths, 27.2% [11.4 to 40.2] for stroke, 9.7% [3.7 to 15.3] for IHD, and 6.5% [3.5 to 9.2] for respiratory diseases) than extremely cold, moderately hot, and extremely hot temperatures. In 2015, 197,000 (121,000 to 259,000) deaths from stroke, IHD, and respiratory diseases at ages 30-69 years were attributable to moderately cold temperature, which was 12- and 42-fold higher than totals from extremely cold and extremely hot temperature, respectively. The main limitation of this study was the coarse spatial resolution of the temperature data, which may mask microclimate effects. CONCLUSIONS: Public health interventions to mitigate temperature effects need to focus not only on extremely hot temperatures but also moderately cold temperatures. Future absolute totals of temperature-related deaths are likely to depend on the large absolute numbers of people exposed to both extremely hot and moderately cold temperatures. Similar large-scale and nationally representative studies are required in other low- and middle-income countries to better understand the impact of future temperature changes on cause-specific mortality.

Renal failure deaths and their risk factors in India 2001-13: nationally representative estimates from the Million Death Study.

BACKGROUND: Renal failure represents a growing but mostly undocumented cause of premature mortality in low-income and middle-income countries. We investigated changes in adult renal failure mortality and its key risk factors in India using the nationally representative Million Death Study. METHODS: In this cross-sectional analysis of population-based data, two trained physicians independently assigned underlying causes to 150 018 deaths at ages 15-69 years from a nationally-representative mortality survey in India for 2001-03 and 2010-13, using the International Classification of Diseases, 10th version (ICD-10). We applied the age-specific proportion of renal failure deaths for the 2010-13 period to the 2015 UN estimates of total deaths in India and calculated age-standardised death rates for renal failure by rural or urban residence, state, and age group. We used proportional mortality of renal deaths (cases) to injuries (controls) to calculate the odds of renal death in the presence of different comorbidities and stratified risks by decade of birth. FINDINGS: In 2001-03, 2·1% of total deaths among 15-69 year olds were from renal failure (1266 [2·2%] of 58 871; unweighted). By 2010-13, the proportion of deaths from renal failure had risen to 2·9% (2943 [3·2%] of 91 147; unweighted) of total deaths and corresponding to 136 000 renal failure deaths (range 108 000-150 000) of 4 688 000 total deaths nationally in 2015. Age-standardised renal death rates were highest in the southern and eastern states, particularly among adults aged 45-69 years in 2010-13. Diabetes, hypertension, and cardiovascular disease were all significantly associated with increased renal failure deaths, with diabetes the strongest predictor-odds ratio (OR) vs control 9·2 (95% CI 6·7-12·7) in 2001-03, rising to 15·1 (12·6-18·1) in 2010-13. In the 2010-13 study population, the diabetes to non-diabetes OR was twice as large in adults born in the 1970s (25·5, 95% CI 17·6-37·1) as in those individuals born during or before the 1950s (11·7, 9·1-14·9). INTERPRETATION: Renal failure is a growing cause of premature death in India. Poorly treated diabetes is the most probable reason for this increase. Strategies aimed at diabetes prevention, and early detection and treatment are urgently needed in India, as well as greater access to renal replacement therapy. FUNDING: US National Institutes of Health, International Development Research Centre, Centre for Global Health Research, University of Toronto.

Deaths from acute abdominal conditions and geographical access to surgical care in India: a nationally representative spatial analysis.

BACKGROUND: Few population-based studies quantify mortality from surgical conditions and relate mortality to access to surgical care in low-income and middle-income countries. METHODS: We linked deaths from acute abdominal conditions within a nationally representative, population-based mortality survey of 1·1 million households in India to nationally representative facility data. We calculated total and age-standardised death rates for acute abdominal conditions. Using 4064 postal codes, we undertook a spatial clustering analysis to compare geographical access to well-resourced government district hospitals (24 h surgical and anaesthesia services, blood bank, critical care beds, basic laboratory, and radiology) in high-mortality or low-mortality clusters from acute abdominal conditions. FINDINGS: 923 (1·1%) of 86,806 study deaths at ages 0-69 years were identified as deaths from acute abdominal conditions, corresponding to 72,000 deaths nationally in 2010 in India. Most deaths occurred at home (71%) and in rural areas (87%). Compared with 567 low-mortality geographical clusters, the 393 high-mortality clusters had a nine times higher age-standardised acute abdominal mortality rate and significantly greater distance to a well-resourced hospital. The odds ratio (OR) of being a high-mortality cluster was 4·4 (99% CI 3·2-6·0) for living 50 km or more from well-resourced district hospitals (rising to an OR of 16·1 [95% CI 7·9-32·8] for >100 km). No such relation was seen for deaths from non-acute surgical conditions (ie, oral, breast, and uterine cancer). INTERPRETATION: Improvements in human and physical resources at existing government hospitals are needed to reduce deaths from acute abdominal conditions in India. Full access to well-resourced hospitals within 50 km by all of India's population could have avoided about 50,000 deaths from acute abdominal conditions, and probably more from other emergency surgical conditions. FUNDING: Bill & Melinda Gates Foundation, Dalla Lana School of Public Health, Canadian Institute of Health Research.

Deaths from acute abdominal conditions and geographic access to surgical care in India: a nationally representative population-based spatial analysis.

BACKGROUND: Acute abdominal conditions have high case-fatality rates in the absence of timely surgical care. In India, and many other low-income and middle-income countries, few population-based studies have quantified mortality from surgical conditions and related mortality to access to surgical care. We aimed to describe the spatial and socioeconomic distributions of deaths from acute abdomen (DAA) in India and to quantify potential access to surgical facilities in relation to such deaths. METHODS: We examined deaths from acute abdominal conditions within a nationally representative, population-based mortality survey of 1·1 million Indian households and linked these to nationally representative facility data. Spatial clustering of deaths from acute abdominal conditions was calculated with the Getis-Ord Gi* statistic from about 4000 postal codes. We compared high or low acute abdominal mortality clusters for their geographic access to well-resourced surgical care (24 h surgical and anaesthesia services, blood bank, critical care beds, basic laboratory, and radiology). FINDINGS: 923 (1·1%) of 86 806 study deaths in those aged 0-69 years were identified as deaths from acute abdominal conditions, corresponding to an estimated 72 000 deaths nationally in India in 2010. Most deaths occurred at home (71%), in rural areas (87%), and were caused by peptic ulcer disease (79%). There was wide variation in rates of deaths from acute abdominal conditions. We identified 393 high-mortality geographic clusters and 567 low-mortality clusters. High-mortality clusters of acute abdominal conditions were located significantly further from well-resourced hospitals than were low-mortality clusters. The odds ratio of a postal code area being a high-mortality cluster was 4·4 (99% CI 3·2-6·0) for living 50 km or more from well-resourced district hospitals (rising to an OR of 16·1 for >100 km), after adjustment for socioeconomic status and caste. INTERPRETATION: Improvements in human and physical resources at existing public hospitals are required to reduce deaths from acute abdominal conditions in India. Had all of the Indian population had access to well-resourced hospitals within 50 km, more than 50 000 deaths from acute abdominal conditions could have been averted in 2010, and likely more from other emergency surgical conditions. Our geocoded facility data were limited to public district hospitals. However, noting the high rate of catastrophic health expenditures in India, we chose to focus on publicly provided services which are the only option usually available to the poor. FUNDING: The Bill & Melinda Gates Foundation, Dalla Lana School of Public Health, and Canadian Institute of Health Research.

Geospatial analysis on the distributions of tobacco smoking and alcohol drinking in India.

BACKGROUND: Tobacco smoking and binge alcohol drinking are two of the leading risk factors for premature mortality worldwide. In India, studies have examined the geographic distributions of tobacco smoking and alcohol drinking only at the state-level; sub-state variations and the spatial association between the two consumptions are poorly understood. METHODOLOGY: We used data from the Special Fertility and Mortality Survey conducted in 1998 to examine the geographic distributions of tobacco smoking and alcohol drinking at the district and postal code levels. We used kriging interpolation to generate smoking and drinking distributions at the postal code level. We also examined spatial autocorrelations and identified spatial clusters of high and low prevalence of smoking and drinking. Finally, we used bivariate analyses to examine the spatial correlations between smoking and drinking, and between cigarette and bidi smoking. RESULTS: There was a high prevalence of any smoking in the central and northeastern states, and a high prevalence of any drinking in Himachal Pradesh, Arunachal Pradesh, and eastern Madhya Pradesh. Spatial clusters of early smoking (started smoking before age 20) were identified in the central states. Cigarette and bidi smoking showed distinctly different geographic patterns, with high levels of cigarette smoking in the northeastern states and high levels of bidi smoking in the central states. The geographic pattern of bidi smoking was similar to early smoking. Cigarette smoking was spatially associated with any drinking. Smoking prevalences in 1998 were correlated with prevalences in 2004 at the district level and 2010 at the state level. CONCLUSION: These results along with earlier evidence on the complementarities between tobacco smoking and alcohol drinking suggest that local public health action on smoking might also help to reduce alcohol consumption, and vice versa. Surveys that properly represent tobacco and alcohol consumptions at the district level are recommended.