Next Article in Journal
An IoT Framework for Screening of COVID-19 Using Real-Time Data from Wearable Sensors
Previous Article in Journal
Association between Triglyceride Glucose Index and Corrected QT Prolongation in Chinese Male Steelworkers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of Vulnerable Populations and Areas at Higher Risk of COVID-19-Related Mortality during the Early Stage of the Epidemic in the United States

1
Department of Geography and Geographic Information Science, University of Cincinnati, Cincinnati, OH 45220, USA
2
Health Geography and Disease Modeling Laboratory, University of Cincinnati, Cincinnati, OH 45220, USA
3
Division of Asthma Research, Cincinnati Children’s Hospital Medical Center, Department of Pediatrics, University of Cincinnati, Cincinnati, OH 45220, USA
4
Department of Biostatistics, College of Public Health and Human Sciences, Oregon State University, Corvallis, OR 97331, USA
5
Geospatial Health Advising Group, University of Cincinnati, Cincinnati, OH 45220, USA
6
Medical College of Georgia, Augusta University, Augusta, GA 30912, USA
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2021, 18(8), 4021; https://doi.org/10.3390/ijerph18084021
Submission received: 24 February 2021 / Revised: 5 April 2021 / Accepted: 6 April 2021 / Published: 12 April 2021

Abstract

:
We characterized vulnerable populations located in areas at higher risk of COVID-19-related mortality and low critical healthcare capacity during the early stage of the epidemic in the United States. We analyze data obtained from a Johns Hopkins University COVID-19 database to assess the county-level spatial variation of COVID-19-related mortality risk during the early stage of the epidemic in relation to health determinants and health infrastructure. Overall, we identified highly populated and polluted areas, regional air hub areas, race minorities (non-white population), and Hispanic or Latino population with an increased risk of COVID-19-related death during the first phase of the epidemic. The 10 highest COVID-19 mortality risk areas in highly populated counties had on average a lower proportion of white population (48.0%) and higher proportions of black population (18.7%) and other races (33.3%) compared to the national averages of 83.0%, 9.1%, and 7.9%, respectively. The Hispanic and Latino population proportion was higher in these 10 counties (29.3%, compared to the national average of 9.3%). Counties with major air hubs had a 31% increase in mortality risk compared to counties with no airport connectivity. Sixty-eight percent of the counties with high COVID-19-related mortality risk also had lower critical care capacity than the national average. The disparity in health and environmental risk factors might have exacerbated the COVID-19-related mortality risk in vulnerable groups during the early stage of the epidemic.

1. Introduction

Early Coronavirus Disease 2019 (COVID-19) data from Europe and Asia suggested high transmission and death rates. In late March of 2020, the United States (U.S.) had the fastest growing population impact across developed countries, leading to about 6,000,000 confirmed COVID-19 cases and 180,000 deaths by September 1. However, there was substantial spatial and temporal variation in the dynamics of the epidemic within the U.S., [1,2] with urban areas in New York, New Jersey, and Maryland experiencing the highest burden of infection early in the epidemic followed by a shift of the epicenter of the disease to rural areas in the South and Midwest late in 2020. Several countries including the U.S. have reported higher mortality rates (MR) for older people with comorbidities, including chronic lower respiratory diseases, diabetes, hypertension, and ischemic diseases among others [3,4,5,6].
Preliminary studies had started to create baseline population characteristics of COVID-19-related deaths during the early stage of the epidemic [7,8,9]. However, the roles of geospatial disparities, including sociodemographic and environmental exposures, and critical care capacity during the early stage of the epidemic are still not well understood. In general, identifying the population groups and areas at higher risk of infection and mortality during the early stage of an epidemic based on underlying health disparities and low critical healthcare capacity is a logical step to developing effective strategies for reducing the risk of infection and mortality throughout a country. Therefore, this study aims to: (a) assess sociodemographic and environmental drivers of COVID-19-related deaths and (b) identify vulnerable areas at higher risk of COVID-19 mortality but with low healthcare capacity during the early stage of the COVID-19 epidemic in the U.S. We use counts of confirmed deaths caused or contributed by COVID-19 to assess the impact of sociodemographic and environmental variables on the epidemic [10,11]. We hypothesized that during the early stage of the COVID-19 epidemic in the U.S., counties with predominantly minority populations, poor air quality, and low critical healthcare capacity were adversely affected by the epidemic.

2. Materials and Methods

2.1. Study Area and Data Sources

Counts of COVID-19-related deaths in the U.S. were obtained from a Johns Hopkins University database [12] for 3142 counties from all 50 states and the District of Columbia for the time period of 22 January to 1 September of 2020. This time period is the early stage of the epidemic in the U.S. that includes the initiation of intervention measures against the pandemic (e.g., reduced interactions and business closures) in March. Sociodemographic data were derived from the 2014–2018 American Community Survey (ACS) [13] and the Center for Disease Control (CDC) Social Vulnerability Index [14]. County-level comorbidities, including chronic lower respiratory disease (CLRD), diabetes mellitus, hypertensive diseases (HTA), and ischemic heart disease were obtained between 2010 and 2018 from U.S. Centers for Disease Control and Prevention [15]. Air pollution was quantified using surface annual particulate matter satellite images from 2000 to 2018 [16]. All data used in this study were recorded and analyzed at the county level.

2.2. Study Variables

