Abstract
Introduction Social determinants of health (SDoH), defined as nonmedical factors that impact health outcomes, have been associated with fatal crash occurrences. Road users who live in communities with negative SDoH may be at increased risk of crash-related mortality, and the risks may be further heightened among geriatric road users and in rural areas. We evaluated the relationship between the county-level measure of SDoH and county-level fatal crash counts among geriatric and non-geriatric road users living in rural, suburban, and urban areas.
Methods For this ecological study, we pooled data from Fatality Analysis Reporting System (2018 to 2020) and the U.S. Census Bureau (2019 data) and limited our analyses to the 3,108 contiguous US counties. The outcome measures were county-level fatal crash counts involving (1) geriatric (65 years and older) road users (2) non-geriatric road users, and (2) the general population. The predictor variable was the Multidimensional Deprivation Index (MDI), a score that measures the five domains of SDoH - economic quality, healthcare access, education, community, and neighborhood quality. We defined the MDI as a three-level categorical variable: at or below the national average, within two-fold of the national average, and higher than two-fold of the national average. We controlled for county-level demographics and crash characteristics. We performed a Bayesian spatial Poisson regression analysis using Integrated Nested Laplace Approximations and reported the crash fatality rate ratios (plus 95% Credible Intervals (CrI)).
Results The median (Q1, Q3) standardized mortality rate ratios among geriatric and non-geriatric road users were 1.3 (0.6, 2.5) and 1.6 (0.9, 2.7), respectively. A total of 283 (9.1%) and 806 (15.9%) counties were classified as very highly deprived and highly deprived, respectively. Clusters of counties with high deprivation rates were identified in the Southern states. Counties classified as very highly deprived and highly deprived had 40% (95% CrI: 1.24 – 1.57) and 25% (95% CrI: 1.17 – 1.34) increased geriatric fatality crash rate ratios and this pattern of association persisted in suburban and urban areas. Also, counties classified as very highly deprived and highly deprived had 42% (95% CI: 1.27 – 1.58) and 32% (95% CI: 1.23 – 1.38) increased fatality crash rate ratios among all road users and this pattern persisted in suburban and urban areas. Counties with more than four-fold increased fatality rate ratios were located commonly in Texas, Oklahoma, Nevada, and Utah.
Conclusion Despite older adults being less frequent road users, county-level deprivation measures of the SDoH are equally associated with geriatric and non-geriatric crash-related fatal rate ratios. Policies that improve county-level SDoH may reduce the county-level fatal rate ratios equally among geriatric and non-geriatric road users.
Introduction
Social determinants of health are conditions located where people live that affect people’s quality of life and health outcomes.1-3 These conditions are commonly grouped into five interrelated domains: (1) economic stability, (2) education access and quality, (3) healthcare access and quality, (4) neighborhood and built environment, and (5) community and social context.3 Each of these wide-ranging conditions has been associated with individual and community-level health disparity and inequities and may serve as predictors of increased susceptibility to acute and chronic disease morbidity and mortality. Poverty, an exemplar of economic stability, predicts substance use, injury, and adverse events from chronic diseases.4 Similarly, rurality, an exemplar of neighborhood and built environment is associated with increased hospital closure,5-7 longer emergency medical service response and hospital transit time,8-11, and increased deaths without resuscitation.10
Understanding the relationship between SDoH and fatal crash counts will further highlight the importance of improving measures of SDoH. The literature on geriatric fatal crash injuries is sparse and very few studies assess the relationship between SDoH and crash injuries.12-15 Yet, SDoH may provide insight into how the conditions surrounding the place we live influence the occurrence of crash injuries. This study, therefore, fills the gap in the SDoH and crash injury literature using three sets of objectives. First, we identified counties with high clusters of negative SDoH and fatal crash rates in the geriatric and non-geriatric populations. We hypothesize that there will be a heterogenous distribution of SDoH and fatal crash rates across the US. Secondly, we assessed the relationship between SDoH and fatal crash counts across the geriatric and non-geriatric populations. We hypothesized that county-level measure of SDoH will be associated with county-level fatal crash counts and the pattern of association will not differ among the geriatric and non-geriatric populations. Lastly, we hypothesize that the effect size of the association between SDoH and geriatric and non-geriatric fatal crash counts will be greater in rural areas compared to urban areas.
Methods
Study Design and Population
For this ecological study, we employed a cross-sectional design to assess the association between county-level measures of fatal crash counts and SDoH. Our unit of analysis was at the county-level
Outcome Variable
The outcome variable is county-level fatal crash counts. We extracted data on fatal crash counts from the Fatality Analysis Reporting System (FARS), one of the national crash databases of the National Highway Traffic Safety Administration.16, 17 The FARS dataset is a census of crash events during which one or more persons died either at the scene or within 30 days of the index crash event.16 The FARS data consists of multiple linkable files that provide information on the accident, person, vehicle, and risk factor-related information. Data in the FARS database is provided either at the individual or crash event levels.
For this study, we defined three sets of county-level fatal crash deaths - across all road users, geriatric (aged 65 years and older), and non-geriatric road users. We defined these three sets of county-level fatal crash counts in three steps. First, from the person file of the FARS data, we created three datasets, restricting each by age (dataset 1: all ages, dataset 2: 65 years and older, dataset 3: less than 65 years).
Next, we restricted each dataset to persons whose injury status was recorded as dead. Lastly, we aggregated the death counts in each dataset by county using the five-digit Federal Information Processing System (FIPS) code. We limited the counties to the 3,108 contiguous US counties, excluding Alaska and the Virgin Islands.
Predictor Variable
The predictor variable is the multidimensional deprivation index (MDI). The MDI is a measure of SDoH developed by the US Census Bureau.18 It consists of six dimensions: standard of living, health, education, economic security, housing quality, and neighborhood quality (Report, A.C.S., Multidimensional Deprivation in the United States: 2017. 2019.).19 Each of the six dimensions is measured at the individual level and has specific defining criteria.19 An individual is considered deprived if at least two criteria are met. The ACS published the national and county-level deprivation index,18 constructed from individual-level measures using the Alkire Foster method.20 Additionally, the ACS defined MDI was defined as a five-level categorical variable defined as county deprivation - very low (i.e. less than half of the national rate), low (ranged from just above half of the national rate and below the lower limit of the 90% confidence interval), at national level (within 90% confidence interval of the national rate), high (above the upper limit of the 90% confidence interval and just below twice the national rate), and very high (at least two times the national rate).18 For this study, we defined MDI as a continuous measure and a three-level categorical variable - very low, low, at the national average, or higher.
Control Variable
We selected county-level sociodemographic and crash characteristics based on prior literature. The county-level sociodemographic and health characteristics included proportions of males, Blacks, Hispanics, the poor, those with bachelor’s degrees, and excessive alcohol intake and emergency department utilization rate. These variables, excluding excessive alcohol intake and emergency department utilization rate, were extracted from American Community Survey.21 Excessive alcohol intake, defined as the proportion of adults with self-reported binge or heavy drinking, was extracted from the County Health Rankings and Roadmaps.22, 23 The emergency department utilization per 1000 Medicare beneficiaries was extracted from the Agency for Healthcare Research and Quality.24 The county-level crash characteristics include the proportions of fatal crash events that occurred at night, during the rush hour period, and the proportion of crash victims screened for alcohol and drugs. Night driving was defined as crash events occurring between 12 midnight and 6 am. Rush hour driving was defined as crash events occurring between 6 and 9 am, and 3 and 7 pm.25 Persons with positive alcohol use were defined as those with a positive blood alcohol level of 0.08% or higher or whose alcohol use status was either self-reported or officer-reported. Persons with positive drug use were defined as those with positive drug tests or whose drug status was either self-reported or officer-reported. All crash characteristics were extracted from FARS.
Stratification
We stratified the counties by rurality/urbanicity using the 2010 Rural-Urban Commuting Area (RUCA) codes (ref).26 The RUCA code is a 10-level classification based on density, urbanization, and daily commuting. We defined rurality/urbanicity in three levels: metropolitan-urban (RUCA codes 1-3), micropolitan-urban (RUCA codes 4-7), and rural and small towns (RUCA codes 8-10).
Spatial Weight Matrix
Using the Euclidean distance across the 3,108 contiguous US counties, we created a spatial weight matrix to map the spatial relationship across the counties. We used the inverse-weighted interpolation, which defines weight as a decreasing function of distance, to generate the spatial weight and account for the presence of spatial autocorrelation in crash occurrences.27
Data Analysis
We reported the proportions of county-level fatal crash injury metrics (counts, case fatality rates (CFR), and standardized mortality ratio (SMR)) across all road users and among geriatric and non-geriatric road users in metropolitan-urban, micropolitan-urban, and rural areas. We reported the frequency distributions, mean (and standard deviation (SD)), and median (and first and third quartiles (Q1, Q3)) for all the predictor and control variables as appropriate. We visualized the distribution of the county-level crude fatality rates and MDI values on choropleths. We assessed the construct validity of the MDI across the crash, health, and population characteristics using correlation analysis. MDI was adjudged to correlate with variables across the domains of SDoH if the correlation coefficient was >0.3. For the regression models, variables with correlation coefficients of >0.3 were deemed to exhibit substantial correlation with MDI and were removed from the models.
Since our outcome variable is a count measure, we performed model testing to compare the Poisson, negative binomial regression, and zero-inflated Poisson and negative binomial regression models, using a spatial Bayesian paradigm with non-informative priors. The offset variable was the county-level expected count, defined as county population multiplied by the ratio of the sum of fatal counts and the sum of the population across the 3108 counties. Consistent with our definition of the outcome variable, we generated three different county population sizes (all ages, geriatric, and non-geriatric), and computed county-level expected counts for all road users, and geriatric and non-geriatric road users. We reported the Akaike Information Criteria and Bayesian Information Criteria.
After performing model fitting and selecting the Poisson model, we performed a spatial Bayesian Poisson regression analysis. For the actual fatal crash counts Y(si),we assumed a Poisson distribution in the data. The mortality rate ratio {λ(Si),i = 1,…,n} over the n = 3, 108 counties were modeled by county-level K independent variables Xk(Si),k = 1, …, K,, a set of regression coefficients βk,K = 0,1, …, K,, spatial random effects terms {V(Si),i = 1, …, n}, and a error components {ε(Si), i = 1, …, n} that follows independent and identical zero-mean Gaussian distributions with variance parameter σ2. The regression equation takes the form:
Thereafter, we reported the unadjusted and adjusted mortality rate ratios (incidence rate ratios) and reported the 95% credible interval. For each model, we computed the mortality rate ratios for all road users as well as geriatric and non-geriatric road users. We computed the predicted rate ratios for each county and visualized the distribution on choropleths. Additionally, we performed stratified analysis and for each of the three outcome measures, we generated mortality rate ratios by rurality/urbanicity.
Results
Across the 3,108 counties, the median (Q1, Q3) MDI rate was 13.0 (9.4, 18.1) (Table 1). The proportion of counties classified as high and very highly deprived were 9% and 26%, respectively. The mean (SD) crash case fatality rate across all age groups was 44.7% (15.1), with the rates lower among the geriatric population (9.1% (8.8)) compared to the non-geriatric population (35.6% (1.0)). The median (Q1, Q3) county SMR across all age groups was 1.6 (1.0, 2.7), with the ratio marginally lower among the geriatric 1.3 (0.6, 2.5) population compared to the non-geriatric population 1.6 (0.9, 2.7). The mean county proportion of males was 50% (2.2) while the median county proportions of non-Hispanic Blacks and Hispanics were 2.6% and 4.5%, respectively. The mean county proportion of the poor classified as poor was 14.5% while those with bachelor’s degree was 19.6%. Also, the mean county proportions of alcohol intake and ED utilization rates were 0.2% and 5.6%, respectively. Approximately a quarter of crashes occurred at night and 43% occur during the rush hour period. Furthermore, mean alcohol and drug screening rates across the counties were 43% and 39%, respectively. Counties with very high deprivation index were found in South Dakota, Arizona, New Mexico, and in some southern region states - Texas, Louisiana, Mississippi, Alabama, Georgia, South Carolina, and Kentucky (Figure 1a). County clusters of high deprivation index were similarly identified in South Dakota, Arizona, New Mexico, and in all southern US states (Figure 1b). Counties with standardized mortality rates of 10 or higher were identified in Montana, North and South Dakota, Wyoming, Colorado, New Mexico, Texas, Nevada, California, and Oregon (Figure 2a). This pattern of high standardized mortality rate was consistent among geriatric and non-geriatric road users (Figure 2b and 2c).
County-level frequency distribution and summary statistics of multidimensional deprivation index, crash, and sociodemographic characteristics by rurality/urbanicity (N=3,108)
Pairwise correlation coefficient assessing the relationship between county-level multidimensional deprivation index, sociodemographic characteristics and fatal crash events
A. U.S. County distribution of the Multidimensional Deprivation Index B. Cluster identification of multidimensional deprivation index across U.S. counties
County-level distribution of the motor vehicle crash standardized mortality rates among (A) road users of all ages (B) geriatric road users, and (C) non-geriatric road users.
We evaluated the construct validity of the MDI across the five constructs of the SDoH. In the economic stability domain, MDI exhibited a statistically significant strong positive correlation with the county proportions of those classified as poor (r=0.72, p<0.01). In the education access and quality domain, MDI exhibited statistically significant moderate negative correlations with county proportions of those with bachelor’s degrees (r= -0.38, p<0.01). In the healthcare access and quality domain, MDI exhibited moderate positive correlations with county-level ED utilization rate (r= 0.38, p<0.01). In the neighborhood and built environment, MDI exhibited a statistically significant mild positive correlation with rurality (r=0.14; p<0.01). In the community and social context domain, MDI exhibited a statistically significant moderate positive correlation with county proportions of Blacks (r=0.49, p<0.01) and a moderate negative correlation with county proportions of excessive alcohol intake (r=-0.48, p<0.01).
The result of the model testing showed that the spatial Poisson regression was the most parsimonious model compared to the negative binomial, zero-inflated negative binomial, and zero-inflated Poisson models (Table 3). Across the entire population and separately among the geriatric and non-geriatric populations, the spatial Poisson model had the lowest AIC and DIC.
Summary of Model Testing
Across all age groups, a unit increase in MDI was associated with a 1% increase in crash fatality rates (Table 4). Also, those living in very highly deprived and highly deprived counties had 33% (95% CrI: 1.20 - 1.48) and 25% (95% CrI: 1.17 - 1.32) increased crash mortality rate ratio compared to those living in counties with normal MDI rate. This pattern of association was consistent among geriatric and non-geriatric road users. Several states in the West and South of the US had four-fold crash mortality rate ratios (Figure 3a). Among geriatric road users, states with several counties with a four-fold fatal mortality rate ratio were found in Oklahoma, Kansas, Texas, Utah, Nevada, California, and Oregon (Figure 3b). Among non-geriatric road users, states with higher than four-fold mortality rate ratio were found in almost all states in the West and South US regions as well as in North and South Dakota, Nebraska, and Kansas states (Figure 3c).
Univariate association between county-level crash fatality and county multidimensional deprivation index, socioeconomic and crash characteristics
Incidence rate ratios of crash fatalities involving older adults, teenagers, and all age groups across U.S. counties with varying levels of deprivation (N=3,108)
County-level distribution of adjusted crash mortality rates among (A) road users of all ages (B) geriatric road users, and (C) non-geriatric road users
Discussion
Earlier studies have reported that social determinants of health are associated with cardiovascular,28-30 diabetes,31, 32 cancer,33, 34 and mental health-related morbidities and mortality.35, 36 Our study adds to the extant literature by reporting that social determinants of health are associated with fatal crash injuries among geriatric and non-geriatric road users. Our study validated that the multidimensional deprivation index is a composite measure that can be used to rank social determinants of health at the county level and we report that about a third of US counties had high to very high multidimensional deprivation indices. Rural-urban disparities exist in both the multidimensional deprivation index and fatal crash rates with rural areas having a larger proportion of counties with high and very high deprivation indices. Also, the geriatric and non-geriatric case fatality rates and SMR were highest in rural areas and lowest in metropolitan-urban areas. Counties with high and very high deprivation indices have significantly elevated fatal crash rate ratios and the ratios are not substantially different among geriatric and non-geriatric road users. Among geriatric road users, rural and micropolitan urban counties with high and very high deprivation indices have significantly elevated fatal crash rate ratios. However, among the non-geriatric road users, micropolitan-urban and metropolitan-urban counties with high and very high deprivation indices have significantly elevated fatal crash rate ratios. Lastly, we identified county hotspots of high deprivation indices as well as counties with more than four-fold increased fatality rate ratios in Texas, Oklahoma, Nevada, and Utah.
Although geriatric road users have higher case fatality rates and standardized mortality rates compared to non-geriatric road users, county-level measures of social determinants of health affect geriatric and non-geriatric road users equally. The elevated geriatric case fatality rate and standardized mortality rates may be a reflection of several factors such as the presence of co-morbidities,37-40, frailty,41-46, and under-appreciation and under-triage of the injury severity in the geriatric trauma population.47-51 However, the comparable effect of the multidimensional deprivation index on fatal crash injuries on geriatric and non-geriatric road users underscores the extent to which social determinants of health affect all road users, irrespective of age. This finding indicates that interventions targeting social determinants could be beneficial for all road users, regardless of age. Efforts to improve socioeconomic conditions, enhance educational opportunities, ensure equitable access to healthcare services, and create safe and supportive neighborhoods will, therefore, have a positive effect on the health and safety outcomes of both older adults and younger individuals. These findings emphasize the need for comprehensive approaches that address social determinants of health as a means of promoting road safety and well-being for all age groups.
We report that geriatric road users in rural counties have significantly higher fatality rate ratios while those geriatric road users in metropolitan-urban counties are at no significantly elevated crash fatality rate ratio. The significantly higher fatality rate ratios observed in rural counties suggest that there may be specific challenges unique to rural areas. Earlier studies have reported that residents in rural areas have longer emergency medical service (EMS) response times,8, 10, 11 limited access to level I or II trauma centers,52-54 and a higher proportion of deaths at the crash scene. In contrast, the lack of a significantly elevated crash fatality rate ratio among geriatric road users in metropolitan-urban counties suggests that factors such as better access to emergency medical services and more developed transportation networks, including helicopter EMS access, may contribute to improved injury outcomes for older adults in these areas. Addressing social determinants of health in rural areas may reduce fatal crash occurrences, especially among geriatric road users. Strategies to improve access to healthcare services, enhance transportation infrastructure, and create age-friendly environments can play a vital role in promoting road safety and protecting the well-being of geriatric road users. There is a need for collaborative efforts among policymakers, healthcare professionals, transportation agencies, and community stakeholders to develop and implement comprehensive interventions that address social determinants of health and improve road safety outcomes in rural areas.
Non-geriatric road users in metropolitan counties have significantly higher fatality rate ratios while those non-geriatric road users in rural counties are at no significantly elevated crash fatality rate ratio. This observed finding may suggest the disproportionate distribution of non-geriatric road users in metropolitan urban areas and the impact of risky driving behaviors associated with young and middle-aged drivers. Earlier studies have reported that non-geriatric road users are more likely to drive while under the influence of drugs and alcohol,9, 55, 56 and engage in phone-related distracted driving.57-60 These unique characteristics may account for the observed elevated fatal crash rate ratio in non-geriatric road users.
We identified fatality rate ratios in several counties located in Texas, Oklahoma, Nevada, and Utah. These states also have county hotspots of very high multidimensional deprivation index. Earlier studies have reported that these states have a high proportion of counties with prolonged EMS response time, county hotspots of rush hour-related fatal crash rates, and fatal crash rates due to non-use of seatbelt.8, 14, 61, 62 Identification of counties with high crash fatality rates and with high multidimensional deprivation index allows for focused interventions. Using spatial models and cluster analyses for fatal crash injury assessment can complement current efforts to address predictors of fatal crash injuries using the Motor Vehicle Prioritizing Interventions and Cost Calculator for States (MV PICCS).63, 64 Using the MV PICCS to identify which specific domain of social determinants of health will be the most cost-effective in reducing fatal crash injury at the state level and using spatial and spatiotemporal models to identify counties within the states that should receive the highest priority may be an effective and efficient way of addressing the social determinants of fatal crash injuries at the state and county-levels.
This study has its limitations. The ecological nature of the study makes it impossible to make causal inferences. We pooled our data from FARS, a data repository that relies on state-level reporting of crash events. Data entry and processing errors, as well as inconsistent reporting of crash and crash-related events across the US, cannot be eliminated. The MDI was computed from multiple data sources. Data errors intrinsic to each of the data are inherently transferred into the MDI computation and misclassification bias cannot be eliminated. Despite these limitations, this study is one of the few studies that assessed the association between social determinants of health and fatal crash rates among geriatric and non-geriatric road users. By identifying counties with high rates of fatal crash injury occurrences as those with a high level of deprivation, this study provides information that can guide targeted interventions to counties and states in greatest need. Additionally, our study provides information that can inform policy and resource allocation on improving the social determinants of health and preventing fatal crash occurrences in the presence of other competing community needs.
Conclusion
Despite older adults being less frequent road users, geriatric road users have higher case fatality rates and standardized mortality rates. County-level deprivation measures of the SDoH are equally associated with geriatric and non-geriatric crash-related fatal rate ratios. Policies and interventions that improve county-level SDoH may reduce the county-level fatal rate ratios equally among geriatric and non-geriatric road users.
Data Availability
All data produced in the present study are available upon reasonable request to the authors