ORIGINAL RESEARCH article

Front. Public Health, 18 June 2025

Sec. Life-Course Epidemiology and Social Inequalities in Health

Volume 13 - 2025 | https://doi.org/10.3389/fpubh.2025.1562610

This article is part of the Research TopicPublic Health Outcomes: The Role of Social Security Systems in Improving Residents' Health WelfareView all 63 articles

Social determinants of health and rehabilitation service areas: an urban and rural mediation analysis

  • 1Department of Occupational Therapy, The University of Texas Health Science Center at San Antonio, San Antonio, TX, United States
  • 2Department of Population Health Sciences, The University of Texas Health Science Center at San Antonio, San Antonio, TX, United States
  • 3Department of Health Studies, Sul Ross State University, Alpine, TX, United States
  • 4University of Texas Medical Branch at Galveston, Galveston, TX, United States

Introduction: Social determinants of health significantly shape community discharge (CD) rates in post-acute rehabilitation settings. Additionally, healthcare disparities between urban and rural regions in the United States can affect these discharge rates. These disparities underscore the critical need to understand how social, economic, educational, and healthcare-related factors influence community discharge outcomes to guide equitable healthcare strategies.

Methods: This observational, cross-sectional study analyzed 40,476 ZIP code tabulation area (ZCTA)-level data points linked to rehabilitation service areas and Agency for Healthcare Research and Quality datasets. Exploratory and confirmatory factor analyses identified five social determinants of health domains—social, economic, educational, physical infrastructure, and healthcare—which were assessed using structural equation modeling to evaluate their direct and mediated effects on community discharge rates.

Results: Significant disparities in community discharge rates were observed across urban and rural areas. Urban areas exhibited lower community discharge rates, influenced by higher social and economic deprivation and limited English proficiency. Conversely, rural areas demonstrated higher rates, attributed to areal social, economic, and education characteristics. Key factors affecting community discharge outcomes included economic inequities, limited healthcare access, and transportation barriers.

Conclusion: Targeted interventions addressing economic inequities, healthcare access, and transportation challenges are essential to improving community discharge outcomes. These findings inform policy and healthcare practices aimed at fostering equitable rehabilitation services and optimizing community reintegration.

Introduction

Social determinants of health (SDH) are main drivers of inequities in healthcare. Defined as the conditions and contexts where people live and work, these determinants represent a complex interaction of social structures, economic systems, and physical environments. Federal and international policymakers and agencies have reported on the importance of SDH as powerful contributors to healthcare utilization and quality of care outcomes (1, 2). SDH are reported to account for 30%−55% of health outcomes (3) and are associated with higher healthcare expenditures (4).

Even though SDH represent a large portion of health outcomes and account for considerable costs, there is no agreed upon theoretical model for SDH. The Dahlgren and Whitehead model of population health visually shows the inter-related layers of social determinant domains on health which has been utilized internationally (5). A common classification scheme by policy makers characterizes factors into categories or domains of social determinants. For example, the Agency for Healthcare Research and Quality (AHRQ) breaks SDH into five domains: economic, education, physical environment, social context, and healthcare (6). Healthy People 2030 domains are similar: economic stability, education access and quality, neighborhood and built environment, social and community context, and healthcare access and quality (7). Additionally, multiple indices exist that combine different SDH measures. Some of these are the Area Deprivation Index (8), the Distressed Community Index (9), the Social Deprivation Index (10), and the Social Vulnerability Index (11). These indices are calculated in various ways with different purposes in mind resulting in the healthcare context domain being largely absent. A review of existing indices is beyond the scope of this paper. These challenges have resulted in limitations in SDH research as it relates to the healthcare domain and healthcare outcomes.

There is growing recognition that both individual and community level social determinates impact health outcomes across healthcare sectors and systems (1216). One healthcare sphere where community context plays a role in outcomes is post-acute rehabilitation. Post-acute rehabilitation is an important component of the healthcare system providing services for a myriad of health-related issues. In 2022, post-acute care accounted for $57.8 billion in health care spending across rehabilitation service venues of long-term care hospitals, inpatient rehabilitation facilities, skilled nursing facilities, and home health agencies which each operate under different admission and payment structures (17). Moreover, substantial geographic variation in the distribution of these service venues exists which influences rehabilitation outcomes (1821).

The National Institutes of Health report, Moving Rehabilitation Research Forward discusses the complex interaction between the environment and health outcomes suggesting the need for refined approaches and policies to improve care (22). Like other sectors of healthcare which have noted disparities linked to social determinants (23), rehabilitation research has found disparities in functional status and hospital readmission connected to the SDH domains: socioeconomic, education, and neighborhood quality (2429).

Rehabilitation service areas (RSA) were constructed to aid in the examination and understanding of geographic differences in care and subsequent rehabilitation outcomes. Like hospital and primary care service areas, RSAs are post-acute care utilization areas for rehabilitation service use in the US. There are 1,711 RSAs nationally built from ZIP code tabulation areas (ZCTA) based on patterns of rehabilitation care (30, 31). RSAs are a tool to describe and analyze resource availability utilization, care patterns, and gaps in rehabilitation care.

A prominent rehabilitation goal is to help return someone to their home with the ability to function independently in the community. This goal has driven the development of successful community discharge (CD) as a quality measure for post-acute care rehabilitation. The Centers for Medicare and Medicaid Services defined community discharge as remaining in the community for 31 days without an unplanned hospital readmission. The recognition of successful community discharge as a quality measure for post-acute rehabilitation is inherently linked to environmental factors (32, 33) and likely area level social determinants of health. Moreover rural and urban differences in rehabilitation use suggest potential disparities (34) and the influence of SDH. In a study of rehabilitation use between urban and rural areas, researchers found significant differences in the use of skilled nursing and home health services (35). The impact of area level social determinant on rural and urban rehabilitation use remains understudied. Although there is growing rehabilitation research on SDH, many of these studies evaluate select variables, are often limited by diagnosis, and do not evaluate the panoply of area characteristics and resources.

This study aims to bridge critical gaps in understanding the influence of area-level SDH on community discharge outcomes within post-acute care rehabilitation settings. By leveraging comprehensive datasets at the ZCTA level and integrating RSA-specific service patterns, this research identifies how community-level factors intersect with healthcare access and utilization. Such analyses are pivotal for developing tailored interventions to enhance post-acute care outcomes and reduce disparities in rehabilitation efficacy. Ultimately, this work seeks to inform policymakers and healthcare stakeholders about strategies to optimize community reintegration, advancing equity and quality in healthcare delivery systems.

Methods

Study design, setting, and data sources

This was an observational, cross-sectional secondary data analysis of US data from two main sources. First, we used the publicly available Social Determinants of Health Database from the Agency for Healthcare Research and Quality (AHRQ) (6). This database was created to facilitate research and ultimately improve health outcomes. It provides a single file with a range of variables across several domains that can easily be linked to other data (6). We used the 2016 ZIP Code Data, which includes 5-year estimates from the 2012–2016 American Community Survey (ACS) at the ZCTA level.