The primary outcome variable in this study was COVID-19-related deaths. The projected expected number of deaths was calculated using an indirect standardization with the ACS standard population. We used the ratio between the number of COVID-19 deaths in a county and the number of deaths that would have been expected to occur based on a standardized population size. Demographic variables included the percent of the population in different age groups (<25, 25–34, 35–44, 45–59, 60–74, and ≥75 years), the percent of the population who self-identified as White, African American, or other races, and, following the Census Bureau definition of ethnicity [17], the percent of the Hispanic or Latino population. The percent of the population in poverty was obtained using the CDC’s Vulnerability Index [14]. We used 2000 to 2018 annual images of ground-level fine particle matter <2.5 μm (PM2.5) to assess the association between COVID-19-related mortality and long-term exposure of a county to several years of air pollutants [16]. These calibrated images are estimated at a 0.01° × 0.01° grid resolution combining satellite and monitoring stations data sources using a geographically weighted regression. We computed long-term exposure by averaging PM2.5 from 2010 to 2018 within each county (See Figure S1). We also calculated a regional air hub and road Connectivity Index for each county to examine the association between deaths and county-level airport hubs and main roads (See Figure S2). A categorical covariate with four geospatial connectivity levels was generated as follows: counties with an airport with more than 50,000 passengers per year (has an airport), counties next to a county with an airport (next to airport), counties crossed by a main road (crossed by a highway), and counties not neighboring a county with an airport and not crossed by a main road (no airport/highway). Due to the reported connection between COVID-19 mortality and underlying health problems [3,4,5,6,18,19], we included four chronic conditions as predictor variables for COVID-19 mortality. Crude county-level mortality rates per 100,000 people between 2010 and 2018 were included, obtained from the CDC Wonder database [15], including CLRD, diabetes mellitus, HTA, and ischemic heart disease. This study follows the guidelines of the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) [20]. Additional information about variables and STROBE checklist is included in the Supplementary Materials (See Table S1).

2.3. Multivariate Analyses of Risk Factors for COVID-19-Related Death

All covariates were selected according to an evidence synthesis process of preliminary reports [7,18,21,22,23,24]. A directed acyclic graph was built to infer causal effects to the observational data. We then removed open paths and checked for colliders and over control in the implied graph (See Figure S3) [25]. An exploratory data analysis of COVID-19 mortality rates was performed using linear regression for each covariate. For disease mapping, a Bayesian spatial Poisson model was used to estimate relative risks of COVID-19-related mortality [26,27]. Spatially correlated random effects were modeled as conditional autoregressive (CAR) priors using an adjacency county matrix (See Figure S4). We calculated a Moran’s I statistic to evaluate for spatial clustering of the residuals of the final model (See Figure S5). Quantile population classification was used to identify the 10 highest COVID-19 mortality risk areas in highly populated counties (4th quartile). A bivariate map combining COVID-19-related mortality risk and number of intensive care unit (ICU) beds was generated to identify areas with low critical healthcare capacity. ICU beds per 100,000 people were included as an index of critical healthcare capacity for each county [28]. All models were adjusted for age group due to the strong association between age and COVID-19 death. All numeric covariates were centered by subtracting the sample mean. The R language with the INLA, SpatialEpi, and raster packages was used for data analysis [29,30,31,32]. Detailed information is included in the Supplementary Materials 2 and 3.

3. Results

3.1. Sociodemographic and Socio-Environmental Variables

Table 1 presents descriptive statistics of COVID-19 deaths included in the study. The total number of deaths reported was 181,937, corresponding to 3.1% of the 5,958,655 COVID-19 confirmed cases. The highest standardized mortality ratios were found in NJ (SMR = 3.19; 15,950 deaths), NY (SMR = 2.98; 32,942 deaths), MA (SMR = 2.35; 9058 deaths), CT (SMR = 2.21; 4466 deaths), and LA (SMR = 1.83; 4821 deaths). Of the 3142 counties included in the study, 2460 had at least one confirmed COVID-19 death. The overall percent estimate of poverty was 15.6 (standard deviation (SD) 6.5) for the entire country. On average, the county-level percent of the population who were white was 83.0 (SD 16.7), while the average percent of African American was 9.1 (SD 14.6) and for other races, it was 7.9% (SD 10.2). For ethnicity, the average percent of Hispanic or Latino population was 9.3 (SD 13.9). The national average PM2.5 exposure was 8.0 µg/m (SD 2.4). For the Connectivity Index, 232 counties had an airport with more than 50,000 passengers per year, 1047 counties were next to airport hubs, 629 counties had a highway or main road, and 1200 counties are categorized as low transportation connectivity. The overall ICU beds capacity was 28.4 per 100,000 (SD 34.6).
The COVID-19-related mortality rate per 100,000 people was higher in counties with high poverty (51.42 per 100,000), PM2.5 greater than 9.5 µg/m (50.75 per 100,000), and counties with an airport (42.23 per 100,000) (Figure S2). Counties with a high percentage of African Americans (60.43 per 100,000) and other races (54.95 per 100,000) had higher average mortality rates, whereas counties with a high proportion of white population had lower rates (15.82 per 100,000). Latino population also had a higher COVID-19-related mortality rate (40.55 per 100,000). Finally, counties with a high proportion of CLRD, HTA, and ischemic heart diseases had high COVID-19-related mortality rates of 51.42, 45.55, and 27.58 per 100,000 people, respectively. Detailed information is included in the Supplementary Materials 2 and 3 (See Figure S6 and Table S3)

3.2. Multivariate Analyses of Risk Factors for COVID-19-Related Death

