Abstract
The coronavirus disease 2019 (COVID-19) outbreak in North, Central, and South America has become the epicenter of the current pandemic. We have suggested previously that the infection rate of this virus might be lower in people living at high altitude (over 2,500 m) compared to that in the lowlands. Based on data from official sources, we performed a new epidemiological analysis of the development of the pandemic in 23 countries on the American continent as of May 23, 2020. Our results confirm our previous finding, further showing that the incidence of COVID-19 on the American continent decreases significantly starting at 1,000 m above sea level (masl). Moreover, epidemiological modeling indicates that the virus transmission rate capacity is lower in the highlands (>1,000 masl) than in the lowlands (<1,000 masl). Finally, evaluating the differences in the recovery percentage of patients, the death-to-case ratio, and the theoretical fraction of undiagnosed cases, we found that the severity of COVID-19 is also decreased above 1,000 m. We conclude that the impact of the COVID-19 decreases significantly with altitude.
Highlights
There is a negative correlation between altitude and COVID-19 incidence on the American Continent starting from 1,000 m above sea level.
The transmission rate of SARS-CoV-2 is lower in the highlands than in the lowlands.
The severity of COVID-19 decreases significantly with increased altitude.
1. Introduction
On March 11, 2020, coronavirus disease 2019 (COVID-19) was declared a pandemic by the World Health Organization (Joshua Berlinger, 2020; World Health Organization, 2020b). In late April, the health crisis began to ease in Asia and Europe (France 24, 2020; Langton, 2020; NPR, 2020), whereas case numbers began to rise in American countries. The first cases of COVID-19 on the American continent were reported in Canada on January 15th and in the United States on January 20th. Before February 25th, Brazil was the first affected country in Latin America. From that date until July 7th, the Pan-American Health Organization (PAHO) (Pan-American Health Organization, 2020) reported 6,004,685 confirmed cases and 268,828 deaths from COVID-19 in the American continent (Pan-American Health Organization, 2020). The spread of the virus was so fast in the region, that it has now become the epicenter of the disease, with United States (1st), Brazil (2nd), Peru (5th), Chile (6th), Mexico (9th), Colombia (19th) and Canada (20th) among the top twenty countries with the highest number of confirmed cases in the world on July 7th (Worldometer, 2020b).
Although American countries had more time and information to prepare for the pandemic than did Europe and Asia, with few exceptions, weaker public health systems, late political responses, and complex cultural and social conditions (i.e., poorly respected quarantines in most countries) have led to a major public health crisis on the continent. A crucial aspect for the spread of the virus is overcrowding in large cities. Clear examples of this fact are the cities of New York (8.1 million people - 10,194 people/km2 - 214,371 cases) (N. Y. C. Health Department, 2020), Montreal (1,8 million people - 4,517 people/km2 - 27,438 cases) (Montréal, 2020), and Rio de Janeiro (6, 7 million people - 5,597 people/km2 - 33,695 cases) (Ministerio de Saude de Brasil, 2020) (data as of July 7). Furthermore, this scenario is particularly critical in the poorer districts of those cities, where large numbers of people occupy the same housing units (Rede Brasil Atual, 2020). Remarkably, however, other densely populated metropolises with large slum areas, but located above 1,000 meters above sea level (masl), such as México city (8 million people - 26,000 people/km2 - 53,423 cases) (Gobierno de la Ciudad de México, 2020), Bogotá (7,4 million people - 4,310 people/km2 - 36,554 cases)(Alcaldía de Bogotá, 2020) and La Paz (2,3 million people - 2,676 people/km2 - 4,413 cases)(SEDES La Paz, 2020) (data as of July 7), seem to show lower incidences of COVID-19. This observation is of crucial importance since in American countries, more than 120 million people live at an altitude higher than 1,500 masl (defined as moderate altitude - MA), and more than 35 million people live at an altitude higher than 2,500 masl (defined as high altitude) (Hornbein & Schoene, 2001).
The impact of altitude on a potentially decreased virulence of SARS-CoV-2 was previously reported by our team. In fact, the global data analysis, which included detailed information from the Tibetan Autonomous Region of China, Bolivia, and Ecuador, suggested that the infection rates for the SARS-CoV-2 virus decrease significantly from 2,500 masl (Arias-Reyes, et al 2020). Subsequent reports from other research groups supported this observation (Huamaní, Velásquez, Montes, & Miranda-Solis, 2020; Ortiz-Prado et al., 2020; Rivero & Montoya, 2020). However, being aware that the course of the pandemic changes from day to day and that more detailed statistical analyzes are required, in this new study, we analyzed the epidemiological data from 23 countries in the American continent as of the 23rd of May. Our results show that the incidence of COVID-19, the virus transmission rate, and the severity of COVID-19 decrease significantly starting at 1,000 masl.
2. Methods
2.1. Data sources
Supplementary Table S1 shows the list of data sources used to collect information from COVID-19 cases through May 23rd, 2020. For Bolivia, Brazil, Canada, Colombia, French Guyana, Panama, and USA, we used information about the number of cases per city/county. For Argentina, Belize, Chile, Costa Rica, Cuba, Ecuador, El Salvador, Haiti, Honduras, Mexico, Paraguay, Peru, Puerto Rico, Dominican Republic, Uruguay, and Venezuela, we used information about the number of cases per state/province/departamento.
All locations with reported positive COVID-19 cases were associated with their respective geographic coordinates (latitude and longitude) using the OpenCage Geocoding API (GmbH, 2020). The altitude information for each location was extracted from the WORLDCLIM digital elevation model (Fick & Hijmans, 2017), and the population density data were assigned from the dataset of the CIESIN (Center for International Earth Science Information Network - CIESIN - Columbia University, 2018). In cases where the value of population density was zero, the information was retrieved from the national statistical institute of the corresponding country. The full list of these locations can be found at: https://doi.org/10.6084/m9.figshare.12685610.v1.
2.2. Assessment of incidence versus altitude
The number of COVID-19 cases by location was normalized by population density (inhabitants per square kilometer). These data were then organized by intervals of 100 meters of altitude. Finally, the natural logarithm (ln) of each normalized value was calculated. The correlation between the number of COVID-19 cases normalized and logarithmized (ln [cases/population density]) and the altitude was analyzed using a Pearson correlation analysis (n= 51). The difference in COVID-19 incidence below and above 1,000 masl at the continental level was calculated using a bidirectional random block ANOVA. The analyzed variable was the number of COVID-19 cases normalized and logarithmized (ln [cases/population density]); the grouping variable was the altitude (> 1,000 masl or < 1,000 masl) and the block variable was the country (n>1,000 masl= 827; n<1,000= 3,659). An identical correlation analysis was performed for each of the 23 countries of the American continent (nArgentina= 12; nBelize= 3; nBolivia= 27; nBrazil= 17; nCanada= 13; nChile= 14; nColombia= 36; nCosta Rica= 20; nCuba= 4; nDominican Republic= 9; nEcuador= 19; nEl Salvador= 8; nFrench Guyana= 1; nHaiti= 8; nHonduras= 12; nMexico= 15; nPanama= 13; nParaguay= 5; nPeru= 17; nPuerto Rico= 8; nUruguay= 2; nUSA= 34; nVenezuela= 12).
2.3. Evaluation of the virus transmission rate
The evaluation of the SARS-CoV-2 virus transmission rate was performed only for Argentina, Bolivia, Colombia, Ecuador, and Peru, as these countries applied similar strong early quarantines and provided daily epidemiological data at state/province/departamento level.
The COVID-19 data retrieved daily by state/province/departmento in each country was classified into two groups: highlands (>1,000 masl) and lowlands (<1,000 masl). The SEIR model (Susceptible - Exposed - Infectious - Removed) was used to calculate the estimated number of susceptible, exposed, infected and removed (recovered + deceased) individuals (Abadie, Bertolotti, & Arnab, 2020). The number of susceptible individuals was calculated as the total population minus the number of infected and exposed subjects. The initial number of infected was set to 1, and the number of exposed subjects was arbitrarily set to 30 (according to (Birbuet & López, 2020). The contact rate (β), recovery rate (γ), and the rate at which exposed individuals become infected (ε) were calculated as conducted elsewhere (Birbuet & López, 2020). The contact rate was estimated as the product of the “interaction frequency” and the “probability of transmission of the disease” (transmission rate) (Birbuet & López, 2020). Finally, the epidemiological behavior of the disease was modeled, calculating the best transmission rate for highlands (>1,000 m) and for lowlands (<1,000 m).
2.4. Assessment of COVID-19 severity
The severity of the disease was estimated based on the percentage of recovered patients and the death-to-case ratio. Indeed, a lower mortality rate per case and a higher percentage of recovered patients suggest a lower severity of the disease (and the opposite). Accordingly, the severity of the disease was estimated for Argentina, Bolivia, Colombia, Ecuador, and Peru, as these countries applied similar containment measures and provided daily epidemiological data at the state/province/departamento level.
The death-to-case ratio and the percentage of recovered patients ([recovered patients/reported cases] * 100) for each country (except Ecuador) were calculated using the data from the last 10 days evaluated (from May 13th to 23rd). The underlying rationale is that these parameters stabilize as the pandemic progresses. As this information was not available for Ecuador on the mentioned dates, these values were calculated only from the data of April 29 (the latest data available for this country).
The Wilcoxon signed-rank test was used to determine if the percentage of recovered patients (npairs= 5) and the death-to-case ratio (npairs= 5) were different between the highlands and the lowlands.
2.5. Assessment of undiagnosed cases
COVID-19 patients can be asymptomatic or symptomatic. According to the severity of the disease, symptomatic patients are classified as follows: 1) Mild: patients with very mild symptoms and without evidence of viral pneumonia or hypoxia; 2) Moderate: patients with signs of pneumonia (fever, cough, dyspnea, rapid breathing but with regular arterial O2 saturation values); 3) Severe: patients with established lung damage, symptoms of pneumonia, cough, fever, dyspnea and hypoxia; and 4) Critical: patients with systemic extrapulmonary inflammation and the presence of large amounts of proinflammatory cytokines (World Health Organization, 2020a). Based on the fact that in most American countries, health policies restrict the diagnosis of COVID-19 to suspects with clear symptoms of infection and/or a proven history of contact with infected people (Aguirre, 2020; Americas Society Council of Administration, 2020; Galindo, 2020), in this study, we assumed that the cases observed and reported officially include mainly symptomatic (mild + moderate + severe + critical). We performed a theoretical calculation to determine the number of undiagnosed cases (asymptomatic + undiagnosed symptomatic). To do so, we used the infection mortality rates (IFR) of New York (IFRNY = 1.4 - considered so far to be the most reliable calculated value available for this parameter (Worldometer, 2020a)) to calculate the “Estimated Total cases” (= observed deaths*IFRNY) (Jombart et al., 2020) and the “Estimated Fraction of Undiagnosed cases” (= [calculated total cases - observed cases]/calculated total cases) (Worldometer, 2020a); Verity et al. (2020)). Data on the observed deaths and recovered patients were obtained from each country’s official government website on May 23rd.
2.6. Evaluation of the impact of quarantines and self-isolation on the highland populations of Argentina, Bolivia, Colombia, Ecuador, and Peru
The impact of quarantines and self-isolation procedures was estimated for Argentina, Bolivia, Colombia, Ecuador, and Peru, as these countries applied early strict quarantines and self-isolation procedures. To do so, the number of people infected in a quarantine-free environment was estimated with SEIR epidemiological models (described above), using a “frequency of interaction” value equal to 30 (according to (Birbuet & López, 2020)). Next, these values were compared to the number of people infected obtained by using “frequency of interaction” values for each country calculated according to (Birbuet & López, 2020): Argentina (8.18), Bolivia (8.1), Colombia (9.03), Ecuador (8.7) and Peru (11.3).
3. Results
3.1. The incidence of COVID-19 is dependent on the population density in the American continent
The correlation between the number of positive cases and the population density for each location was evaluated at the continental level (23 countries). Our results show a significant positive correlation (Pearson’s correlation, p <0.0001; R2 = 0.004) between these two variables (S2). These results are in line with previous studies showing that the contagion rate and, therefore, the number of COVID-19 cases are positively correlated with population density. (Engle, Stromme, & Zhou, 2020; Jason Barr, 2020; Rocklöv & Sjödin, 2020; Rogers, 2020). As such, the number of COVID-19 cases normalized by population density (number of cases/population density) was used in the following analyses.
3.2. The incidence of COVID-19 decreases above 1,000 masl in the American countries
The correlation between the number of COVID-19 cases (normalized, logarithmized, and organized by intervals of 100 meters of altitude) and the altitude at the continental level (data from 23 countries together) was analyzed using a Pearson correlation analysis. Our results revealed a strong negative correlation (p<0.0001; R2= 0.603) between these variables, underlining a decrease in the incidence of COVID-19 cases with increasing altitude (see Fig. 1a). Furthermore, we remarked that a significant decrease in the number of normalized and logarithmized (but not yet organized by intervals of 100 meters of altitude) COVID-19 cases occurs already approximately 1,000 masl, with a further significant decrease at higher altitudes (Fig. 1b). Repeated correlation analysis performed at altitudes above 800, 1,000, 1,500 and 2,500 masl confirmed this observation. Nonsignificant correlation was found for data below 1,000 masl (Pearson correlation, p=0.568; R2= 0.042) (Fig. 1c), while a strongly significant correlation between COVID-19 incidence (ln [cases/population density]) and altitude was obtained above 1,000 masl (Pearson correlation, p<0.0001; R2= 0.427) (Fig. 1d). Furthermore, considering that various factors, such as public health policies, diagnostic strategies, confinement rules, and cultural aspects, may influence the number of reported cases of COVID-19, we performed a randomized block design ANOVA test to compare the number of COVID-19 cases (normalized, logarithmized, and organized by intervals of 100 meters of altitude) above and below 1,000 m. The advantage of this type of statistical analysis is that it considers the internal variability of each country in the incidence analysis. In line with our previous results, this test revealed a significant difference between the number of COVID-19 cases observed above versus below 1,000 masl (RBD-ANOVA F = 5,273; df= 45; p = 0.022;).
In the next step, we wanted to test whether the negative correlation found between the incidence of COVID-19 and the altitude for the American continent can be independently reproduced in each of the 14 American countries (of 23 analyzed) that reported cases at more than 1,000 masl. Our results showed significant negative correlations in 9 of 14 countries (Table 2) (S3): Argentina (p=0.003; R2= 0.597), Brazil (p<0.0001; R2= 0.789), Canada (p<0.012; R2= 0.012), Colombia (p<0.0001; R2= 0.376), Costa Rica (p=0.004; R2= 0.381), Ecuador (p<0.0001; R2= 0.647), Mexico (p=0.002; R2= 0.286), Peru (p=0.027; R2= 0.285) and USA (p<0.0001; R2= 0.577).
Taken together, these data confirm our previous observations (Arias-Reyes et al., 2020), suggesting that the virulence of SARS-CoV-2 gradually decreases with altitude. However, although we previously suggested that the incidence of COVID-19 decreases at high altitude (2,500 masl), the data of the present study show a significant decrease in the disease incidence starting at 1,000 m of altitude.
3.3. The transmission rate of SARS-CoV-2 is lower in the highlands compared to lowlands
To investigate whether the transmission rate of SARS-CoV-2 differs between highlands (>1,000 masl) and lowlands (<1,000 masl), we used SEIR epidemiological models only in countries that applied similar strong early quarantines and provided daily epidemiological data at state/province/departamento level: Argentina, Bolivia, Colombia, Ecuador, and Peru. To do so, the “transmission rate” values that allowed the most faithful theoretical reproduction of the epidemiological curves observed in the highlands and lowlands were obtained (Fig. 2). Our results show that compared to lowlands (Argentina=4.2%, Bolivia=4%, Colombia 3.7%, Ecuador=3.8%, Peru=4%), lower “transmission rates‥ at highlands lead to better modeling data of the real epidemiological curves (Argentina=2.5%, Bolivia=3%, Colombia 3.6%, Ecuador=3.6%, Peru=3%). These results strongly support the hypothesis of decreased SARS-CoV-2 virulence in highlands compared to lowlands.
Of note, the “transmission rate” values between the highlands and the lowlands in Colombia show a small difference. This result is explained by the fact that 60% of its population (∼ 30 million inhabitants) is settled above 1,000 m altitude, which are also the most densely populated areas. This causes the requirement of a higher infection probability value in the SEIR model to fit the real data, thus approaching the values of the transmission rate between highlands and lowlands.
3.4 The severity of COVID-19 is reduced in highlands compared to lowlands
Classically, the estimation of the severity of a disease is performed by a comparison between recovery rates and infection mortality rates (IFR). However, since this information is still limited for COVID-19, we used a theoretical rationalization approach to estimate this parameter. We evaluated the differences in the percentage of recovered patients (from the total reported cases) and in the death-to-case ratio (deaths/total reported cases) as indicators of the recovery rate and the IFR for Argentina, Bolivia, Colombia, Ecuador, and Perú (countries that applied similar strong early quarantines and provided daily epidemiological data at the state/province/departamento level). We found a significantly higher percentage of recovered patients in the highlands (Wilcoxon signed-rank test p= 0.031) versus lowlands, suggesting higher a recovery rate in the highlands versus lowlands. On the other hand, our results did not show significant differences of death-to-case ratio between the highlands and the lowlands (Table 3). This may occur only in the case in which the number of undiagnosed cases (asymptomatic + undiagnosed symptomatic) is not similar between the highlands and lowlands. Indeed, our theoretical approach for the assessment of undiagnosed cases showed approximately 76% of undiagnosed cases in the highlands and 73% of undiagnosed cases in the lowlands (1.04 ± 0.12-fold higher) (Table 4).
Although, these values are theoretical (and should be interpreted carefully), they suggest that the severity of COVID-19 is lower at highlands compared to lowlands.
3.5 Strict quarantine and self-isolation measures at altitude significantly decreased the rapid progression of COVID-19
To close our analysis, we wanted to assess the impact of quarantine and social distancing measures in countries that applied early strict measurements, such as Argentina, Bolivia, Colombia, Ecuador, and Peru. By using SEIR epidemiological models, we compared the number of people infected under conditions in which the “frequency of interaction” value is associated with a quarantine-free environment (= 30) to the number of infected people under conditions in which the “frequency of interaction” value is associated with a strict quarantine environment: 8.18 for Argentina, 8.1 for Bolivia, 9.03 for Colombia, 8.7 for Ecuador, and 11.3 for Peru (11.3). Our results show that the application of containment measures drastically delayed the number of COVID-19 infections (Figure 4).
4. Discussion
By expanding our analysis to 23 countries of the American continent, this work clearly demonstrates our previous observations suggesting that the virulence of SARS-CoV-2 decreases significantly with altitude. Indeed, in a previous study we reported that global data analysis (that included detailed information from the Tibetan Autonomous Region of China, Bolivia, and Ecuador) suggested that COVID-19 decreases significantly at high altitude (2,500 masl) (Arias-Reyes, et al 2020). While subsequent reports from other teams confirmed this observation (Accinelli & Leon-Abarca, 2020; Huamaní et al., 2020; Ortiz-Prado et al., 2020; Rivero & Montoya, 2020), in this new study we appropriately corrected the COVID-19 data by population density. Remarkably, our results show that a clear turning point in the incidence of the disease occurs at 1,000 masl. Moreover, independent examination from all the American countries that reported cases at more than 1,000 masl (14 out of 23) confirmed this result, except for four nations: Bolivia, in which the low number of cases reported for moderate-altitude regions (1,000 – 2,500 masl) masked a possible correlation; Chile, in which its long and narrow geography (flanked by the Andean mountains on the east and the Pacific Ocean) does not allow obtaining precise altitude values for COVID-19 cases; and Dominican Republic, Haiti and Panamá, in which the highest altitude with reported cases is just over 1,000 masl. Therefore, these data may have a low weight in the statistical evaluation.
To better understand this result, we performed additional theoretical approaches to determine the capability of the virus to move from one host to another (the transmission rate of the virus). Our results show that this parameter decreases significantly in the highlands compared to lowlands, suggesting that environmental factors influence the virulence of SARS-CoV-2 at above 1000 m. Indeed, as altitude increases, the environment is characterized by more drastic changes in temperature between night and day, higher air dryness, and higher levels of ultraviolet (UV) light radiation. (United-States-Environmental-Protection-Agency, 2017). In particular, UV light radiation was suggested to be an important natural sanitizer at altitude that may shorten the half-life of any given virus (Andrade, 2020; G. Zubieta-Calleja, 2020; G. R. Zubieta-Calleja & Zubieta-DeUrioste, 2017). In addition, the solar radiation is also more intense at altitude, and a recent study reported that this factor may also be a key factor leading to the deactivation of the virus (Asyary & Veruswati, 2020; Li, Li, Zhang, & Liu, 2020). Taken together, these factors may lead to a gradual reduction of the “survival” and “virulence” capacity of the virus as altitude progresses. Finally, the size of the virus inoculum in the air should gradually decrease as the barometric pressure decreases and the distance between the air molecules increases.
Apart from environmental factors that may decrease the transmission capacity of the virus, our results also suggest that physiological mechanisms associated with a prolonged exposure to barometrical hypoxia help to decrease the severity of the infection. Thus far, two factors have been suggested that may be involved in this phenomenon: A decreased expression of the virus’s gateway to the body, angiotensin converting enzyme 2 (ACE2) (Arias-Reyes, et al, 2020), and an increased level of erythropoietin (EPO) (Soliz, et al., 2020). Stabilization of the hypoxia inducible factor 1 alpha (HIF1-a - a master regulator of the response to hypoxia) may lead the regulation of both parameters. Indeed, it was shown that exposure of human pulmonary artery smooth muscle cells (hPASMC) to hypoxia markedly decreases the expression of ACE2 (Zhang et al., 2009). Similar results were obtained in heart cells of Sprague Dawley rats after 28 days exposure to conditions equivalent to 10% O2 hypoxia (Dang et al., 2020). Furthermore, it was shown that changes in oxygen availability as small as 2% are sufficient to induce HIF activation in many tissues (Jiang, Semenza, Bauer, & Marti, 1996; Stroka et al., 2001). In the same direction, HIF is the main booster of EPO production in the kidney and other tissues, and it was observed that the HIF-related oxygen sensing mechanism precisely senses altitude differences of 300 m even in low to moderate altitudes (Gassmann et al., 2019). While EPO is the central factor leading the stimulation of red blood cells (Martin & Windsor, 2008; Storz & Moriyama, 2008), EPO also promotes adaptive cellular responses to hypoxic challenges and tissue-damaging insults in nonhematopoietic tissues (Juul & Felderhoff-Mueser, 2007). As such, in the infectious context induced by the SARS-CoV-2 virus, since EPO may help improve oxygen delivery to the tissues, it may also protect other tissues from a multiple organ dysfunction and inflammation (Chen et al., 2020; Yang et al., 2020).
Although the results of this work clearly indicate that the virulence of SARS-CoV-2 decreases with altitude, these results should not be misinterpreted either by citizens or by government systems responsible for imposing sanitary restrictions. Indeed, to make this point clear, we performed calculations to determine the impact of quarantine and self-isolation measures in highlands of countries that applied early strict measurements, such as Argentina, Bolivia, Colombia, Ecuador, and Peru. The results of this analysis show that without relevant restriction measures, the number of currently reported cases of COVID-19 would be increased by 10,000 times (Figure 4).
Finally, the limitations of the study include possible failures and delays in reporting cases as well as that our analysis does not consider the effect of different risk groups such as age, sex, comorbidities, and the possibility of reinfection.
In conclusion, epidemiological analyses of the present work strongly suggest that the incidence of COVID-19 significantly decreases with altitude with a turning point at 1,000 masl. This effect seems to be related both to a decrease in the transmission capacity of the virus and to the physiological characteristics of altitude residents. Finally, our results suggest that knowledge of the mechanisms of physiological acclimatization to hypoxia may help to better understand the viral nature of SARS-CoV-2 as well as facilitate the discovery of new treatment strategies.
Data Availability
Raw, normalized, and adjusted data of COVID-19 cases are available and registered at figshare. Raw epidemiological daily data of COVID-19 of Argentina, Bolivia, Colombia, Ecuador, and Peru are available and registered at figshare.
DATA AVAILABILITY
Raw, normalized, and adjusted data of COVID-19 cases are available at https://doi.org/10.6084/m9.figshare.12685478
Raw epidemiological daily data of COVID-19 of Argentina, Bolivia, Colombia, Ecuador, and Peru are available at https://doi.org/10.6084/m9.figshare.12685523
AUTHOR CONTRIBUTIONS
C.A.-R., L.P.-M., F.A.-R., and F.C.-R. gathered the epidemiological and geographical data. F.C.- R. processed the geographical data and created the distribution maps. C.A.-R. and J.S. analyzed the data. J.S., C.A.-R., L.P.-M., F.A.-R, D.A.M., N. Z.D., R.A.-T., E.M.S.-G., G.Z.-C., and M.D. wrote this manuscript.
COMPETING INTERESTS
The authors have no conflict of interests to declare.
Supplementary material
S1. Table of sources of epidemiological data
S2. Effect of population density on the incidence of COVID-19 cases in the American continent. Epidemiological data were retrieved on May 23. Data on population density were extracted from the dataset created by CIESIN (Center for International Earth Science Information Network - CIESIN - Columbia University, 2018) or the corresponding country’s national statistics institute on May 23.
S3. Effect of altitude on the incidence of COVID-19 cases in American countries. Epidemiological data were retrieved on May 23. Data on population density were extracted from the dataset created by CIESIN (Center for International Earth Science Information Network - CIESIN - Columbia University, 2018) or the corresponding country’s national statistics institute on May 23. Data were normalized by the population density of the same location, adjusted by calculating the natural logarithm (ln) of each value, and summed in intervals of 100 m of elevation. Raw, normalized, and adjusted data are available at https://doi.org/10.6084/m9.figshare.12685478
ACKNOWLEDGMENTS
The authors wish to acknowledge Dr Gonzalo Leonardini for his advice and assistance in the mathematical evaluation of our models, the platform CovidBot Peru for giving us access to their database of systematized daily epidemiological information from Peru, and Teresa Roncal for her inputs and thoughts on the redaction of the manuscript.
Footnotes
The name of one of the authors was updated