Second, we used a crosswalk dataset that contained three variables: (1) rehabilitation service areas (RSAs) (30), (2) ZCTA, and (3) community discharge rates per RSA. This dataset contains 1,711 observations, corresponding to the 1,711 RSAs across the United States and is publicly available for use by all researchers (36). Whereas, RSAs are geographic boundaries similar in concept to Primary Care Service Areas and Hospital Service Areas, RSAs represent the unique patterns of post-acute care services that other small-area boundaries do not capture. We merged the crosswalk dataset and the AHRQ dataset through ZCTA variables and extracted SDH and community discharge rate from the AHRQ dataset. The RSA and community discharge data were developed using 2013–2015 national Medicare claims data, which temporally aligns with the 2012–2016 American Community Survey period used in the AHRQ SDH dataset.

Variables

Social determinants of health

The AHRQ SDH Database 2016 ZCTA Data contains 313 total variables that includes the following contexts form the ACS: social, economic, education, physical infrastructure, and healthcare. The 2016 AHRQ Codebook indicates the following number of variables per domain: social (123 variables), economic (80 variables), education (10 variables), physical infrastructure (69 variables), and healthcare (31 variables) (6).

From these variables, an expert panel experienced with older adult health outcomes selected variables thought to affect successful community discharge. This expert panel included professionals from demography, population health, data science, rehabilitation science, and a rehabilitation clinician. The panel arrived at 34 variables present in the 2016 data at the ZCTA level. The 34 ACS variables finally selected represent the percentage at the ZCTA levels (Supplementary Appendix S1). By domain, the following variables remained: social (10 variables), economic (six variables), education (three variables), physical infrastructure (nine variables), and healthcare (six variables).

A multidisciplinary expert panel consisting of professionals in demography, rehabilitation science, data science, and public health convened weekly over a 6-week period to review and select variables from the AHRQ SDH dataset. Variables were chosen based on two primary criteria: (1) percentage-based metrics to ensure comparability across ZCTAs, and (2) conceptual relevance to one of the five major SDH domains—social, economic, education, physical infrastructure, and healthcare. Although a formal Delphi process was not used to identify variables the expert panel engaged in structured discussions of SDH domains to ensure content validity and representation.

Community discharge

The community discharge (CD) rate is a metric designed to assess the rate at which patients are discharged to and remain in the community for 31 days without an unplanned readmission to the hospital. This metric is consistent with CMS measurement specifications for post-acute care rehabilitation (37, 38). In this study, the CD rate variable in the crosswalk dataset was merged into the AHRQ dataset by the Zip Code Tabulation Area (ZCTA) variable.

Rural-Urban commuting area (RUCA)

The Rural-Urban commuting area (RUCA) system classifies areas based on commuting patterns and population density. The RUCA system divides areas into four categories based on proposed guidelines (39). Urban core areas are highly urbanized, characterized by high population density and excellent access to resources and services. Suburban areas, located on the outskirts of cities, have dense populations and strong commuting connections to urban cores. Large rural areas are rural regions with relatively higher populations and moderate commuting ties to urban centers, whereas small town/rural areas include small communities and remote rural regions with limited commuting flows and restricted access to resources.

These classifications are well-suited for explaining variations in health outcomes and service accessibility across geographic areas (40). In this study, RUCA was particularly applied to analyze the mediating effect of service accessibility on the relationship between the five domains of SDH and successful community discharge. By leveraging RUCA, the study identified how regional characteristics indirectly influence rehabilitation outcomes and successful community reintegration, mediated by SDH domains such as economic, social, educational, physical, and healthcare contexts.

Study size

We utilized a merged dataset derived from the crosswalk dataset and the AHRQ dataset, using the ZCTA variable as the key for merging. The original AHRQ dataset consisted of 40,824 observations. Through the ZCTA variable, the RSA, and CD rate variables from the crosswalk dataset were incorporated into the AHRQ dataset. During this process, 348 ZCTAs were excluded due to missing a RUCA classification or values for one or more of the 34 selected SDH variables. In total, 348 responses (0.85% of the original dataset) were excluded, resulting in 40,476 ZCTAs included in the final analysis. Moreover, we conducted a sensitivity analysis using inclusion status as the outcome and SDH variables as predictors. None of the variables were significantly associated with inclusion, suggesting that the exclusion process did not introduce meaningful bias across the SDH characteristics examined.

Statistical analysis

Exploratory factor analysis (EFA)

In this study, we conducted an exploratory factor analysis (EFA) using 34 ACS variables initially selected to explore the latent factor structure among variables. All analyses were conducted using Mplus software version 8.4, applying Maximum Likelihood estimation with Robust standard errors (MLR) to address data non-normality (41). The MLR estimator also accommodates missing data through full information maximum likelihood under the assumption of data missing at random, allowing all available information to be utilized in the estimation process. The scree plot test was used to determine the number of variables by identifying the elbow point on the graph for the factors with eigenvalues >1. These methods are widely used to determine the number of factors, and to select the factors that explain the most variance in the data. In the analysis process, Mplus's default Geomin rotation (oblique rotation) was applied to allow for correlations between factors. The Geomin rotation was chosen over alternatives like Promax to align with expected interdependence among the latent variables in the dataset. Additionally, only items with factor loadings of 0.3 or higher were included in the subsequent confirmatory factor analysis (CFA), based on previous research suggesting that this threshold is sufficient. While some recommend a higher threshold such as 0.40 for large samples, the 0.30 threshold was deemed appropriate to capture conceptually relevant items given the multidimensional nature of the SDH constructs (42, 43).

Confirmatory factor analysis (CFA)

In this study, CFA was conducted to validate the factor structure derived from the EFA. The variance of the latent variables was fixed at one to ensure model identification, adding constraints for proper model specification. Correlations between factors were freely estimated to interpret factor loadings by fixing the relative size of each factor and to evaluate the relationships between factors. Similar to the EFA, only items with factor loadings of 0.3 or higher were included in the final analysis (43). Modification indices were reviewed to identify potential areas of misfit, but modifications were only made when they were theoretically justified to maintain the integrity of the model. In addition, standardized residuals were inspected to confirm that no large discrepancies existed between observed and model-implied covariances.

Structural equation modeling (SEM) regression

In this study, structural equation modeling (SEM) regression was conducted using the final variables selected through CFA to investigate the relationships between variables. The SEM analysis employed a multiple-factor model to evaluate the relationships among latent variables. The key latent variables included in the analysis were social factors, economic factors, education factors, physical infrastructure factors, and healthcare factors. We analyzed their effects on the dependent variable, CD rate (44). Bootstrapping methods were applied to estimate the confidence intervals for indirect effects, providing more robust inferences regarding mediation effects. Additionally, modification indices were reviewed to address potential areas of misfit, and model refinements were made iteratively where theoretically justified.

SEM mediation