Estimated relative risks (RR) of COVID-19 mortality from a Bayesian spatial Poisson regression analysis are presented in Table 2. Population over 75 years old or more was associated to 5% higher risk of mortality compared to population under 25 (RR = 1.05, 95% credible interval (CI): 1.01–1.08). For race, black and other races exhibit higher risk of COVID-19-related death compared to white population, 1% (RR = 1.01, CI: 1.01–1.02) and 2% (RR = 1.02, CI: 1.01–1.02) respectively. In ethnicity, the proportion of Hispanic or Latino population if infected with COVID-19 was associated with a 2% higher risk of COVID-19-related death (RR = 1.02, CI: 1.02–1.03). We found no association between the proportion of poverty and comorbidities and the risk of COVID-19-related death at county-level. For the long-term exposure to air pollution, we found that every additional unit of PM2.5 (1.0 µg/m) increased the risk of COVID-related death by 14% (RR = 1.14, CI: 1.08–1.20). Lastly, counties with an airport (RR = 1.31, CI: 1.14–1.51), and counties near to airports (RR = 1.13, CI: 1.03–1.24), had a higher risk of COVID-19-related death compared to counties with low transportation connectivity. More detailed information, including model formulation is included in the Supplementary Materials 3.

3.3. COVID-19 Disease Mapping

Figure 1 and Table S4 illustrate the RR for COVID-19 at the state level. Seven of the 51 territories (50 states and D.C.) had a risk higher than average, namely, IA (RR = 1.87, 95% credible interval (CI): 1.13–3.11), IN (RR = 1.92, CI: 1.18–3.15), MA (RR = 3.05, CI: 1.51–6.08), LA (RR = 2.19, CI: 1.27–3.83), MS (RR = 1.76, CI: 1.05–2.97), TX (RR = 1.80, CI: 1.07–3.04), and AZ (RR = 3.08, CI: 1.52–6.34). Conversely, only four of the 50 states (RI, AK, HI, and WY) had an RR lower than the national average.
Maps in Figure 2 illustrate the COVID-19 relative risk of death for all counties included in the study. Overall, we found that 315 counties from 33 states and D.C. had a higher risk of COVID-19-related mortality. Five states had 21 or more counties with high mortality risk, including GA (39), LA (42), MS (44), NJ (21), and TX (27). About 46% of the overall variance was explained by states, and nearly 44% was spatially structured indicating a spatial correlation of COVID-19-relative risk of death. This spatial pattern is present in Northeast and Southern states (Moran’s I p-value < 0.001). Table 3 shows the 10 highest COVID-19 mortality areas in highly populated counties (4th quartile). NY (4) and NJ (5) had nine of the top-10 highest risk locations. Six of 10 counties had higher poverty levels than the national county average (15.6%). These counties had on average a lower proportion of white population (48.0%) and higher proportions of black population (18.7%) and other races (33.3%) compared to the national averages of 83.0%, 9.1%, and 7.9%, respectively. Similarly, Hispanic and Latino population proportions were higher in these counties (29.3% vs. 9.3%). Nine of these 10 counties had a long-term PM2.5 exposure of at least 1.6 µg/m above the national average (8.0 µg/m). Nine of these 10 counties had an airport or were next to a county with an airport and 9 out of 10 had lower ICU beds availability than the national average of 28.4 ICU beds per 100,000 inhabitants.
Figure 3 illustrates a bivariate map of the COVID-19-related mortality risk and ICU beds availability for all counties in the conterminous U.S. There were counties with moderate to high mortality risk and low ICU availability in the Northeastern and southern regions. About 53% of these counties were in GA, LA, MS, NJ, and TX. In addition, areas with low risk of COVID-19-related mortality risk but high ICU availability were observed in KS, MT, ND, NE, and SD. Additional details can be found in the Supplementary Materials 4.

4. Discussion

