Spatial Analysis of Intra-Urban Air Pollution Disparities through an Environmental Justice Lens: A Case Study of Philadelphia, PA

.


Introduction
Globally, air pollution is a major environmental health hazard, especially in urban areas [1].The World Health Organization (WHO) estimates that outdoor air pollution (PM2.5 and ozone) caused 4.2 million premature deaths around the world in 2019 [2].Currently, over half of the world's population resides in cities, but this is predicted to swell to 70% by 2050 due to the economic pressures of industrialization and globalization [3].People living in cities face greater exposure to air pollution because of continued urban population growth, and consequently, these populations also face a greater risk of the negative effects of air pollution on human health and well-being [4].
PM2.5, or particulate matter with an aerodynamic diameter less than or equal to 2.5 μm, is of great concern for urban areas such as Philadelphia, as it is primarily generated from anthropogenic actions such as industrial emissions, fossil fuel combustion, and gasfueled vehicles [5].Because of the small particle size, PM2.5 can penetrate deep into the lungs, leading to increased incidence of asthma and other respiratory illnesses, as well as all-cause mortality rates [6,7].As of 2024, the U.S. Environmental Protection Agency (EPA) has set their 24-h PM2.5 standard at 35 μg/m 3 and their annual PM2.5 standard at 9.0 μg/m 3 [8].The WHO has set their 24-h PM2.5 standard at 15 μg/m 3 and their annual PM2.5 standard at 5 μg/m 3 [9].As of 2021, 37.5 million Americans were living in counties with ambient concentrations of PM2.5 above the acceptable levels set by the National Ambient Air Quality Standards (NAAQS).The American Lung Association reports fine particulate pollution (i.e., PM2.5) as one of the most widespread and dangerous air pollutant along with ozone in its 2022 "State of the Air" report [10].Black carbon (BC), also known as soot, is a subset of particles, mainly PM2.5, that is especially detrimental to the environment and human health because of its warming capability and small size.BC is generated from incomplete fossil fuel combustion, which, in the United States, primarily is a consequence of vehicle emissions and household energy consumption [11].Currently, there are no set standards on ambient BC concentrations in the U.S., although it is regulated indirectly through the PM2.5 standards.
Air pollution is a pertinent issue for environmental justice advocates because it disproportionately burdens vulnerable populations based on race, ethnicity, and socioeconomic status [12,13].Nationwide, people of color are 61% more likely than White people to live in a country with a failing grade for at least one air pollutant [10].Some groups, such as children, the elderly, and low-socioeconomic status populations, have a greater risk for negative health outcomes from air pollution exposure [14][15][16] due to their limited capacity to cope with the effects of exposure.The social determinants of health, non-medical factors at the individual, community, and systemic level that shape the conditions of daily life, can compound health risk from air pollution exposure as people can face multiple layers of vulnerability.Vulnerability refers to both the negative socioeconomic and health consequences of exposure to an environmental hazard, as well as the ability of an individual or a group to cope and recover from the exposure [17].In the case of air pollution, exposure is often ongoing, which means affected populations have limited strategies for coping with the air pollution on an individual scale, which can lead to feelings of helplessness [18].Improving the understanding of the intersections of vulnerabilities like race, socioeconomic status, and health with air quality is crucial for fostering collective action to address environmental inequities and ensure justice.
While air pollution and environmental justice have been widely studied, research at the intra-city scale has been stymied due to insufficient fine-scale air pollution data.Several previous studies have relied on proxies, such as health outcomes [19] or proximity analysis [20], and spatially limited monitor-based air pollution observations [21].However, air pollution in an urban context is highly variable due to the dynamic physical landscape of the city, meaning concentrations can vary dramatically over a short distance [22,23].For example, Miller et al. found that differences in exposure to particulate matter across individual cities showed greater differences than overall exposure compared between cities [24].Thus, fine-scale, citywide air pollution data are critical to illuminating the intra-urban implications of air pollution for environmental justice.Our group recently used a mobile monitoring approach to measure PM2.5 and BC concentrations in the city of Philadelphia with Structure of Urban Landscapes (STURLA) classifications to interpolate citywide air pollution prediction models for PM2.5 and BC [25].This air pollution prediction model assigned air pollutant concentrations to 120 m 2 grid cells across the city, enabling citywide research and analysis.The main objective of this study is to use these finescale data to investigate the granular relationships between urban air pollution and environmental justice.
Studies show that accounting for spatial dependence in the data is important both in research on air pollution and environmental justice.Use of spatially autoregressive [26], geographically weighted [27] and land-use regression models [28] have been shown to improve the accuracy of air pollution estimates and reduce uncertainty.In the case of Shen et al., their model of air pollution across the European continent using spatially varying linear regression even outperformed machine learning methods [27].There is also demonstrated utility of considering spatial dependence in models investigating air pollution exposure through an environmental justice lens [16,[29][30][31].Park et al. [32] suggests that the lack of consideration for spatial autocorrelation may have biased previous assessments of air pollution's influence on asthma risk and shows that dealing with spatial autocorrelation can create stronger asthma risk models.
This study aims to use spatially autoregressive modeling techniques to understand the intra-urban relationship between air pollution, specifically PM2.5 and BC, and socioeconomic and health data across the city of Philadelphia, PA.Specifically, poverty, non-White populations, and asthma will be used to demonstrate the social, economic, and health dimensions of vulnerability, although we recognize that these variables do not encapsulate all of the complexities that make an individual or group vulnerable.This study utilizes 120-m city wide air pollution predictions in the City of Philadelphia from Cummings et al. [25].We aggregated air pollution data to the block group scale to (a) assess the presence of spatial dependence in air pollution, socio-economic, and health data; (b) compare spatially autoregressive models to more commonly used ordinary least squares (OLS) regression; and (c) highlight intra urban disparities in air pollution exposure.Access to air pollution data at such a fine spatial scale is rare, and it creates a unique opportunity to assess the distributional injustice of air pollution within a city.

Study Area
The study area for this analysis is the city of Philadelphia, Pennsylvania, located on the east coast of the United States between the Schuylkill and Delaware Rivers.Philadelphia is Pennsylvania's largest city, encompassing 367 km 2 of land with a population of 16,037,909 residents as of 2020 [33].Philadelphia consists primarily of high-density development with a population density of 4609 people per square kilometer [33] and, as the population continues to grow, projections predict that the city will gain more dense urban development in the coming years [34].The city's heavy industry is primarily located on the southern and eastern borders along the banks of the Delaware and Schuylkill rivers.In contrast, the northern and western areas contain many large parks and an overall greater presence of green space [35].Motor vehicle emissions, construction, and industrial activities are major sources of air pollution for the city [36].
Philadelphia contains a significant presence of marginalized communities.23.1% of the population of Philadelphia lives below the poverty line, making it one of the poorest cities in the United States [33].As of 2020, the population comprised 41.4% Black, 39.3% White, 15.1% Hispanic or Latino, and 7.4% Asian peoples [33].Policies such as redlining and urban renewal have made Philadelphia into one of the most segregated cities in the US [37].Philadelphia has been identified as having the seventh-highest level of segregation between Black and White populations of all US cities [37].White populations are most concentrated in South Philadelphia, Center City, and the surrounding suburbs, while Black communities are located in the Northern and Western areas of the city (Figure 1).
The American Lung Association (ALA) ranked Philadelphia 18th in year-round particle pollution in 2022 [10].The ALA highlighted people of color and people experiencing poverty as groups facing significant air quality risk in Philadelphia County [38].The city's poor air quality has translated into negative health outcomes for its residents.Childhood asthma prevalence in the city of Philadelphia is 21%, which is almost three times greater than the national prevalence of 5.8% as of 2021 [39].The burden of disease caused by air pollution has disparate impacts on younger and older populations, as well as communities of color.For example, Black and Hispanic children are over four times more likely to be hospitalized for asthma-related complications than White children [40].

Data Source
This study focused on five variables in total with two representing air pollution exposure (predicted PM2.5 and predicted BC) and three representing social, economic, and health vulnerabilities (percent non-White population, percent of population in poverty, and asthma prevalence in the population).This study employs fine-scale air pollution prediction from the work of Cummings et al. [25], which used a mobile monitoring method to collect PM2.5 and BC concentrations [41] and Structure of Urban Landscapes (STURLA) classifications to interpolate citywide air pollution prediction models for PM2.5 and BC [25].The air pollution data were averaged for each block group within the study area to allow for spatially consistent analysis with the social, economic, and health variables (Figure 1).
Race and poverty data were downloaded from the U.S. Census Bureau [42,43].Race data were obtained from the 2020 Decennial Census while the poverty data were sourced from the 2020 American Community Survey.Asthma data were downloaded from the 2019 Centers for Disease Control 500 Cities Project dataset for Philadelphia on the census tract level [44].Block groups with no population were excluded from our analysis.Further, not all block groups contained data for asthma prevalence, so block groups without this data were excluded from the health models.

Geospatial Analysis
Socioeconomic, health, and air pollution data were mapped using ArcGIS Pro 3.2 to visualize their spatial distribution within the city (Figure 1).These data were imported into R Studio and assessed for spatial dependence using the Moran's I test.Spatial dependence refers to the relationship between a measured variable at neighboring locations.The Moran's I test measures spatial dependence through spatial autocorrelation.The test calculates a Moran's I value measuring how similar or dissimilar the data are at neighboring locations (Equation ( 1)) by comparing the standard deviations of each observation with the mean standard deviation of its neighbors weighted by their spatial relationships [45].The Moran's I value ranges from −1 to 1 with negative values indicating negative spatial autocorrelations and dissimilarity among neighbors, and positive values indicating positive spatial autocorrelation and the presence of clustering.
where n is the number of observed locations,   and   are observations of the variable of interest at neighboring locations, ̅ is the mean of the variable of interest at neighboring locations, W is the spatial weights matrix, and   is the weight between I and j in the spatial weights matrix.

Statistical Analysis
To understand the relationships between the variables, univariate Ordinary Least Squares (OLS) regressions and univariate spatial autoregressive (spatial error, spatial lag) regressions were conducted in R Program ver 4.3.OLS is often used to quantify the relationship between a dependent variable and one or more independent variables by estimating the coefficients of a linear regression equation that captures the underlying relationship between the variables.As spatial dependence is present in the datasets, the foundational assumption of independence for linear regression is violated.Thus, spatial lag and spatial error models are employed, since they account for the spatial autocorrelation of the data.OLS analysis is conducted to provide a reference for the importance of considering spatial autocorrelation.
Univariate regression was purposely chosen to isolate the relationships between air pollution and vulnerability variables.Multicollinearity, particularly of the vulnerability variables, was of concern as it can lead to model overfitting and limits the interpretability of multiple regression models.Previous research has established the correlations between populations of color, poverty, and asthma [40,[46][47][48].For instance, in a nationwide study, Keet et al. [48] found that individual and neighborhood-level poverty, as well as Black racial identity, were associated with increased risk of asthma.Thus, multiple regression would be inappropriate, as we know from empirical research that the vulnerability variables are highly correlated.
In the Spatial Lag model, a new spatially lagged dependent variable is introduced to represent the spatial effects of the dependent variable in neighboring observations [49].The model is expressed mathematically as: where y is the dependent variable, W is the spatial weights matrix, p is the coefficient to the spatial autocorrelation term which represents the strength of the spatial dependence, X is a matrix of observations on the explanatory variables, β is a vector of regression coefficients, and ϵ is the error term.The pWy term is the spatial lag term which captures spatial autocorrelation based on how much the value of y at a location is influenced by the neighboring location's y values.
In contrast to the Spatial Lag model's approach to capturing spatial autocorrelation as a dependent variable, the Spatial Error model accounts for spatial dependence through the model's error term [50].In the error term, the error from each observation is correlated with the error from its neighboring observations as shown in Equation ( 3): where y is the dependent variable, X is a matrix of observations on the explanatory variables, β is a vector of regression coefficients,  is the vector of error terms, spatially weighted using the weights matric (W),  is the Spatial Error coefficient, and  is a vector of independent identically distributed error terms.This model is most useful when the spatial dependence in the data is caused by factors not included in the model.For the socio-economic factors, air pollution served as the dependent variable, while the socioeconomic data served as independent variables.For the relationship with asthma prevalence, the air pollutants were independent variables, since outdoor air pollution has a well-documented causal relationship with asthma [51,52].

Results
The results of Moran's I tests (Table 1) indicate that spatial dependence is present at some magnitude in all the datasets used in this analysis.The test shows high levels of spatial autocorrelation in PM2.5 and BC predictions, as well as % non-White and asthma prevalence.Spatial autocorrelation was assessed at the block group scale for all variables.Comparisons of the results of OLS and spatially autoregressive models are shown in Tables 2-4 for the socioeconomic and health variables.In all cases, the spatial autoregressive models had higher r-squared values than the OLS model, indicating the importance of accounting for spatial dependence in terms of model performance.
Results differed for PM2.5 and BC models.Specifically, the Spatial Error model performed best in terms of both r-squared and AIC for evaluating the relationships between PM2.5 and percent in poverty and percent non-White.Given that the distribution of air pollution is related to but not completely caused by socio-economic variables, this model's strong performance, especially in comparison to the Spatial Lag model, is to be expected.However, in the case of BC, although the Spatial Error models had higher r-squared and lower AIC than the Spatial Lag models, the p-values of the regression coefficient in the spatial error models were not significant for the poverty and asthma models.This indicates that there is not sufficient evidence to reject the null hypothesis that the coefficient is equal to zero.Thus, there is not definitive evidence of an association with BC for these specific Spatial Error models and the model results are insignificant.This may be due to spatial heterogeneity in relationship, meaning that the relationship between the vulnerability variables and BC varies across space, thus requiring more complex models to accurately discern this.However, the p-value of the Spatial Error model assessing the relationship between BC and non-White population was less than 0.05, indicating a significant relationship between the two variables.
When accounting for spatial autocorrelation in the residuals, there is a strong relationship between PM2.5 concentrations and the socioeconomic variables (Tables 2 and 3).The r-squared values for the percent non-White and percent in poverty Spatial Error models were both 0.85.This is much greater than the r-squared values calculated for the OLS (percent non-White = 0.11 and percent in poverty = 0.13) and Spatial Lag (percent non-White = 0.60 and percent in poverty = 0.57) models.Further, both variables' models showed low AIC values (−313.232for non-White and −306.48 for poverty) in comparison with OLS AIC values (1687.62 for non-White and 1668.99 for poverty), demonstrating the strong fit of the Spatial Error model.The Spatial Error model results show there is a small positive relationship between PM2.5 and both percent non-White (regression coefficient = 0.002) and percent in poverty (regression coefficient = 0.001) for p ≤ 0.05.Overall, all coefficients reflect a small magnitude positive association between socioeconomic variables and PM2.5.The Spatial Error coefficient was equal to 0.94 in both variables' models, indicating a high level of spatial dependence in their residuals.The Breusch-Pagan test shows the presence of heteroskedasticity in all models, which is to be expected given the spatial dependence present in the residuals.The Spatial Lag models outperformed other models when assessing the relationship between asthma and the air pollutants (Table 4).The Spatial Lag model had an r-squared value of 0.88 and an AIC of 2788.12 for PM2.5 and an r-squared value of 0.88 and an AIC of 2796.51 for BC.This is far greater than the r-squared of 0.18 and far less than the AIC of 5065.48 for the PM2.5 OLS model and the r-squared of 0.09 and far less than the AIC of 5179.95 for the BC OLS model, indicating a stronger fit for the spatially autoregressive model.While the Spatial Error models had high r-squared and AIC values, the p-values of the regression coefficients were not significant, and there is insufficient evidence to conclude that the air pollution variables have a significant effect on asthma prevalence.The Spatial Error models produced smaller regression coefficients than the other models, which may mean its effect is hard to detect (Supplemental Figures S1-S3).Nevertheless, in both Spatial Lag models, there is a positive relationship between the air pollutants and asthma prevalence, although the coefficients were relatively small.The regression coefficient is 0.45 for the PM2.5 model and 0.69 for the BC model for p ≤ 0.05.In contrast with the PM2.5, the p-values of the Breusch-Pagan test for the PC spatial models were greater than 0.05, indicating there is not a significant level of heteroskedasticity in the models.This test reveals that there is relatively consistent variability in the residuals of these models, an important consideration for regression performance and reliability.

Discussion
This research aims to investigate how socio-economic and health vulnerabilities intersect with air pollution in Philadelphia.Overall, accounting for spatial dependence in any capacity through the use of either Spatial Lag or Spatial Error models improved model performance when compared to the OLS model.This finding is consistent with previous studies' findings that also emphasize the importance of spatial autocorrelation in environmental justice analyses [16,30,31].Positive relationships were found between PM2.5 and the socioeconomic variables tested when accounting for spatial dependence in the residuals using the Spatial Error model.In effect, higher concentrations of PM2.5 are associated with higher percentages of the population living in poverty and being non-White.Higher asthma prevalence was associated with greater PM2.5 and BC concentrations for the Spatial Lag models (r-squared = 0.88), while the regression coefficients of the Spatial Error models (r-squared = 0.88) did not have significant p-values, despite also having a high r squared value.It should be noted, however, that in both the socioeconomic and health models, the significant regression coefficients were small.
Previous work has shown air pollution exposure varies across the United States for socioeconomic population subgroups [53][54][55].Our results are consistent with these studies in finding that populations living in census block groups with greater proportions of people of color and living in poverty face greater exposure to PM2.5 [53][54][55][56].Modeled air quality data are commonly used to provide the spatially continuous air quality data needed to investigate spatial inequities [20,55,56].For example, Bravo et al. and Collins et al. both conducted their analysis using EPA air pollution data and the Downscaler methods to create air quality data for the United States.These methods yield 12 × 12 km predictions, which are too coarse for meaningful urban analysis [55,56].Maroko used air dispersion modeling to model PM2.5 for the 733,517 populated tax lots across New York City, NY, to show heightened exposure for Latinx and poorer populations within the Bronx and Brooklyn [20].Although their results did not show significant city-wide exposure inequities, the air dispersion modeling used by Maroko only accounts for PM2.5 from local stationary sources, meaning their results do not include transportation emissions, which account for a large portion of urban air pollution [20,57].
Our study showed statistically significant relationships between BC and percent non-White as well as asthma prevalence, although the regression coefficients were small.While the OLS and Spatial Lag models for percent poverty had statistically significant regression coefficients, their low r-squared values indicate that the models do not fully explain BC distribution.The Spatial Error model had a much higher r-squared value of 0.70 although the regression coefficient was not statistically significant.BC emissions are a product of the incomplete combustion of carbon-based fuels most often used to power diesel engines in transportation vehicles [58] and heavy machinery [59].BC can also originate from other anthropogenic and natural sources such as waste burning [60] and wildfires [61].As seen in Figure 1, concentrations of BC are greatest along the eastern border of the city where Interstate-95 runs, a major highway that connects much of the east coast of the United States.I-95 serves as a critical transportation linkage between Philadelphia and New York City, and the portion that runs through Center City Philadelphia sees on average 160,000 vehicles a day, including many large trucks that run on diesel [62].Most of the block groups adjacent to the highway are sparsely populated, since they are primarily used for industrial purposes which may have impacted the BC models.
While there are many environmental justice studies that look at PM2.5 as a whole, there are few that specifically investigate ambient BC despite its deleterious effects on public health and environmental quality [63][64][65].Our study found a positive relationship between the percentage of people of color living in a census block group and BC concentrations.This is in line with previous studies that show how living in a neighborhood with a greater proportion of non-White population is positively correlated with increased BC exposure [66][67][68].For example, Northcross et al. conducted a case study that measured higher concentrations of BC and PM2.5 in Ivy City, a historically Black neighborhood of Washington D.C., when compared with the rest of the city [67].Their study presents evidence that the city's current air monitoring network was not accurately capturing the burden of air pollution faced by the neighborhood with vulnerable populations and prevented a previously approved parking lot from being constructed.Our and Northcross et al.'s studies show the power of accurate pollution data to reveal inequalities and reparative policies and actions.
In cities with documented segregations such as Philadelphia [69], where variations in the urban landscape can be scars of discriminatory development practices, it is especially critical to understand the social implications of intra-city variability in hazardous air pollution exposure.While the distributional injustice of air pollution has been widely studied [70,71], there have been difficulties finding empirical evidence of this pattern at the intra-city scale due to a lack of fine-scale air pollution data [20,72].This research is unique to previous work in that we used spatially autoregressive models and finer-scale air pollution data (census block group).Air pollution is not spatially independent due to the complex and overlapping interactions of emission sources, atmospheric processes, topography, meteorologic conditions, and other geographic factors that affect its distribution.Thus, using OLS to study air pollution is inappropriate, as the foundational assumption of independence is violated.In this study, spatial autoregressive models were used to account for air pollution's spatial dependence and improve model strength.It is also important to study air pollution at the finest scale possible, since concentrations can vary dramatically over a small geographic area due to many of the same temporal and geographic factors that induce its spatial dependence.These small-scale variations become more prevalent in an urban environment where the composition and topography of the landscape are highly varied and there are many possible sources of pollution present [22,73].
Association of air pollution exposure with non-White and poor populations highlights a significant environmental justice problem for the city of Philadelphia.This association represents an instance of distributive injustice where neighborhoods that have higher concentrations of poverty and people of color are saddled with a greater air pollution burden than the rest of the city.The regression coefficients of our models were, in almost all cases, very small.Nevertheless, over the course of a lifetime, even small magnitude differences in air pollution exposure can have negative consequences for health and well-being.Weichenthal et al. found that PM2.5 still has a significant impact on mortality even at very low concentrations (<5 μg/m 3 ), estimating that around 1.5 million deaths annually can be attributed to low level PM2.5 [74].As PM2.5 concentrations used in this study ranged from 11.11 to 14.08 μg/m 3 , small-scale differences of concentrations in tracks with varying levels of vulnerability can have significant impacts on the population's quality of life and overall well-being.
Previous studies in the city of Philadelphia have documented many other cases of distributive environmental injustices from heat exposure [75] to tree canopy cover [76] to access to high-quality parks [77].Taken together, these injustices have a compounding effect, causing the continued reproduction of systemic inequities and negative consequences for the people living in these communities.For example, proximity to a point source pollution can affect home prices in an area, preventing members from being able to sell their homes and relocate away from harm and making the community less desirable for outside investment [78].In effect, the relationships identified in this study serve to further demonstrate the spatially dependent environmental inequities that exist within the city and provide quantifiable evidence of the need for transformative actions.

Uncertainties and Limitations
The citywide air pollution dataset in this study is a modeled dataset carrying inherent limitations.The modeled dataset from Cummings et al. [25] is based on air pollution data that were collected through mobile monitoring over 12 days (6 replications) and driven on a 483 km (300 mile) route with the sensors mounted to a vehicle roof.The authors acknowledge that limiting their data collection to roadways may have introduced a layer bias into their results.Although their air pollution observations were greater than those collected at stationary EPA stations, both measurements captured similar patterns [25].Further limitations to this data include the fact that air quality was only measured during the summer, that there were only six sampling repetitions at each location, and that the measurements were collected at different times throughout the day, rather than at the same time.PM2.5 and BC concentrations vary seasonally and diurnally based on changes in meteorological conditions and anthropogenic activity [79].A dense network of qualitycontrolled air pollution monitoring stations may be needed to create a more accurate representation of fine-scale air pollution distribution across the city [80].
This study was limited in that it only investigated two types of air pollution (PM2.5 and BC) and three vulnerability variables (percent non-White, percent poverty, and asthma prevalence).Philadelphia also experiences elevated levels of other pollutants, such as ozone.The publicly available air quality data (airnow.gov(accessed on 11 June 2024)) show several days exceeding the national ambient ozone standard level in Philadelphia [81].Vulnerability is nuanced, and we recognize that these three variables do not fully capture the many intersecting physical, social, economic, and health interactions that contribute to an individual or group's overall vulnerability.Due to concerns over multicollinearity, we refrained from conducting multiple regression analysis inhibiting our ability to understand the interaction effects of the chosen vulnerability variables.Future research should seek to expand the breadth and amount of vulnerability-related variables they investigate, and attempt to understand their interacting effects.

Conclusions
This study investigates the intersection of socio-economic and health vulnerabilities with air pollution in Philadelphia.Through the use of spatial autoregressive models and finer-scale air pollution data, we reveal spatial patterns of environmental disparity within the city.Our findings highlight the disproportionate burden of air pollution borne by vulnerable communities, particularly those characterized by higher percentages of poverty and non-White populations.The significance of spatial dependence in our models underscores the complexity of air pollution distribution and emphasizes the importance of accounting for spatial autocorrelation in environmental justice analyses.By revealing the spatially dependent nature of environmental inequities, our research contributes to a deeper understanding of intra-city variability in hazardous air pollution exposure.Despite the small magnitude of regression coefficients, the health implications of even minor disparities in air pollution exposure cannot be overlooked.These findings underscore the urgent need for transformative actions to address environmental injustices in Philadelphia and beyond.Implementing policies and interventions that prioritize the well-being of vulnerable communities is essential for promoting environmental equity and ensuring a healthier future for all residents.Future research should aim to expand the breadth of vulnerability-related variables investigated and explore their interacting effects.Additionally, efforts to improve the accuracy of air pollution models and coverage of air pollution monitoring networks to capture intra-urban variability are critical for developing a more comprehensive understanding of fine-scale air pollution distribution and its impacts on public health and environmental quality in urban areas, especially in vulnerable communities.

Supplementary Materials:
The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/atmos15070755/s1, Figure S1: Graphs comparing observed PM2.5 and BC concentrations with those modeled using OLS, Spatial Lag (SLag) and Spatial Error (SpErr) models and percent poverty, Figure S2: Graphs comparing observed PM2.5 and BC concentrations with those modeled using OLS, Spatial Lag (SLag) and Spatial Error (SpErr) models and percent nonwhite, Figure S3: Graphs comparing observed asthma prevalence with those modeled using OLS, Spatial Lag (SLag) and Spatial Error (SpErr) models and PM2.5 or BC concentrations.

Author Contributions:
Conceptualization, P.K. and K.M.S.; Methodology, P.K. and K.M.S.; Validation, M.S., C.B., P.K. and K.M.S.; Formal Analysis, M.S. and C.B.; Investigation, M.S., C.B., P.K. and K.M.S.; Resources, P.K. and K.M.S.; Data Curation, M.S., C.B., P.K. and K.M.S.; Writing-Original Draft Preparation, M.S. and C.B.; Writing-Review & Editing, M.S., C.B., P.K. and K.M.S.; Visualization, M.S.; Supervision, P.K. and K.M.S.; Project Administration, P.K. and K.M.S.; Funding Acquisition, P.K. and K.M.S.All authors have read and agreed to the published version of the manuscript.Funding: This research was supported by the NSF HEGS grant #1832407.Institutional Review Board Statement: Not applicable.Informed Consent Statement: Not applicable.Data Availability Statement:The data used in this paper comes from previous research from Cummings et al. 2022, which is cited in the "Data Source" portion of the Method section, or publicly available sources like the U.S. Census Bureau.

Table 1 .
Summary of global Moran's I statistic for each variable.

Table 2 .
The percent of non-White population correlated with PM2.5 and BC concentrations at the census block group level.

Table 3 .
The percent of population in poverty correlated with PM2.5 and BC concentrations at the census block group level.

Table 4 .
PM2.5 and BC concentrations correlated with asthma prevalence at the census block group level.