In this study, mediation analysis using SEM was conducted to evaluate the indirect effects between the independent and dependent variables, employing the final variables selected through CFA. The SEM mediation analysis assessed both the direct and indirect effects of social factors, economic factors, education factors, physical infrastructure factors, and healthcare factors on the dependent variable, CD rate, through the mediator variable, RUCA. Specifically, RUCA served as a mediator in the relationship between the independent variables and the dependent variable, analyzing how the independent variables indirectly influenced CD rates. In the mediation model, RUCA classification was treated as an observed categorical mediator representing four distinct geographic contexts (urban core, suburban, large rural, and small town/rural).

The use of RUCA classification as a mediator is grounded in theories of urban health ecology and spatial health inequality, which posits that geographic context—particularly the degree of urbanization—influences how social and structural determinants of health translating into healthcare access and outcomes (45). Urban–rural distinctions shape service availability, care-seeking behaviors, and environmental constraints, positioning the RUCA classification as a conceptually grounded mediator through which area-level social determinants of health influence community discharge outcomes —consistent with theories of urban health ecology and spatial health inequality (46).

The significance of the mediation effects was tested using the bootstrapping method, which involved repeatedly resampling the data to estimate the confidence intervals of the indirect effects. This method enhances reliability by accounting for the uncertainty in the sample distribution (41). The significance of the indirect paths was confirmed through bootstrapping, and the mediation effects were evaluated based on the significant p-values derived from the estimated standard errors (47).

We calculated the indirect effects as the product of the a-path (SDH to RUCA) and the b-path (RUCA to CD). Additionally, to determine the total effect, we summed the direct and indirect effects. Finally, to assess the mediation effect, we divided the indirect effect by the total effect and expressed the extent of mediation as a percentage.

Model fit

In this study, the model fit of each model was assessed using the following indices: root mean square error of approximation (RMSEA), comparative fit index (CFI), Tucker–Lewis index (TLI), and standardized root mean square residual (SRMR). The criteria for each fit index are as follows. First, RMSEA is an index for evaluating model fit, where a value of 0.06 or lower indicates good fit, and values below 0.08 are considered acceptable (48). Second, CFI and TLI values above 0.95 indicate good fit, while values above 0.90 are deemed acceptable (48, 49). Third, SRMR values below 0.08 represent good fit, and this criterion was used to evaluate the SRMR values in this study (48). It is important to note that these fit indices may vary slightly depending on the complexity of the model and the sample size, which was considered in the interpretation of the results.

Results

Exploratory factor analysis (EFA)

As the first step of this study, an EFA was conducted using 34 ACS variables categorized in the initial classification. The results are as follows (Supplementary Appendix S2). The Social factor was divided into two components: the first factor included S5: does not speak English well or at all (λ = 0.860), S9 (λ = 0.913), and S10 (λ = 0.989), while the second factor included S4: single parent families (λ = 0.490), S6 (λ = 0.893), and S7 (λ = 0.969; RMSEA = 0.068, CFI/TLI = 0.909/0.843, and SRMR = 0.053). The economic factor was initially divided into three components, but after excluding duplicated items, it was ultimately classified into two factors. The first factor included EC1: under 3.99 of the poverty threshold (λ = 0.556) and EC3: not in labor force (λ = 0.701), while the second factor included EC2: household income < $24,999 (λ = 0.513), EC4: Asian population below poverty level (λ = 0.384), EC5: Black or African American population below poverty level (λ = 0.432), and EC6: some other race population below poverty level (λ = 0.433; RMSEA = 0.000, CFI/TLI = 1.000/1.000, and SRMR = 0.000). The education factor was grouped into a single factor, including ED1: a bachelor's degree (λ = 0.825), ED2: master's or doctorate (λ = 0.784), and ED3: only high school diploma (λ = −0.840; RMSEA = 0.000, CFI/TLI = 1.000/1.000, and SRMR = 0.000). The physical infrastructure factor was divided into two components: the first factor included PI2: housing units lacking compete kitchen facilities (λ = 0.867) and PI3: housing units lacking plumbing facilities (λ = 0.947), and the second factor included PI5: housing units with no vehicle available (λ = 0.916), PI6: workers taking public transportation, excluding taxicab (λ = 0.664), PI7: works taking taxicab, motorcycle or other means to work (λ = 0.365), PI8: workers walking to work (λ = 0.492), and PI9: workers in households with no vehicle available (λ = 0.960; RMSEA = 0.039, CFI/TLI = 0.950/0.905, and SRMR = 0.026). Finally, the Healthcare factor was also divided into two components: the first factor included H3: Medicare, Medicaid, TRICARE/military, U.S Department of Veterans Affairs coverage (λ = 0.902) and H6: no health insurance coverage (λ = 0.339), while the second factor included H1: other private-only health insurance combinations (λ = 0.611) and H5: TRICARE/military or VA health insurance coverage only (λ = 0.969; RMSE = 0.035, CFI/TLI = 0.937/0.763, and SRMR = 0.017).

Confirmatory factor analysis (CFA)

In this study, CFA was conducted to confirm the validity of the factor structure identified in the EFA, allowing for correlations between the nine factors. The results of the CFA analysis are as follows (Supplementary Appendix S3) For the social factor, the first factor showed factor loadings ranging from 0.859 to 0.985, indicating strong associations, while the second factor had moderate to strong loadings ranging from 0.459 to 0.977. For the economic factor, the first factor displayed loadings between 0.445 and 0.834, and the second factor showed relatively lower loadings between 0.380 and 0.442, suggesting a weaker association. The education factor was represented by a single factor, with loadings of 0.801 and a negative factor loading of −0.794, which may indicate reverse coding or a differing pattern in the data. For the physical infrastructure factor, the first factor exhibited strong loadings between 0.884 and 0.931, and the second factor showed loadings ranging from 0.378 to 0.945. For the healthcare factor, the first factor had loadings between 0.572 and 0.989, indicating variability in factor strength, while the second factor displayed lower loadings from 0.428 to 0.599. Nonetheless, all factor loadings exceeded the threshold of 0.300, and thus all factors were retained for further analysis. The model fit indices for the CFA were as follows: RMSE = 0.058, CFI/TLI = 0.803/0.754, and SRMR = 0.075.

Table 1 presents descriptive statistics for the final items selected through CFA, showing their distribution across different residential regions. The most notable variables by region are as follows: in urban core areas, the most notable characteristic was the higher percentage of foreign-born individuals (S9: 10.66%) compared to other regions. Additionally, the use of public transportation was highest here (PI6: 4.21%). In Suburban areas, the highest proportion of the population under 3.99 of the poverty thresholds was observed (EC1: 72.30%). Furthermore, a significant percentage of the population had only a high school diploma (ED4: 36.92%). In large rural areas, the population without labor force participation was notably high (EC3: 44.01%). This region also had the highest percentage of households lacking kitchen facilities (PI2: 5.82%). In small town/rural areas, the highest percentage of workers walked to work (PI8: 5.85%), and the region also had the greatest proportion of households lacking plumbing facilities (PI3: 7.59%).

Table 1
www.frontiersin.org

Table 1. Social determinants of health demographic characteristics by region variables (N = 40,476).

SEM regression