This study provides state- and county-level characterization of the COVID-19-related mortality risk, including sociodemographic and socio-environmental factors during the early stage of the epidemic across the U.S. In addition, our study assessed the spatial link between COVID-19-related mortality risk and current critical healthcare capacity. Overall, we identified highly populated and polluted areas, regional air hub areas, race minorities (non-white population), and Hispanic or Latino population with an increased risk of COVID-19-related death during the first phase of the epidemic. The 10 most populated counties with the highest mortality risk had about fourfold higher risk than the national average with higher proportions of minorities residing in these counties. Finally, our spatial analysis highlights high-risk and low ICU availability areas that are worth noting for geographically targeted interventions, particularly during an early phase of an epidemic.
We found that the 10 highest COVID-19 mortality areas in highly populated counties (4th quartile) had about fourfold higher mortality risk than the national average during the early stage of the epidemic. Notably, these counties had on average, a lower white population and higher proportions of population in poverty, African American, other minorities, and Hispanic or Latino populations compared to their national averages. States and counties with a historically lower proportion of white population such as LA, NJ, and TX were at greater COVID-19-related mortality risk than other states. These demographic disparities in COVID-19-related mortality have been recognized in preliminary results of several major cities in other countries, including London [21,33,34]. Moreover, non-pharmaceutical interventions (e.g., school closing, physical distancing, lockdowns, and additional sanitation) are difficult to implement in these groups [35]. As a result, the effectiveness and benefit of these non-pharmaceutical interventions can be diluted by the work activities that involve person to person interaction and are more common in these low-income groups exposing them to a higher risk of infection and thus, higher mortality risk.
Air pollution was directly associated with higher COVID-19-related mortality risk. Of note, 8 of the top 10 counties with the highest mortality risk were very close to or above the cutoff for healthy air quality levels (PM2.5 of 12.0 µg/m), according to the U.S. national ambient air quality standards for particle pollution [36]. Air pollution is one of the leading risk factors for respiratory-related death globally [34], and this factor could be playing a key role in exacerbating the numbers of COVID-19-related deaths in highly polluted areas. Air pollution has an indirect impact on most of the organs and systems of the human body and indirectly comorbidities. Although we did not find significant associations between COVID-19 and comorbidities at the county level, air pollution has been identified as a contributing factor for many respiratory diseases like chronic obstructive pulmonary disease (COPD) [33,37], asthma [37,38,39,40,41,42], and lung cancer [43,44,45,46], which are concomitant comorbidities that had a strong association with COVID-19-related deaths at the individual level [3]. The health effects of air pollution depend on the components and sources of pollutants, which can vary among counties, seasons, and times. Initial evidence of incidence and mortality with comorbidities have been reported in Italy, with regional differences between the northern and southern region [8,23]. Although we found a strong positive association between air pollution and the risk of COVID-19-related death, the role of long-term exposure to poor air quality on COVID-19-related deaths in the U.S. is still not well understood, and thus more studies are needed, including for major U.S. cities taking into account long-term exposure to outdoor and indoor pollution. Furthermore, our results suggest that counties with airports have a higher COVID-19-related mortality risk than those with less connectivity. The lift of travel restrictions and high connectivity in these counties generated by airports can produce a high influx of imported infections that boost the local transmission of the virus in the county and consequently the number of COVID-19-related deaths.
According to our spatial analysis, 315 counties (10.0% of the total number of counties included in the study) from 34 states had a higher risk of COVID-19-related mortality than the national average. About 55% (173 out of 315) of these counties were located in five states: GA, LA, MS, NJ, and TX. Reasons for the increased mortality risk in these areas could be numerous. First, NJ shares borders with NY, which had experienced an intensive COVID-19 outbreak with about 5% of worldwide cases [47] and about 18% of U.S. deaths. Many people that work in the New York metropolitan area reside in NJ border counties, which might have contributed to the initial spread of the disease. Second, most of these states had a lower amount of white population and higher proportions of at least one minority group (African American and Latino) than the national average (See Table S4). State-average poverty was larger than the 15.6% study average for GA (20.7%), LA (22.0%), MS (24.1%), and TX (16.1%), but not NJ (10.2%). In terms of air pollution, four of these five states had long-term PM2.5 averages above the national average (8.0 µg/m), ranging from 9.4 µg/m (LA) to 10.6 µg/m (GA). For high populated and high COVID-19-related mortality risk counties, 9 of the 10 most populated counties were located in the New York metropolitan area, which is ranked 12 out of the 25 most ozone-polluted cities in the U.S. [48].
Our bivariate analysis showed that GA (22), LA (23), MS (21), NJ (18), and TX (17) represent 53% of the counties with moderate to high COVID-19-related mortality risk but lower critical care capacity than the national average of 28.4 ICU per 100,000 inhabitants. Most of the counties (79 out of 192) were also in areas with lower critical care capacity than the national average of 28.4 ICU per 100,000 inhabitants. With the onset of COVID-19 and the current lift of lockdown measures across the U.S., critical healthcare capacity might be potentially overwhelmed in several of these counties, not only with respect to ICU beds capacity but also in mechanical ventilators and staffing. Therefore, counties with high COVID-19-related mortality but low healthcare capacity identified in our study should be prioritized for strategies aimed at diminishing the overall number of COVID-19-related deaths, including patient relocation, strengthening of critical healthcare infrastructure and supply chains, and staff step-up [49].
Our study has several limitations worth noting. First, COVID-19 data that include comorbidity data are not available at the unit of analysis (county level). This issue might hamper precise comparisons of the real epidemic burden in specific groups and comorbidities. Second, we analyzed air pollution based on PM2.5 measures at county-level resolution. Smaller geographic areas and other pollutants including sulfate (SO4), nitrate (NO3), ammonium (NH4), organic matter (OM), black carbon (BC), mineral dust (DUST), and sea-salt (SS) would produce more refined inference [50]. Third, our Connectivity Index is based in a simplified geospatial attribute. Only airports and main roads were used as a Connectivity Index for each county. Other potential geospatial and population characteristics might influence the influx of travelers, and, therefore, the transmission dynamic. A further limitation relates to the challenges in translating cross-sectional associations into conclusions on causation of COVID-19-related deaths at county level. In addition, an ecological study is an approach for examining the association between factors and diseases by conducting analysis at the population level in specific areal units. In ecological studies, due to the lack of individual data, it is difficult to adjust for all potential confounding factors even if multivariate analysis is performed, leading to a possible ecological fallacy. Hence, our results should be interpreted with caution. Finally, this analysis only focuses on the early stage of this epidemic and results might not be comparable through time. Additional analysis focused on subsequent time periods and epidemic phases might reveal other disparities and dominant factors.

5. Conclusions

This study used geospatial approaches to examine population risk determinants of COVID-19-related deaths at the county level and identify vulnerable areas and populations at higher risk of COVID-19-related mortality during the early stage of the epidemic. The results have significant public health implications with respect to the critical healthcare infrastructure for an effective response to a pandemic, particularly during an epidemic’s early phase. The social gradient of health and environment in which underserved groups are highly vulnerable to more severe health outcomes can also be an important driver of the geographical and social disparity observed in the early stage of the COVID-19 epidemic in the U.S. Moreover, the healthcare capacity’s substantial regional disparities might increase the vulnerability of areas already at higher risk of disease spread during the early stage of the epidemic. Therefore, the results from this study can be used to guide the development of strategies for identifying locations and groups for targeted prevention efforts in vulnerable communities at high risk of disease spread and mortality during the early stage of an epidemic in the U.S.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ijerph18084021/s1, Figure S1: U.S. 2000 to 2018 Long-Term Mean PM2.5 Concentrations by County, mean = 7.98 µg/m (range is 1.42–13.30), Figure S2: U.S. Connectivity index by county, Figure S3: Directed acyclic graph for Coronavirus Disease 2019 (COVID-19) mortality, Figure S4: U.S. counties adjacency matrix for the intrinsic CAR model, Figure S5: Bayesian spatial random effects (σ), Moran’s I statistic standard deviate = −4.61, p-value = 1, Figure S6: Exploratory data analysis covariates vs COVID-19 mortality rate. Table S1: Strobe checklist, Table S2: State Abbreviations List, Table S3: State summary for sociodemographic factors, Table S4: Relative risk by state.

