key: cord-1024905-7zrzl77x authors: Lee, Kyung-Shin; Min, Hye Sook; Jeon, Jae-Hyun; Choi, Yoon-Jung; Bang, Ji Hwan; Sung, Ho Kyung title: The association between greenness exposure and COVID-19 incidence in South Korea: An ecological study date: 2022-04-02 journal: Sci Total Environ DOI: 10.1016/j.scitotenv.2022.154981 sha: 7f75e520bef62618f035719f8029464e76302caf doc_id: 1024905 cord_uid: 7zrzl77x Background The rapid spread of COVID-19 has caused an emergency situation worldwide. Investigating the association between environmental characteristics and COVID-19 incidence can be of the occurrence and transmission. The objective of this study was to evaluate the association between greenness exposure and COVID-19 cases at the district levels in South Korea. We also explored this association by considering several environmental indicators. Methods District-level data from across South Korea were used to model the cumulative count of COVID-19 cases per 100,000 persons between January 20, 2020, and February 25, 2021. Greenness exposure data were derived from the Environmental Geographic Information Service of the Korean Ministry of Environment. A negative binomial mixed model evaluated the association between greenness exposure and COVID-19 incidence rate at the district level. Furthermore, we assessed this association between demographic, socioeconomic, environmental statuses, and COVID-19 incidence. Results Data from 239 of 250 districts (95.6%) were included in the analyses, resulting in 127.89 COVID-19 cases per 100,000 persons between January 20, 2020 and February 25, 2021. Several demographic and socioeconomic variables, districts with a higher rate of natural greenness exposure, were significantly associated with lower COVID-19 incidence rates (incidence rate ratio (IRR), 0.70; 95% confidence interval (CI), 0.54–0.90; P-value = 0.008) after adjusting covariates, but no evidence for the association between built greenness and COVID-19 incidence rates was found. Conclusion In this ecological study of South Korea, we found that higher rates of exposure to natural greenness were associated with lower rates of COVID-19 cases. The rapid spread of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has caused an emergency situation worldwide (Sohrabi et al., 2020) . SARS-CoV-2 is a highly pathogenic and transmissible virus which causes coronavirus disease 2019 (COVID-19) (Russette et al., 2021) . Since vaccines were not available until the end of 2020 and due to the scarcity of effective treatments, population-level public health measures were put in place to control the outbreak, including social distancing, stay-at-home orders, nationwide lockdowns, and public education initiatives (Alfano and Ercolano, 2020; Tammes, 2020) . In South Korea, the first patient was diagnosed on January 20, 2020. On May 2, a total of 10,799 people were diagnosed with COVID-19. The fatality rate was 2.4%, which translates to 4.94 fatalities per million (Park et al., 2020c) . During the COVID-19 pandemic in South Korea, public health authorities conducted rapid and strict contact tracing, isolation of confirmed cases, and quarantined the suspected COVID-19 cases (Park et al., 2020a) . The government implemented strict nationwide social distancing guidelines, however, school shutdowns and strict closures were not implemented (Shim et al., 2021) . South Korea experienced three waves of the COVID-19 until December 2020. The first wave from February till March was related to religious group gatherings in the local area. The daily new cases in the second wave, comprising 441 cases on August 26, was related to worship services in the densely populated Seoul areas. The third wave in December 2020, however, had more than 1000 daily new cases. The third wave also had a fatality rate higher than that in the second wave (1.26% vs. 0.91%) (Seong et al., 2021) . Regarding non-pharmaceutical interventions (NPI) to contain the spread of COVID-19, South Korea is one of the few countries that had not implemented restrictions such as large-scale lockdowns or stay-at-home orders (Dighe et al., 2020; Jeong et al., 2020; Yoo et al., 2021) . Instead of stringent NPI, the South Korean government focused on an aggressive -test, trace, isolate‖ strategy (Dighe et al., 2020) . Although social etiquette, as well as avoiding crowds, not using public transportation, and maintaining at least 6 feet distance from others. A recent epidemiological study showed that outdoor physical activities increased during the pandemic compared to time before the implementation of social distancing policy (Park et al., 2020b) . The study also suggested if physical distancing measures are maintained carefully during the pandemic, outdoor physical activity can be safe for people and the spread of infection can be stabilized (Park et al., 2020b) . Environmental characteristics are indivisibly linked with the occurrence and spread of diseases (Jian et al., 2014) . Previous studies have shown that several environmental factors, including air pollution (Wu et al., 2020a) , greenness exposure (Klompmaker et al., 2021; Russette et al., 2021) , temperature , and humidity (Wang et al., 2020) , were associated with the spread, transmission, and reduction of COVID-19 due to the stability of the virus as well as host susceptibility. Spotswood et al. suggested that nature plays an important role in building defenses against viruses through boosting natural killer cells. Their result showed that an increase of 0.1 in the Normalized Difference Vegetation Index was associated with a 4.1% decrease in COVID-19 incidence rate (Spotswood et al., 2021) . A recent emerging and growing body of data has led many to suggest that greenness exposure could benefit health outcomes of various diseases including infectious diseases (Cunha et al., 2021; Eldeirawi et al., 2019; Ferrante et al., 2020; Huang et al., 2018; Squillacioti et al., 2020; Sun et al., 2020) . Links between greenness and health elaborating the potential pathways was noted which included harm reduction and restoring and building capacities (Markevych et al., 2017) . The mediating factors of this association were also reported to include reduced stress (McCormick, 2017) , increased physical activity (Marquet et al., 2020) , filtered out air pollution (Nowak et al., 2006) , decreased traffic-related noise (Dzhambov et al., 2019) , and enhanced social interaction (Liu et al., 2019) . Several studies have provided the importance of public green spaces because exposure to more green during COVID-19 seemed to reduce air pollutants and increased time in nature has been linked to positive behavioral changes (Diener and Mudu, 2021; Labib et al., 2021) . Scientific evidence supports the beneficial role of environmental factors such as greenness and biodiversity in regulating the human immune system (Copat et al., 2020; Marselle et al., 2021; Rook, 2013) , however, only a few studies have considered greenness exposure in relation to COVID-19 incidence rate. Four US studies suggested that counties with higher exposure to greenness had better COVID-19 outcome measures, including incidence and mortality rates (Klompmaker et al., 2021; Lee et al., 2021c; Russette et al., 2021; Spotswood et al., 2021) . However, the impact Moreover, it is still unclear whether greenness exposure and the incidence of COVID-19 cases are related and to what extent sociodemographic factors, such as age and socioeconomic status, play a role in explaining the number of COVID-19 cases in the previous studies. Therefore, this study aimed to evaluate the association between greenness exposure and COVID-19 cases across 239 districts in South Korea using an ecological study design. We also explored this association by considering several environmental indicators. South Korea is composed of 250 districts, the basic administrative units of Korea, in eight metropolitan cities and nine provinces. We excluded from analysis eight districts near the border with North Korea for which greenness exposure data were not available due to security reasons, one district in Daegu city as its outbreak pattern differed from that of other places because of its association with a religious institute, and two districts due to the lack of exact information on COVID-19 cases. Finally, we included 239 districts in this study. We obtained COVID-19 counts for the included South Korean districts from the Central Disaster Management Headquarters (http://ncov.mohw.go.kr/en/). This source provides the most comprehensive districtlevel COVID-19 data to date and is updated daily by the Korean government. We collected the cumulative number of confirmed cases for each district up to February 25, 2021, the period before vaccinations started in South Korea. District-level COVID-19 incidence rates were defined for this analysis as the number of COVID-image data with a shapefile, then we combined it to the administrative district map from the open data (http://www.gisdeveloper.co.kr/) using a geometric intersection tool of the ArcGIS desktop version 10.5 software (ESRI Inc., Redlands, CA, USA). We calculated the ratio of the quotient of green area divided by total area (square meter divided by square meter). We separately analyzed two types of greenness: natural greenness, which included forests and natural grasslands, and built greenness, which included golf courses, cemeteries, artificial grassland, and street trees. The definition of each type of greenness was presented in table S1(Lee et al., 2021b). Based on previous studies and considering the goodness of fit test in the model, we selected ten potential confounding factors for each district that included demographic, socioeconomic, environmental, comorbid or behavioral variables (Klompmaker et al., 2021; Lee et al., 2021c; Russette et al., 2021) . We obtained from Statistics Korea (http://kostat.go.kr/portal/eng/index.action) information on the percentage of older people (aged 65 years or more) in the population, and this is a demographic variable. We also obtained data on four socioeconomic variables, house prices (median apartment price), percentage of people educated to the high school level or above in the adult population (aged >20 years), deprivation index (range, 0-500), and the proportion of apartments residents. The data on house prices were collected from the Ministry of Land, Infrastructure and Transport and those on the education level from Statistics Korea. The deprivation index was calculated as the sum of the scores of five domains, including unemployment, poverty, housing, labor, and social networks. Higher deprivation index score indicated greater deprivation (Shin et al., 2009 ). Since COVID-19 spreads quickly through close contact in an enclosed space, we included the proportion of apartment residents per administrative unit as a covariate. The district-level environmental variables included NO 2 concentration and temperature based on previous studies about the relationship between exposure to air pollution, temperature, district-level information on the number of urbanity, and the number of confirmed cases less than 10 (Klompmaker et al., 2021; Raines et al., 2021; Wu et al., 2020b) . The NO 2 concentration was computed by annual average concentration in 2019 from the mean hourly data from the Air Korea website (http://www.airkorea.or.kr/eng/). The hourly NO 2 for the 364 Air Korea monitoring stations were measured using the chemiluminescent methods, which provide the data J o u r n a l P r e -p r o o f Journal Pre-proof regarding the outdoor air quality and is managed by the Korean Ministry of Environment(Lee et al., 2021a). Because the monitoring stations were not covered for all districts, we estimated that through kriging interpolation for each district (Lee et al., 2020a; Oliver and Webster, 1990) . We also computed the average summer (June-August) and winter (December-February) temperatures from the hourly ambient temperature data in 2019, which was obtained from the Meteorological Administration. To show the different transmission dynamics between urban and rural areas, as well as a proxy for the stage of the COVID-19 outbreak, we included a dummy variable for urban areas (i.e., eight metropolitan cities) or rural areas (i.e., nine provinces), and a dummy variable for a district with confirmed cases less than 10 (Yes) or 10 cases and more (No). District-level comorbid or behavioral variables in the adult population included the percentages of people with diabetes mellitus, hypertension, and moderate physical activity. We collected these district-level variables from the Community Health Survey (Kang et al., 2015) . Detailed information on how these data were collected is described in table S2. We tested a variance inflation factor(VIF) as a diagnostic tool for multi-collinearity to assess whether the covariates were inflated when the VIF was large (VIF >10) (Table S3 ) (Tu et al., 2005) . We also displayed a cross-correlation to estimate how the variables correlated with each other ( Figure S1 ). The ArcGIS desktop version 10.5. software (ESRI Inc., Redlands, CA, USA) was used to visualize the spatial distribution of greenness exposure and COVID-19 incidence rate ratio (IRR) across South Korea. We examined the descriptive statistics of the 239 districts. We estimated the significance of the association between greenness and covariates, and the association between covariates and COVID-19 IRR to select the modifier or pure confounder (Blakely and Woodward, 2000) . We also examined the association between greenness and covariates using Pearson's correlation analysis for continuous covariates, and Student's t-test for categorical covariates (Table S4) . We estimated the association between the covariates and COVID-19 IRR in univariate models on COVID-19 IRR. Subsequently, we used a negative binomial generalized linear mixed model, with districts as random intercepts and a population size(log transformed) as offset variable, to account for the overdispersion of the data, we performed model fitting to estimate the linearity for the association of greenness with COVID-19 IRR using the glmmTMB function in R package (Magnusson et al., 2017) . We also plotted the relationship of greenness with each covariate to test this linearity ( Figure S2 ). For the main model, the associations of J o u r n a l P r e -p r o o f Journal Pre-proof greenness exposure per interquartile range (IQR) increase were expressed in terms of the increase in COVID-19 incidence rate using the crude, multivariate models adjusting for pure confounders and multivariate models adjusting for the overall covariates. These analyses were conducted using a negative binomial mixed model, with population size offset and random intercept by districts, to account for potential correlation between districts within the same local area for estimating the association between exposure to greenness and COVID-19 IRR; the 95% confidence interval (CI) per IQR greenness increase was determined using PROC GLIMMIX in SAS (version 9.4; SAS Institute, Cary, NC, USA). Compound symmetry was assumed for the covariance structure, considering the spatial proximity of residents in each district in previous studies (Lee et al., 2009; Lee et al., 2019a; Lee et al., 2020a; Lee et al., 2019b) . In sensitivity analyses, we conducted the association between greenness and COVID-19 IRR by stratifying into quantiles of the population density, and deprivation index to estimate whether there were different effects depending on these variables as a moderator. The first case of COVID-19 in South Korea was reported on January 20, 2020. As of February 25, 2021, the day before vaccination against COVID-19 started in the country, 88,922 confirmed cases of COVID-19 were reported in South Korea. Table 1 shows the mean COVID-19 incidence rates in eight cities and nine provinces in South Korea until February 25, 2021. Seoul City, Daegu City, and Jeju Province had the highest mean incidence rates, reaching 264, 303, and 297 cases per 100,000 persons, respectively. Our study included 239 out of 250 districts of South Korea. The highest percentages of greenness were seen in districts along the east coast and the southern areas (Figure 1) . At the district level, the top three highest COVID-19 incidence rates were seen in Cheorwon-gun, Sunchang-gun, and Cheongdo-gun (Figure 2 ). On the other hand, the highest percentages of built greenness were seen in districts with metropolitan cities. The mean percentage ± standard deviation (SD) of adults older than 65 years and adults educated to the high school level or above were 19.89 ± 8.23% and 66.55 ± 14.39%, respectively. The mean natural, and built greenness percentages in all 239 districts were 46.99%, and 2.21%, respectively. The mean COVID-19 incidence was 127.9 cases per 100,000 persons, calculated for the period from January 20, 2020, to February 25, 2021 (Table J o The explored associations between various covariates and COVID-19 IRR are presented in Table 3 . In the univariate models, districts with higher house prices (IRR, 1.69; 95% CI, 1.46-1.96), a higher proportion of apartments (IRR, 1.00; 95% CI, 1.00-1.01), districts with urbanity (IRR, 1.65; 95% CI, 1.32-2.06 compared to the rural area), a higher percentage of adults educated to the high school level or above (IRR, 1.02; 95% CI, 1.01-1.02), higher deprivation index (IRR, 1.32; 95% CI, 1.09-1.60), and higher annual NO 2 concentrations (IRR, 1.04; 95% CI, 1.03-1.05), a higher average temperature for the summer season (IRR, 1.16; 9%% CI, 1.03-1.31) had a significantly higher risk of high COVID-19 IRR, while districts with a higher prevalence of people aged 65 years (IRR, 0.96; 95% CI, 0.95-0.98), with a lower average temperature for the winter season (IRR, 0.91; 95% CI, 0.86-0.95) had a significantly lower risk. The overall exposure-response curve for COVID-19 incidence showed a linear negative relationship with natural greenness but a positive relationship with built greenness, after adjusting for covariates using a negative binomial mixed model (Figure 3 ). In the multivariate model adjusting for pure confounders, we found no association between natural or built greenness, and COVID-19 IRR. However, in the full models, we found IRRs of 0.70 (95% CI, 0.54-0.90) when one-IQR natural greenness increased. These associations between built greenness and COVID-19 still remained no significant even after further adjustment for overall covariates in the full model (IRR, 1.00; 95% CI, 0.88-1.14) ( Table 4 ). In subgroup analysis, we found that there are no statistically significant associations between greenness and COVID-19 incidence by stratifying into the quantiles of population density, and deprivation index, and we also did not find any significant interaction effects by these strata in Table S5 . Journal Pre-proof Figure 3 . The association between greenness and COVID-19 incidence rate among 239 districts. Abbreviations: COVID-19, coronavirus 2019. Models included a population size(log transformed) offset and a random intercept by district and adjusted for percentage of people aged ≥65 years, house price (log transformed), deprivation index (log transformed), annual NO 2 concentration, average summer temperature, average winter temperature, the proportion of apartments residents, Districts with urbanity, Districts with confirmed cases less than 10, and percentages of adults with high school education status, adults with diabetes mellitus, adults with hypertension, and adults with moderate physical activity using negative binomial mixed model. We found that greenness was significantly associated with the COVID-19 IRR. Interestingly, an IQR increase in natural greenness was associated with a COVID-19 IRR decrease after adjusting for potential district-level confounders. We also found that South Korean districts with higher house prices, higher annual NO 2 concentrations, higher deprivation index, higher proportion of apartments, greater urbanization, higher average temperature during the summer season, lower average temperature during winter season, higher education status, lower population of the elderly, were at a higher risk of having a high COVID-19 IRR. Understanding the relationship that environmental and socioeconomic factors have with infectious diseases could lead to better disease prevention strategies (Huang et al., 2018) ; however, few studies previously have reported on the impact of greenness on infectious diseases including COVID-19. Huang et al. reported the association between greenness and a dengue fever (DF) outbreak in Taiwan using a spatial autocorrelation method and generalized linear mixed models (Huang et al., 2018) . They suggested a high aggregation of DF hotspots in the metropolitan areas of Taiwan, a positive correlation between parks and DF, and a negative association between other land types and DF. However, a study in South America reported that there was no significant association between greenness and DF. This study suggested that socioeconomic vulnerability was a strong modifiable confounder for this association (Cunha et al., 2021) . Recently, four US ecological studies have shown an association between greenness exposure and low COVID-19 mortality and incidence rates (Klompmaker et al., 2021; Lee et al., 2021c; Russette et al., 2021; Spotswood et al., 2021) . Russette et al. found that the top three greenspace deciles were significantly associated with a reduced risk of COVID-19 mortality (Russette et al., 2021) . They suggested that counties with a higher prevalence of health coverage for very low-income earners(Medicaid) had a significantly lower risk of COVID-19 mortality. Klompmaker et al. reported stronger associations between greenness and COVID-19 incidence than mortality (Klompmaker et al., 2021) , possibly because greenness affected contact rates during the COVID-19 pandemic. COVID-19 mortality was related to treatments and host susceptibility due to factors such as age and the presence of chronic diseases (Klompmaker et al., 2021) . Lee et al. showed that population density was associated with worse COVID-19 outcomes, including transmission and mortality (Lee et al., 2021c); thus, areas with higher greenness had lower COVID-19 incidence, transmission, and mortality rates. They also suggested that higher greenness was strongly associated with lower population density and that this provided many opportunities for outdoor activities, J o u r n a l P r e -p r o o f Journal Pre-proof resulting in a lower transmission risk than being involved in indoor activities (Lee et al., 2021c) . Several studies reported that people spent significantly more time in greenness and visited nature more often during the pandemic than before (Robinson et al., 2021; Soga et al., 2021; Venter et al., 2020) . The study by Robinson et al. estimated the individual mental well-being index using a digital questionnaire and well-being instruments. They suggested that increased greenness surrounding a participant's residence was associated with better mental well-being during the COVID-19 pandemic (Robinson et al., 2021) . In line with this study, several previous studies showed that greenness was beneficial to mental and physical health outcomes, such as mental well-being scale and self-reported general health, during the COVID-19 pandemic (Slater et al., 2020; Ugolini et al., 2020) . Potential mechanisms for this association could be that neighborhoods with higher greenness could affect the host susceptibility (Klompmaker et al., 2021) , lower the rates of person-to-person contact (Klompmaker et al., 2021) , improved the air quality (Diener and Mudu, 2021; Ji et al., 2020) , and present a wide diversity of microorganisms that could regulate the human immune system (Roslund et al., 2020) . However, further studies should investigate the role of biological mechanisms in the association between greenness exposure and infectious diseases. We found that the districts with higher natural greenness were associated with a lower COVID-19 incidence rate in South Korea, whereas, higher built greenness had no association with COVID-19 incidence rate. The result showed that districts with higher urbanization such as higher population, a high proportion of apartment residents, and high percentages of people aged at least 65 years had a lower rate of natural greenness. However, after adjusting for these covariates, there was a significant association noted between natural greenness and COVID-19 IRR. It suggested the probability that natural and built greenness has different effects on health outcomes. However, only a few previous studies have considered the differences between types of greenness. Thus, more research using the type of greenness as an indicator should be conducted in the future. Our univariate model results suggested that districts with higher house prices, higher rates of people with high school education or above, higher annual NO 2 concentrations, and lower rates of older people were associated with an increased risk of COVID-19 incidence. Consistent with previous studies, people living in these environments usually have a more active social life and, therefore, have higher COVID-19 transmission and incidence rates (Boehmer et al., 2020; Lee et al., 2020b) . Our finding of districts with high deprivation index was also consistent with the results of previous studies. These people may have been at a higher risk due to poor 2 Adjusted Model1 included a population size offset and a random intercept by district and adjusted for pure confounders such as percentage of people aged ≥65 years, house price (log transformed), deprivation index (log transformed), annual NO 2 concentration, average summer temperature, the proportion of apartments residents, Districts with urbanity, Districts with confirmed cases less than 10, and percentages of adults with high school education status, using negative binomial mixed model. The efficacy of lockdown against COVID-19: a cross-country panel analysis. Appl. Health Econ. Health Policy Ecological effects in multi-level studies Changing age distribution of the COVID-19 pandemic-United States Applicability of Supervised Classification for Subdividing Forested Areas Using SPOT-5 and KOMPSAT-2 Data Control CfD, Prevention. 2020. Health equity considerations and racial and ethnic minority groups The role of air pollution (PM and NO2) in COVID-19 spread and lethality: a systematic review Disentangling associations between vegetation greenness and dengue in a Latin American city: Findings and challenges How can vegetation protect us from air pollution? A critical review on green spaces' mitigation abilities for air-borne particles from a public health perspective-with implications for urban planning Use of Databases Available on the Web to Describe Covid-19 Morbidity and Mortality Trends Estimation of Risk Factor's Contribution to mortality from COVID-19 in Highly Populated European Countries The effect of residential urban greenness on allergic respiratory diseases in youth: A narrative review Accepting the limits of ecologic studies: Drs. Greenland and Robins reply to Drs Spatial clustering of dengue fever incidence and its association with surrounding greenness Understanding South Korea's response to the COVID-19 outbreak: a real-time analysis Interaction between residential greenness and air pollution mortality: analysis of the Chinese longitudinal healthy longevity survey Environmental forcing and densitydependent controls of Culex pipiens abundance in a temperate climate (Northeastern Italy) Korea community health survey data profiles County-level exposures to greenness and associations with COVID-19 incidence and mortality in the United States Pathways linking biodiversity to human health: A conceptual framework Does access to green space impact the mental well-being of children: A systematic review Air pollution removal by urban trees and shrubs in the United States. Urban For. Urban Green Kriging: a method of interpolation for geographical information systems Contact transmission of Covid-19 in South Korea: Novel investigation techniques for tracing contacts Social Distancing and Outdoor Physical Activity During the COVID-19 Outbreak in South Korea: Implications for Physical Distancing Strategies Application of testing-tracing-treatment strategy in response to the COVID-19 outbreak in Seoul The transmission of SARS-CoV-2 is likely comodulated by temperature and by relative humidity Nature's role in supporting health during the COVID-19 pandemic: A geospatial and socioecological study Regulation of the immune system by biodiversity from the natural environment: an ecosystem service essential to health Biodiversity intervention enhances immune regulation and health-associated commensal microbiota among daycare children Greenspace exposure and J o u r n a l P r e -p r o o f COVID-19 mortality in the United States Comparison of the second and third waves of the COVID-19 pandemic in South Korea: Importance of early public health intervention Impact of temperature on the dynamics of the COVID-19 outbreak in China Spatial variability in reproduction number and doubling time across two waves of the COVID-19 pandemic in South Korea Development of composite deprivation index for Korea: the correlation with standardized mortality ratio Peer Reviewed: Recommendations for keeping parks and green space accessible for mental and physical health during COVID-19 and other pandemics Impacts of the COVID-19 pandemic on human-nature interactions: Pathways, evidence and implications World Health Organization declares global emergency: A review of the 2019 novel coronavirus (COVID-19) Nature inequity and higher COVID-19 case rates in less-green neighbourhoods in the United States Greenness availability and respiratory health in a population of urbanised children in North-Western Italy Air pollution associated respiratory mortality risk alleviated by residential greenness in the Chinese Elderly Health Service Cohort Social distancing, population density, and spread of COVID-19 in England: a longitudinal study Problems of correlations between explanatory variables in multiple regression analyses in the dental literature Effects of the COVID-19 pandemic on the use and perceptions of urban green space: An international exploratory study. Urban For. Urban Green Urban nature in a time of crisis: recreational use of green space increases during the COVID-19 outbreak in Oslo High temperature and high humidity reduce the transmission of COVID-19 Exposure to air pollution and COVID-19 mortality in the United States Air pollution and COVID-19 mortality in the United States: Strengths and limitations of an ecological regression analysis Systematic assessment of South Korea's capabilities to control COVID-19