Table 2 analyzes the relationship between the nine factors and CD rate based on the SEM regression results. First, with the exception of the first factor of physical infrastructure (p = 0.204) and the first factor of healthcare (p = 0.111), all factors showed significant differences in their effects on CD rate. For the social factor, both factors were negatively associated with CD rate, with estimates of −0.052 and −0.029, respectively. This refers to that as the proportions of social variables such as single-parent families, foreign-born individuals, and minority populations increase, the CD rate decreases. The economic factor had contrasting results. The first factor exhibited a strong negative relationship with CD rate (β = −0.261), indicating that increases in variables like the percentage of the population below the poverty threshold and the percentage of the population not in the labor force are associated with a decrease in CD rate. Conversely, the second factor had a positive association (β = 0.093), indicating that as other economic indicators, such as the percentage of the population below the poverty level for specific racial groups, increase, CD rate may increase. The educational factor showed a substantial negative effect on CD rate (β = −0.444). This result points to that as the proportion of the population with higher education levels increases, CD rate tends to decrease. For the physical infrastructure factor, the second component showed a small negative relationship with CD rate (β = −0.031), indicating that increases in variables like the percentage of workers walking to work are associated with a slight reduction in CD rate. Lastly, the second component of the healthcare factor had a positive relationship with CD rate (β = 0.044), suggesting that as the percentage of the population covered by specific health insurance types, such as TRICARE or VA coverage, increases. The model fit indices for the SEM regression analysis were as follows: RMSE = 0.057, CFI/TLI = 0.806/0.754, and SRMR = 0.073.

Table 2
www.frontiersin.org

Table 2. SEM regression model results for effects of social, economic, educational, physical infrastructure, and healthcare factors on community discharge rate.

SEM mediation

Table 3 presents the results of the SEM mediation analysis, which investigated the mediating effect of RUCA on the relationship between the nine factors and CD rate. The diagrams in Figures 1, 2 provide a visual representation of these effects, with the left panels depicting the structural pathways for each factor's mediating impact through RUCA. Complementing these, the right-side graphs illustrate distinct patterns in the relationship outcomes: Figure 1 highlights contrasting values across factors, showcasing opposing trends, while Figure 2 underscores the aligned values, indicating consistent directional relationships among the factors. The study found that, with the exception of the second social factor (p = 0.174), all direct effects between the factors and CD rate showed significant differences. Additionally, significant differences were observed in the relationships between all factors and region, as well as between all regions and CD rate. This indicates that region acts as a partial mediator for all factors except the second social factor. For the first social factor, the indirect effect was −0.005, with a total effect of −0.145, indicating that region accounted for 3.45% of the partial mediation effect (RMSE = 0.059, CFI/TLI = 0.929/0.875, and SRMR = 0.029). For the first economic factor, the indirect effect was −0.007, with a total effect of −0.123, showing that region accounted for 5.38% of the partial mediation effect. The second economic factor showed an indirect effect of 0.007 and a total effect of 0.221, with region mediating 3.07% of the effect (RMSE = 0.024, CFI/TLI = 0.941/0.876, and SRMR = 0.015). For the education factor, the indirect effect was −0.006, with a total effect of −0.276, indicating that region accounted for 2.17% of the partial mediation effect (RMSE = 0.030, CFI/TLI = 0.982/0.956, and SRMR = 0.008). For the first physical infrastructure factor, the indirect effect was 0.011, with a total effect of 0.094, showing that region had the highest partial mediation effect at 11.70%. The second physical infrastructure factor had an indirect effect of −0.067 and a total effect of −0.116, with region accounting for 3.33% of the partial mediation effect (RMSE = 0.041, CFI/TLI = 0.916/0.869, and SRMR = 0.051). In the case of the two healthcare factors, the initial analysis aimed to examine both factors together. However, due to a negative residual variance estimate in the H3 item of the first factor, which caused model fit issues, the first factor was excluded. The final analysis proceeded with only the second healthcare factor, which showed an indirect effect of 0.008 and a total effect of 0.085, with region mediating 9.41% of the effect (RMSE = 0.024, CFI/TLI = 0.973/0.840, SRMR = 0.006).

Table 3
www.frontiersin.org

Table 3. Mediation analysis results of residence area between social determinants of health factor and community discharge.

Figure 1
www.frontiersin.org

Figure 1. Diagram of mediation analysis results (factors with opposite outcomes). p > 0.05*, p > 0.001**, p > 0.0001***.

Figure 2
www.frontiersin.org

Figure 2. Diagram of mediation analysis results (factors with similar outcomes). p > 0.05*, p > 0.001**, p > 0.0001***.

Discussion

The primary aim of this study was to explore the relationship between community-level SDH and CD rates. Through EFA and CFA, we identified five major SDH contexts: social, economic, education, physical infrastructure, and healthcare similar to those used by the WHO and healthy people 2030 (3, 7). These factors were analyzed for their impact on CD rates across different urban/rural categories: urban core area, suburban area, large rural area, and small town/rural area. Our results indicate significant variability in CD rates based on community-level SDH, with factors such as economic deprivation, social inequality, and lack of healthcare access emerging as critical barriers to CD.

The findings from the SEM regression analysis underscore the significant impact of community social structures on CD outcomes. Higher percentages of single-parent families, minority populations, and foreign-born residents were linked to lower CD rates. This aligns with existing research that highlights how social fragmentation and minority status contribute to poorer health outcomes, often through increased stress, reduced access to healthcare, and limited social support (50). Furthermore, the mediation analysis revealed that urban/rural categories partially mediated the relationship between social factors and CD rates. Urban core areas with higher percentages of foreign-born residents and non-citizens showed lower CD rates. This is consistent with prior research demonstrating that geographic disparities in healthcare access disproportionately impact minority communities, exacerbating health disparities (51). Additionally, the complex relationship between language proficiency and health outcomes was evident in the data. Urban areas with a higher proportion of individuals with limited English proficiency exhibited lower CD rates. This supports previous studies showing that language barriers can obstruct effective communication between patients and healthcare providers, ultimately affecting the quality of care and discharge planning (52, 53).

Economic factors play a complex role in influencing CD rates, particularly in the context of poverty and labor force non-participation. According to our SEM regression analysis, increased levels of poverty and labor force non-participation are associated with a significant decrease in CD rates. This finding underscores the connection between economic hardships and limited access to essential healthcare services in general (54, 55), particularly following rehabilitation. This finding also aligns with existing research at the individual level highlighting that populations experiencing economic instability often face barriers to healthcare which leads to poor outcomes (56). In contrast, the mediation analysis indicates distinct trends across urban and rural areas. Specifically, higher poverty rates among minority populations correlate with a decrease in CD rates and are more prevalent in urban settings, suggesting that urban areas face unique challenges related to economic hardship, where healthcare infrastructure and resources are often strained (57, 58). This urban context likely exacerbates the negative impact of economic disadvantages on rehabilitation outcomes, as these communities may struggle with an overcrowded healthcare systems and/or limited access to essential services. Conversely, increased poverty and labor force non-participation in rural areas appear to be associated with an increase in CD rates, suggesting that rural areas may benefit from stronger community networks or targeted support initiatives that help counterbalance economic challenges (59, 60).