Author Contributions

Conceptualization: E.C.-A., T.B.M., and D.F.C.; methodology: E.C.-A., T.B.M., and D.F.C.; software: E.C.-A.; validation: E.C.-A., D.F.C., A.J.B., and N.J.M.; writing—original draft preparation: E.C.-A. and D.F.C., writing—review and editing: E.C.-A., T.B.M., A.J.B., and N.J.M. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by NIH NHLBI award R01HL132344 and NIEHS award P30ES006096 from Center for Environmental Genetics (CEG), Department of Environmental and Public Health Sciences, University of Cincinnati.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors thank to Centers for Disease Control and Prevention (CDC) for releasing the national surveys in the service of science and the United States Agencies and other donors supporting these initiatives.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. CDC COVID-19 Response Team. Geographic Differences in COVID-19 Cases, Deaths, and Incidence—United States, February 12–April 7, 2020. Morb. Mortal. Wkly. Rep. 2020, 69, 465–471. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Miller, I.F.; Becker, A.D.; Grenfell, B.T.; Metcalf, C.J.E. Disease and healthcare burden of COVID-19 in the United States. Nat. Med. 2020, 26, 1212–1217. [Google Scholar] [CrossRef]
  3. Richardson, S.; Hirsch, J.S.; Narasimhan, M.; Crawford, J.M.; McGinn, T.; Davidson, K.W. The Northwell COVID-19 Research Consortium. Presenting Characteristics, Comorbidities, and Outcomes among 5700 Patients Hospitalized with COVID-19 in the New York City Area. JAMA 2020. Available online: https://jamanetwork.com/journals/jama/fullarticle/2765184 (accessed on 20 April 2020). [CrossRef] [PubMed]
  4. Shi, S.; Qin, M.; Shen, B.; Cai, Y.; Liu, T.; Yang, F.; Gong, W.; Liu, X.; Liang, J.; Zhao, Q.; et al. Association of Cardiac Injury with Mortality in Hospitalized Patients with COVID-19 in Wuhan, China. JAMA Cardiol. 2020. Available online: https://jamanetwork.com/journals/jamacardiology/fullarticle/2763524 (accessed on 20 April 2020). [CrossRef] [Green Version]
  5. Grasselli, G.; Zangrillo, A.; Zanella, A.; Antonelli, M.; Cabrini, L.; Castelli, A.; Cereda, D.; Coluccello, A.; Foti, G.; Fumagalli, R.; et al. Baseline Characteristics and Outcomes of 1591 Patients Infected with SARS-CoV-2 Admitted to ICUs of the Lombardy Region, Italy. JAMA 2020, 323, 1574–1581. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Wu, M.; Huang, S.; Liu, J.; Shu, Y.; Luo, Y.; Wang, L.; Li, M.; Wang, Y. Association between comorbidities and the risk of death in patients with COVID-19: Sex-specific differences. medRxiv 2020. Available online: https://www.medrxiv.org/content/10.1101/2020.05.22.20109579v1 (accessed on 29 May 2020). [CrossRef]
  7. Williamson, E.; Walker, A.J.; Bhaskaran, K.J.; Bacon, S.; Bates, C.; Morton, C.E.; Curtis, H.J.; Mehrkar, A.; Evans, D.; Inglesby, P.; et al. OpenSAFELY: Factors associated with COVID-19-related hospital death in the linked electronic health records of 17 million adult NHS patients. medRxiv 2020. Available online: https://www.medrxiv.org/content/10.1101/2020.05.06.20092999v1 (accessed on 10 May 2020). [CrossRef]
  8. Remuzzi, A.; Remuzzi, G. COVID-19 and Italy: What next? Lancet 2020, 395, 1225–1228. [Google Scholar] [CrossRef]
  9. Li, R.; Rivers, C.; Tan, Q.; Murray, M.B.; Toner, E.; Lipsitch, M. Estimated Demand for US Hospital Inpatient and Intensive Care Unit Beds for Patients with COVID-19 Based on Comparisons with Wuhan and Guangzhou, China. JAMA Netw. Open 2020, 3, e208297. [Google Scholar] [CrossRef] [PubMed]
  10. Schneider, E.C. Failing the Test—The Tragic Data Gap Undermining the U.S. Pandemic Response. N. Engl. J. Med. 2020, 383, 299–302. [Google Scholar] [CrossRef]
  11. CDC. COVID-19 Death Data and Resources. Available online: https://www.cdc.gov/nchs/nvss/covid-19.htm (accessed on 2 April 2021).
  12. Johns Hopkins University. 2019 Novel Coronavirus COVID-19 (2019-nCoV) Data Repository by Johns Hopkins CSSE. Available online: https://github.com/CSSEGISandData/COVID-19 (accessed on 24 April 2020).
  13. Bureau, U.S.C. American Community Survey 2014–2018 5-Year Estimates. Available online: https://www.census.gov/data/developers/data-sets/acs-5year.2018.html (accessed on 20 April 2020).
  14. Gregory, E.W.; Hallisey, E.; Heitgerd, J.L.; Lewis, B. A Social Vulnerability Index for Disaster Management. J. Homel. Secur. Emerg. Manag. 2011, 8, 3. [Google Scholar]
  15. Centers for Disease Control and Prevention, National Center for Health Statistics. Underlying Cause of Death 1999–2018 on CDC WONDER Online Database, Released in 2020. Available online: http://wonder.cdc.gov/ucd-icd10.html (accessed on 26 April 2020).
  16. van Donkelaar, A.; Martin, R.V.; Li, C.; Burnett, R.T. Regional Estimates of Chemical Composition of Fine Particulate Matter Using a Combined Geoscience-Statistical Method with Information from Satellites, Models, and Monitors. Environ. Sci. Technol. 2019, 53, 2595–2611. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Bureau, C. Race & Ethnicity. Available online: https://bber.unm.edu/media/files/race-ethnicity-onepager.pdf (accessed on 20 April 2020).
  18. CDC COVID-19 Response Team. Preliminary Estimates of the Prevalence of Selected Underlying Health Conditions among Patients with Coronavirus Disease 2019. Morb. Mortal. Wkly. Rep. 2020, 69, 382–386. [Google Scholar]
  19. Hsu, H.E.; Ashe, E.M.; Silverstein, M.; Hofman, M.; Lange, S.J.; Razzaghi, H.; Mishuris, R.G.; Davidoff, R.; Parker, E.M.; Penman-Aguilar, A.; et al. Race/Ethnicity, Underlying Medical Conditions, Homelessness, and Hospitalization Status of Adult Patients with COVID-19 at an Urban Safety-Net Medical Center—Boston, Massachusetts, 2020. Morb. Mortal. Wkly. Rep. 2020, 69, 864–869. [Google Scholar] [CrossRef] [PubMed]
  20. Von Elm, E.; Altman, D.G.; Egger, M.; Pocock, S.J.; Gotzche, P.C.; Vandenbroucke, J.P. The Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) Statement: Guidelines for Reporting Observational Studies. Epidemiology 2007, 18, 800–804. [Google Scholar] [CrossRef] [Green Version]
  21. Niedzwiedz, C.L.; Donnell, C.A.; Jani, B.D.; Demou, E.; Ho, F.K.; Celis-Morales, C.; Nicholl, B.I.; Mair, F.; Welsh, P.; Sattar, N.; et al. Ethnic and socioeconomic differences in SARS-CoV-2 infection: Prospective cohort study using UK Biobank. medRxiv 2020, 18, 1–14. [Google Scholar] [CrossRef]
  22. Halpin, D.M.G.; Faner, R.; Sibila, O.; Badia, J.R.; Agusti, A. Do chronic respiratory diseases or their treatment affect the risk of SARS-CoV-2 infection? Lancet Respir. Med. 2020, 8, 436–438. [Google Scholar] [CrossRef]
  23. Pansini, R.; Fornacca, D. Initial evidence of higher morbidity and mortality due to SARS-CoV-2 in regions with lower air quality. medRxiv 2020. Available online: https://www.medrxiv.org/content/medrxiv/early/2020/04/07/2020.04.04.20053595.full.pdf (accessed on 10 May 2020). [CrossRef]
  24. Fattorini, D.; Regoli, F. Role of the atmospheric pollution in the Covid-19 outbreak risk in Italy. medRxiv 2020, 264, 114732. [Google Scholar] [CrossRef] [Green Version]
  25. Hernán, M.A.; Hernández-Díaz, S.; Robins, J.M. A Structural Approach to Selection Bias. Epidemiology 2004, 15, 615–625. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Besag, J.; York, J.; Mollié, A. Bayesian image restoration, with two applications in spatial statistics. Ann. Inst. Stat. Math. 1991, 43, 1–20. [Google Scholar] [CrossRef]
  27. Lu, Z. Handbook of Spatial Statistics; Gelfand, A.E., Diggle, P.J., Montserrat Fuentes and Peter Guttorp, Eds.; Chapman & Hall/CRC: Boca Raton, FL, USA, 2010; p. xii+607. ISBN 978-1-4200-7287-7. [Google Scholar] [CrossRef]
  28. Healthcare, D. USA Hospital Beds. Available online: https://coronavirus-resources.esri.com/datasets/1044bb19da8d4dbfb6a96eb1b4ebf629_0?geometry=80.507%2C-16.820%2C-105.469%2C72.123 (accessed on 5 May 2020).
  29. Hijmans, R.J. Raster: Geographic Data Analysis and Modeling. 2019. Available online: https://rdrr.io/cran/raster/#vignettes (accessed on 1 April 2020).
  30. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2018. [Google Scholar]
  31. Bürkner, P.-C. brms: An R Package for Bayesian Multilevel Models Using Stan. J. Stat. Softw. 2017, 80, 1–28. [Google Scholar] [CrossRef] [Green Version]
  32. Kim, A.Y.; Jon, W. SpatialEpi: Methods and Data for Spatial Epidemiology; R Package Version 1.2.3; University of Washington: Seatle, WA, USA, 2018. [Google Scholar]
  33. Owen, W.F., Jr.; Carmona, R.; Pomeroy, C. Failing Another National Stress Test on Health Disparities. JAMA 2020, 323, 1905–1906. [Google Scholar] [CrossRef] [Green Version]
  34. NYC Department of Health and Mental Hygiene. COVID-19: Data. Available online: https://www1.nyc.gov/site/doh/covid/covid-19-data.page (accessed on 17 May 2020).
  35. Hernandez, A.; Correa-Agudelo, E.; Kim, H.; Branscum, A.J.; Miller, F.D.; MacKinnon, N.; Cuadros, D.F. On the impact of early non-pharmaceutical interventions as containment strategies against the COVID-19 pandemic. medRxiv 2020. Available online: https://www.medrxiv.org/content/10.1101/2020.05.05.20092304v1 (accessed on 10 May 2020). [CrossRef]
  36. United States Environmental Protection Agency. Revised Air Quality Standards for Particle Pollution and Updates to the Air Quality Index (AQI). Available online: https://www.epa.gov/sites/production/files/2016-04/documents/2012_aqi_factsheet.pdf (accessed on 28 March 2021).
  37. Kelly, F.J.; Fussell, J.C. Air pollution and airway disease. Clin. Exp. Allergy 2011, 41, 1059–1071. [Google Scholar] [CrossRef] [PubMed]
  38. Wu, S.; Ni, Y.; Li, H.; Pan, L.; Yang, D.; Baccarelli, A.A.; Deng, F.; Chen, Y.; Shima, M.; Guo, X. Short-term exposure to high ambient air pollution increases airway inflammation and respiratory symptoms in chronic obstructive pulmonary disease patients in Beijing, China. Environ. Int. 2016, 94, 76–82. [Google Scholar] [CrossRef] [PubMed]
  39. Minelli, C.; Wei, I.; Sagoo, G.; Jarvis, D.; Shaheen, S.; Burney, P. Interactive effects of antioxidant genes and air pollution on respiratory function and airway disease: A HuGE review. Am. J. Epidemiol. 2011, 173, 603–620. [Google Scholar] [CrossRef] [PubMed]
  40. Viegi, G.; Simoni, M.; Scognamiglio, A.; Baldacci, S.; Pistelli, F.; Carrozzi, L.; Annesi-Maesano, I. Indoor air pollution and airway disease. Int. J. Tuberc. Lung Dis. 2004, 8, 1401–1415. [Google Scholar] [PubMed]
  41. Ulmer, W.T. Air pollution and airway disease. Der Internist 1985, 26, 233–240. [Google Scholar] [PubMed]
  42. Karakatsani, A.; Analitis, A.; Perifanou, D.; Ayres, J.G.; Harrison, R.M.; Kotronarou, A.; Kavouras, I.G.; Pekkanen, J.; Hameri, K.; Kos, G.P.; et al. Particulate matter air pollution and respiratory symptoms in individuals having either asthma or chronic obstructive pulmonary disease: A European multicentre panel study. Environ. Health 2012, 11, 75. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Vineis, P.; Hoek, G.; Krzyzanowski, M.; Vigna-Taglianti, F.; Veglia, F.; Airoldi, L.; Autrup, H.; Dunning, A.; Garte, S.; Hainaut, P.; et al. Air pollution and risk of lung cancer in a prospective study in Europe. Int. J. Cancer 2006, 119, 169–174. [Google Scholar] [CrossRef] [PubMed]
  44. Raaschou-Nielsen, O.; Vineis, P.; Brunekreef, B.; Nieuwenhuijsen, M.; Hoffmann, B.; Forastiere, F.; Oudin, A.; Hoek, G.; Katsouyanni, K.; Schwarze, P.; et al. Air pollution and lung cancer in Europe—Authors’ reply. Lancet Oncol. 2013, 14, e440. [Google Scholar] [CrossRef]
  45. Sax, S.N.; Zu, K.; Goodman, J.E. Air pollution and lung cancer in Europe. Lancet Oncol. 2013, 14, e439–e440. [Google Scholar] [CrossRef]
  46. Nawrot, T.S.; Nackaerts, K.; Hoet, P.H.; Nemery, B. Lung cancer mortality and fine particulate air pollution in Europe. Int. J. Cancer 2007, 120, 1825–1826. [Google Scholar] [CrossRef] [PubMed]
  47. Nature. Coronavirus: The First Three Months As It Happened. Available online: https://www.nature.com/articles/d41586-020-00154-w (accessed on 19 May 2020).
  48. American Lung Association. State of the Air 2020; American Lung Association: Chicago, IL, USA, 2020. [Google Scholar]
  49. Institute for Health Metrics and Evaluation. New COVID-19 Forecasts: US Hospitals Could Be Overwhelmed in the Second Week of April by Demand for ICU Beds, and US Deaths Could Total 81,000 by July. Available online: www.healthdata.org/news-release/new-covid-19-forecasts-us-hospitals-could-be-overwhelmed-second-week-april-demand-icu (accessed on 10 May 2020).
  50. Achakulwisut, P.; Brauer, M.; Hystad, P.; Anenberg, S.C. Global, national, and urban burdens of paediatric asthma incidence attributable to ambient NO2 pollution: Estimates from global datasets. Lancet Planet. Health 2019, 3, e166–e178. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Relative risk for Coronavirus Disease 2019 (COVID-19) mortality rates by state.
