Geospatial Variability in Excess Death Rates during the COVID-19 Pandemic in Mexico: Examining Socio Demographic and Population Health Characteristics ====================================================================================================================================================== * Sushma Dahal * Ruiyan Luo * Monica H. Swahn * Gerardo Chowell ## Abstract **Background** Mexico has suffered one of the highest COVID-19 mortality rates in the world. In this study we examined how socio demographic and population health characteristics shape the geospatial variability in excess mortality patterns during the COVID-19 pandemic in Mexico. **Methods** Weekly all-cause mortality time series for all 32 Mexican states, from January 4, 2015 to April 10, 2021, were analyzed to estimate the excess mortality rates using Serfling regression models. The association between socio-demographic, health indicators and excess mortality rates were determined using multiple linear regression analyses. Finally, we used functional data analysis to characterize clusters of states with distinct mortality growth rate curves. **Results** The overall all-cause excess deaths rate during the COVID-19 pandemic in Mexico until April 10, 2021 was estimated at 39.66 per 10 000 population. The lowest excess death rates were observed in southeastern states including Chiapas (12.72), Oaxaca (13.42) and Quintana Roo (19.41) whereas Mexico City had the highest excess death rate (106.17), followed by Tlaxcala (51.99) and Morelos (45.90). We found a positive association of excess mortality rates with aging index (P value<.0001), marginalization index (P value<.0001), and average household size (P value=0.0003) in the final adjusted model (Model R2=76%). We identified four distinct clusters with qualitatively similar excess mortality curves. **Conclusion** Central states exhibited the highest excess mortality rates whereas the distribution of aging index, marginalization index, and average household size explained the variability in excess mortality rates across Mexico. Our findings can help tailor interventions to mitigate the mortality impact of the pandemic. **Key message** * This study quantified and examined spatial patterns of excess mortality across states of Mexico, with lower rates of excess mortality in southeastern states and higher rates in central states. * Mexico City recorded 1 in 5 of all excess deaths in reported in Mexico, which accounted for 19% of total excess deaths across the country. * Findinds indicate that aging index, marginalization index, and average household size played a significant role on excess death rates across Mexican states during the COVID-19 pandemic. * Four distinct clusters characterized the excess mortality curves across Mexican states. Key words * excess mortality * COVID-19 pandemic * Mexico * states * socio-demographic factors * spatial variation ## Introduction Monitoring all-cause excess mortality, above an expected level of total deaths, as a pandemic unfolds is one of the key ways to evaluate its mortality impact (1). All-cause excess mortality estimates include deaths that are directly or indirectly attributed to the pandemic (2, 3). While the excess deaths can be directly attributed to COVID-19, other excess deaths attributed to the pandemic such as those that could be related to denied or delayed care for acute emergencies (4) (5) or other chronic conditions (6), the disruption of routine health care services in an overburdened health care system(7), unaddressed mental health concerns including suicide and self-harm (8, 9), and drug overdoses (10). Detailed analyses of excess mortality can help determine where the mortality impact of the pandemic has been most significant. Mexico is one of the countries in Latin America that is bearing the brunt of COVID-19 pandemic with the fourth-highest number of COVID-19 deaths in the world, after the USA, Brazil, and India, as of late June 2021 (11). In fact, Mexico has reported a total of 2 487 747 (1.38% of global cases) confirmed cases of COVID-19, including 231 847 deaths (5.96% of global deaths), as of June 25, 2021 (12). A previous study reported a high all-cause excess death rate of 26.10 per 10 000 population in Mexico in 2020, with COVID-19 deaths accounting for only 38.64% of the estimated excess deaths (13). While this relatively low proportion of COVID-19 deaths out of all excess deaths could be the result of low testing rates, misclassification of COVID-19 deaths, and delays in reporting COVID-19 deaths (14), a substantial number of deaths during the pandemic could be due to the indirect causes (3) and need to be examined in more depth. The distribution of indirect causes of deaths depends on several factors such as sociodemographic characteristics, population health and the selection, timing and intensity of any public health interventions, in addition to the efficiency and reach of the health and social care system (15). In Mexico, pandemic control measures have varied widely (16). Therefore, a more detailed understanding of the mortality burden of the pandemic can be obtained by quantifying spatial heterogeneity in excess deaths at a state level and by examining the influence of underlying sociodemographic, economic, and health system related factors. In this study, we pose the question whether in a country such as Mexico, with very high covid-19 mortality, potential spatial variability in the excess deaths can be explained by underlying sociodemographic and population health indicators. To answer that research question, we first estimated the all-cause excess deaths during the COVID-19 pandemic in Mexico comprising 31 states and Mexico City. Next, we evaluated the potential associations between different socio-demographic factors and excess mortality patterns at the state level in Mexico. Furthermore, we also conducted a cluster analysis to characterize the shapes of the excess mortality curves into different groups that describe the potential geospatial variability in excess mortality. Analyses such as these are critically important for understanding excess mortality and for guiding intervention strategies. ## Methods ### Data We obtained weekly all-cause death counts updated on May 25, 2021 for Mexico at the state level and for Mexico City, based on epidemiological weeks from January 2020 until April 10, 2021 and for the preceding 5□years (2015–2019) to establish a baseline mortality level(17). We accessed publicly available weekly mortality data from the National Institute of Statistics and Geography (INEGI) for the years from 2015 to 2018, and data from National Population Registry (RENAPO) for the years 2020 and 2021(17). For the year 2019, we chose either INEGI or RENAPO as the data source, based on the value of the last week of 2018 and the first week of 2019 for each state. We obtained the national and state-level population size estimates from the National Population Council (CONAPO) of Mexico(18). Mortality data was not available for the state of Tlaxcala for the last six weeks of the study period. For this reason, this state was excluded from our regression and functional cluster analyses. For each state, including Mexico City, we obtained data on six variables: population density (2020), aging index (2020), average household size (2020), marginalization index (2020), rate of new cases of depression per 100 000 population (2019), and public spending on health as percent of GDP (2019). Data on population density, aging index, average household size, and rate of new case of depression were obtained from INEGI (19), data on public spending on health was obtained from the subsystem of health accounts at the federal and state level (SICUENTAS)(20), and the data on the marginalization index was available from CONAPO(21). Summary statistics of these variables from 31 states and Mexico city are provided in Table 1. View this table: [Table 1.](http://medrxiv.org/content/early/2021/08/12/2021.08.11.21261930/T1) Table 1. Descriptive statistics for six study variables included in the multiple regression analysis of Excess Mortality in Mexico (*n*=32) ### Pandemic period and excess deaths For both the national data and the data for each state, we separately estimated the baseline mortality level by fitting cyclical Serfling regression models to all-cause deaths in the non-COVID-19 period, after excluding data from March 2020 to April 2021 by employing established methodology (2, 22-25). Details on the model equation that was used can be found in ref (13). After establishing a weekly baseline and the corresponding 95% CI at the national level, we defined the periods of COVID-19 pandemic as the weeks in 2020 and 2021 where the observed all-cause mortality rate at the national level in Mexico exceeded the upper 95% confidence limit of the national baseline mortality rate. The excess mortality rate was estimated at the state level and for Mexico City for the same defined period of COVID-19 pandemic. Excess all-cause mortality rate was estimated as the difference between the observed and model adjusted baseline mortality rates for each week constituting the pandemic period. The overall pandemic excess mortality in 2020 and 2021 was calculated by summing the excess death rates across the pandemic weeks in the given year (13, 22, 24). Negative excess mortality estimates were replaced by zeros in our analyses to account for underreporting due to reporting delay (3, 26). ### Multiple regression analysis After estimating the total excess mortality rate for each state, we explored the association between the total excess mortality rate and the predictor variables. Because the population density and rate of new case of depression distributions were skewed, we transformed these variables to log base 10. Since we identified Mexico City as a potential influential point, we performed sensitivity analysis by comparing the results of different models, including and excluding Mexico City. Since there was no significant change in the statistical inference of the parameters, we included Mexico City in the multiple regression analysis. ### Cluster analysis We followed the analytic methods described in (27) to pre-process the weekly cumulative all-cause excess deaths for 30 states and Mexico city (excluding Tlaxcala, refer to study setting for details). Then, we analyzed the shapes of the excess all-cause death rate curves to compare, cluster, and summarize growth rates. We employed the following steps to smooth and normalize the weekly all-cause excess death data: 1. Smoothing: Cumulative excess deaths curves were smoothed using smooth function in Matlab which uses a moving average filter over a 10-week span. 2. Time differencing: If *f**i*(*t*) denotes the given cumulative number of excess deaths for state i on week t, then per week growth rate at time t is given by *g**i*(*t*) = *f**i*(*t*) − *f**i*(*t* − 1). 3. Re-scaling: We rescaled each curve by dividing each *g**i*(*t*) by the total excess number of deaths for a given state i, which is equivalent to computing *h**i*(*t*) = *g**i*(*t*) / *r**i* where ![Graphic][1] and K is the number of weeks in the period. 4. Smoothing: We then smoothed the normalized curves over a 5-weeks span, using the smooth function in Matlab. To identify the clusters by comparing the curves, we used a simple metric. For any two rate curves, *h**i* and *h**j*, we compute the norm ∥hi −hj∥, where the double bars denote the *L**2* norm of the difference function, i.e., ![Graphic][2], which is approximated by ![Graphic][3], where K is the number of weeks in the period. To perform clustering of thirty-one curves into smaller groups, we applied the dendrogram function in Matlab using the “ward” linkage as explained in ref. (27). The number of clusters was decided empirically by inspecting the overall clustering results. After clustering the states into different groups, we derived the average curve for each cluster using a time wrapping algorithm (27, 28). ## Results Observed death rate was greater than the upper 95% confidence interval of the baseline in 52 weeks starting from week of April 12-18, 2020, out of total 58 weeks from March 1, 2020, to April 10, 2021 (Figure 1). The excess mortality rate first peaked during the week of July 12-18, 2020, with the excess death rate of 1.04 per 10 000 population, then declined slightly for a few weeks and then increased again from the week of December 20-26, 2020, and reached a peak with an all-cause excess death rate of 1.99 per 10 000 population on the week of January 17-23, 2021. The excess death rate remained below 0.5 since the week of 28 February 2021 until the end of the study period. ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/08/12/2021.08.11.21261930/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2021/08/12/2021.08.11.21261930/F1) Figure 1. Mortality rate per 10 000 population, Mexico, January 2015–March 2021. The black curve is the observed weekly death rate. The grey curve is the predicted baseline death rate. Square dotted curves indicate the upper and lower 95% confidence intervals of the baseline death rate. The long dashed line indicates the COVID-19 pandemic period. Table 1 shows the total all-cause excess death rate for the national level, Mexico City, and 31 states of Mexico. The map showing the estimates at the state level is given in Figure 2. While the total excess death rate in Mexico was at 39.66 per 10 000 population, equivalent to a total of ∼508 289 excess deaths, the excess mortality rate in Mexico City was the highest and estimated at 106.17 per 10 000 population (∼95 690 total number of excess deaths). Among 31 states, Tlaxcala (51.99), Morelos (45.90), Puebla (45.12), and Mexico (44.43) were among the states with the highest excess mortality rates. The states with the lowest death rates included Chiapas (12.72), Oaxaca (13.42), Quintana Roo (19.41), and Yucatan (21.11) (Table 2). Only one state, Chiapas, had no excess deaths in 2021. COVID-19 accounted for only 42.16% of total exces deaths at the national level ranging from 20.97% in Chiapas to 76.05% in Quintana Roo. View this table: [Table 2.](http://medrxiv.org/content/early/2021/08/12/2021.08.11.21261930/T2) Table 2. Estimates for all cause excess mortality rate by state per 10 000 population during COVID-19 pandemic in Mexico, March 1, 2020-April 10, 2021 ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/08/12/2021.08.11.21261930/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2021/08/12/2021.08.11.21261930/F2) Figure 2. Map showing the excess death rate per 10 000 population at state level in Mexico. Table 3 shows the results of fitting a taxonomy of multiple regression models of excess mortality rate at the state level in Mexico. Since R2 was slightly different for Models 4, 5, 6, and 7, we performed a multiple partial F-test to select the final model. We failed to find a significant contribution of adding population density, depression rate, and public expenditure on health on predicting excess mortality rate after accounting for the contribution of aging index, marginalization index, and average household size (F-value3,24=0.39, P-value=0.7631). Therefore we selected Model 4 as the final model. Our final model was able to explain 76% of the observed variance in the excess mortality rate (Coefficient of determination (R2)=0.76). View this table: [Table 3.](http://medrxiv.org/content/early/2021/08/12/2021.08.11.21261930/T3) Table 3. Results of fitting a taxonomy of multiple regression models of excess mortality rate at the state level in Mexico (*n*=31) As shown in table 4, we found a positive association of excess mortality rate with aging index, marginalization index, and average household size in the adjusted model at 0.05 level of significance. View this table: [Table 4.](http://medrxiv.org/content/early/2021/08/12/2021.08.11.21261930/T4) Table 4. Results for the Final Regression Model 4 of excess mortality rate at the state level in Mexico (*n*=31) The results of our clustering analyses is displayed in a dendrogram plot (Supplemental Figure 1). Specifically, we identified the following four prominent clusters based on the shapes of excess growth rate curves at state level: *Cluster 1:* Baja California, Coahuila, Guanajuato, Hidalgo, Jalisco, Mexico, Mexico City, Michoacan, Morelos, Nayarit, Nuevo Leon, San Luis Potosi *Cluster 2*: Aguascalientes, Chihuahua, Durango, Queretaro, Zacatecas *Cluster 3*: Baja California Sur, Colima, Guerrero, Oaxaca, Puebla, Quintana Roo, Sinaloa, Sonora, Tabasco, Tamaulipas, Veracruz, Yucatan *Cluster 4*: Campeche and Chiapas ![Supplemental Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/08/12/2021.08.11.21261930/F3.medium.gif) [Supplemental Figure 1.](http://medrxiv.org/content/early/2021/08/12/2021.08.11.21261930/F3) Supplemental Figure 1. A: Dendrogram plot, B: map of Mexico showing the states in four different clusters Figure 3 shows the average growth rate curves and one standard deviation band around it. The growth patterns in each cluster are very distinct. For cluster 1, we see two different peaks in growth rate, first small peak in July 2020 and the second big peak in January 2021. For cluster 2, there is a rapid increase in growth rate since July 2021 that peak on around December 2020. Unlike cluster 2, in cluster 3, the first big peak in July is followed by a small peak in January. Finally, in cluster 4, the growth rate rapidly increases from April to July followed by a rapid fall and a small rise in January 2021. In overall, first peak in most of the states occurred in around July, 2020 and the second peak occurred in around January, 2021. ![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/08/12/2021.08.11.21261930/F4.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2021/08/12/2021.08.11.21261930/F4) Figure 3. Average growth rate in each cluster, the dotted blue line are the one standard deviation band around the average growth rate ## Discussion In this study we investigated the excess mortality patterns during COVID-19 pandemic at the national and subnational level in Mexico from March 1, 2020, to April 10, 2021. We estimated an excess all-cause mortality rate of 39.66 per 10 000 population at the national level (a total of ∼508 289 excess deaths), indicating a devastating mortality impact of COVID-19 pandemic in Mexico. We found that the excess mortality rate has continuously declined after the second COVID-19 peak during the week of January 17-23, 2021. Mexico City alone accounted for about 19% of total excess deaths in Mexico, with an excess mortality rate of 106.17 per 10 000 population. Interestingly, we found that the states with the highest excess death rate (i.e., Mexico City, Tlaxcala, Morelos, Puebla, Mexico) were the central states in Mexico, while the lowest excess death rates were observed in the southern states (i.e., Chiapas, Oaxaca, Quintana Roo, Yucatan, Campeche). In Mexico, the majority of the indigenous population are located in the southern states (29). In 2015, Oaxaca had the highest proportion of native population (65.73%), followed by Yucatan (65.40%), Campeche (44.54%), and Quintana Roo (44.44%) (30). Compared to non-native groups, the indigenous populations across continents have suffered health disparities and a greater burden of diseases, including higher infant mortality, and lower life expectancy (31). During the pandemic, the indigenous populations have remained at higher risk of infection and death (32, 33). A previous study demonstrated higher excess deaths in U.S. states with higher concentration of native Americans (34). In contrast, in this study we found a lower excess mortality rate during the pandemic in Mexican states with a higher proportion of the native population. This is an intriguing finding that warrant further inquiry and examination as it may provide great insight to factors that may buffer against the impact of the pandemic and other adverse health events including natural disasters. In our analyses of data form Mexico, we found that COVID-19 specifc deaths accounted for only 42.16% of total excess death at the national level, lowest in Chiapas (20.97%) and highest in Quintana Roo (76.05%). At the state level in Mexico, the timing and the rigor of implementation of public policies to contain the virus has varied widely (16). For example, some of the states such as Veracruz, Yucatan, Nuevo Leon, and Tamaulipas, established policies to promote social distancing before the federal government enacted those policies (16). While states such as Chiapas, Tabasco, San Luis Potoshi, and Zacatecas underperformed in implementation of public policy measures (16), some other states, with a relatively low excess mortality rate, such as Baja California Sur and Nayarit implemented public information campaigns and international travel restrictions for longer periods, despite the potential adverse impact on tourism, which is a major economic activity (16). We found a positive association between the aging index and excess mortality, in the adjusted model confirming previous studies linking older age and COVID-19. The aging index is defined as the number of older adults (60 years of age and older) for every 100 children and youth (0 to 14 years of age)(35), and it increases as the population ages. Older age is a significant predictor of COVID-19 mortality as well as mortality from other causes(36) (37) (38). Previous studies that analysed excess mortality pattern during the first wave of COVID-19 in 21 industrialized countries has shown that those aged 65 years and above comprised 94% of all excess deaths, indicating a very high risk of death among older aged population specifically due to COVID-19 (15). Similarly, our finding of a positive association between the marginalization index and and excess mortality support previous research underscoring the close link between social disadvantage and COVID mortality and the overall increased burden of the pandemic in marginalized populations. The marginalization index that we used is an indicator of the inequities in quality of housing, access to basic public services like electricity and drinking water, schooling, proportion of poorly paid population and other sociodemographic and population health characteristics (21, 39). There may be several explanations for our findings. For example, public health measures such as social distancing and sheltering in place to combat the COVID-19 pandemic resulted in a disproportionate burden to vulnerable and marginalized populations (40-42). Marginalized groups are also more likely to be infected by the coronavirus due to their the context of their living arrangements and which may limit the ability to self-isolate and socially distance. Similarly, it is well demonstrated that marginalized populations tend to have a higher prevalence of chronic conditions such as obesity, hypertension and diabetes which are all strong risk factors associated with poor prognostic outcomes among those infected with COVID-19 (42). To complicate matters further, marginalized population are also often at greater risk of dying due to other indirect causes such as limited access to already-stressed health care system, poor mental health outcomes, food insecurity, abuse and violence among other social ills. (41, 42). Interestingly, Chiapas, a Southern Mexico with a high marginalization index(21), a higher concentration of indigenous population(30), and lower average performance in implementation of public policies to combat COVID-19 (16), had the lowest excess mortality rate among the states examined in these analyses. While there is no clear explanation for these finding, the varying climate across the states examined may be a contributing factor. For example, the southern states in Mexico have a weak seasonality and a tropical climate throughout the year (43). According to previous studies conducted in Mexico, tropical climate delayed the local transmission of SARS-CoV-2 at regional level. As such, the temperate climate regions like Tlaxcala and Jalisco may have been more vulnerable for local transmission than the tropical climate regions such as Chiapas and Veracruz. These findings should be replicated in other setitngs that comprise multiple climate regions to determine the impact of seasonality in the transmission spread and impact of excess mortality patterns. Additionally, more research is needed to elucidate the factors associated with lower all-cause excess death rate in the relatively marginalized southern states as observed in this study. Such insight may provide mitigation strategies for other regions with higher impact. We also found a positive association between average household size and excess death rate in the adjusted model. Although the links between average family size and excess death rate at the state level have not been reported elsewhere, the average family size could have interacted with other social determinants of health such as poverty, food insecurity, and lack of access to health care. Further studies are needed to understand the potential mechanism underlying this association and to more specifically consider family size as a potential population-level indicator of communities at risk for increased impact. Our classification of excess deaths growth rate curves at the state level reflects four distinct categories of Mexican states. In all of the clusters the first peak of excess deaths growth rate curve occurred in around July, 2020 which happened after the phased reopening of non-essential services in June, 2020 in Mexico. The reopening of the country coincided with both increase in driving and walking trends, and the highest levels of COVID-19 deaths that remained at a high level during June and July 2020 (44). The visual analysis of the growth rate curve indicates that the coastal tropical southeastern states were most affected during the first few months of the pandemic compared to other states. However, these states exhibited the lower overall excess death rate that could indicate the effect of temperature and other environmental factors. Moreover information on the growth rate curves can be utilized at the state level to guide the implementation of medical and public health measures. Besides, learning from the public health measures implemented in states of one cluster (for example, cluster 4) can be helpful to the other states (for example, states in cluster 1). To our knowledge, this is the first study that assesses the growth rate curves of excess deaths at the state level in Mexico. In our study, the estimates of excess death rate, as well as the proportion of COVID-19 attributed deaths, could be underestimated due to factors such as low testing rates in Mexico, misclassification of COVID-19 deaths, and delay in reporting COVID-19 deaths. This limitation should be considered when interpreting these findings. Additionally, other potential confounders that were not measures may explain the patterns of excess mortality across states. ## Conclusion Our estimate of all-cause excess death rate in Mexico was 39.66 per 10 000 with central states exhibiting higher rates and southern states exhibiting lower rates. Our study highlights that several population measures including the aging index, marginalization index, and average household size were significantly associated with the all-cause excess mortality rates across Mexican states during the COVID-19 pandemic. Our excess mortality estimates can help tailor state specific medical and public health interventions to prevent excess mortality in vulnerable areas but targeting specific regions and socio-economic indicators. We also recommend further studies that investigate the lower excess death rate in southern states, and studies that explore the role of environmental factors, particularly the social determinants of health, in spatial variation in excess death rate in Mexico and other regions heavily impacted by COVID-19. ## Data Availability All the data used in the study are publicly available. The sources of those data are provided under the methods section. ## Declarations ### Ethics approval Not applicable ### Author contribution GC and SD conceptualized and designed the Study. GC directed the study’s implementation. SD conducted the formal analysis and drafted the manuscript. RL, MS, GC reviewed and edited the manuscript. All the authors have read and approved the final version of this manuscript. ### Data availability All the data used in the study are publicly available. The source of those data are provided under the methds section. ### Funding SD was funded by 2CI Doctoral Fellowship at Georgia State University. ### Conflict of interest None declared * Received August 11, 2021. * Revision received August 11, 2021. * Accepted August 12, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NoDerivs 4.0 International), CC BY-ND 4.0, as described at [http://creativecommons.org/licenses/by-nd/4.0/](http://creativecommons.org/licenses/by-nd/4.0/) ## References 1. 1.Weinberger DM, Chen J, Cohen T, Crawford FW, Mostashari F, Olson D, et al. Estimation of excess deaths associated with the COVID-19 pandemic in the United States, March to May 2020. JAMA Internal Medicine. 2020;180(10):1336–44. 2. 2.Serfling RE. Methods for current statistical analysis of excess pneumonia-influenza deaths. Public health reports. 1963;78(6):494. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2307/4591848&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19316455&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F12%2F2021.08.11.21261930.atom) 3. 3.CDC. Excess Deaths Associated with COVID-19 [Available from: [https://www.cdc.gov/nchs/nvss/vsrr/covid19/excess\_deaths.htm](https://www.cdc.gov/nchs/nvss/vsrr/covid19/excess_deaths.htm) 4. 4.Schirmer CM, Ringer AJ, Arthur AS, Binning MJ, Fox WC, James RF, et al. Delayed presentation of acute ischemic strokes during the COVID-19 crisis. Journal of neurointerventional surgery. 2020;12(7):639–42. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTE6Im5ldXJpbnRzdXJnIjtzOjU6InJlc2lkIjtzOjg6IjEyLzcvNjM5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDgvMTIvMjAyMS4wOC4xMS4yMTI2MTkzMC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 5. 5.Maringe C, Spicer J, Morris M, Purushotham A, Nolte E, Sullivan R, et al. The impact of the COVID-19 pandemic on cancer deaths due to delays in diagnosis in England, UK: a national, population-based, modelling study. The lancet oncology. 2020;21(8):1023–34. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1470-2045(20)30388-0&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F12%2F2021.08.11.21261930.atom) 6. 6.Douglas M, Katikireddi SV, Taulbut M, McKee M, McCartney G. Mitigating the wider health effects of covid-19 pandemic response. Bmj. 2020;369. 7. 7.Roberton T, Carter ED, Chou VB, Stegmuller AR, Jackson BD, Tam Y, et al. Early estimates of the indirect effects of the COVID-19 pandemic on maternal and child mortality in low-income and middleincome countries: a modelling study. The Lancet Global Health. 2020;8(7):e901–e8. 8. 8.Kawohl W, Nordt C. COVID-19, unemployment, and suicide. The Lancet Psychiatry. 2020;7(5):389–90. 9. 9.Sahoo S, Bharadwaj S, Parveen S, Singh AP, Tandup C, Mehra A, et al. Self-harm and COVID-19 Pandemic: An emerging concern–A report of 2 cases from India. Asian journal of psychiatry. 2020. 10. 10.Currie JM, Schnell MK, Schwandt H, Zhang J. Trends in drug overdose mortality in Ohio during the first 7 months of the COVID-19 pandemic. JAMA network open. 2021;4(4):e217112–e. 11. 11.statistica. Number of novel coronavirus (COVID-19) deaths worldwide as of June 24, 2021, by country: statistica; 2021 [Available from: [https://www.statista.com/statistics/1093256/novel-coronavirus-2019ncov-deaths-worldwide-by-country/](https://www.statista.com/statistics/1093256/novel-coronavirus-2019ncov-deaths-worldwide-by-country/) 12. 12.WHO. Mexico situation 2021 [Available from: [https://covid19.who.int/region/amro/country/mx](https://covid19.who.int/region/amro/country/mx) 13. 13.Dahal S, Banda JM, Bento AI, Mizumoto K, Chowell G. Characterizing all-cause excess mortality patterns during COVID-19 pandemic in Mexico. BMC Infectious Diseases. 2021;21(1):1–10. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s12879-021-05950-x&link_type=DOI) 14. 14.Gutierrez E, Rubli A, Tavares T. Delays in death reports and their implications for tracking the evolution of COVID-19. Available at SSRN 3645304. 2020. 15. 15.Kontis V, Bennett JE, Rashid T, Parks RM, Pearson-Stuttard J, Guillot M, et al. Magnitude, demographics and dynamics of the effect of the first wave of the COVID-19 pandemic on all-cause mortality in 21 industrialized countries. Nature medicine. 2020;26(12):1919–28. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-020-1112-0&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F12%2F2021.08.11.21261930.atom) 16. 16.Knaul F, Arreola-Ornelas H, Porteny T, Touchton M, Sánchez-Talanquer M, Méndez Ó, et al. Not far enough: Public health policies to combat COVID-19 in Mexico’s states. Plos one. 2021;16(6):e0251722. 17. 17.Government of Mexico. Excess mortality in Mexico [Available from: [https://coronavirus.gob.mx/exceso-de-mortalidad-en-mexico/](https://coronavirus.gob.mx/exceso-de-mortalidad-en-mexico/) 18. 18.CONAPO. Projections of the population of the municipalities of Mexico 2015-2030[ 19. 19.INEGI. Data [Available from: [https://www.inegi.org.mx/datos/](https://www.inegi.org.mx/datos/). 20. 20.General Directorate of Health Information. Subsystem of Health Accounts at the Federal and State Level (SICUENTAS): Government of Mexico; 2021 [Available from: [http://www.dgis.salud.gob.mx/contenidos/sinais/s_sicuentas.html](http://www.dgis.salud.gob.mx/contenidos/sinais/s_sicuentas.html) 21. 21.CONAPO. Marginalization index by state 2020 2020 [Available from: [https://datos.gob.mx/busca/dataset/indice-de-marginacion-carencias-poblacionales-por-localidad-municipio-y-entidad/resource/b0d4c947-ac1d-438b-8f96-1661ad021aeb](https://datos.gob.mx/busca/dataset/indice-de-marginacion-carencias-poblacionales-por-localidad-municipio-y-entidad/resource/b0d4c947-ac1d-438b-8f96-1661ad021aeb). 22. 22.Chowell G, Simonsen L, Flores J, Miller MA, Viboud C. Death patterns during the 1918 influenza pandemic in Chile. Emerging infectious diseases. 2014;20(11):1803. 23. 23.Chowell G, Viboud C, Simonsen L, Miller MA, Acuna-Soto R, Díaz JMO, et al. The 1918–19 influenza pandemic in Boyaca, Colombia. Emerging infectious diseases. 2012;18(1):48. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22257780&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F12%2F2021.08.11.21261930.atom) 24. 24.Dahal S, Jenner M, Dinh L, Mizumoto K, Viboud C, Chowell G. Excess mortality patterns during 1918–1921 influenza pandemic in the state of Arizona, USA. Annals of epidemiology. 2018;28(5):273–80. 25. 25.Viboud C, Eisenstein J, Reid AH, Janczewski TA, Morens DM, Taubenberger JK. Age-and sex-specific mortality associated with the 1918–1919 influenza pandemic in Kentucky. The Journal of infectious diseases. 2013;207(5):721–9. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jis745&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23230061&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F12%2F2021.08.11.21261930.atom) 26. 26.Aron J, Muellbauer J. Transatlantic excess mortality comparisons in the pandemic. nine. 2020. 27. 27.Srivastava A, Chowell G. Understanding spatial heterogeneity of COVID-19 pandemic using shape analysis of growth rate curves. medRxiv. 2020. 28. 28.Srivastava A, Klassen EP. Functional and shape data analysis: Springer; 2016. 29. 29.Britannica. Mexico [6/29/2021]. Available from: [https://www.britannica.com/place/Mexico](https://www.britannica.com/place/Mexico). 30. 30.Wikipedia. Indigenous peoples of Mexico. 31. 31.Curtice K, Choo E. Indigenous populations: left behind in the COVID-19 response. The Lancet. 2020;395(10239):1753. 32. 32.Power T, Wilson D, Best O, Brockie T, Bourque Bearskin L, Millender E, et al. COVID-19 and Indigenous Peoples: An imperative for action. Journal of Clinical Nursing. 2020. 33. 33.CIDRAP. Reports detail high COVID-19 burden in Native Americans: University of Minnesota [Available from: [https://www.cidrap.umn.edu/news-perspective/2021/04/reports-detail-high-covid-19-burden-native-americans](https://www.cidrap.umn.edu/news-perspective/2021/04/reports-detail-high-covid-19-burden-native-americans) 34. 34.Dahal S, Mizumoto K, Bolin B, Viboud C, Chowell G. Natality decline and spatial variation in excess death rates during the 1918–1920 influenza pandemic in Arizona, United States. American journal of epidemiology. 2018;187(12):2577–84. 35. 35.Instituto Nacional de Estadística y Geografía (INEGI). índice de envejecimiento por entidad federativa, serie de años censales de 1990 a 2020 [Available from: [https://www.inegi.org.mx/app/tabulados/interactivos/?pxq=Poblacion\_Poblacion\_05\_da611615-0bef-4433-933d-d6746c924ad4](https://www.inegi.org.mx/app/tabulados/interactivos/?pxq=Poblacion_Poblacion_05_da611615-0bef-4433-933d-d6746c924ad4) 36. 36.Bello-Chavolla OY, Bahena-López JP, Antonio-Villa NE, Vargas-Vázquez A, González-Díaz A, Márquez-Salinas A, et al. Predicting mortality due to SARS-CoV-2: a mechanistic score relating obesity and diabetes to COVID-19 outcomes in Mexico. The Journal of Clinical Endocrinology & Metabolism. 2020;105(8):2752–61. 37. 37.Ho FK, Petermann-Rocha F, Gray SR, Jani BD, Katikireddi SV, Niedzwiedz CL, et al. Is older age associated with COVID-19 mortality in the absence of other risk factors? General population cohort study of 470,034 participants. PloS one. 2020;15(11):e0241824. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F12%2F2021.08.11.21261930.atom) 38. 38.CDC. Older Adults [6/29/2021]. Available from: [https://www.cdc.gov/coronavirus/2019-ncov/need-extra-precautions/older-adults.html](https://www.cdc.gov/coronavirus/2019-ncov/need-extra-precautions/older-adults.html). 39. 39.World Bank. World Development Report 2009: Reshaping Economic Geography. Washington DC. 40. 40.Kantamneni N. The impact of the COVID-19 pandemic on marginalized populations in the United States: A research agenda. Journal of Vocational Behavior. 2020;119. 41. 41.Benfer EA, Wiley LF. Health Justice Strategies To Combat COVID-19: Protecting Vulnerable Communities During A Pandemic. Health Affairs. March 19 2020. 42. 42.Anderson G, Frank JW, Naylor CD, Wodchis W, Feng P. Using socioeconomics to counter health disparities arising from the covid-19 pandemic. Bmj. 2020;369. 43. 43.Burton T. Mexico’s seven climate regions 2013 [Available from: [https://geo-mexico.com/?p=9512](https://geo-mexico.com/?p=9512) 44. 44.Tariq A, Banda JM, Skums P, Dahal S, Castillo-Garsow C, Espinoza B, et al. Transmission dynamics and forecasts of the COVID-19 pandemic in Mexico, March-December 2020. PloS one. 2021;16(7):e0254826. [1]: /embed/inline-graphic-1.gif [2]: /embed/inline-graphic-2.gif [3]: /embed/inline-graphic-3.gif