The education context also had a significant impact on CD rates. This study found that higher proportions of individuals with advanced education were associated with lower CD rates. Individuals with higher education levels generally have greater access to social and economic resources, and this finding suggests that they may be more likely to utilize healthcare following rehabilitation potentially due to their access to greater resources. Notably, the negative association between higher education levels and CD rates was more pronounced in urban areas. Although higher educational attainment is generally associated with improved health outcomes, prior studies have shown that individuals with greater health literacy and more complex care expectations may be more likely to seek additional services, remain within the healthcare system longer, or delay discharge planning, especially in urban contexts (61, 62). In such settings, well-educated individuals may advocate for extended rehabilitation or pursue more specialized care pathways, potentially leading to lower community discharge rates.

The physical infrastructure context was associated with CD rates. This suggests that communities where residents are reliant on public transportation, regardless of urban/rural status, may see significant challenges in achieving successful reintegration into the community. These findings are consistent with previous research demonstrating the negative impact of limited transportation accessibility on health outcomes (63). For those dependent on public transit, accessing healthcare services post-discharge can be particularly constrained, posing a major barrier to recovery during the rehabilitation process (64). Additionally, the lack of transportation accessibility can diminish patients' ability to live independently, further hindering their reintegration into the community. In the mediation analysis, areas with a higher proportion of individuals who use public transportation or walk to work, typically urban regions, were associated with lower CD rates. This may reflect the challenges of navigating healthcare in urban environments, where access to transportation services plays a crucial role in post-discharge recovery (63). These findings underscore the importance of improving transportation access to promote better health outcomes, particularly in underserved communities (65). Additionally, we found that areas with housing lacking complete kitchen and plumbing in rural areas were associated with an increase in CD rates in the mediation analysis, meaning higher successfully discharged to the community in areas with limited facilities. This is a counterintuitive finding as one might expect that housing unit challenges in rural areas would hinder CD. However, this finding may be explained by the strength of informal caregiving networks and supports that are more prevalent in rural areas. In these rural areas, despite inconsistent housing characteristics, family members or neighbors may provide essential support that facilitates successful reintegration into the community (66). These findings suggest that rural areas may have unique social and healthcare dynamics that enable successful CD even under challenging circumstances.

The healthcare context also played a crucial role in determining CD rates. In this study, a higher proportion of individuals covered by TRICARE or VA health insurance was associated with increased CD rates. This suggests that individuals with military or veteran' insurance may receive more comprehensive and coordinated care that supports successful reintegration into the community. This finding aligns with existing research showing that access to comprehensive healthcare coverage improves health outcomes, particularly in vulnerable populations (67). The mediation analysis further revealed that areas with higher proportions of individuals covered by TRICARE or VA insurance are more likely to be in rural areas, where CD rates tend to be higher. This highlights the importance of healthcare coverage in enhancing discharge outcomes, particularly in rural settings where TRICARE or VA insurance is more prevalent (68). Rural areas may benefit from stronger community networks and tailored healthcare resources, which can positively influence CD rates. These results underscore the necessity of expanding access to comprehensive healthcare coverage, especially in rural regions, to improve rehabilitation outcomes across diverse populations (69).

Although this study is based on U.S.-specific geographic areas, such as rehabilitation service areas (RSAs) and ZIP code tabulation areas (ZCTAs), the underlying methodological approach is transferable. Our framework links area-level social determinants of health to rehabilitation outcomes through geographic mediators. This approach can be applied in other health systems that utilize small-area geographies and population-based data. Countries like the United Kingdom, Canada, and Australia maintain similar administrative structures in healthcare planning and service delivery. In these settings, comparable methods could be used to explore equity in post-acute care access and outcomes.

Limitations

There are certain limitations in this study that warrant attention. First, the merging of the RSA crosswalk dataset with its 1,711 records with over 40,000 AHRQ records, linked using the ZCTA contextual variables, resulted in a nested data structure. However, because the CD rate variable was calculated based on RSAs, we were unable to apply multilevel modeling for CD as the outcome variable. To address this, we applied SEM regression and mediation analysis to ensure analytical rigor despite the limitation. Additionally, our model incorporated five SDH domains and nine factors derived from EFA and CFA, making the analysis inherently complex. This complexity, combined with the inclusion of CD rate as the outcome and RUCA categories as a mediator, led to some model fit indices falling below standard thresholds. However, we mitigated this by refining factor loadings and iterating the model to optimize its fit. Given the theoretical importance of preserving distinct SDH domains, we retained the full nine-factor structure despite moderate CFI and TLI values. Previous research has shown that complex models with large sample sizes and multiple latent constructs often yield lower comparative fit indices, yet still provide valid insights when supported by conceptual coherence and acceptable RMSEA and SRMR values. Since this analysis was conducted on the contextual/ecological level, it is important to not interpret these findings on the individual level. Additionally, although exploration and confirmatory factor analyses were conducted at the ZCTA level (n = 40,476) to leverage the variability and maintain statistical power, the outcome variable, community discharge rate, was aggregated at the RSA level (n = 1,711). This analytical misalignment between the level of measurement for independent and dependent variables may limit the direct interpretability of latent factors at the outcome level.

Conclusion and future directions

Our analysis demonstrated that area-level social determinants of health relate to community discharge outcomes following post-acute rehabilitation. Moreover, these relationships have differing impacts in rural and urban areas. Urban/rural areas mediated the relationship between social, economic, educational, physical infrastructure and healthcare factors and successful community discharge rates. This highlights how geographic context shapes the impact of these factors on rehabilitation outcomes. These finding confirm the multifaceted connection between area-level social determinants of health and community discharge following rehabilitation. Therefore, careful consideration of drivers of community integration and healthcare use in rural and urban areas is warranted. Future studies of rural and urban outcome differences across RSAs could benefit from exploring more advanced analytical techniques, such as multilevel SEM, to improve model precision and better inform public health approaches. Despite these challenges, our study provides valuable insights into how area-level contexts (social, economic, educational, physical infrastructure, and healthcare) influence community discharge outcomes following post-acute care rehabilitation.

Data availability statement

Publicly available datasets were analyzed in this study. This data can be found here: Agency for Healthcare Research and Quality available at: https://www.ahrq.gov/sdoh/data-analytics/sdoh-data.html; Rehabilitation Service Areas (for post-acute care), doi: https://doi.org/10.5281/zenodo.14291764.

Ethics statement

The studies involving humans were approved by University of Texas Medical Branch, and University of Texas Health Science Center at San Antonio. The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation was not required from the participants or the participants' legal guardians/next of kin in accordance with the national legislation and institutional requirements.

Author contributions

SN: Writing – original draft, Writing – review & editing, Data curation, Formal analysis, Methodology, Visualization. SS: Writing – review & editing. JD: Writing – review & editing. AB: Writing – review & editing, Data curation. TR: Writing – review & editing, Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Writing – original draft.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This work was supported by the Agency for Healthcare Research and Quality (R01-HS024711, R01-HS029891), the National Center for Advancing Translational Sciences, the National Institutes of Health through Grant UM1 TR004538, and the National Institute on Aging, San Antonio Claude D. Pepper Older Americans Independence Center (P30AG044271). The authors acknowledge DUA 51261 connected to the data used in this project.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Gen AI was used in the creation of this manuscript.

Publisher's note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpubh.2025.1562610/full#supplementary-material

References

1. Centers for Medicare & Medicaid Services: Office of Minority Health. The CMS Framework for Health Equity (2022-2023) (2022).

Google Scholar

2. Office of the Assistant Secretary for Planning and Evaluation Department of Health Human Services. Second Report to Congress on Social Risk Factors and Performance in Medicare's Value-Based Purchasing Program (2020).

Google Scholar

3. World Health Organization. Social Determinants of Health 2024. Available online at: https://www.who.int/health-topics/universal-health-coverage/social-determinants-of-health#tab=tab_1 (accessed November 4, 2024).

Google Scholar

4. Mohan G, Gaskin DJ. Social determinants of health and US health care expenditures by insurer. JAMA Netw Open. (2024) 7:e2440467-e. doi: 10.1001/jamanetworkopen.2024.40467

PubMed Abstract | Crossref Full Text | Google Scholar

5. Dahlgren G, Whitehead M. The Dahlgren-Whitehead model of health determinants: 30 years on and still chasing rainbows. Public Health. (2021) 199:20–4. doi: 10.1016/j.puhe.2021.08.009

PubMed Abstract | Crossref Full Text | Google Scholar

6. Agency for Healthcare Research and Quality. Social Determinants of Health Database. Rockville, MD: Agency for Healthcare Research and Quality (2023). Available online at: https://www.ahrq.gov/sdoh/data-analytics/sdoh-data.html (accessed November 20, 2024).

Google Scholar

7. Gómez CA, Kleinman DV, Pronk N, Wrenn Gordon GL, Ochiai E, Blakey C, et al. Addressing health equity and social determinants of health through healthy people 2030. J Pub Health Manage Pract. (2021) 27(Supplement 6):S249–57. doi: 10.1097/PHH.0000000000001297

PubMed Abstract | Crossref Full Text | Google Scholar

8. Hu J, Kind AJ, Nerenz D. Area deprivation index predicts readmission risk at an urban teaching hospital. Am J Med Qual. (2018) 33:493–501. doi: 10.1177/1062860617753063

PubMed Abstract | Crossref Full Text | Google Scholar

9. Mehaffey JH, Hawkins RB, Charles EJ, Turrentine FE, Hallowell PT, Friel C, et al. Socioeconomic “Distressed Communities Index” improves surgical risk-adjustment. Ann Surg. (2020) 271:470–4. doi: 10.1097/SLA.0000000000002997

PubMed Abstract | Crossref Full Text | Google Scholar

10. Butler DC, Petterson S, Phillips RL, Bazemore AW. Measures of social deprivation that predict health care access and need within a rational area of primary care service delivery. Health Serv Res. (2013) 48:539–59. doi: 10.1111/j.1475-6773.2012.01449.x

PubMed Abstract | Crossref Full Text | Google Scholar

11. Schmidtlein MC, Deutsch RC, Piegorsch WW, Cutter SL. A sensitivity analysis of the social vulnerability index. Risk Anal. (2008) 28:1099–114. doi: 10.1111/j.1539-6924.2008.01072.x

PubMed Abstract | Crossref Full Text | Google Scholar

12. Reyes A, Prabhakaran D, Banegas MP, Shih JJ, Iragui-Madoz VJ, Almane DN, et al. Individual-and community-level social determinants of health are associated with cognition in older adults with focal epilepsy. Epilepsy Behav. (2024) 158:109927. doi: 10.1016/j.yebeh.2024.109927

PubMed Abstract | Crossref Full Text | Google Scholar

13. Horwitz LI, Chang C, Arcilla HN, Knickman JR. Quantifying health systems' investment in social determinants of health, by sector, 2017–19: study analyzes the extent to which US health systems are directly investing in community programs to address social determinants of health. Health Aff. (2020) 39:192–8. doi: 10.1377/hlthaff.2019.01246

PubMed Abstract | Crossref Full Text | Google Scholar

14. Espinoza JC, Sehgal S, Phuong J, Bahroos N, Starren J, Wilcox A, et al. Development of a social and environmental determinants of health informatics maturity model. J Clin Transl Sci. (2023) 7:e266. doi: 10.1017/cts.2023.691

PubMed Abstract | Crossref Full Text | Google Scholar

15. Ehrman RR, Malik AN, Haber BD, Glassman SR, Bowen CA, Korzeniewski SJ, et al. The role of place-based factors and other social determinants of health on adverse post-sepsis outcomes: a review of the literature. Front Disaster Emerg Med. (2024) 2:1357806. doi: 10.3389/femer.2024.1357806

PubMed Abstract | Crossref Full Text | Google Scholar

16. Nichols LM, Taylor LA. Social determinants as public goods: a new approach to financing key investments in healthy communities. Health Aff. (2018) 37:1223–30. doi: 10.1377/hlthaff.2018.0039

PubMed Abstract | Crossref Full Text | Google Scholar

17. Medicare Payment Advisory Commission. Health Care Spending and the Medicare Program: A Data Book. Washington, DC: Medicare Payment Advisory Commission (2024).

Google Scholar

18. Burney IL, Schieber GJ, Blaxall MO, Gabel JR. Geographic variation in physicians' fees. Payments to physicians under medicare and medicaid. JAMA. (1978) 240:1368–71. doi: 10.1001/jama.1978.03290130062023

Crossref Full Text | Google Scholar

19. Freburger JK, Holmes GM, Ku LJE, Cutchin MP, Heatwole-Shank K, Edwards LJ. Disparities in post-acute rehabilitation care for joint replacement. Arthritis Care Res. (2011) 63:1020–30. doi: 10.1002/acr.20477

PubMed Abstract | Crossref Full Text | Google Scholar

20. Middleton A, Graham JE, Bettger JP, Haas A, Ottenbacher KJ. Facility and geographic variation in rates of successful community discharge after inpatient rehabilitation among Medicare fee-for-service beneficiaries. JAMA Netw Open. (2018) 1:e184332-e. doi: 10.1001/jamanetworkopen.2018.4332

PubMed Abstract | Crossref Full Text | Google Scholar

21. Reistetter TA, Kuo Y-F, Karmarkar AM, Eschbach K, Teppala S, Freeman JL, et al. Geographic and facility variation in inpatient stroke rehabilitation: multilevel analysis of functional status. Arch Phys Med Rehabil. (2015) 96:1248–54. doi: 10.1016/j.apmr.2015.02.020

PubMed Abstract | Crossref Full Text | Google Scholar

22. Frontera WR, Bean JF, Damiano D, Ehrlich-Jones L, Fried-Oken M, Jette A, et al. Rehabilitation research at the National Institutes of Health moving the field forward (executive summary). Rehabil Psychol. (2017) 62:387–96. doi: 10.1037/rep0000164

PubMed Abstract | Crossref Full Text | Google Scholar