Figure 1. Relative risk for Coronavirus Disease 2019 (COVID-19) mortality rates by state.
Ijerph 18 04021 g001
Figure 2. U.S. Relative Risk for COVID-19 mortality by county, mean = 0.63 (CrI range is 0.00–14.8) (a). U.S. COVID-19-related relative risk (RR) of death and (b) COVID-19 spatial effect.
Figure 2. U.S. Relative Risk for COVID-19 mortality by county, mean = 0.63 (CrI range is 0.00–14.8) (a). U.S. COVID-19-related relative risk (RR) of death and (b) COVID-19 spatial effect.
Ijerph 18 04021 g002
Figure 3. Intensive care units (ICU) bed availability per 100,000 (ICU information was not available in AK and HI). Dark purple indicates counties with high ICU availability and low mortality risk, whereas areas in darker green-blue indicate counties with high mortality risk but low ICU availability. Both variables were classified with a tertile scheme as follows: COVID-19-related RR (0–1 lower risk, 1–3 medium risk, and 3 > high risk) and ICU beds per 100,000 (<28.4 low availability, 28.4–100 medium availability, and >100 high availability).
Figure 3. Intensive care units (ICU) bed availability per 100,000 (ICU information was not available in AK and HI). Dark purple indicates counties with high ICU availability and low mortality risk, whereas areas in darker green-blue indicate counties with high mortality risk but low ICU availability. Both variables were classified with a tertile scheme as follows: COVID-19-related RR (0–1 lower risk, 1–3 medium risk, and 3 > high risk) and ICU beds per 100,000 (<28.4 low availability, 28.4–100 medium availability, and >100 high availability).
Ijerph 18 04021 g003
Table 1. Baseline characteristics.
Table 1. Baseline characteristics.
COVID-19n
  Confirmed Cases5,958,655
  Confirmed Deaths181,937
