key: cord-0933221-3dv7dkxc authors: Sarria-Guzmán, Yohanna; Bernal, Jaime; De Biase, Michele; Muñoz-Arenas, Ligia C.; González-Jiménez, Francisco Erik; Mosso, Clemente; De León-Lorenzana, Arit; Fusaro, Carmine title: Using demographic data to understand the distribution of H1N1 and COVID-19 pandemics cases among federal entities and municipalities of Mexico date: 2021-03-24 journal: PeerJ DOI: 10.7717/peerj.11144 sha: 5eee70488df497ba7a768e3e53f5d35cb113619f doc_id: 933221 cord_uid: 3dv7dkxc BACKGROUND: The novel coronavirus disease (COVID-19) pandemic is the second global health emergency the world has faced in less than two decades, after the H1N1 Influenza pandemic in 2009–2010. Spread of pandemics is frequently associated with increased population size and population density. The geographical scales (national, regional or local scale) are key elements in determining the correlation between demographic factors and the spread of outbreaks. The aims of this study were: (a) to collect the Mexican data related to the two pandemics; (b) to create thematic maps using federal and municipal geographic scales; (c) to investigate the correlations between the pandemics indicators (numbers of contagious and deaths) and demographic patterns (population size and density). METHODS: The demographic patterns of all Mexican Federal Entities and all municipalities were taken from the database of “Instituto Nacional de Estadística y Geografía” (INEGI). The data of “Centro Nacional de Programas Preventivos y Control de Enfermedades” (CENAPRECE) and the geoportal of Mexico Government were also used in our analysis. The results are presented by means of tables, graphs and thematic maps. A Spearman correlation was used to assess the associations between the pandemics indicators and the demographic patterns. Correlations with a p value < 0.05 were considered significant. RESULTS: The confirmed cases (ccH1N1) and deaths (dH1N1) registered during the H1N1 Influenza pandemic were 72.4 thousand and 1.2 thousand respectively. Mexico City (CDMX) was the most affected area by the pandemic with 8,502 ccH1N1 and 152 dH1N1. The ccH1N1 and dH1N1 were positively correlated to demographic patterns; p-values higher than the level of marginal significance were found analyzing the % ccH1N1 and the % dH1N1 vs the population density. The COVID-19 pandemic data indicated 75.0 million confirmed cases (ccCOVID-19) and 1.6 million deaths (dCOVID-19) worldwide, as of date. The CDMX, where 264,330 infections were recorded, is the national epicenter of the pandemic. The federal scale did not allow to observe the correlation between demographic data and pandemic indicators; hence the next step was to choose a more detailed geographical scale (municipal basis). The ccCOVID-19 and dCOVID-19 (municipal basis) were highly correlated with demographic patterns; also the % ccCOVID-19 and % dCOVID-19 were moderately correlated with demographic patterns. CONCLUSION: The magnitude of COVID-19 pandemic is much greater than the H1N1 Influenza pandemic. The CDMX was the national epicenter in both pandemics. The federal scale did not allow to evaluate the correlation between exanimated demographic variables and the spread of infections, but the municipal basis allowed the identification of local variations and “red zones” such as the delegation of Iztapalapa and Gustavo A. Madero in CDMX. Disease outbreaks have frequently occurred in human history and have had significant impacts on many sectors such as public health, political and economic systems at a local and global scale. The novel coronavirus disease pandemic is the second global health emergency that governments, health organizations and citizens worldwide have faced over the last two decades (WHO, 2020). The widely known as "Swine flu", the first influenza pandemic of the 21 century, occurred in (CDC, 2020 caused by the virus A (H1N1) pdm09; the novel virus presented genetic material from human, pig and bird flu viruses. It was very different from H1N1 viruses formerly identified and, therefore, vaccination of seasonal flu offered little cross-protection (Fraser et al., 2009; Lim & Mahmood, 2011) . Most young people did not show immunity to this virus (Crum-Cianflone et al., 2009 ), but around one-third of the elderly manifested antibodies probably acquired from exposure to other H1N1 virus in their youth (CDC, 2020) . The symptoms of the Influenza A (H1N1) 2009 in people have ranged from mild to severe and include fever, cough, sore throat, nose, body aches, headache, chills, fatigue, diarrhea and vomiting (Witkop et al., 2010; Petrova & Russell, 2018; Saunders-Hastings et al., 2017) . There were approx. 2.8 million Influenza cases and between 151.7 thousand to 575.4 thousand deaths worldwide (April 2009 -April 2010 mostly people younger than 65 years (CDC, 2020) . Many countries, especially in Africa and Southeast Asia, lacked the ability to perform detection tests and still today, there is no certainty about factual pandemic deaths (WHO, 2020) . Mexico was among the first countries to notify the World Health Organization (WHO) of viral respiratory disease "Swine flu" and was, also, one of the most affected nations by this pandemic (Kain & Fowler, 2019; CDC, 2020) . The H1N1 Influenza confirmed cases (ccH1N1) were approx. 72.4 thousand while the deaths (dH1N1) were 1.2 thousand (data updated to 10th April, 2010) (http://www.cenaprece.salud.gob.mx/). On August 10th of 2010, the WHO declared the end of sanitary emergency related to Influenza A (H1N1) 2009 pandemic. However, the A (H1N1) pdm09 virus continues to circulate as a seasonal flu virus causing several hospitalizations and deaths in many developing countries every year (WHO, 2020) . At the end of 2019, cases of pneumonia associated with the novel coronavirus (SARS-CoV-2) were reported in the Chinese province of Wuhan. The SARS-CoV-2 is the seventh member of the coronavirus family to infect humans. It is characterized by prolonged incubation time and shedding from asymptomatic patients (Widders, Broom & Broom, 2020) . The COVID-19 outbreak spread quickly worldwide avoiding medical control and it was declared pandemic by the WHO on the 11th of March 2020 (WHO, 2020) . The SARS-CoV-2 infected patients are mostly asymptomatic or experience mild symptoms such as fever, dry cough, and sore throat . However, some patients, principally advanced age people or individuals with comorbities, develop severe and even fatal medical complications (Mehta et al., 2020; Sohrabi et al., 2020) . Risk factors and health conditions that could aggravate the clinical picture of patients with COVID-19 include hypertension, diabetes, cardiovascular disease, chronic respiratory disease, chronic kidney disease, immune compromised status, neurological disorders, cancer, smoking and obesity (Alqahtani et al., 2020; Nieman & Wentz, 2019; Roncon et al., 2020) . People with diabetes increased the odds to develop severe symptoms or die from COVID-19 by 2.7 (Williamson et al., 2020) . Patients that suffer from hypertension, cardiovascular or cerebrovascular disease were up to three times more likely to have severe COVID-19 symptoms than people in good health . Williamson et al. (2020) suggested that the risk of mortality from COVID-19 increase in patients with asthma or respiratory disease. The risk of severe symptoms or death from COVID-19 significantly increase in people with cancer, moreover, among patients with blood cancers the mortality rate is twice that of patients with solid tumors (Meng et al., 2020) . According to Simonnet et al. (2020) , people that suffer from obesity increase seven times the odds of developing severe COVID-19. Environmental issues such as air pollution could compromise lung function, and consequently increase the vulnerability to respiratory infections such as COVID-19. It has been proved that other local factors such us overcrowded areas, frequent social interactions and full public transport can increase the number of contagious by SARS-CoV-2 due its mechanism of transmission; the virus can persist in aerosols and surfaces. With this in mind, the proximity and contact among citizens result key determinants (Meyerowitz et al., 2021; Jordan, Adab & Cheng, 2020) . The incidence and morbidity of COVID-19 were different among countries and even regions of the same country, probably due to the different implementation and timing of mitigation strategies adopted by governments (social distancing, travel restrictions and quarantine) (Gerli et al., 2020) and social demographic indicators (population density, poverty and unemployment) (Ramírez & Lee, 2020) . The SARS-CoV-2 has infected more than 75.0 million people and we have approx. 1.6 million deaths worldwide by COVID-19 as of 18th December, 2020 (Coronavirus Research Center, 2020 . The WHO indicates Latin America as "red zone" of coronavirus transmissions in the world (3rd June, 2020). More specifically, in Mexico the ccCOVID-19 (1.2 million) and dCOVID-19 (115.1 thousand) are inexorably increasing day by day (data updated to 16th December, 2020) (COVID-19 México, 2020) . Spread of infectious diseases is frequently associated with increase in population size, population density, aged population density and construction land area (Vazquez-Prokopec et al., 2010; Barreto et al., 2011; Fang et al., 2012; Pequeno et al., 2020; Poole, 2020; Rodrigue, Luke & Osterholm, 2020; Verity et al., 2020) . Epidemics in smaller cities are generally focused on a shorter period (Dalziel et al., 2018) , while the incidence is diffuse in the bigger urban areas leading to more interaction among local residents, workers and tourist which make them potential hotspots for pandemics spread; therefore, communication routes such as national highways or freeways, and public transport that is, subway, buses and trams could magnify the impact of the outbreaks. Understanding the spatial distribution of the COVID-19 confirmed cases (ccCOVID-19), and their association with demographic factors can benefit to isolate the areas of highest risk of infection and the most vulnerable groups of the population at local scale. The geographical scales (national, regional or local scale) are key elements in determining the correlations between demographic factors and spread of outbreaks. The aims of this article were: (a) to collect the official Mexican data (counts of confirmed cases and deaths) related to the H1N1 Influenza pandemic (CENAPRECE, 2020) and the official governmental data related to the COVID-19 pandemic (COVID-19 México, 2020) both on a federal and municipal scale; (b) to produce thematic maps (population size-confirmed cases, population size-deaths); (c) to investigate the correlations between the pandemics indicators and demographic patterns. Mexico is constituted by 32 Federal Entities (FEs) and covers a geographical area of 1,964,375 km 2 . It is the second most populated country in Latin America after Brazil. Its population, that has increased in the last 10 years (+13.72%), passing from 112.34 to 127.79 million inhabitants, is mainly grouped in urban areas. All of FEs have reported an upward trend, although the increase in population has been slower than in the first decade of this century. The population density (2020) varied significantly among the FEs and ranged between 11 inhabitants per km 2 (hab./km 2 ) in Baja California Sur and 6,033 hab./km 2 in Mexico City (CDMX). The capital is far more crowded than any other city in Mexico and is one of the most populated cities in North America, with 9.02 million inhabitants distributed in the "relatively small" surface area of 1,495 km 2 . Its population density is higher than the national average (65 hab./km 2 ); in particular, the delegations of Iztapalapa (1.82 million hab.) and Gustavo A. Madero (1.18 million hab.) are among the most overpopulated areas in the world. The state of Mexico (22,351 km 2 ), in the central part of the country, neighboring with CDMX, is the most populated among the FEs (17.43 million hab.-780 hab./km 2 ), its inhabitants are concentrated mainly in the cities of Ecatepec de Morelos (1.71 million hab.) and Nezahualcóyotl (1.13 million hab.). Tlaxcala, although is the second smallest Mexican State in size (4,016 km 2 ), is in the fourth position in the ranking of the most densely populated places of the republic (344 hab./km 2 ). Its residents live mostly in Tlaxcala city (103,435 hab.) . Guerrero (63,596 km 2 -3.66 million hab.) in the Pacific coast and Yucatan (39,524 km 2 -2.26 million hab.) in the Gulf of Mexico reported a population density concordant with the national average. Acapulco de Juárez (840,795 hab.) and Merida (963,861 hab.) are the urban areas more populated in addition to being among the significant national and international tourist destinations. Baja California Sur (73,909 km 2 -804,708 hab.) in the Pacific coast and Durango (123,317 km 2 -1.87 million hab.) in the central part of the country are among the least populated FEs, with population density about five times lower than the national average. The demographic patterns that is, population size (hab.), population density (hab./km 2 ) and surface (km 2 ) characteristics of each one of thirty-two FEs and all Mexican municipalities were taken from the database of "Instituto Nacional de Estadística y Geografía" of Mexico (INEGI, 2020) . The topographic map (scale 1:250,000) updated to 2018 was adopted as cartographic base (CONABIO, 2020). As stated before, risk factors that could aggravate the symptoms derived by COVID-19 include hypertension, diabetes, chronic respiratory or cardiovascular disease, immune compromised status, neurological disorders, cancer, smoking and obesity (Alqahtani et al., 2020; Nieman & Wentz, 2019; Roncon et al., 2020) . In addition, various factors such overcrowded areas, frequent social interactions, full public transport, poor attention to preventive measures in particular disrespect the social distancing could be involved in the growth of the number of contagious during pandemics, which cause the collapse of the national health system and the drastic reduction in the number of ICUs available in Mexican hospitals. Using the official Mexican data, we have chosen to analyze only the correlation between demographic patterns that is, population size and density with the number of contagious and victims in the two pandemics (Pequeno et al., 2020) . The official data of CENAPRECE (http://www.cenaprece.salud.gob.mx/) related to ccH1N1 and dH1N1 were used in this analysis while the geoportal of Mexico's government (COVID-19 México, 2020 ) was adopted to find the data of ccCOVID-19 and deaths dCOVID-19 related to the pandemic. The results are presented as tables, graphs and thematic maps. The thematic maps were obtained using ArcMap 10.5 (ESRI, Redlands, CA, USA); two attributes have been drawn (population size-counts of confirmed cases or counts of deaths) in each map. A special focus was directed to CDMX. A Spearman correlation was performed to assess the associations between the pandemics data and demographic patterns. Correlations with p values < 0.05 were considered significant. The Spearman coefficient values range between −1 and 1and indicate negative or positive correlation, the values close to zero indicate poor or no correlation between the variables. The nonparametric method locally weighted linear regression (LOESS), developed thought the programing language Python, have been used according to Austin & Steyerberg (2014) for the estimation of relationships between the demographic and pandemics indicators; the results were presented through XY graphs. The H1N1 influenza pandemic The first ccH1N1 were recorded in the small rural town of La Gloria, Veracruz (central part of Mexico) at the end of March, 2009. The national ccH1N1 and dH1N1 registered by CENAPRECE (http://www.cenaprece.salud.gob.mx/) during the sanitary emergency were 72.4 thousand and 1.2 thousand respectively. The thematic maps (Figs. 1A and 1B) describe the spatial distribution of ccH1N1 and dH1N1 in Mexico and their association with population density. The CDMX was the most affected area by the "Swine flu" pandemic with 8,502 ccH1N1 and 152 dH1N1. Other four FEs passed the threshold of 4,000 ccH1N1 that is, State of Mexico (4,682 ccH1N1), San Luis Potosí (4,446 ccH1N1), Nuevo León (4,358 ccH1N1) and Jalisco (4,333 ccN1N1). The ccH1N1 were less than a thousand people in six FEs that is, Baja California Sur, Morelos, Quintana Roo, Sinaloa, Coahuila and Campeche. The ccH1N1 ranged between 1,000 and 4,000 in the rest of the FEs (Table 1) . The Spearman correlation at federal scale ( Fig. 2A) indicated different results. The ccH1N1 and dH1N1 (in absolute terms) were positively correlated (high correlation) with the demographic patterns. The % ccH1N1 showed a negative correlation with the population size (Spearman rho: −0.43, p: 0.011). A p-value higher than the level of marginal significance (i.e., p > 0.05) was found analyzing the correlation between the % ccH1N1 and the population density, also the correlation between the % dH1N1 and population density was not significant. It is probable than the federal scale was not suitable for describing the correlation between the variables under investigation. The data at municipal scale is not available. Mexican health authorities confirmed the first SARS-CoV-2 positive case in the country at the end of February 2020, a young man who was in CDMX, hospitalized at the National Institute of Respiratory Diseases (INER) that had previously traveled to Italy. Starting from March 2020 the ccCOVID-19 and consequently the victims have increased considerably, placing Mexico among the first places in the unpleasant ranking of countries most affected by COVID-19 pandemic. The thematic maps (Figs. 3A and 3B) describe the spatial distribution of ccCOVID-19 and dCOVID-19 in Mexico and their association with population density, while the Figs. 4A and 4B show a special focus on CDMX. The official data of national government (COVID-19 México, 2020) indicated 1.2 million ccCOVID-19 and 115.1 thousand dCOVID-19 in Mexico, as of date (18th December, 2020). The CDMX, where 264,330 infections were registered, is the national epicenter of the pandemic ( Table 2 ). The more crowed districts of the capital such as Iztapalapa (37,332 ccCOVID-19-2,812 dCOVID-19), Gustavo A. Madero (31,778 ccCOVID-19-2,570 dCOVID-19) and Tlalpan (24,598 ccCOVID-19-765 dCOVID-19) are the principal "red zones" of the country where the pressure on the health system is greatest. Other FEs that have exceeded the threshold of 40,000 ccCOVID-19 are mainly located in the north-central area of the country or in the Gulf of Mexico, that is, State of Mexico (15, and Nezahualcóyotl (12,929 ccCOVID-19-1,695 dCOVID-19) connected through a dense public transport network with CDMX are the most exposed areas to the COVID-19 health emergency. In Tlaxcala (9,984 ccCOVID-19-1,314 dCOVID-19), another State close to the capital, the majority of the ccCOVID-19 and dCOVID-19 were recorded in urban areas of Tlaxcala city (1, . Approx. two thirds of the ccCOVID-19 in Guerrero (24,723 ccCOVID-19-2,581 dCOVID-19) were found in the Acapulco de Juárez (10,781 ccCOVID-19-1,218 dCOVID-19) and in the neighboring city of Chilpancingo de los Bravo (4,572 ccCOVID-19-276 dCOVID-19), the municipalities of the interior have, at least so far, clearly lower numbers of ccCOVID-19. Yucatan (25, showed similar numbers to Guerrero, approx. half of the infections were registered in Merida (15, . Baja California Sur (15, , , Durango (23, and Colima (7, were among the FEs of Mexico with less incidence of the COVID-19 pandemic. The disease fatality rate (quotient between dCOVID-19 and ccCOVID-19) at federal scale (Fig. 5A ) ranged between 4.49% of Baja California Sur to 16.01% of Sinaloa. The Spearman correlation, at federal scale (Fig. 2B) , showed different results. On one hand the ccCOVID-19 and dCOVID-19 (in absolute terms) were positively correlated (high correlation) with the population size. On the other hand p-values higher than the level of marginal significance (i.e., p > 0.05) were found analyzing the correlations between the % ccCOVID-19 and the demographic patterns; also the correlations between the % dH1N1 and the demographic patterns were not significant. These results lead us to suppose that the federal scale do not allow to observe the correlation between the pandemics indicators and the two demographic patterns, therefore the next step was to choose a more detailed geographical scale, that is, on a municipal basis. The Spearman analysis, at municipal scale ( Fig. 2C) (2,455 municipalities), showed high correlation between the data. The ccCOVID-19 and dCOVID-19 (in absolute terms) were high positively correlated with the demographic patterns; also the % ccCOVID-19 and % dCOVID-19 showed moderate high correlations with demographic patterns. The municipal scale resulted more suitable than the federal scale to describe the spread of COVID-19 pandemic. The R 2 of LOESS ( Fig. 6A ) was equal to 0.76 indicating medium high correlation between ccCOVID-19 and municipal population, while the R 2 of 0.91 indicated high correlation between dCOVID-19 and the municipal population (Fig. 6B) . However, it is important to emphasize that correlation is not causation and it is not right to say that high population sizes caused more dCOVID-19 despite the analysis regression showed that the variables were related to each other's. The urban areas, that is, federal capitals and principal municipalities, were particularly affected by the COVID-19 pandemic; some example: in the city of La Paz, where approximately the 37.52% of the total Baja California Sur' population lives, the 50.22% of federal ccCOVID-19 were registered; Ciudad Juárez, the biggest city of Chihuahua with the 44.06% of federal population, showed the 53.81% of federal ccCOVID-19; also in the port city of Veracruz, one of the most important trading centers in the Gulf of Mexico where the 7.34% of federal population lives, the federal ccCOVID-19 were the 21.56%. The discrepancy between the % of the federal population and the % of federal ccCOVID-19 was visible in most cases throughout all of Mexico (Table 3) The disease fatality rate, at municipal basis ranged until approx. 20.00% (Baja California); values of approximately 14.00% were found in many cities of Mexico including Cuernavaca (Morelos), Culiacan (Sinaloa), Ecatepec de Morelos and Nezahualcoyotl (State of Mexico), Veracruz (Veracruz de Ignacio de la Llave). In CDMX the disease fatality rate (Fig. 5B ) ranged between 2.50% of Milpa Alta and approx. 7.50% of the Iztapalapa and Iztacalco. The COVID-19 pandemic is increasing worldwide, with the American continent as one of its epicenters and Mexico is among the ten worst countries in terms of infections and deaths. The magnitude of COVID-19 pandemic is much higher than the "Swine flu" pandemic; nevertheless it is possible to identify common aspects and differences between the two health emergencies. As in 2009 the most populated FEs of Mexico are the "red spots of infection" and report the greatest number of confirmed cases and deaths. The authorities have been trying to reduce the spread of the COVID-19 by mitigation activities and strategies such as social distancing, travel restrictions, closing schools, shutting down nonessential activities and quarantine (Méndez-Arriaga, 2020), as they had previously done during the 2009 health emergency (Chowell et al., 2008) . The major trouble has been identifying the root and the new cluster of infection in the urban areas. The SARS-CoV-2 is transmitted by human-to-human, symptomatic people are the most frequent source of contagious; the virus can also be spread by contaminated objects and aerosol inhalation (Helmy et al., 2020) . The metropolitan areas, involving districts tightly linked through socio-economic relationships and permanent commuting of citizens, workers and tourists are more exposed to the pandemic outbreaks compared to rural areas (Hamidi, Sabouri & Ewing, 2020) . Several causes could explain the correlation, observed at municipal scale, between demographic factors and the COVID-19 spread. The urban areas or those with high population size and density imply that more people live in per unit of area increasing the direct and indirect contact among the citizens, workers and tourist (Macharia, Joseph & Okiro, 2020) . Chhikara et al. (2020) indicate that high population density increased the possibility of COVID-19 diffusion, the urban-rural segmentation reduced the spread of outbreaks between urban and rural areas. Méndez-Arriaga (2020), investigating environmental climate patterns of COVID-19 spread in Mexico, showed that the highest national local transmission in CDMX and State of Mexico were not only the result of the environmental factors (temperature and regional climate) but also the demographic characteristics (population size and density). According to Chowell et al. (2011) , the three H1N1 pandemic waves across Mexican FEs were driven by the most populated cities, the earlier onset occurred in the overcrowded FEs located in the central part of the country, as it is nowadays during the COVID-19 emergency; however, in their publication the authors indicated that incidence rates observed in large population centers was lowest if compared with rural zone. The rural zones, as well as the degraded peripheral areas, in most cases, have structural gaps that is, small and poorly efficient hospitals, precarious sanitary conditions that could accelerate the spread of infectious diseases. In regard to COVID-19 pandemic, several countries are facing a second wave of infections, however, Mexico has not been the case yet and this study does not cover that scenario. The CDMX is the national epicenter point of the actual outbreak; all its citizens, that live closely in a "reduced space", appear to have accelerated the virus spread through multiple contacts and interactions in packed subway trains, busy playgrounds and crowded residential buildings creating "red spots" of infections. In particular, the districts of Iztapalapa and Gustavo A. Madero had the highest incidence of COVID-19 confirmed cases at national scale. According to Zepeda-Lopez et al. (2010) , also during the outbreak of the novel influenza A (H1N1)v virus in 2009, CDMX showed the higher number of ccH1N1 principally in the crowed districts of Iztapalapa and Gustavo A. Madero. According to Ponnambalam et al. (2012) and , population density, size if population and urban land are key factors in the infectious diseases transmission and, consequently, the urban areas present a significant challenge for national public health. The underground (metro), with approximately 4.5 million commuters every day (1.6 billion ridership per year), is the most crowded urban transportation system of CDMX. It brings people in close contact and facilitates the virus spread by human-to-human transmission. During the daily massive movements, the passengers, sharing the same elements (surfaces, objects) and air, promoted a significant interchange of human or environmental microbiota causing the wider spread of pathogens (Vargas-Robles et al., 2020) . The terminal stations such as Indios Verdes in the delegation of Gustavo A. Madero (north of CDMX) and Pantitlan in Venustiano Carranza (east of CDMX) that connect the inner city to suburban areas of Ecatepec de Morelos and Nezahualcóyotl (State of Mexico) were considered among the main "red spots" of COVID-19 contagious in CDMX by national government. In contrast to the policies implemented by most of countries worldwide, Mexico in March 2020 preferred not to apply restrictions for air travel keeping the Benito Juárez International Airport of CDMX opened, that is the nation's main gateway with approximately 50 million commuters every year. Foreign tourists or Mexican citizens returning may have been vectors of the virus. Massive religious events, such as the representation of the crucifixion of Jesus in Iztapalapa (Iztapalapa Passion Play 2009), which was held few days before the identification of the virus A (H1N1) pdm09 and was attended by more than 2 million people, were indeed dangerous events for the spread of the virus at that time. According to Zepeda-Lopez et al. (2010) , the congregation of people in the Iztapalapa Passion Play has contributed to the roll out of the virus A (H1N1) pdm09 during the initial outbreak phase throughout CDMX and beyond. Even more, while music events such as Estéreo Picnic, Lollapalooza Chile, Argentina y Brazil were reprogramed for next year, the "Vive Latino", one of the most important Rock festivals in Latin America took place in the midst of COVID-19 pandemic, on 14th and 15th March at the Foro Sol Stadium in CDMX, gathering more of 40,000 people. The government and festival organization, instead of suspending the musical show and postpone it on another occasion, have chosen to take health preventive measures installing eight medical tents, filters with infrared thermometers at the entrances, and a team of 92 paramedics, 10 doctors and 8 ambulances (Salud, 2020) . The government of CDMX identified seven market places including "Central de Abasto" in Iztapalapa and "La Merced" in Cuauhtémoc as significant sources of infection. They recommended the citizens to avoid buying groceries in those places. In particular, the "Central de Abasto", that supplies the capital and the center of the country and receives perishable food from all of Mexico, is an agglomeration point in the city representing high risk infection. The entry at this place was prohibited to pregnant women, children and older adults and the use of facemask mandatory. Misleading statements by government representatives and the lack of official restrictions that is, other than the closure of schools and suspension of classes in all academic levels, had left many of the preventative measures to citizen's "good judgment", since self isolation was voluntary, in the poorest areas of the CDMX such as Iztapalapa where several thousands of informal workers that have no other means of subsistence, and live on a "day to day" basis, need to work daily to meet basic family needs putting their health and their families at risk. Hence the decision to close the capital's "Zócalo" starting from the last week of March with almost 270 metal fences that protected the place in order to avoid crowds, potential meetings or social outings. Public and private universities in CDMX have taken mitigation measures against the strain of COVID-19, suspending face-to-face activities and starting online programs. Chowell et al. (2011) analyzing the epidemiological patterns of "Swine flu" pandemic highlighted the importance of school closure and other restriction measures to mitigate the impact of "Swine flu" pandemic and future sanitary emergencies in urban and rural areas. High densely populated cities within the same country have shown a grater out brake of COVID-19 than less densely populated cities (Tira, 2020) . However, it is important to mention the limitations of the observations presented here since the analysis were made using phew parameters that is, population, km 2 , population density, ccCOVID-19 and dCOVID-19; several other studies have shown the need to focus also on weather differences among regions (Rashed et al., 2020) , seasonality (Smit et al., 2020; Carlson et al., 2020) and social parameters such as income (Bonacini, Gallo & Scicchitano, 2020) . The geographical scale results a key factor to accurately describe the spread of viral diseases in Mexico; an adequate geographical scale of reference is crucial for designing more efficient control or mitigation measures in such a heterogeneous and socially complex population during pandemic emergencies. The municipality scale that used of 2,455 data against the 32 data of federal scale decreasing the uncertainty of the analysis and confirmed the correlation between demographic parameters and pandemic indicators; in our research a detailed geographic scale has produced better correlates with available data and highlighted the areas most exposed to the health emergency such us the delegation of Iztapalapa and Gustavo A. Madero in CDMX. 2020. Prevalence, severity and mortality associated with COPD and smoking in patients with COVID-19: a rapid systematic review and meta-analysis Graphical assessment of internal and external calibration of logistic regression models by using loess smoothers Successes and failures in the control of infectious diseases in Brazil: social and environmental context, policies, interventions, and research needs Working from home and income inequality: risks of a 'new normal' with COVID-19 Misconceptions about weather and seasonality must not misguide COVID-19 response Centers for disease control and prevention Corona virus SARS-CoV-2 disease COVID-19: infection, prevention and clinical advances of the prospective chemical drug therapeutics The 1918-1919 influenza pandemic in England and Wales: spatial patterns in transmissibility and mortality impact Characterizing the epidemiology of the 2009 influenza A/H1N1 pandemic in Mexico Comisión Nacional para el Conocimiento y Uso de la Biodiversidad CONACYT-CentroGeo-GeoInt-DataLab Clinical and epidemiologic characteristics of an outbreak of novel H1N1 (swine origin) influenza A virus among United States military beneficiaries Urbanization and humidity shape the intensity of influenza epidemics in U.S. cities Distribution and risk factors of 2009 pandemic influenza A (H1N1) in mainland China The WHO Rapid Pandemic Assessment Collaboration COVID-19 mortality rates in the European Union, Switzerland, and the UK: effect of timeliness, lockdown rigidity, and population density Gobierno de México-Secretaría de Does density aggravate the COVID-19 pandemic The COVID-19 pandemic: a comprehensive review of taxonomy, genetics, epidemiology, diagnosis, treatment, and control Instituto Nacional de Estadística, Geografía e Informática Covid-19: risk factors for severe disease and death Preparing intensive care for the next pandemic influenza Influenza A H1N1 2009 (Swine Flu) and pregnancy Characteristics of asymptomatic patients with SARS-CoV-2 infection in Jinan A vulnerability index for COVID-19: spatial analysis at the subnational level in Kenya COVID-19: consider cytokine storm syndromes and immunosuppression The temperature and regional climate effects on communitarian COVID-19 contagion in Mexico throughout phase 1 Cancer history is an independent risk factor for mortality in hospitalized COVID-19 patients: a propensity score-matched analysis Transmission of SARS-CoV-2: a review of viral, host, and environmental factors The compelling link between physical activity and the body's defense system Air transportation, population density and temperature predict the spread of COVID-19 in Brazil The evolution of seasonal influenza viruses Understanding the socioeconomic heterogeneity in healthcare in US counties: the effect of population density, education and poverty on H1N1 pandemic mortality Seasonal influences on the spread of SARS-CoV-2 (COVID19), causality, and forecastabililty COVID-19 emergence and social and health determinants in Colorado: a rapid spatial analysis Influence of absolute humidity, temperature and population density on COVID-19 spread and decay durations: multi-prefecture study in Japan Transportation and pandemics: the geography of transport systems. Fifth Edition Diabetic patients with COVID-19 infection are at higher risk of ICU admission and poor short-term outcome Effectiveness of personal protective measures in reducing pandemic influenza transmission: a systematic review and meta-analysis LICORN and the Lille COVID-19 and obesity study group, high prevalence of obesity in severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) requiring invasive mechanical ventilation Winter is coming: a southern hemisphere perspective of the environmental drivers of SARS-CoV-2 and the potential seasonality of COVID-19 World Health Organization declares global emergency: a review of the 2019 novel coronavirus (COVID-19) About the sustainability of urban settlements Passenger-surface microbiome interactions in the subway of Mexico City Quantifying the spatial dimension of dengue virus epidemic spread within a tropical urban environment Estimates of the severity of coronavirus disease 2019: a model-based analysis Does comorbidity increase the risk of patients with COVID-19: evidence from meta-analysis SARS-CoV-2: the viral shedding vs infectivity dilemma Factors associated with COVID-19-related death using OpenSAFELY Novel Influenza A (H1N1) outbreak at the U.S. Air Force Academy: epidemiology and viral shedding duration World Health Organization Distribution of COVID-19 morbidity rate in association with social and economic factors in Wuhan, China: implications for urban development Inside the outbreak of the 2009 influenza A (H1N1)v virus in Mexico The authors thank one anonymous reviewer for helpful comments on a previous version of the manuscript. The authors received no funding for this work. The authors declare that they have no competing interests. Yohanna Sarria-Guzmán conceived and designed the experiments, performed the experiments, prepared figures and/or tables, authored or reviewed drafts of the paper, searched the official data in the specific sources, and approved the final draft. Jaime Bernal performed the experiments, authored or reviewed drafts of the paper, and approved the final draft. Michele De Biase performed the experiments, analyzed the data, authored or reviewed drafts of the paper, and approved the final draft. Ligia C. Muñoz-Arenas analyzed the data, authored or reviewed drafts of the paper, and approved the final draft. Francisco Erik González-Jiménez performed the experiments, prepared figures and/or tables, authored or reviewed drafts of the paper, and approved the final draft. Clemente Mosso analyzed the data, authored or reviewed drafts of the paper, and approved the final draft. Arit De León-Lorenzana analyzed the data, prepared figures and/or tables, authored or reviewed drafts of the paper, and approved the final draft. Carmine Fusaro conceived and designed the experiments, performed the experiments, analyzed the data, prepared figures and/or tables, authored or reviewed drafts of the paper, and approved the final draft. The following information was supplied regarding data availability:Raw data are available as a Supplemental File. Supplemental information for this article can be found online at http://dx.doi.org/10.7717/ peerj.11144#supplemental-information.