23. Palmer RC, Ismond D, Rodriquez EJ, Kaufman JS. Social determinants of health: future directions for health disparities research. Am J Public Health. (2019). doi: 10.2105/AJPH.2019.304964

PubMed Abstract | Crossref Full Text | Google Scholar

24. Botticello AL, Boninger M, Charlifue S, Chen Y, Fyffe D, Heinemann A, et al. To what extent do neighborhood differences mediate racial disparities in participation after spinal cord injury? Arch Phys Med Rehabil. (2016) 97:1735–44. doi: 10.1016/j.apmr.2016.04.007

PubMed Abstract | Crossref Full Text | Google Scholar

25. Cimarolli VR, Falzarano F, Hicks S. Predictors of rehospitalization in post-acute rehabilitation among different ethnic groups. J Am Med Dir Assoc. (2020) 21:513–8. doi: 10.1016/j.jamda.2019.10.008

PubMed Abstract | Crossref Full Text | Google Scholar

26. Flores LE, Verduzco-Gutierrez M, Molinares D, Silver JK. Disparities in health care for Hispanic patients in physical medicine and rehabilitation in the United States: a narrative review. Am J Phys Med Rehabil. (2020) 99:338–47. doi: 10.1097/PHM.0000000000001342

PubMed Abstract | Crossref Full Text | Google Scholar

27. Freburger JK, Holmes GM, Ku LJE. Postacute rehabilitation care for hip fracture: who gets the most care? J Am Geriatr Soc. (2012) 60:1929–35. doi: 10.1111/j.1532-5415.2012.04149.x

PubMed Abstract | Crossref Full Text | Google Scholar

28. Inneh IA, Clair AJ, Slover JD, Iorio R. Disparities in discharge destination after lower extremity joint arthroplasty: analysis of 7924 patients in an urban setting. J Arthroplasty. (2016) 31:2700–4. doi: 10.1016/j.arth.2016.05.027

PubMed Abstract | Crossref Full Text | Google Scholar

29. Zuccarino S, Fattore G, Vitali S, Antronaco G, Frigerio S, Colombo M. The association between education and rehabilitation outcomes: a population retrospective observational study. Arch Gerontol Geriatr. (2020) 91:104218. doi: 10.1016/j.archger.2020.104218

PubMed Abstract | Crossref Full Text | Google Scholar

30. Reistetter TA, Dean JM, Haas AM, Prochaska JD, Jupiter DC, Eschbach K, et al. Development and evaluation of rehabilitation service areas for the United States. BMC Health Serv Res. (2023) 23:204. doi: 10.1186/s12913-023-09184-2

PubMed Abstract | Crossref Full Text | Google Scholar

31. Reistetter TA, Eschbach K, Prochaska J, Jupiter DC, Hong I, Haas AM, et al. Understanding variation in postacute care: developing rehabilitation service areas through geographic mapping. Am J Phys Med Rehabil. (2021) 100:465–72. doi: 10.1097/PHM.0000000000001577

PubMed Abstract | Crossref Full Text | Google Scholar

32. Cary Jr MP, Prvu Bettger J, Jarvis JM, Ottenbacher KJ, Graham JE. Successful community discharge following postacute rehabilitation for medicare beneficiaries: analysis of a patient-centered quality measure. Health Serv Res. (2018) 53:2470–82. doi: 10.1111/1475-6773.12796

PubMed Abstract | Crossref Full Text | Google Scholar

33. Leland NE, Roberts P, De Souza R, Hwa Chang S, Shah K, Robinson M. Care transition processes to achieve a successful community discharge after postacute care: a scoping review. Am J Occup Ther. (2019) 73:7301205140p1-p9. doi: 10.5014/ajot.2019.005157

PubMed Abstract | Crossref Full Text | Google Scholar

34. Burke RE, Jones CD, Coleman EA, Falvey JR, Stevens-Lapsley JE, Ginde AA. Use of post-acute care after hospital discharge in urban and rural hospitals. Am J Accountable Care. (2017) 5:16–22.

PubMed Abstract | Google Scholar

35. Kosar CM, Loomer L, Ferdows NB, Trivedi AN, Panagiotou OA, Rahman M. Assessment of rural-urban differences in postacute care utilization and outcomes among older US adults. JAMA Netw Open. (2020) 3:e1918738-e. doi: 10.1001/jamanetworkopen.2019.18738

PubMed Abstract | Crossref Full Text | Google Scholar

36. Reistetter T, Bokov, A. Rehabilitation Service Areas (for post-acute care) (v0.1.4-beta). Zenodo (2024).

Google Scholar

37. RTI International. Discharge to Community-Post Acute Care (PAC) LongTerm Care Hospital (LTCH) Quality Reporting Program (QRP). Centers for Medicare & Medicaid Services (2021). Report No.: NQF ID 3480.

Google Scholar

38. RTI International. Discharge to Community(DTC)-Post Acute Care (PAC) Skilled Nursing Facility (SNF) Quality Reporting Program (QRP). Available online at: https://cmit.cms.gov/CMIT_public/ViewMeasure?MeasureId=2846 (accessed November 1, 2021).

Google Scholar

39. Hailu A, Wasserman C. Guidelines for Using Rural-Urban Classification Systems for Community Health Assessment. Washington State Department of Health (2016).

Google Scholar

40. Surbhi S, Tolley EA, Cossman RE, Dashputre AA, Bailey JE. Refining a traditional urban-rural classification approach to better assess heterogeneity of treatment effects in patient-centered outcomes research. MethodsX. (2021) 8:101299. doi: 10.1016/j.mex.2021.101299

PubMed Abstract | Crossref Full Text | Google Scholar

41. Muthén LK, Muthén BO. Mplus User's Guide, 8th edn. Los Angeles, CA: Muthén & Muthén (2017).

Google Scholar

42. Fabrigar LR, Wegener DT, MacCallum RC, Strahan EJ. Evaluating the use of exploratory factor analysis in psychological research. Psychol Methods. (1999) 4:272. doi: 10.1037//1082-989X.4.3.272

Crossref Full Text | Google Scholar

43. Hair JF, Black WC, Babin BJ, Anderson RE. Multivariate Data Analysis, 7th edn. Upper Saddle River, NJ: Pearson Prentice Hall (2010).

Google Scholar

44. Byrne BM. Structural Equation Modeling with Mplus: Basic Concepts, Applications, and Programming. New York, NY: Routledge (2013). doi: 10.4324/9780203807644

Crossref Full Text | Google Scholar

45. Diez Roux AV. Investigating neighborhood and area effects on health. Am J Pub Health. (2001) 91:1783–9. doi: 10.2105/AJPH.91.11.1783

PubMed Abstract | Crossref Full Text | Google Scholar

46. Cummins S, Curtis S, Diez-Roux AV, Macintyre S. Understanding and representing ‘place' in health research: a relational approach. Soc Sci Med. (2007) 65:1825–38. doi: 10.1016/j.socscimed.2007.05.036

PubMed Abstract | Crossref Full Text | Google Scholar