SociodemographicMean(SD)
  Age (years)
    Percent under 2531.2(4.8)
    Percent 25–3411.8(2.3)
    Percent 35–4411.6(1.5)
    Percent 45–5920.2(2.2)
    Percent 60–7417.4(3.7)
    Percent 75+7.8(2.4)
Percent of population in poverty15.6(6.5)
Race
Percent of white population83.1(16.9)
Percent of black population9.1(14.5)
Percent of other races7.9(10.2)
Ethnicity
Percent of not Hispanic or Latino population90.7(13.8)
Percent of Hispanic or Latino population9.3(13.8)
Crude mortality ratesMean(SD)
Chronic lower respiratory disease 69.9(26.0)
Diabetes mellitus33.5(14.7)
Hypertension27.1(16.9)
Ischemic heart disease 151.2(57.2)
EnvironmentMean(SD)
Long-term PM2.5 exposure8.0(2.4)
Connectivity Index (n)
Counties with no airport/highway1200
Counties crossed by a highway629
Counties next to airport1047
Counties with an airport232
Table 2. Estimated relative risks and 95% credible intervals from a Bayesian spatial Poisson regression analysis.
Table 2. Estimated relative risks and 95% credible intervals from a Bayesian spatial Poisson regression analysis.
County-Level CovariatesRRCrI: [2.5%, 97.5%]
Sociodemographic
Age
Under 25RefRef
25–340.98(0.961.01)
35–441.01(0.971.04)
45–591.02(1.001.05)
60–740.98(0.961.00)
75+1.05(1.011.08)
Percentage of population in poverty1.01(1.001.02)
Race
Percent of white populationRefRef
Percent of black population1.01(1.011.02)
Percent of other races1.02(1.011.02)
Ethnicity
Percent of non-Hispanic or Latino populationRefRef
Percent of Hispanic or Latino population1.02(1.021.03)
Crude mortality rates
Chronic lower respiratory disease1.00(1.001.00)
Diabetes mellitus1.00(1.001.00)
Hypertension1.00(1.001.01)
Ischemic heart disease1.00(1.001.00)
Environment
Long-term exposure to PM2.51.14(1.081.20)
Connectivity Index
Counties with no airport/highwayRefRef
Counties crossed by a highway1.10(1.001.20)
Counties next to airport1.13(1.031.24)
Counties with an airport1.31(1.141.51)
Table 3. Ten highest Coronavirus Disease 2019 (COVID-19) mortality risk areas in highly populated counties (4th quartile).
Table 3. Ten highest Coronavirus Disease 2019 (COVID-19) mortality risk areas in highly populated counties (4th quartile).
LocationObserved CountsExpected CountsConnectivity
Index
PM25 (u/gml)Poverty (%)
Bronx, NY4912810Next to airport11.729.1
McKinley, NM24341Crossed by a highway3.036
Queens, NY72241295Airport11.213
Kings, NY72901465Next to airport11.521.1
Essex, NJ2116447Airport11.216.4
Passaic, NJ1245284Next to airport9.616.7
Union, NJ1351312Next to airport11.49.8
Richmond, NY1083267Next to airport11.312.8
Hudson, NJ1508377Next to airport12.316.3
Bergen, NJ2035524Next to airport11.37
LocationICU Per 100,000White (%)Black (%)Other Races (%)Latino (%)RR CI: [2.5%, 97.5%]
Bronx, NY19.121.334.144.655.96.07 [5.90, 6.24]
McKinley, NM35.7150.784.314.35.85 [5.13, 6.61]
Queens, NY6.43918.342.7285.58 [5.45, 5.71]
Kings, NY10.843.532.623.919.24.98 [4.86, 5.09]
Essex, NJ28.542.139.818.122.74.74 [4.54, 4.94]
Passaic, NJ10.562.211.426.440.94.39 [4.15, 4.63]
Union, NJ13.956.221.222.631.14.33 [4.11, 4.57]
Richmond, NY15.274.310.215.518.34.06 [3.82, 4.3]
Hudson, NJ13.355.112.432.543.24.0 [3.80, 4.21]
Bergen, NJ13.171.46.022.619.43.88 [3.72, 4.05]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Correa-Agudelo, E.; Mersha, T.B.; Branscum, A.J.; MacKinnon, N.J.; Cuadros, D.F. Identification of Vulnerable Populations and Areas at Higher Risk of COVID-19-Related Mortality during the Early Stage of the Epidemic in the United States. Int. J. Environ. Res. Public Health 2021, 18, 4021. https://doi.org/10.3390/ijerph18084021

AMA Style

Correa-Agudelo E, Mersha TB, Branscum AJ, MacKinnon NJ, Cuadros DF. Identification of Vulnerable Populations and Areas at Higher Risk of COVID-19-Related Mortality during the Early Stage of the Epidemic in the United States. International Journal of Environmental Research and Public Health. 2021; 18(8):4021. https://doi.org/10.3390/ijerph18084021

Chicago/Turabian Style

Correa-Agudelo, Esteban, Tesfaye B. Mersha, Adam J. Branscum, Neil J. MacKinnon, and Diego F. Cuadros. 2021. "Identification of Vulnerable Populations and Areas at Higher Risk of COVID-19-Related Mortality during the Early Stage of the Epidemic in the United States" International Journal of Environmental Research and Public Health 18, no. 8: 4021. https://doi.org/10.3390/ijerph18084021

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop