Influence of socio-ecological factors on COVID-19 risk: a cross-sectional study based on 178 countries/regions worldwide ======================================================================================================================== * Dai Su * Yingchun Chen * Kevin He * Tao Zhang * Min Tan * Yunfan Zhang * Xingyu Zhang ## Abstract **Background** The initial outbreak of COVID-19 caused by SARS-CoV-2 in China in 2019 has been severely tested in other countries worldwide. We aimed to describe the spatial distribution of the COVID-19 pandemic worldwide and assess the effects of various socio-ecological factors on COVID-19 risk. **Methods** We collected COVID-19 pandemic infection data and social-ecological data of 178 countries/regions worldwide from three database. We used spatial econometrics method to assess the global and local correlation of COVID-19 risk indicators for COVID-19. To estimate the adjusted incidence rate ratio (IRR), we modelled negative binomial regression analysis with spatial information and socio-ecological factors. **Findings** The study indicated that 37, 29 and 39 countries/regions were strongly opposite from the IR, CMR and DCI index “spatial autocorrelation hypothesis”, respectively. The IRs were significantly positively associated with GDP per capita, the use of at least basic sanitation services and social insurance program coverage, and were significantly negatively associated with the proportion of the population spending more than 25% of household consumption or income on out-of-pocket health care expenses and the poverty headcount ratio at the national poverty lines. The CMR was significantly positively associated with urban populations, GDP per capita and current health expenditure, and was significantly negatively associated with the number of hospital beds, number of nurses and midwives, and poverty headcount ratio at the national poverty lines. The DCI was significantly positively associated with urban populations, population density and researchers in R&D, and was significantly negatively associated with the number of hospital beds, number of nurses and midwives and poverty headcount ratio at the national poverty lines. We also found that climatic factors were not significantly associated with COVID-19 risk. **Conclusion** Countries/regions should pay more attention to controlling population flow, improving diagnosis and treatment capacity, and improving public welfare policies. Keywords * socio-ecological factors * COVID-19 risk * cross-sectional study * 178 countries/regions worldwide ## 1. Introduction The novel coronavirus disease (COVID-19) that has spread to more than one hundred countries and killed hundreds of thousands of people has officially been categorized as a pandemic by the World Health Organization. The initial outbreak of COVID-19 caused by SARS-CoV-2 in China in 2019 has been severely tested in other countries worldwide. As of April 6, 2020, COVID-19 had infected 1,345,048 patients in 184 countries/regions and caused 74,565 deaths, as countries worldwide responded to a human-to-human respiratory disease pandemic caused by COVID-19. Emerging infectious diseases (EIDs), such as SARS and COVID-19, pose a vast economic and public health burden worldwide [1,2]. COVID-19 not only seriously endangers people’s life safety and health but also greatly affects economic globalization. To address the challenges posed by COVID-19, the links among the transmission of COVID-19, socio-economic factors and climatic factors must be understood to suggest better strategies for predicting, preventing, coping with and mitigating the associated challenges. Simultaneously, given that the climate and socio-economic context are unlikely to change in the short term, it is easier to intervene accordingly [3]. The spread of many EIDs has been reported to be influenced by socio-ecological factors, including socio-economic and climate factors [1,2,4–7]. Previous studies have found that climatic conditions limit the geographical and seasonal distribution of EIDs, and weather affects the timing and intensity of outbreaks [8–12]. In addition, whereas climate patterns may control the potential global distribution of EIDs, the actual size and spatial scope of a region may be controlled by several non-climatic factors associated with transmission, including epidemiological, socio-economic and demographic factors [13–18]. However, research on the climatic and socio-economic drivers of COVID-19 transmission remains lacking, especially regarding the effects of socioeconomic factors and the total effects of socio-ecological factors. Ignoring important non-climatic factors or other confounding factors (such as urban development, economic growth, poverty, health, infrastructure, science and technology, social security and labor) would overestimate the effects of climate change. Therefore, studying the influence of socio-ecological factors on the transmission risk of EIDs is highly important. For most EIDs, three elements are essential: an agent (or pathogen), a host (or vector) and the environment of transmission [19]. Appropriate climatic and weather conditions are necessary for the survival, reproduction, distribution and transmission of disease pathogens, vectors and hosts. Therefore, changes in climate or weather conditions may affect EIDs by affecting pathogens, vectors, hosts and their living environments [19–21]. Although many climate variables may influence the transmission of EIDs, some studies have shown that changes in the four main variables have the greatest effects on infected diseases with strong environmental components (temperature, precipitation, relative humidity (RH) and wind) [22–26]. In recent studies, although the severity of some cases of COVID-19 has mimicked that of SARS-CoV cases [27–30], the reproductive number (average R0=3.28) of COVID-19 is higher than that of SARS-CoV; therefore, considering the climate and environment may improve understanding of the pathogen’s vectorial capacity and basic reproduction number, and the risk of transmission of COVID-19 [31,32]. In recent decades, many rapid and pronounced changes in human social ecology have altered the likelihood of the emergence and spread of infectious diseases [33–35]. These changes include increases in population size and density; urbanization; persistent poverty (especially in the expansion of urban slums); the number and movement of political, economic and environmental refugees; differences in infrastructure and science and technology; and poor health awareness [36]. The socio-economic environment contributes significantly to the health of individuals as well as communities [37] and is the root cause of health and health equity. These socio-economic drivers have contributed to the shifting global ecology of vector transmission that enabled COVID-19 to emerge worldwide, by dangerously uniting the human hosts, vectors and pathogen. Socioeconomic changes interact with environmental changes in promoting EID spread and increase the harm of EIDs to humans. The purpose of this study was to describe the spatial distribution of the COVID-19 pandemic worldwide, and assess the effects of different socio-ecological factors, including climate and socio-economic factors, on COVID-19 risk in 178 countries/regions worldwide, including incidence rate (IR), cumulative mortality rate (CMR) and daily cumulative index (DCI). In addition, this study analyzed intervention policies in different countries and regions to establish early warning and decision support systems and provide guidance for COVID-19 management in different countries/regions. ## 2. Methods ### 2.1. Concept model According to previous research [38–41], we established the Potential Risk Assessment Framework for COVID-19 (Figure 1). The influence of global socio-ecological factors (climate and socio-economic factors) on the risk of COVID-19 can be tested by its influence on the following three disease components: agent (or pathogen), host (or vector) and the environment of transmission. A combination of natural and human influences led to the COVID-19 pandemic. We used three main variables to assess the potential risk of COVID-19: IR, CMR and DCI. ![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/04/2020.04.23.20077545/F1.medium.gif) [Figure 1.](http://medrxiv.org/content/early/2020/05/04/2020.04.23.20077545/F1) Figure 1. Climate, socio-economic conditions and COVID-19 transmission ### 2.2 Definitions of different cases for COVID-19 A confirmed case of COVID-19 infection was defined by laboratory confirmation of the virus causing COVID-19 infection, regardless of clinical signs and symptoms [42–44]. However, some reported case numbers from China have included people with symptoms of COVID-19 without laboratory confirmation. The definitions of COVID-19 related deaths differ across countries. In Italy, any death of a person with positive reverse transcriptase–polymerase chain reaction (RTPCR) testing for SARS-CoV-2 is considered COVID-19 related. ### 2.3 Data collection #### 2.3.1 Outcome variables A dashboard published and hosted by researchers at the Center for Systems Science and Engineering, Johns Hopkins University (JHU-CSSE) [45] shows the numbers and locations of confirmed COVID-19 cases, deaths and recoveries in all affected countries. All collected data on COVID-19 from the Johns Hopkins University are made freely available by the researchers through a GitHub repository. All manual updates (for countries and regions outside mainland China) are coordinated by a team at Johns Hopkins University. We extracted the global time series data of confirmed and recovered cases and deaths due to COVID-19 from the JHU-CSSE GitHub repository. The data were recorded from January 22, 2020 and were updated once daily around 23:59 (UTC). We selected the cross-sectional data from April 6, 2020. On the basis of the availability of the data, we extracted 178 countries from the database (excluding countries/regions without COVID-19 cases and some unmatched countries/regions, such as Taiwan, China). The first-level geographical unit of the dataset is the country/region, and the second-level geographical unit is the province/state. We uniformly selected first-level geographical units (countries/regions). In addition, we further classified the countries/regions in the data set according to UN geographical divisions and divided the countries with epidemic COVID-19 into 20 regions. As previously mentioned, for outcome variables, we selected IR, CMR and DCI as indicators to measure COVID-19 risk. The specific calculation process is shown in Table 1. View this table: [Table 1](http://medrxiv.org/content/early/2020/05/04/2020.04.23.20077545/T1) Table 1 The calculation process of outcome variables for COVID-19 IR was used to describe the distribution of COVID-19, explore the etiological factors, propose an etiological hypothesis, and evaluate the efficacy of detection and prevention measures. CMR reflects the total deaths due to COVID-19 and is an indicator of the risk of death from COVID-19. DCI mainly describes the growth rate of COVID-19 in different countries/regions and is a measure of the risk of disease transmission. The World Health Organization, on March 11, 2020, declared the COVID-19 outbreak a global pandemic, thus indicating that COVID-19 had broadly spread worldwide. Therefore, when we considered IR, CMR and DCI in different countries/regions, these measures reflected not only the rapid growth in the number of people infected with COVID-19 but also the detection level in the entire country/region, which was used to identify more people infected with COVID-19. #### 2.3.2 Climate data We obtained daily meteorological observation values from the Global Surface Summary of the Day (GSOD) via The Integrated Surface Hourly (ISH) dataset. The ISH dataset includes global data obtained from the USAF Climatology Center, which is located in the Federal Climate Complex with NCDC. GSOD comprises 12 daily averages computed from global hourly station data. Except in United States stations, 24-hour periods are based on UTC times. The latest daily summary data are normally available 1–2 days after the date-time of the observations used in the daily summaries. More than 9,000 stations’ data worldwide are typically available. Daily weather elements include mean values of temperature, dew point temperature, sea level pressure, station pressure, visibility, wind speed, maximum and minimum temperature, maximum sustained wind speed and maximum gust, precipitation amount, snow depth and weather indicators. However, we chose the climate data from April 6, 2020 and selected four variables from the GSOD dataset that significantly affected COVID-19 risk: (1) mean temperature (.1 Fahrenheit); (2) mean dew point (.1 Fahrenheit); (3) mean wind speed (.1 knots); and (4) precipitation amount (.01 inches). The reason for extracting the average dew point variable was to calculate the RH value by using this variable and the temperature variable. The temperature and dew point in Celsius were used to calculate the RH according to the temperature and dew point at each time point [46]: ![Formula][1] where *C* is the temperature in Celsius, *F* is the temperature in Celsius, *D* is the mean dew point for the day in Celsius, *T* is the mean temperature for the day in Celsius, and *e* is the base of the natural log. #### 2.3.3 Socio-economic data Indicators of socio-economic factors affecting the spread of COVID-19 were derived from the World Development Indicators dataset, the primary World Bank collection of development indicators, which compiles relevant, high-quality and internationally comparable statistics about global development and the fight against poverty. The database contains 1,600 time series indicators for 217 economies and more than 40 country groups, and data for many indicators cover a period of more than 50 years. As shown in Table 2, we selected 32 indicators affecting COVID-19 risk in seven dimensions in 178 countries/regions. The index value for 2019 was taken as the priority for each indicator. If the index was missing in 2019, the index value of the most recent year was selected as a substitute. View this table: [Table 2.](http://medrxiv.org/content/early/2020/05/04/2020.04.23.20077545/T2) Table 2. Socioeconomic indicators influencing the spread of COVID-19 ### 2.4 Statistical analysis #### 2.4.1. Spatial econometrics method First, we used Moran’s I to measure the global correlation of COVID-19 risk indicators [47]. Global Moran’s I is a measure of global spatial autocorrelation, and the value of Moran’s I usually ranges from −1 to +1. Values significantly below -1/(N-1) indicate negative spatial autocorrelation, and values significantly above -1/(N-1) indicate positive spatial autocorrelation. If significant global spatial autocorrelation was found, we then used local indicators of spatial autocorrelation (LISA) to evaluate the locations of COVID-19 clusters. The meaning of local Moran’s Ii is similar to that of global Moran’s I. A positive Ii indicates that the high (or low) value of region i is surrounded by the surrounding high (or low) value; A negative Ii indicates that the high (or low) value of region i is surrounded by the surrounding low (or high) value. The general models are described in Eq. 1–2. ![Formula][2] ![Formula][3] where *n* is the number of spatial units indexed by *i* and *j*, *x* is the variable of social ecology factors, ![Graphic][4] is the mean of *X*, and *wij* is a matrix of spatial weights with zeroes on a diagonal (i.e., *wii*=0). Second, to better approximate the real infectious disease spatial spread process, we fit a one- order spatial autoregressive regression model comprising spatial lags. Because we believed that COVID-19 risk transmission in a certain country/region would be different for neighboring countries/regions, we sought to reflect this difference in the model. A one-order spatial autoregressive process takes the form (Eq. 3) [48]: ![Formula][5] where *δ* is the spatial autoregressive coefficient, *W* is the *i*, *j*-th element of the exogenous,non-negative *N*×*N* spatial weight matrix with zero diagonal elements that describes the arrangement of the spatial units in the countries/regions, and *ε* is i.i.d. innovations with zero mean and finite variance *σ*2. For simplicity, in this paper, we assumed that the spatial weight matrix *W* was non-standardized and also used a queen spatial weight matrix. #### 2.4.2. Processing of missing values Before negative binominal regression, the k-nearest neighbors (k-NN) approach was used to impute missing data for some socio-ecological variables. For a given patient with missing values, the k-NN method identifies the k-nearest countries/regions on the basis of Euclidean distance. Using these countries/regions, we then replaced missing values by using a majority vote for discrete variables and weighted means for continuous features. One advantage of using this method is that missing values in all features are imputed simultaneously without the need to treat features individually. #### 2.4.3. Negative binomial regression First, we established the correlation matrix of socio-economic factors to check for multicollinearity. If there was a strong correlation (> 0.8) among socio-economic factors, then we removed the factor with strong correlation with other variables. Then, the incidence rate ratio (IRR) of each socio-ecological factor was calculated with single factor negative binomial regression analysis, that is, the effect of each socio-ecological factor on COVID-19 risk by changing the average COVID-19 risk value by a specific unit quantity. The spatial autoregressive models comprising spatial lags, which were a weighted average of observations on the diseases over neighboring units, were input into the model to adjust for spatial variation in COVID-19 risk. Modeled values of climate factors were centered on the mean values for each station in every country/region [49]. The factors with P < 0.05 were included in the multi-factor negative binomial regression analysis with spatial information to calculate the adjusted IRR (aIRR). The general model is described in Eq. 4. ![Formula][6] where *yi* denotes the daily counts of COVID-19 risk indicators in county/region *i*; *WY* represents spatial lags, and *W*; is the spatial weight; *Sn* represents socio-economic factors (all variables are in Table 2); *Tt* is the mean temperature in county/region *i*; *Ht* is the RH in county/region *i*; *Mt* is the wind speed in county/region *i*; *Pt* is the precipitation amount in county/region *i*; and *εi* is a random intercept. Sensitivity analyses with maximum and minimum temperatures instead of average temperatures were also conducted with the same procedures, in which we used the same non-informative priors for the minimum and maximum temperatures [49, 50]. All statistical analyses were performed in Stata statistical software Version 15, and p-values were two-tailed, with statistical significance set at.05. ArcMap 10.7 and Geoda software were used to process basic geographic information. Data visualization mainly used RStudio software Version 1.2.5033. ## 3. Results ### 3.1 Characteristics of 178 countries/regions with reported cases of COVID-19 As of April 6, 2020, a total of 178 countries/regions worldwide had reported data and were included in this study (Table S1). The three countries/regions with the highest IR worldwide were Andorra (Southern Europe, IR=313.80), Iceland (Northern Europe, IR=215.90) and Gibraltar (United Kingdom) (Southern Europe, IR=178.52). The three countries/regions with the highest CMR worldwide were San Marino (Southern Europe, CMR=947.17), Spain (Southern Europe, IR=285.53) and Italy (Southern Europe, IR=273.42). The three countries/regions with the highest DCI worldwide were the United States (North America, DCI=4823.87), Spain (Southern Europe, DCI=2070.83) and Italy (Southern Europe, DCI=1978.31). View this table: [Table S1](http://medrxiv.org/content/early/2020/05/04/2020.04.23.20077545/T5) Table S1 Characteristics of 178 countries/regions with reported cases of COVID-19 as of 6 April 2020 ### 3.2 Spatial clustering evaluation for COVID-19 #### 3.2.1 Test results for global spatial correlation The number and distribution of first-order neighbors in different countries/regions are shown in Figure 2. The number of neighboring countries/regions was mainly concentrated in 0-6, accounting for 87.08% of the total number of neighboring countries. Among them, China and Russia had the largest number of neighboring countries/regions. ![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/04/2020.04.23.20077545/F2.medium.gif) [Figure 2.](http://medrxiv.org/content/early/2020/05/04/2020.04.23.20077545/F2) Figure 2. Number of first-order neighbors in different countries/regions Table S2 shows the global spatial Moran’s I indexes of the IR, CMR and DCI of 178 countries/regions worldwide according to the one-order spatial contiguity matrix. The Moran’s I indexes of IR, CMR and DCI were all positive, and all index values were significant at the level of 1%, thus indicating that the IR, CMR and DCI of 178 countries/regions worldwide had strong spatial aggregation effects. Meanwhile, the Moran’s I index of different indicators showed significant differences, thus indicating to some extent that IR, CMR and DCI have different aggregation effects in different countries/regions. View this table: [Table S2.](http://medrxiv.org/content/early/2020/05/04/2020.04.23.20077545/T6) Table S2. The Global Moran’s I index of IR, CMR and DCI for COVID-19 #### 3.2.2 Test results for local spatial correlation A total of 37 countries/regions were strongly opposite from the IR index “spatial autocorrelation hypothesis,” including 11 countries/regions with high–high patterns, mainly concentrated in Western Europe, southern Europe and Canada; 24 countries/regions with low–low patterns, mainly concentrated in parts of Africa, parts of Asia (China, India, Laos); Cuba with a low–high pattern; and Djibouti with a high–low pattern. Simultaneously, 29 countries/regions were strongly opposite from the CMR index “spatial autocorrelation hypothesis,” which was suitable for 10 countries with high–high patterns, mainly in Western Europe and southern Europe (e.g., France, Italy and Spain); 16 countries/regions with low–low patterns, mainly concentrated in parts of Africa and China; and three countries with low–high patterns, including Morocco and Slovenia. In addition, 39 countries/regions strongly did not support the hypothesis of “no spatial autocorrelation” of the DCI index, among which six countries/regions had high–high patterns (Canada, France, Portugal, Belgium, the Netherlands, and Switzerland), and 22 countries/regions had low–low patterns, mainly in parts of Africa and Honduras. Eleven countries were in the low–high pattern category, including Mexico, Cuba, Morocco, Denmark and Luxembourg, etc. The above results are consistent with the global spatial autocorrelation test results, thus indicating that IR, CMR and DCI indicators in some countries/regions may be affected by the COVID-19 epidemic in neighboring countries/regions and may show clear geographical characteristics. To directly reflect the local spatial characteristics of IR, CMR and DCI, LISA scatter diagrams of the three indexes are shown in Figure 3. Most of the three indicators fell into the third quadrant (low–low), but the countries/regions whose IR and DCI index fell into the first quadrant (high– high) and the second quadrant (low–high) had indicator values exceeding the CMR. Thus, among the 178 countries/regions worldwide, the countries/regions with low IR, CMR and DCI indicators showed a spatial agglomeration effect, as did the countries/regions with high IR, CMR and DCI indicators. In addition, some neighboring countries/regions showed some differences in IR, CMR and DCI (high–low and low–high). ![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/05/04/2020.04.23.20077545/F3.medium.gif) [Figure 3.](http://medrxiv.org/content/early/2020/05/04/2020.04.23.20077545/F3) Figure 3. LISA scatter diagram of IR, CMR and DCI indexes for COVID-19 in 178 countries/regions (z is the value of the variable, and Wz is the local Moran’s Ii value of the variable.) ### 3.3 Analysis of the influence of socio-ecological factors on COVID-19 risk #### 3.3.1 Correlation analysis of socio-economic factors To eliminate the influence of the collinearity between the socio-economic indicators on the estimation effect of the model, we established a correlation matrix of the socio-economic indicators (Table S3). The indexes with strong correlation (> 0.8) were screened, and one of the effective indexes was reserved for model analysis. We excluded eight socio-economic indicators in Table 2, numbered 7 (current health expenditure per capita), 12 (total life expectancy at birth), 13 (maternal mortality ratio), 14 (infant mortality rate), 17 (access to basic handwashing facilities including soap and water), 20 (population growth), 21 (proportion of the population spending more than 10% of household consumption or income on out-of-pocket health care expenditure) and 26 (technicians in R&D), and we retained 20 socio-economic indicators. View this table: [Table S3.](http://medrxiv.org/content/early/2020/05/04/2020.04.23.20077545/T7) Table S3. Correlation matrix between socio-economic factors for the risk of COVID-19 #### 3.3.2 Negative binomial regression analysis of socio-ecological factors on COVID-19 risk We analyzed the effects of socio-ecological factors on COVID-19 risk in 178 countries. The results of single-factor and multi-factor negative binomial regression analysis are shown in Table 4. View this table: [Table 3](http://medrxiv.org/content/early/2020/05/04/2020.04.23.20077545/T3) Table 3 Characteristics of UN geographical divisions with reported cases of COVID-19 as of 6 April 2020 View this table: [Table 4.](http://medrxiv.org/content/early/2020/05/04/2020.04.23.20077545/T4) Table 4. The results of single-factor and multi-factor negative binomial regression analysis for COVID-19 in 178 countries/regions The IR was significantly positively associated with GDP per capita (aIRR=1.029, 95%CI: 1.013–1.045), use of at least basic sanitation services (aIRR=1.022, 95%CI: 1.005–1.039) and coverage of social insurance programs (aIRR=1.047, 95%CI: 1.009–1.086), and was significantly negatively associated with the proportion of the population spending more than 25% of household consumption or income on out-of-pocket health care expenses (aIRR=0.846, 95%CI: 0.750–0.955) and the poverty headcount ratio at national poverty lines (aIRR=0.970, 95%CI: 0.948–0.993). The CMR was significantly positively associated with urban populations (aIRR=1.027, 95%CI: 1.010–1.044), GDP per capita (aIRR=1.031, 95%CI: 1.021–1.041) and current health expenditure (aIRR=1.211, 95%CI: 1.040–1.410), and was significantly negatively associated with the number of hospital beds (aIRR=0.799, 95%CI: 0.696–0.916), number of nurses and midwives (aIRR=0.837, 95%CI: 0.749–0.936) and poverty headcount ratio at the national poverty lines (aIRR=0.960, 95%CI: 0.940–0.982). The DCI was significantly positively associated with urban populations (aIRR=1.021, 95%CI: 1.009–1.034), population density (aIRR=1.000, 95%CI: 1.000–1.000) and researchers in R&D (aIRR=1.000, 95%CI: 1.000–1.001), and was significantly negatively associated with the number of hospital beds (aIRR=0.731, 95%CI: 0.641–0.833), number of nurses and midwives (aIRR=0.904, 95%CI: 0.820–0.997) and poverty headcount ratio at the national poverty lines (aIRR=0.963, 95%CI: 0.945–0.982). The results of the sensitivity analysis are reported in Table S4 and Table S5: we used the maximum and minimum temperatures instead of the average temperature, and then incorporated the two climate factors into the single-factor and multi-factor negative binomial regression. The results showed that the significance of different socio-ecological factors was essentially consistent. We found that only the variable of poverty headcount ratio at the national poverty lines (percentage of population) became significant after sensitivity analysis on IR, thus indicating that the analysis results were relatively reliable. View this table: [Table S4.](http://medrxiv.org/content/early/2020/05/04/2020.04.23.20077545/T8) Table S4. The sensitive analysis of single-factor negative binominal regression for the risk of COVID-19 View this table: [Table S5.](http://medrxiv.org/content/early/2020/05/04/2020.04.23.20077545/T9) Table S5. The sensitive analysis of multiple-factor negative binominal regression for the risk of COVID-19 ## 4. Discussion By evaluating the spatial aggregation characteristics of three indicators—IR, CMR and DCI—on COVID-19 risk in 178 countries, Western Europe, Southern Europe, East Asia and some African countries, we found that all showed relatively large spatial correlations, thus indicating that COVID-19 broadly affects these countries/regions. Because COVID-19 is highly contagious, after an outbreak occurs in a country/region, the virus tends to spread rapidly in surrounding countries. Italy was the first country in Europe to have a large outbreak of COVID-19, but the Italian health system adopted several control measures, such as timely intervention and containment measures brought about by decentralization, flexible financing mechanisms, private and public sector partnerships, and human resources mobilization, so that the IR and DCI could be effectively controlled [51]. However, because some European countries did not perform effective prevention and control measures, such as blockading countries or cities, in early stages of the outbreak, the epidemic gradually broke out in countries including France, Germany, Spain and Portugal. In North America, the development of COVID-19 presents progressive characteristics (high–high and low–high mode), and Canada is also significantly affect clearly the United States by the COVID–19 outbreak, but the DCI in Mexico and other countries/regions in Central America remained relatively low. The United States also recently closed its border with Canada and Mexico and decreased the flow of people across the border. Because, blockading and quarantining provide very good protection, taking these measures in countries or cities is very important to decrease the risk of COVID-19 multinational spread. The potential transmission of COVID-19 in South America must not be ignored. Second, both IR and CMR presented low–low patterns in China, thus indicating that China and some neighboring countries (such as South Korea and Singapore) have effectively reduced the risk to neighboring countries by implementing strong prevention and control measures against COVID-19 transmission and have also acquired valuable experience useful to other countries in fighting COVID-19 virus [52–56]. However, notably, the IR of COVID-19 in India presents a low– low model, thus indicating that India is at low risk of COVID-19 transmission from surrounding countries, and consequently has a low DCI. However, as the world’s second most populous country, India may have a high risk of COVID-19 transmission because of inadequate medical conditions and detection levels. Simultaneously, China, South Korea and other countries must strengthen screening of imported cases from other countries, reduce social contact among travelers and prevent the possible secondary transmission of COVID-19 [52,57,58]. Third, most countries/regions in Africa remain in a low–low mode, but the short distance and frequent contacts between Western and Southern Europe and North Africa may place North Africa at high risk of COVID-19 spread; Morocco currently has a low–high mode representing an early warning, and COVID-19 viruses must further be prevented from entering other parts of Africa. Although African countries took measures to prevent the Ebola outbreak in 2014, Africa remains one of the poorest countries worldwide, and it has a shortage of health resources to quickly control the outbreak. Studies have shown that the current spread of COVID-19 in West Africa urgently requires action to control the further spread of COVID-19 and improve the response capacity of affected countries in West Africa [59]. Although most parts of Africa are in the low–low mode, they may also face threats. Many COVID-19 cases may be undetected, thus potentially explaining the current low-level indicators (IR, CMR and DCI). We found that, in terms of urban development, both CMR and DCI were significantly positive associated with the urban population (percentage home of total population), and people per sq. km of land area had significant positive effects on the DCI. Non-drug intervention measures have already been implemented, and if traffic restrictions, social isolation and family measures are not ensured, the increase in population density and urbanization may result in many problems, such as public traffic, rural population health inequities, poor housing conditions, inadequate freshwater supply, and poor sanitation and ventilation systems, thus accelerating the spread of the COVID-19 virus, in agreement with previous research [60,61]. The higher the urban population (percentage of total population), the faster the urbanization process of the country/region; consequently, aging and young people participating in social activities become more likely to aggravate the spread of the virus and increase the burden on the health system, in agreement with the results of one study [62]. In addition, studies have shown that, with urbanization, the risk of infection and the chances of survival after COVID-19 infection among older individuals with complications is greatly increased, thus resulting in a significant increase in the CMR in the country/region [63–64]. In terms of the economy and growth, we found that GDP per capita (current 1,000 US$) was significantly positive associated with the IR and CMR of COVID-19, possibly because the GDP per capita tends to reflect a country’s economic development level: with higher GDP per capita, governments can invest more in screening and treatment of patients with mild and severe cases of COVID-19. Consequently, with more confirmed cases and deaths, classification strategies can be considered for COVID-19 in low-income groups. Especially in economically underdeveloped areas such as Africa, similar symptoms can be used as a basis to implement a series of diagnostic tests [65]. This method of raising clinical diagnostic standards was used in Wuhan, China. In terms of health, we found that increasing the proportion of residents using at least basic sanitation services was significantly positive associated with the IR of COVID-19. For example, improving basic sanitation services and increasing contact between primary health workers and potential and diagnosed COVID-19 patients is very important. In particular, the government of Wuhan, China implemented nucleic acid testing on each resident via hospitals and primary health workers, thus enabling COVID-19 detection in a larger proportion and facilitating rapid control of COVID-19 risk transmission. Second, the numbers of hospital beds (per 1,000 people), nurses and midwives (per 1,000 people) were significantly negative associated with the IR and CMR of COVID-19, thus suggesting that COVID-19 risk should be controlled, and the number of hospital beds and nurses should be increased in a short period of time. The increase in the numbers of hospital beds and nurses can help achieve standardized management of patients and allow more medical resources to be concentrated on the treatment of severe cases. Some research has shown that some countries, such as Italy, China and the United States, have established Fangcang shelter hospitals or field hospitals and increased the numbers of regular hospital beds, intensive care beds and medical workers (by transferring resources from other regions and the military), reopened closed hospitals, and considered use of medical volunteers in the treatment of mild and severe COVID-19 cases; these measures are effective ways to reduce the IR and CMR [66, 67]. Third, among people infected with COVID-19, the proportion of the population spending more than 25% of household consumption or income on out-of-pocket health care expenses were significantly negative associated with the IR, whereas coverage with social insurance plans positively influences the IR. Higher income, enhanced health insurance coverage and decreased burden of medical treatment significantly increase the IR, thus suggesting that governments and health insurance providers should cooperate in the prevention and control of COVID-19. In addition to financial subsidies, the government should also reduce or grant exemptions for patient co-payments, to increase the possibility of COVID-19 patients receiving testing and treatment. In science and technology, the number of researchers in R&D (per million people) was significantly positive associated with the DCI. This improvement includes facilitating health science and technology input, strengthening basic life science research, fostering international cooperation between science and technology (such as in the development and use of effective drugs), providing more convenient testing technology, shortening testing times, expanding the scale of detection, improving treatment technology and performing ongoing vaccine development to reduce present and future COVID-19 transmission. In addition, the poverty headcount ratio at the national poverty lines (percentage of population) was significantly negative associated with the IR, CMR and DCI. Increases in the population in poverty and in racial discrimination greatly diminish accessibility to medical services. Government and society must address these problems through economic stimulus plans, unemployment relief programs, welfare and health safeguarding measures, and plans to decrease health spending by these groups [68]. We also found that climatic factors (temperature, RH, precipitation and wind speed) were not significantly associated with COVID-19 risk, in agreement with the results of some studies [69]. However, previous studies have primarily considered the effects of single climate factors, thus potentially affecting the estimates of the results [70–72]. There is no sufficient evidence indicating that climate factors have specific effects on the spread of COVID-19. This study also shows that in the measurement of COVID-19 risk, the influences of other factors should be considered—such as the constraints of economic development, transportation and other factors—to improve understanding of the mechanisms underlying interrelationships among factors. ## 5. Limitations This study has several limitations. First, we selected cross-sectional data for spatial analysis and regression modeling; therefore, the results may not reflect more changes in time, thus potentially decreasing the statistical ability to detect the relationships among various factors and COVID-19 risk. Second, owing to data matching across databases, some aspects of country/region data may have been lost, thus potentially affecting the spatial weight matrix estimation and regression modeling results. Third, because of the socio-ecological study design, we were unable to access data at the individual level, such as age, sex, occupation, economic and health status, and the actual exposure temperature of each person. However, future studies could adopt hybrid study designs, which use individual-level data from subpopulations to improve ecological extrapolation. ## 6. Conclusion By using the data from 178 countries/regions, we found that socio-economic factors can significantly reduce the risk of COVID-19. As a next step in COVID-19 prevention, different countries/regions should focus on controlling urban populations, providing economic subsidies and medical resource supplies, and taking broad views of social welfare. Strategies may include population isolation, travel restrictions, case screening, cross-regional or national science and technology exchange to promote diagnosis and treatment, public welfare policy improvement, as well as decreasing the burden of low-income groups in obtaining medical treatment. Simultaneously, we must be alert to the COVID-19 risk in some countries in Africa and Asia, and must curb the second wave of COVID-19 transmission. ## Data Availability The data used for the analyses are publicly available from the Johns Hopkins University Center for Systems Science and Engineering (https://github.com/CSSEGISandData/COVID-19), the World Bank (https://datacatalog.worldbank.org/dataset/world-development-indicators) and National Oceanic and Atmospheric Administration, Department of Commerce (https://catalog.data.gov/dataset/global-surface-summary-of-the-day-gsod). [https://github.com/CSSEGISandData/COVID-19](https://github.com/CSSEGISandData/COVID-19) [https://datacatalog.worldbank.org/dataset/world-development-indicators](https://datacatalog.worldbank.org/dataset/world-development-indicators) [https://catalog.data.gov/dataset/global-surface-summary-of-the-day-gsod](https://catalog.data.gov/dataset/global-surface-summary-of-the-day-gsod) ## Contributors DS, XZ contributed to the conception and design of the project; DS, TZ, KH, XZ contributed to the analysis and interpretation of the data; MT, YZ contributed to the data acquisition and provided statistical analysis support; DS drafted the article. DS and XZ are the guarantors. The corresponding author attests that all listed authors meet authorship criteria and that no others meeting the criteria have been omitted. ## Funding This study was funded by National Natural Science Foundation of China (No. 71473096; No. 71673101; No. 71974066) and Michigan Institute for Clinical and Health Research (MICHR No. UL1TR002240). ## Declaration of interests All other authors declare no competing interests. ## Patient and public involvement This research was done without patient involvement. Patients were not invited to comment on the study design and were not consulted to develop patient relevant outcomes or interpret the results. Patients were not invited to contribute to the writing or editing of this document for readability or accuracy. ## Patient consent for publication Not required. ## Ethics approval This study was approved by the Ethics Committee of the Tongji Medical College, Huazhong University of Science and Technology (IORG No: IORG0003571). ## Data sharing statement The data used for the analyses are publicly available from the Johns Hopkins University Center for Systems Science and Engineering ([https://github.com/CSSEGISandData/COVID-19](https://github.com/CSSEGISandData/COVID-19)), the World Bank ([https://datacatalog.worldbank.org/dataset/world-development-indicators](https://datacatalog.worldbank.org/dataset/world-development-indicators)) and National Oceanic and Atmospheric Administration, Department of Commerce ([https://catalog.data.gov/dataset/global-surface-summary-of-the-day-gsod](https://catalog.data.gov/dataset/global-surface-summary-of-the-day-gsod)). ## Acknowledgments The authors would like to thank the National Natural Science Foundation of China and the National School of Development, Peking University, University of Michigan, and other members for their support and cooperation. We would also like to thank the study samples from The Center for Systems Science and Engineering, Johns Hopkins University (JHU-CSSE), Global Surface Summary of the Day (GSOD) from the Integrated Surface Hourly (ISH) dataset and World Development Indicators (WDI) dataset for providing the information in our research. * Received April 23, 2020. * Revision received April 30, 2020. * Accepted May 4, 2020. * © 2020, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. [1].Cleaveland, S., M. K. Laurenson, and L. H. Taylor “Diseases of humans and their domestic mammals: pathogen characteristics, host range and the risk of emergence.” Philosophical Transactions of the Royal Society of London. Series B: Biological Sciences 356.1411 (2001): 991–999. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rstb.2001.0889&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11516377&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F04%2F2020.04.23.20077545.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000170315900003&link_type=ISI) 2. [2].Halliday, Jo EB, et al. “Driving improvements in emerging disease surveillance through locally relevant capacity strengthening.” Science 357.6347 (2017): 146–148. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNTcvNjM0Ny8xNDYiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wNS8wNC8yMDIwLjA0LjIzLjIwMDc3NTQ1LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 3. [3].Schneider, Maria Cristina, and Gustavo Machado. “Environmental and socioeconomic drivers in infectious disease.” The Lancet Planetary Health 2.5 (2018): e198-e199. 4. [4].Daszak, Peter, Andrew A. Cunningham, and Alex D. Hyatt. “Emerging infectious diseases of wildlife--threats to biodiversity and human health.” science 287.5452 (2000): 443–449. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIyODcvNTQ1Mi80NDMiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMC8wNS8wNC8yMDIwLjA0LjIzLjIwMDc3NTQ1LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 5. [5].Weiss, Robin A., and Anthony J. McMichael. “Social and environmental risk factors in the emergence of infectious diseases.” Nature medicine 10.12 (2004): S70-S76. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nm1150&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15577934&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F04%2F2020.04.23.20077545.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000225733900002&link_type=ISI) 6. [6].Woolhouse, Mark EJ, and Sonya Gowtage-Sequeria. “Host range and emerging and reemerging pathogens.” Emerging infectious diseases 11.12 (2005). 7. [7].Patz, Jonathan A., et al. “Unhealthy landscapes: policy recommendations on land use change and infectious disease emergence.” Environmental health perspectives 112.10 (2004): 1092–1098. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1289/ehp.6877&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15238283&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F04%2F2020.04.23.20077545.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000222482300030&link_type=ISI) 8. [8].Kuhn, K., et al. “Using climate to predict infectious disease outbreaks: A review.” World Health Organization (2004). 9. [9].Wu, XiaoXu, et al. “Impact of global change on transmission of human infectious diseases.” Science China Earth Sciences 57.2 (2014): 189–203. 10. [10].Hamrick, Patricia Najera, et al. “Geographic patterns and environmental factors associated with human yellow fever presence in the Americas.” PLoS neglected tropical diseases 11.9 (2017): e0005897. 11. [11].Muñoz-Zanzi, Claudia, et al. “Leptospira contamination in household and environmental water in rural communities in southern Chile.” International journal of environmental research and public health 11.7 (2014): 6666–6680. 12. [12].Semenza, Jan C., et al. “Linking environmental drivers to infectious diseases: the European environment and epidemiology network.” PLoS neglected tropical diseases 7.7 (2013). 13. [13].Froment, A., et al. “Biodiversity and health: the place of parasitic and infectious diseases. Biodiversity change and human health.” (2009): 211–227. 14. [14].Yang, Hyun M., and Marcelo U. Ferreira. “Assessing the effects of global warming and local social and economic conditions on the malaria transmission.” Revista de saude publica 34.3 (2000): 214–222. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10920442&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F04%2F2020.04.23.20077545.atom) 15. [15].Martens, Pim, and Susanne C. Moser. “Health impacts of climate change.” Science 292.5519 (2001): 1065–1066. [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIyOTIvNTUxOS8xMDY1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMDUvMDQvMjAyMC4wNC4yMy4yMDA3NzU0NS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 16. [16].McMichael, Anthony J., and Rosalie E. Woodruff. “Detecting the health effects of environmental change: scientific and political challenge.” (2005): 1–3. 17. [17].Casman, Elizabeth A., and Hadi Dowlatabadi, eds. The contextual determinants of malaria. Resources for the Future, 2002. 18. [18].Tol, Richard SJ, and Hadi Dowlatabadi. “Vector-borne diseases, development & climate change.” Integrated Assessment 2.4 (2001): 173–181. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1023/A:1013390516078&link_type=DOI) 19. [19].Epstein, Paul R. “Climate change and emerging infectious diseases.” Microbes and infection 3.9 (2001): 747–754. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S1286-4579(01)01429-0&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11489423&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F04%2F2020.04.23.20077545.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000170529200008&link_type=ISI) 20. [20].Wu, Xiaoxu, et al. “Impact of climate change on human infectious diseases: Empirical evidence and human adaptation.” Environment international 86 (2016): 14–23. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.envint.2015.09.007&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26479830&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F04%2F2020.04.23.20077545.atom) 21. [21].Wu, XiaoXu, et al. “Impact of global change on transmission of human infectious diseases.” Science China Earth Sciences 57.2 (2014): 189–203. 22. [22].Parham, Paul E., and Edwin Michael. “Outbreak properties of epidemic models: The roles of temporal forcing and stochasticity on pathogen invasion dynamics.” Journal of theoretical biology 271.1 (2011): 1–9. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jtbi.2010.11.015&link_type=DOI) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000286492100001&link_type=ISI) 23. [23].Turell, Michael J., Lee W. Cohnstaedt, and William C. Wilson. “Effect of Environmental Temperature on the Ability of Culex tarsalis and Aedes taeniorhynchus (Diptera: Culicidae) to Transmit Rift Valley Fever Virus.” Vector-Borne and Zoonotic Diseases (2020). 24. [24].Kioutsioukis, Ioannis, and Nikolaos I. Stilianakis. “Assessment of West Nile virus transmission risk from a weather-dependent epidemiological model and a global sensitivity analysis framework.” Acta tropica 193 (2019): 129–141. 25. [25].Peci, Adriana, et al. “Effects of absolute humidity, relative humidity, temperature, and wind speed on influenza activity in Toronto, Ontario, Canada.” Appl. Environ. Microbiol. 85.6 (2019): e02426–18. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F04%2F2020.04.23.20077545.atom) 26. [26].Sung, Minki, et al. “Airflow as a Possible Transmission Route of Middle East Respiratory Syndrome at an Initial Outbreak Hospital in Korea.” International journal of environmental research and public health 15.12 (2018). 27. [27].Huang, Chaolin, et al. “Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China.” The Lancet 395.10223 (2020): 497–506. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F04%2F2020.04.23.20077545.atom) 28. [28].Chen, Nanshan, et al. “Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: a descriptive study.” The Lancet 395.10223 (2020): 507–513. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/s01406736(20)30211-7&link_type=DOI) 29. [29].Wang, Dawei, et al. “Clinical characteristics of 138 hospitalized patients with 2019 novel coronavirus–infected pneumonia in Wuhan, China.” Jama (2020). 30. [30].Guan, Wei-jie, et al. “Clinical characteristics of coronavirus disease 2019 in China.” New England Journal of Medicine (2020). 31. [31].Liu, Ying, et al. “The reproductive number of COVID-19 is higher compared to SARS coronavirus.” Journal of travel medicine (2020). 32. [32].Parham, Paul E., et al. “Understanding and modelling the impact of climate change on infectious diseases–progress and future challenges.” Climate Change—Socioeconomic Effects (2011). 33. [33].Caldwell, John C. “Rethinking the African AIDS epidemic.” Population and development review 26.1 (2000): 117–135. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1728-4457.2000.00117.x&link_type=DOI) 34. [34].Butler, Colin. “HIV and AIDS, poverty, and causation.” The Lancet 356.9239 (2000): 1445–1446. 35. [35].Auvert, Betran, et al. “Ecological and individual level analysis of risk factors for HIV infection in four urban populations in sub-Saharan Africa with different levels of HIV infection.” Aids 15 (2001): S15-S30. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/00002030-200108004-00003&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11686462&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F04%2F2020.04.23.20077545.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000171519300003&link_type=ISI) 36. [36].Weiss, Robin A., and Anthony J. McMichael. “Social and environmental risk factors in the emergence of infectious diseases.” Nature medicine 10.12 (2004): S70-S76. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nm1150&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15577934&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F04%2F2020.04.23.20077545.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000225733900002&link_type=ISI) 37. [37].Houéto, D. “The social determinants of emerging infectious diseases in Africa.” MOJ Public Health 8.2 (2019): 57–63. 38. [38].Wu, Xiaoxu, et al. “Impact of climate change on human infectious diseases: Empirical evidence and human adaptation.” Environment international 86 (2016): 14–23. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.envint.2015.09.007&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26479830&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F04%2F2020.04.23.20077545.atom) 39. [39].Kovats, Sari, et al. Climate change and human health: impact and adaptation. World Health Organization (WHO), 2000. 40. [40].Houéto, D. “The social determinants of emerging infectious diseases in Africa.” MOJ Public Health 8.2 (2019): 57–63. 41. [41].Parham, Paul E., et al. “Climate, environmental and socio-economic change: weighing up the balance in vector-borne disease transmission.” Philosophical Transactions of the Royal Society B: Biological Sciences 370.1665 (2015). 42. [42].Zhang S, Diao M Y, Yu W, et al. Estimation of the reproductive number of Novel Coronavirus (COVID-19) and the probable outbreak size on the Diamond Princess cruise ship: A data-driven analysis[J]. International Journal of Infectious Diseases, 2020, 93: 201–204. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijid.2020.02.033&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32097725&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F05%2F04%2F2020.04.23.20077545.atom) 43. [43].Li, Qun, et al. “Early transmission dynamics in Wuhan, China, of novel coronavirus–infected pneumonia.” New England Journal of Medicine (2020). 44. [44].Xie X, Zhong Z, Zhao W, et al. Chest CT for typical 2019-nCoV pneumonia: relationship to negative RT-PCR testing[J]. Radiology, 2020: 200343. 45. [45].Dong, Ensheng, Hongru Du, and Lauren Gardner. “An interactive web-based dashboard to track COVID-19 in real time.” The Lancet Infectious Diseases (2020). 46. [46].Chen, Biqing, et al. “Roles of meteorological conditions in COVID-19 transmission on a worldwide scale.” medRxiv (2020). 47. [47].Anselin, Luc. “A local indicator of multivariate spatial association: extending Geary’s C.” Geographical Analysis 51.2 (2019): 133–150. 48. [48].Elhorst, J. Paul, Donald J. Lacombe, and Gianfranco Piras. “On model specification and parameter space definitions in higher order spatial econometric models.” Regional Science and Urban Economics 42.1–2 (2012): 211–220. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.regsciurbeco.2011.09.003&link_type=DOI) 49. [49].Phung, Dung, et al. “The effects of socioecological factors on variation of communicable diseases: A multiple-disease study at the national scale of Vietnam.” PloS one 13.3 (2018). 50. [50].Hondula, David M., and Adrian G. Barnett. “Heat-related morbidity in Brisbane, Australia: spatial variation and area-level predictors.” Environmental health perspectives 122.8 (2014): 831–836. 51. [51].Armocida, Benedetta, et al. “The Italian health system and the COVID-19 challenge.” The Lancet Public Health (2020). 52. [52].Leung, Kathy, et al. “First-wave COVID-19 transmissibility and severity in China outside Hubei after control measures, and second-wave scenario planning: a modelling impact assessment.” The Lancet (2020). 53. [53].Kraemer, Moritz UG, et al. “The effect of human mobility and control measures on the COVID-19 epidemic in China.” Science (2020). 54. [54].Chinazzi, Matteo, et al. “The effect of travel restrictions on the spread of the 2019 novel coronavirus (COVID-19) outbreak.” Science (2020). 55. [55].Prem, Kiesha, et al. “The effect of control strategies to reduce social mixing on outcomes of the COVID-19 epidemic in Wuhan, China: a modelling study.” The Lancet Public Health (2020). 56. [56].Koo, Joel R., et al. “Interventions to mitigate early spread of SARS-CoV-2 in Singapore: a modelling study.” The Lancet Infectious Diseases (2020). 57. [57].Xu, Shunqing, and Yuanyuan Li. “Beware of the second wave of COVID-19.” The Lancet (2020). 58. [58].Niehus, Rene, et al. “Using observational data to quantify bias of traveller-derived COVID-19 prevalence estimates in Wuhan, China.” The Lancet Infectious Diseases (2020). 59. [59].Martinez-Alvarez, Melisa, et al. “COVID-19 pandemic in west Africa.” The Lancet Global Health (2020). 60. [60].Lee, Vernon J., et al. “Epidemic preparedness in urban settings: new challenges and opportunities.” The Lancet Infectious Diseases (2020). 61. [61].Dowd, Jennifer Beam, et al. “Demographic science aids in understanding the spread and fatality rates of COVID-19.” medRxiv (2020). 62. [62].Sajadi, Mohammad M. and Habibzadeh, Parham and Vintzileos, Augustin and Shokouhi, Shervin and Miralles-Wilhelm, Fernando and Amoroso, Anthony, Temperature, Humidity and Latitude Analysis to Predict Potential Spread and Seasonality for COVID-19 (March 5, 2020). Available at SSRN: [https://ssrn.com/abstract=3550308](https://ssrn.com/abstract=3550308) or [http://dx.doi.org/10.2139/ssrn.3550308](http://dx.doi.org/10.2139/ssrn.3550308) 63. [63].Ruan, Shigui. “Likelihood of survival of coronavirus disease 2019.” The Lancet Infectious Diseases (2020). 64. [64].Fang, Lei, George Karakiulakis, and Michael Roth. “Are patients with hypertension and diabetes mellitus at increased risk for COVID-19 infection?.” The Lancet. Respiratory Medicine (2020). 65. [65].Ayebare, Rodgers R., et al. “Adoption of COVID-19 triage strategies for low-income settings.” The Lancet Respiratory Medicine (2020). 66. [66].Chen, Simiao, et al. “Fangcang shelter hospitals: a novel concept for responding to public health emergencies.” The Lancet (2020). 67. [67].Miani, Alessandro, et al. “The Italian war-like measures to fight coronavirus spreading: Reopen closed hospitals now.” EClinicalMedicine (2020). 68. [68].McCloskey, Brian, et al. “Mass gathering events and reducing further global spread of COVID-19: a political and public health dilemma.” Lancet (London, England) 395.10230 (2020). 69. [69].Yao, Y., Pan, J., Liu, Z., Meng, X., Wang, W., Kan, H., & Wang, W. (2020). No Association of COVID-19 transmission with temperature or UV radiation in Chinese cities. The European respiratory journal, 2000517. Advance online publication. [https://doi.org/10.1183/13993003.00517-2020](https://doi.org/10.1183/13993003.00517-2020) 70. [70].Islam, Nazrul, Sharmin Shabnam, and A. Mesut Erzurumluoglu. “Meteorological factors and Covid-19 incidence in 310 regions across the world.” medRxiv (2020). 71. [71].Zhu, Yongjian, and Jingui Xie. “Association between ambient temperature and COVID-19 infection in 122 cities from China.” Science of The Total Environment (2020): 138201. 72. [72].Wang, Mao, et al. “Temperature significant change COVID-19 Transmission in 429 cities.” medRxiv (2020). [1]: /embed/graphic-3.gif [2]: /embed/graphic-5.gif [3]: /embed/graphic-6.gif [4]: /embed/inline-graphic-1.gif [5]: /embed/graphic-7.gif [6]: /embed/graphic-8.gif