47. Preacher KJ, Hayes AF. Asymptotic and resampling strategies for assessing and comparing indirect effects in multiple mediator models. Behav Res Methods. (2008) 40:879–91. doi: 10.3758/BRM.40.3.879

PubMed Abstract | Crossref Full Text | Google Scholar

48. Hu Lt, Bentler PM. Cutoff criteria for fit indexes in covariance structure analysis: conventional criteria versus new alternatives. Struct Equ Modeling. (1999) 6:1–55. doi: 10.1080/10705519909540118

PubMed Abstract | Crossref Full Text | Google Scholar

49. Bentler PM. Comparative fit indexes in structural models. Psychol Bull. (1990) 107:238. doi: 10.1037//0033-2909.107.2.238

PubMed Abstract | Crossref Full Text | Google Scholar

50. Adler NE, Cutler DM, Fielding JE, Galea S, Glymour MM, Koh HK, et al. Addressing social determinants of health and health disparities: a vital direction for health and health care. NAM Perspectives. (2016).

Google Scholar

51. Williams DR, Mohammed SA. Discrimination and racial disparities in health: evidence and needed research. J Behav Med. (2009) 32:20–47. doi: 10.1007/s10865-008-9185-0

PubMed Abstract | Crossref Full Text | Google Scholar

52. Sentell T, Shumway M, Snowden L. Access to mental health treatment by English language proficiency and race/ethnicity. J Gen Intern Med. (2007) 22:289–93. doi: 10.1007/s11606-007-0345-7

PubMed Abstract | Crossref Full Text | Google Scholar

53. Carrasquillo O, Orav EJ, Brennan TA, Burstin HR. Impact of language barriers on patient satisfaction in an emergency department. J Gen Intern Med. (1999) 14:82–7. doi: 10.1046/j.1525-1497.1999.00293.x

PubMed Abstract | Crossref Full Text | Google Scholar

54. Dahlgren G, Whitehead M. European Strategies for Tackling Social Inequities in Health: Levelling up Part 2. WHO Regional office for Europe Copenhagen (2006).

Google Scholar

55. Braveman P, Egerter S, Williams DR. The social determinants of health: coming of age. Annu Rev Public Health. (2011) 32:381–98. doi: 10.1146/annurev-publhealth-031210-101218

PubMed Abstract | Crossref Full Text | Google Scholar

56. Adler NE, Stewart J. Health disparities across the lifespan: meaning, methods, and mechanisms. Ann N Y Acad Sci. (2010) 1186:5–23. doi: 10.1111/j.1749-6632.2009.05337.x

PubMed Abstract | Crossref Full Text | Google Scholar

57. Benjamins MR, Silva A, Saiyed NS, De Maio FG. Comparison of all-cause mortality rates and inequities between black and white populations across the 30 most populous US cities. JAMA Netw Open. (2021) 4:e2032086-e. doi: 10.1001/jamanetworkopen.2020.32086

PubMed Abstract | Crossref Full Text | Google Scholar

58. Elo IT. Social class differentials in health and mortality: patterns and explanations in comparative perspective. Annu Rev Sociol. (2009) 35:553–72. doi: 10.1146/annurev-soc-070308-115929

Crossref Full Text | Google Scholar

59. Berkman LF, Glass T, Brissette I, Seeman TE. From social integration to health: Durkheim in the new millennium. Soc Sci Med. (2000) 51:843–57. doi: 10.1016/S0277-9536(00)00065-4

PubMed Abstract | Crossref Full Text | Google Scholar

60. Fiscella K, Franks P, Gold MR, Clancy CM. Inequality in quality: addressing socioeconomic, racial, and ethnic disparities in health care. JAMA. (2000) 283:2579–84. doi: 10.1001/jama.283.19.2579

PubMed Abstract | Crossref Full Text | Google Scholar

61. Berkman ND, Sheridan SL, Donahue KE, Halpern DJ, Crotty K. Low health literacy and health outcomes: an updated systematic review. Ann Intern Med. (2011) 155:97–107. doi: 10.7326/0003-4819-155-2-201107190-00005

PubMed Abstract | Crossref Full Text | Google Scholar

62. Clouston SA, Manganello JA, Richards M. A life course approach to health literacy: the role of gender, educational attainment and lifetime cognitive capability. Age Ageing. (2017) 46:493–9. doi: 10.1093/ageing/afw229

PubMed Abstract | Crossref Full Text | Google Scholar

63. Syed ST, Gerber BS, Sharp LK. Traveling towards disease: transportation barriers to health care access. J Community Health. (2013) 38:976–93. doi: 10.1007/s10900-013-9681-1

PubMed Abstract | Crossref Full Text | Google Scholar

64. Mattson J. Transportation, distance, and health care utilization for older adults in rural and small urban areas. Transp Res Rec. (2011) 2265:192–9. doi: 10.3141/2265-22

PubMed Abstract | Crossref Full Text | Google Scholar

65. Williams DR, Lawrence JA, Davis BA. Racism and health: evidence and needed research. Annu Rev Public Health. (2019) 40:105–25. doi: 10.1146/annurev-publhealth-040218-043750

PubMed Abstract | Crossref Full Text | Google Scholar

66. Henning-Smith C, Ecklund A, Kozhimannil K. Rural-urban differences in social isolation and its relationship to health. Innov Aging. (2018) 2:770. doi: 10.1093/geroni/igy023.2851

Crossref Full Text | Google Scholar

67. Miller DR, Safford MM, Pogach LM. Who has diabetes? Best estimates of diabetes prevalence in the Department of Veterans Affairs based on computerized patient data. Diabetes Care. (2004) 27(suppl_2):b10–21. doi: 10.2337/diacare.27.suppl_2.B10

PubMed Abstract | Crossref Full Text | Google Scholar

68. Woolhandler S, Campbell T, Himmelstein DU. Costs of health care administration in the United States and Canada. N Engl J Med. (2003) 349:768–75. doi: 10.1056/NEJMsa022033

PubMed Abstract | Crossref Full Text | Google Scholar

69. Oster A, Bindman AB. Emergency department visits for ambulatory care sensitive conditions: insights into preventable hospitalizations. Med Care. (2003) 41:198–207. doi: 10.1097/01.MLR.0000045021.70297.9F

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: social determinants of health, rehabilitation, healthcare disparities, urban health, rural health

Citation: Nam S, Schmidt S, Dean JM, Bokov AF and Reistetter TA (2025) Social determinants of health and rehabilitation service areas: an urban and rural mediation analysis. Front. Public Health 13:1562610. doi: 10.3389/fpubh.2025.1562610

Received: 31 January 2025; Accepted: 30 May 2025;
Published: 18 June 2025.

Edited by:

Chao Ma, Southeast University, China

Reviewed by:

Zhuoying Qiu, University of Health and Rehabilitation, China
Jeong Min Choi, Mokpo National University, Republic of Korea

Copyright © 2025 Nam, Schmidt, Dean, Bokov and Reistetter. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Timothy A. Reistetter, cmVpc3RldHRlckB1dGhzY3NhLmVkdQ==

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.