Globally, excess deaths during 2020–21 outnumbered documented COVID-19 deaths by 9.5 million, primarily driven by deaths in low- and middle-income countries (LMICs) with limited vital surveillance. Here we unravel the contributions of probable COVID-19 deaths from other changes in mortality related to pandemic control measures using medically-certified death registrations from Madurai, India—an urban center with well-functioning vital surveillance. Between March, 2020 and July, 2021, all-cause deaths in Madurai exceeded expected levels by 30% (95% confidence interval: 27–33%). Although driven by deaths attributed to cardiovascular or cerebrovascular conditions, diabetes, senility, and other uncategorized causes, increases in these attributions were restricted to medically-unsupervised deaths, and aligned with surges in confirmed or attributed COVID-19 mortality, likely reflecting mortality among unconfirmed COVID-19 cases. Implementation of lockdown measures was associated with a 7% (0–13%) reduction in all-cause mortality, driven by reductions in deaths attributed to injuries, infectious diseases and maternal conditions, and cirrhosis and other liver conditions, respectively, but offset by a doubling in cancer deaths. Our findings help to account for gaps between documented COVID-19 mortality and excess all-cause mortality during the pandemic in an LMIC setting.
A reported 5.4 million people died of COVID-19 in 2020 and 2021 globally. However, this number is believed to be vastly lower than true mortality attributable to SARS-CoV-2. Excess mortality—the difference between the total number of deaths that have occurred and the number of deaths that would have been expected in the absence of the pandemic—has been used extensively to quantify the direct and indirect impacts of the pandemic 1 . According to estimates by the World Health Organization, the number of deaths occurring globally in 2020–21 exceeded pre-pandemic expectations by 14.9 million, with 86% of this burden occurring in low- and middle-income countries (LMICs) that reported smaller proportions of all confirmed COVID-19 cases and deaths 2 . In India alone, roughly 3.2–6.5 million excess deaths are estimated to have occurred through the first two waves of the COVID-19 pandemic 3,4 , a figure unmatched by any other demographic event in the country’s history since Independence in 1947 5,6,7,8 .
While understanding factors contributing to this substantial loss of life remains crucial to ongoing efforts aimed at documenting the burden of COVID-19 globally 9,10,11,12 , reliable primary data on cause-specific mortality in India and other LMIC settings remain lacking. Lessons from the COVID-19 pandemic may inform assessments of policy decisions around the implementation and relaxation of non-pharmaceutical interventions as well as both global and within-country vaccine distribution 13 . An early country-wide lockdown was effective in slowing SARS-CoV-2 transmission and blunting the first wave of COVID-19 cases; 14 however, concerns arose that such measures—when implemented in LMIC settings including India—would lead to disproportionate harm through secondary effects on healthcare provision, food security, accidents, and interpersonal violence 15,16,17,18,19,20,21 . Excess death estimates encompass deaths directly caused by COVID-19 and deaths indirectly resulting from overburdening of healthcare systems with COVID-19 cases. In addition, these estimates include both increases and decreases in deaths due to other causes which may have been altered by implementation of non-pharmaceutical interventions. Whereas studies of cause-specific mortality are thus needed to better distinguish the role of COVID-19 and other factors in excess pandemic-associated mortality 22,23 , a lack of a lack of functioning vital registration systems has precluded nationwide assessments of mortality due to COVID-19 and other causes within India and most other LMIC settings 24 .
While country-wide mortality and cause-of-death surveillance are infeasible in India, subnational settings with well-functioning vital surveillance systems provide a valuable opportunity to probe changes in the frequency and causes of death during the COVID-19 pandemic. The southern state of Tamil Nadu ranks highly among Indian states in its per-capita public health investment and medical workforce, and is well regarded for the effectiveness of its primary healthcare delivery system 25 . It is also one of few Indian states with a well-performing Civil Registration System, estimated to capture 100% of deaths, based on concordance with mortality estimates from India’s parallel sample-based mortality surveillance system 26 . Madurai, a major city and administrative district within Tamil Nadu, was the site of an active SARS-CoV-2 surveillance program during the early COVID-19 pandemic which included expanded case-detection efforts through door-to-door polymerase chain reaction (PCR) testing, serological testing, and syndromic surveillance. Previous analyses of data collected through these efforts have estimated that only one in nine expected COVID-19 deaths during the first wave was captured through case-based surveillance 27 , consistent with experience in other regions of India 3,28,29 . To understand this differential between observed and expected mortality, as well as changes in mortality attributed to different causes, we analyzed Civil Registration System (CRS) records, including causes of death recorded by registered medical providers, among all decedents in Madurai over the period encompassing the initial country-wide lockdown and first two COVID-19 waves from March 2020 to July 2021.
Because implementation of lockdown measures could impede deaths registration 30 , we first aimed to validate the tabulation of new deaths occurring during this transition period in early 2020. Whereas non-pharmaceutical interventions may have impacted cause-specific mortality beginning from March 24, 2020, changes in recorded mortality between 1 and 23 March 2020—before lockdown measures were implemented or SARS-CoV-2 transmission became widely established—would be expected to signify changes in CRS data accuracy based on reporting timelines for deaths in Madurai (see Methods). Comparing observed mortality from the “control” period of 23 March 2020 to expectations based on 2018–19 observations, the incidence rate ratio (IRR) of deaths was 1.06 (95% confidence interval: 0.94–1.20) overall, 1.11 (0.95–1.31) for males, and 1.00 (0.83–1.21) for females (Table S1). No statistically significant change was apparent for deaths occurring in healthcare facilities (IRR = 1.05 [0.84–1.31]) or in the community (IRR = 1.07 [0.93–1.24]). While this analysis does not exclude the possibility that changes in accuracy may have occurred during the later stages of the pandemic, the time immediately surrounding India’s initial country-wide lockdown was expected to be associated with the most acute interruptions in ordinary CRS procedures. Previous studies have also reported that excess mortality estimates from CRS data closely resemble findings from other (e.g., survey-based) sources in Tamil Nadu 3,31 , supporting the use of these data to analyze pandemic-associated changes in mortality within Madurai.
Projecting mortality levels from 2018–19, we expected 15,377 (14,995–15,768) deaths to occur between 1 March 2020 and 31 July 2021 (Table 1). In total, 20,004 deaths were recorded, representing a 30% (27–33%) increase over baseline expectations. Increases occurred in both medically-supervised deaths within healthcare facilities (IRR = 1.14 [1.09–1.19]) and unsupervised deaths in the community (IRR = 1.38 [1.34–1.42]; Table S2; Fig. 1). Sex-stratified analyses yielded similar findings, with deaths increasing 31% (27–36%) among males and 29% (24–34%) among females (Table S3; Table S4). These patterns varied across ages, with the greatest increases in mortality apparent among individuals aged 70–79 years and ≥80 years (44% [37–51%] and 37% [29–45%] increases, respectively), consistent with susceptibility to severe outcomes of SARS-CoV-2 infection within these age groups (Fig. S1). In contrast, fewer deaths were recorded among younger males than expected (37% [11–55%], 38% [10–57%], and 18% [–2–34%] reductions in mortality at ages 0–9 years, 10–19 years, and 20–29 years, respectively). As these changes were not apparent during the control period (Table S1), and were not reflected among females of the same ages (Table S2), interruptions in CRS function are unlikely to explain the observed reduction in mortality among boys and young adult men.
Contrary to expectations that indirect harms associated with India’s lockdown could outweigh benefits associated with blunting SARS-CoV-2 transmission 15,16,17,18,19,20 , the period from 24 March to 31 May, 2020 saw 7% (0–3%) lower-than-expected all-cause mortality (Table 1). This observation was driven by a 33% (23–42%) reduction in healthcare facility deaths, whereas medically-unsupervised community deaths did not differ appreciably from expectations (IRR = 1.06 [0.98–1.15]; Table S2). While our analyses were underpowered for demonstrating statistically-significant changes in mortality during this period within age- and sex-specific strata, apparent reductions in mortality were greatest among boys and young adult men (80% [–3–94%], 52% [–60–82%], and 28% [–49–63%] reductions at ages 0–9 years, 10–19 years, and 20–29 years, respectively; Table S3). Overall, males experienced 14% (5–22%) fewer deaths than expected during the early lockdown, whereas among females, deaths did not differ from expectations (IRR = 1.03 [0.92–1.14]) or show clear age-associated patterns of change (Table S4).
Overall, mortality was 49% (43–56%) and 85% (78–93%) higher than expected during June–September 2020 and March–July 2021, the periods encapsulating the most acute phases of the first and second waves of COVID-19 cases (Table 1). Increases in healthcare facility deaths were modest during the first wave (11% [2–21%]) in comparison to the second wave (95% [82–110%]); this difference was less clearly pronounced for medically-unsupervised deaths (68% [60–78%] and 80% [71–90%] increases during the first and second waves, respectively; Table S2). Increases in mortality during these periods were apparent among both males and females (54% [45 –63%] and 43% [33–54%] increases during the first wave, respectively, and 87% [77–97%] and 82% (71–95%) increases, respectively, during the second wave; Table S3; Table S4). The greatest increases in mortality during the two waves occurred among adults aged ≥60 years (76–118% increases over the two waves). Among children aged 0–9 years, 45% (6–67%) fewer deaths than expected were recorded during the second wave; stratifying by sex, mortality was 58% (6–79%) lower than expected among boys aged 0–9 years, whereas increases among girls of the same ages were smaller and not statistically significant (22% [–71–62%]).
We used cause-of-death data from Medically-Certified Cause of Death (MCCD) reports issued for each death in Madurai to understand how changes in all-cause mortality tracked with changes in attributions of deaths to various immediate causes. Whereas medically-certified causes of death are available for only 43.9% of deaths within Tamil Nadu (in part reflecting rural-urban differences in reporting effort) 26 , all deaths are issued medically-certified causes within Madurai; this process continued through the pandemic until late within the second wave (Fig. 2). Physicians have responsibility for reporting causes of death for patients dying under their care in healthcare facilities, and for patients who die in the community after receiving medical treatment. For deaths in the community not preceded by medical care, the registered medical provider who declares the death has the responsibility for assigning a cause of death based on their own observations or information provided by survivors and authorities who came into contact with decedents. Cause-of-death assignments in MCCD data may be imperfect; attributions to vascular diseases, cancer, respiratory diseases, and non-respiratory infections were estimated to capture 72.9%, 49.9%, 65.7%, and 53.7%, respectively, of deaths attributed to the same conditions by verbal autopsy in historical studies dating to 1995–97 32 . In contrast, a greater proportion of deaths are attributed to senility or other uncategorized causes in MCCD data 33 . Consistent with these prior observations, MCCD data from Madurai during 2018–19 attributed a lower proportion of deaths to infectious diseases, and attributed a higher proportion of deaths to uncategorized causes including senility, in comparison to verbal autopsy-based estimates from all of India 34 (Table S5) and from other low-mortality districts 35 (Table S6). These comparisons should be interpreted with caution due to the lack of a gold-standard reference measure for causes of death in Madurai (e.g., due to the lower burden of malaria 36 and tuberculosis 37 in this setting), but provide important context when interpreting cause-of-death assignments from MCCD data.
While the above factors suggest MCCD data may not reveal the same distribution of causes of death as other methods, imperfect sensitivity and specificity of cause-of-death assignments would not be expected to bias the ability of MCCD data to capture changes in cause-specific mortality over time. We therefore focused our analyses on temporal changes in death attributions. We pair our presentation of changes in cause-specific mortality with earlier findings on the reliability of cause-of-death attributions in MCCD data, and with observations from the “control” period of 1–23 March, 2020 (Table S7).
Comparing observed data to expectations under a continuation of 2018–19 mortality levels, the period of 1 March 2020 to 31 July 2021 saw 50% (49–59%) fewer deaths attributed to infectious and maternal conditions, 18% (7–28%) fewer deaths attributed to injuries, and 21% (18–25%) more deaths attributed to noncommunicable diseases than expected (Table 2). Additionally, deaths attributed to senility and other uncategorized causes were 84% (61–112%) and 41% (33–50%) higher than expected, respectively. These patterns differed over time and across healthcare and community settings for each of the attributed causes (Fig. 2; Table S8–S12), as summarized below.
Wards with greater degrees of deprivation experienced greater excess mortality attributed to infectious diseases and maternal conditions, and greater excess mortality attributable to noncommunicable diseases, throughout the total pandemic period and during the second wave, in particular (Table S15). In contrast, community deprivation was not clearly associated with excess mortality attributed to other causes, although these analyses encountered limited statistical power in comparison to those addressing excess all-cause mortality. Measures aiming to capture the completeness of COVID-19 death reporting, including the ratio of deaths among individuals with documented SARS-CoV-2 infection to all-cause excess deaths, and the ratio of deaths attributed to COVID-19 to all-cause excess deaths, were not clearly associated with community deprivation indicators (Table S16).
Globally, excess mortality during the first two years of the COVID-19 pandemic exceeded deaths attributed to COVID-19 by nearly three-fold (14.9 vs. 5.4 million through 31 December 2021) 2 . The majority of unaccounted-for deaths are believed to have occurred in LMIC settings, where excess mortality over this period amounted to 12.7 million deaths, in contrast to 1.2 million deaths reported among confirmed COVID-19 cases. India alone accounts for 3.2–6.5 million excess deaths through 2021, representing 21–44% of the estimated total globally. Although India is a diverse country with significant variation in health systems, demographics, and completeness of civil registration records, findings from our analyses may provide insight into changes in mortality with implications extending beyond Madurai, where continuity of CRS and MCCD functions uniquely enabled our study. Discrepancies between all-cause and reported COVID-19 mortality in Madurai closely resemble findings in other settings within India, both over time and across demographic strata 3,27,55,56,57,58 .
While all-cause deaths increased 30% overall during the period from March 2020 to July 2021, excess deaths were concentrated during periods of peak SARS-CoV-2 transmission associated with the ancestral and Delta (B.1.617.2) variants. The first and second waves saw 25 and 32% of all deaths occurring during the study period, respectively, and 35 and 63% of all excess deaths estimated to have occurred, respectively; mortality occurring during the intervening periods was slightly lower than expected. Excess deaths attributed to infectious and maternal conditions—most prominently including lower respiratory tract infections—as well as cardiovascular/cerebrovascular conditions, diabetes, senility, and other uncategorized causes, exhibited strong temporal associations with surges in deaths attributed to COVID-19 or occurring among individuals with confirmed SARS-CoV-2 infection. Thus, these findings may reveal attributed causes of death among fatal COVID-19 cases, a majority of which likely went undiagnosed.
In contrast, deaths attributed to several other causes, including unintentional injuries and cirrhosis and related liver conditions, exhibited sustained reductions throughout the pandemic which were not observed during a “control” period from 1–23 March 2020, when lockdown-associated disruptions may have been expected to impact reporting. Non-pharmaceutical interventions may have played a role in reducing individuals’ risk of death due to these causes, consistent with findings in several high-income settings 59 . Implementation of lockdown measures accompanied expanded social welfare programs targeting older adults and other vulnerable populations in Madurai and other regions of India; these included increasing food rations for individuals covered by the National Food Security Act and door-to-door mobilization of social and community healthcare workforces to deliver essential supplies during the lockdown as well as screen individuals for symptoms and infection 27 . These efforts may have helped to mitigate vulnerable individuals’ risk of death due to causes unrelated to COVID-19.
While the 30% increase in deaths in Madurai due to all causes is consistent with observations in other settings across India 3,30 , this total greatly exceeds global average increases of 8% in 2020 and 18% in 2021 4 . Other settings with >25% increases in mortality during the pandemic have included Russia and LMICs within Latin America, Eastern Europe, and central Asia 1,4 . As total excess mortality encompasses both COVID-19 deaths as well as increases and decreases in deaths due to causes unrelated to COVID-19, it is important to note that excess mortality ratios do not directly measure the burden of deaths attributable to COVID-19 in any setting.
While excess medically-supervised deaths attributed to lower respiratory infections varied in association with deaths among confirmed COVID-19 cases and deaths attributed to COVID-19, excess deaths attributed to cardiovascular/cerebrovascular conditions, diabetes, and senility or other uncategorized causes occurred primarily in community settings. Such attributions should be viewed with particular scrutiny, as medical practitioners issuing cause-of-death assignments for individuals who died in the community without prior medical care may have had limited access to clinical information to inform their assessment; attributions to senility and uncategorized causes pose particular risks of inaccuracy 33,60 . Deaths for which these causes were assigned could also have increased in the community as a result of avoidance of healthcare settings, or due to patients’ lack of access to healthcare facilities during periods with substantial COVID-19 caseload.
Our study has several limitations. First, cause-of-death determinations in MCCD data are expected to be imperfect 33,61 . Our findings should thus be interpreted as representing changes in assigned causes of death over time rather than characterizing the true distribution of causes of death within Madurai. Many deaths were not assigned causes at the height of the second wave, when the number of deaths occurring overwhelmed local health and vital surveillance systems. While our study benefits from continuous mortality records in 2018–19 to define baseline expectations, unstable counts for rare death attributions and in small age- or sex-specific strata may limit our ability to reliably predict expected mortality levels in 2020–21. Analyses of ward-level socioeconomic characteristics associated with excess mortality are likewise limited by the lack of updated Census of India data since 2011. As slums or makeshift settlements may be closely interspersed with higher-income communities, analyses undertaken at the level of city wards, and drawing on older data on socioeconomic characteristics of wards, may mask the full extent of socioeconomic variation in excess and cause-specific mortality associated with the COVID-19 pandemic. While the population of Madurai Corporation was estimated at 1,017,865 as of the 2011 census, a lack of reliable up-to-date estimates of population size, overall and within age- and sex-specific strata, and uncertainty about migration which may have occurred during the pandemic in association with lockdown measures, precluded direct estimation of per-capita mortality associated with excess COVID-19 deaths. Last, although it is reassuring that age- and sex-specific observations of all-cause mortality in Madurai broadly reflect those reported in other settings within India 3,28,57 , Madurai is only one city, and observations in this setting may not be uniformly generalizable across regions or to rural contexts. For example, whereas maternal and perinatal conditions are not a major cause of death in Madurai, poorer and rural regions of the country experience markedly higher maternal and infant mortality 62 . Notwithstanding these limitations, Madurai presents an important setting for analyses, as CRS data have limited reliability throughout much of the rest of the country. Data on medically-certified causes of deaths are even more scarcely available within other Indian states and LMICs, where understanding of the impact of the COVID-19 pandemic on cause-specific mortality remains limited.
Correctly accounting for deaths and their causes during acute emergencies carries important societal ramifications 9,10,11,12 . While mitigation of COVID-19 morbidity and mortality during the first wave was largely dependent upon non-pharmaceutical interventions, administration of the ChAdOx1 vaccine in India began 16 January 2021, prior to the Delta variant surge. Effectiveness of ChAdOx1 against Delta variant-associated severe disease has been estimated at ≥80% in multiple postlicensure studies 63,64 , including in India 65 . Thus, a majority of COVID-19 deaths during India’s second wave could have been programmatically preventable under a scenario with improved vaccine coverage. Reduced burden on healthcare systems may have also helped to prevent mortality associated with causes unrelated to COVID-19. While our study identifies that lockdowns may have been associated with increases in deaths due to cancer and other conditions for which treatment was interrupted, it is important to note that all-cause mortality declined during this period in Madurai, consistent with observations elsewhere in India 57 . This finding, together with prior evidence of the effectiveness of India’s nationwide lockdown in delaying widespread SARS-CoV-2 transmission 14 , presents a reassuring contrast to expectations in early 2020 that lockdowns could cause greater harm than public health benefit, although measures of impact besides mortality remain important to consider.
India is one of many settings where excess all-cause mortality during the pandemic vastly exceeded reported COVID-19 deaths 1 . While such gaps have received particular attention in LMICs 66,67,68 , where surveillance and vital registration systems may encounter particular strain, the need to reconcile COVID-19 deaths with all-cause excess mortality has also arisen in high-income countries 69,70,71 . Our findings support ongoing efforts to quantify mortality associated with the COVID-19 pandemic globally, to identify causes of systematic undercounting, and to evaluate the impact of COVID-19 lockdown measures on mortality associated with various causes. Changes in cause-specific mortality during acute phases of the COVID-19 pandemic should inform planning for future public health emergencies necessitating non-pharmaceutical interventions.
Data for this study were generated through routine surveillance via CRS and MCCD functions within the municipality of Madurai (Madurai Corporation), including MCCD reporting. Vital surveillance in India is mandated under the Registration of Births and Deaths Act of 1969, which provides standardized elements for reporting births and deaths. Operating procedures of the CRS across India are decentralized, with the expectation that local registration units (states and districts) develop customized strategies adapted to their unique contexts. In Tamil Nadu, a local Coordinating Committee leads vital surveillance within districts. Heads of affected households, in coordination with executive officers of lower administrative units (e.g., taluks, blocks, wards, or villages) have legal responsibility for notification of all births and deaths occurring in the community to their local registration unit. Police, community healthcare workers, operators of crematories or cemeteries, and other officials who come into contact with the deceased have a responsibility to report deaths that they are the first to observe. Under the MCCD system, attending physicians have responsibility for reporting deaths occurring in healthcare facilities under their supervision, and for assigning medically-certified causes of death to these patients. In addition, physicians who provide care to individuals who die in the community have the responsibility for assigning causes of death; causes of death for individuals who die in the community without prior care are assigned by registered medical providers at the point of declaring each death 72 . Data included unique records for each death within Madurai Corporation during the periods of interest (2018–19 and 2020–21) abstracted from standardized MCCD forms (Table S17) and CRS forms (Table S18) including, for each decedent, their age, sex, date of death, ward of residence, attributed immediate cause of death, and name of the facility where the death occurred (for deaths occurring in healthcare facilities). Recording of fatal outcomes among individuals with SARS-CoV-2 infection has been described previously for this setting 27 . Death registrations in MCCD data for confirmed COVID-19 cases who experienced fatal outcomes within ≤30 days of a positive test outcome were linked via manual record review.
We aimed to compare observed mortality during the pandemic period (2020–21) to expectations under a status-quo scenario of pre-pandemic mortality levels observed in 2018–19. Periods of interest for analysis included the total pandemic period (1 March 2020 to 31 July 2021); early lockdown period (24 March–31 May 2020); first wave (1 June–30 September 2020); and second wave (16 March–15 July 2021). For each period, we summed total deaths across the applicable date range during the pandemic (2020–21) and pre-pandemic (2018–19) periods. We also defined a “control” period from 1–23 March 2020 to assess changes in reporting potentially associated with implementation of non-pharmaceutical interventions. Whereas deaths during 1–23 March 2020 are not likely to have been attributable to COVID-19, and preceded implementation of non-pharmaceutical interventions, reporting of deaths occurring during this period could have been affected by acute lockdown-associated disruptions (for deaths between 1 January and 29 February 2020, median time to registration was 12 days, and 30.8% were registered >21 days after occurring). We, therefore, aimed to compare observed to expected mortality during this period to determine whether implementation of lockdown measures was associated with changes in deaths reporting.
For primary analyses of excess all-cause mortality, we fit Poisson regression models to total mortality counts for each period of interest; data included period-specific total mortality during the applicable date ranges in 2018, 2019, 2020, or 2021. For each period, we defined population offsets accounting for expected changes in age- and sex-specific population sizes (log-transformed) for 2020–21 (Table S19). We used this framework to estimate the IRR of mortality comparing the pandemic period to 2018 and 2019, accounting for interannual variability during these years as well as expected changes in population size based on estimated year-on-year changes in population within age- and sex-specific population strata for Madurai District. To quantify absolute mortality expectations under a continuation of pre-pandemic mortality levels, we sampled from a Poisson distribution with a rate parameter defined as the product of the fitted model intercept (corresponding to averaged mortality rates during 2018 and 2019) and person-time at risk within the pandemic period. We subtracted these results for projected mortality from observed mortality to quantify absolute excess mortality. We used the same framework to compare observed mortality within age- and sex-specific population strata, and mortality attributed to specific causes, against expectations based on pre-pandemic observations.
To understand the potential association of various attributed causes of death with unconfirmed, fatal COVID-19 cases, we next aimed to determine whether cause-specific excess death time series were correlated with two measures of mortality related to COVID-19: (1) the number of deaths attributed to COVID-19 in MCCD data, and (2) the number of deaths occurring among individuals with confirmed SARS-CoV-2 infection. We generated 2-week moving average time series of expected mortality for all causes and for specific causes, in both healthcare and community settings, by sampling from Poisson distributions for total deaths each calendar day; these Poisson distributions were parameterized by mortality rates over the period from 7 days before to 7 days after each calendar date. We defined two-week moving averages of excess mortality by subtracting sampled draws from the distribution of expected mortality over each period 7 days before to 7 days after each period of interest from observed mortality over the same period. We then computed Pearson correlation coefficients measuring the association of the resampled time series of 2-week moving averages of excess mortality with each of the two independent variables measuring COVID-19-related mortality.
We also aimed to measure associations of excess mortality with socioeconomic characteristics of communities within Madurai. We used data from the 2011 Census of India to characterize socioeconomic attributes of wards within Madurai (N = 100); whereas the census is ordinarily carried out on a decennial basis, data collection for the 2021 Census of India was delayed to 2023 due to disruptions from the COVID-19 pandemic. Consistent with prior analyses 57 , extracted variables included the following: household crowding (measured as the mean number of individuals per room within a household); adult illiteracy (as a proportion among all adults); membership in scheduled castes or tribes (measured as a proportion among all adults); household condition (measured as the proportion of households classified as residing within dilapidated structures, within non-permanent structures, within structures with unfinished flooring, within structures without electric lighting, or reliant on solid cooking fuels); water and sanitation access (measured as the proportion of households lacking tapped, treated water sources, without onsite water sources, without onsite latrines, without onsite sewer connections, or reliant on open defecation); access to banking (measured as the proportion of households with bank accounts); and lack of material assets (measured as the proportion of households lacking each of the following index assets: computer, phone, bicycle, scooter/moped/motorcycle, and car/jeep van; Table S14). We generated a single index measuring community disadvantage as the first principal component of all extracted variables. The first principal component accounted for 48% of variation across all indicators and was positively associated with each measure besides household crowding at the two-sided p < 0.05 threshold.
We quantified associations of excess deaths with socioeconomic characteristics of communities via a regression model defining log-transformed mortality IRRs during each period (as estimated using the Poisson regression framework described above), within each ward, as the outcome, and the community deprivation index as the exposure. To propagate uncertainty in ward-level excess mortality, we repeated regression analyses across 10,000 samples from the distribution of excess mortality estimates by ward. As a secondary analysis, we also estimated associations of excess ward-level mortality with each of the deprivation indicators used to generate the principal component-based index individually, applying the same regression framework.
Finally, we aimed to assess whether community socioeconomic attributes were associated with the completeness of COVID-19 mortality reporting. We assessed reporting completeness using case-based and mortality-based surveillance, defining two independent variables as potential measures of completeness: (1) the proportion of all excess deaths occurring among reported cases within each ward, obtained by dividing ward-level deaths among confirmed COVID-19 cases by total excess deaths; (2) the proportion of all excess deaths attributed (on death certificates) to COVID-19, obtained by dividing ward-level deaths attributed to COVID-19 by total excess deaths. Each of these analyses propagated uncertainty in excess death measures according to the framework described above; parameter estimates are pooled across analyses undertaken on unique draws from the distribution of excess mortality by ward.
Local authors (C.M.B., G.K., R.L.) were involved in all components of the research process including study conception in relation to local priorities, study design, and authorship of publications. Analyses of de-identified mortality data generated through routine vital surveillance functions of the Civil Registration System were considered to constitute non-human subjects research and ethical review of the project was not required by local institutions (Madurai Municipal Corporation) as a condition for data sharing, or by foreign institutions with which the authors are affiliated (University of California, Berkeley; Princeton University) as a condition for data analysis (NIH Exemption Category 4 for non-human subjects research). This research project did not involve risks to human or animal subjects or environments, or international transfer of biological materials or cultural artefacts from India. Studies undertaken in India have been cited appropriately.
Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.
