Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Utilizing MicroGenetic Algorithm for Optimal Design of Permanent-Magnet-Assisted WFSM for Traction Machines
Previous Article in Journal
Multi-Modal Emotion Recognition for Online Education Using Emoji Prompts
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

How May Building Morphology Influence Pedestrians’ Exposure to PM2.5?

by
Yogita Karale
and
May Yuan
*
School of Economic, Political and Policy Sciences, The University of Texas at Dallas, Dallas, TX 75080, USA
*
Author to whom correspondence should be addressed.
Appl. Sci. 2024, 14(12), 5149; https://doi.org/10.3390/app14125149
Submission received: 13 May 2024 / Revised: 7 June 2024 / Accepted: 11 June 2024 / Published: 13 June 2024
(This article belongs to the Section Environmental Sciences)

Abstract

:
Due to their sparse distribution and placement in open areas, fixed air-quality-monitoring stations fail to characterize the effect of contextual factors such as buildings on the dispersion of PM2.5. This study evaluated the effects of building morphology on PM2.5 dispersion in a pedestrian-friendly area on the University of Texas at Dallas campus, spanning approximately 0.5 km2. The study collected PM2.5 data along five distinct paths exhibiting varying building morphological characteristics in terms of size, height, density, and spacing at a high spatial resolution. The interquartile range of PM2.5 levels across nine data-collection runs varied from 0.3 µg/m3 to 1.7 µg/m3, indicating relatively uniform PM2.5 levels within the study area. Furthermore, weather-related variables played a dominant role in PM2.5 distribution as temporal variation over-powered spatial variation in the PM2.5 data. The study employed a fixed-effects model to assess the effect of time-invariant morphological characteristics of buildings on PM2.5 and found that the buildings’ morphological characteristics explained 33.22% variation in the fixed effects in the model. Furthermore, openness in the direction of wind elevated the PM2.5 concentration.

1. Introduction

Information about PM2.5 concentration from the nearest air-quality-monitoring station may not represent PM2.5 in a microenvironment. Most monitoring stations are in open areas to avoid the effect of obstructions such as buildings [1]. Also, because of their fixed locations, PM2.5 observations from these stations provide limited opportunities to study spatial variabilities of PM2.5 in places where most people live, and the health effects are most impactful.
In contrast to data from fixed-ground monitoring stations, which are often sparsely located, mobile monitoring platforms provide flexibility to gather data at high spatial resolutions and temporal frequencies. Recently, several studies used mobile monitoring platforms to measure pollutant concentration [2,3,4,5,6]. These studies involved data collection across roads in urban areas using vehicles as a platform for PM2.5 monitoring [2,3,4,5,6]. Some studies used aerial platforms to study the horizontal and vertical profiles of PM2.5 [7,8].
Spatial obstructions like buildings and trees affect the dispersion and concentration of PM2.5 [3,9,10,11]. Land use regression models incorporating the buildings’ morphological characteristics to estimate PM2.5 and PM10 during stable meteorological conditions improved the model’s R2 by 10% in the high-density area of downtown Hong Kong [3]. Other studies examined how urban morphology influences the flow of air pollutants in the high-density environment of Hong Kong. Some of the key variables affecting air pollution in these studies were site coverage, average building height, distance between buildings, and degree of enclosure [12,13].
In addition to physical structures or trees that might influence dispersion, PM2.5 distributions have complex relationships with geographic features. PM2.5 refers to all particles with an aerodynamic diameter smaller than 2.5 μm. The chemical components of PM2.5 vary depending on the sources of emissions in a region. Vehicle exhaustion is one of the major sources of black carbon, and black carbon correlates better with traffic patterns than with PM2.5 concentration [14]. PM2.5 observations from near-road sites appeared to be strongly correlated with nearby sites in terms of background pollution [11], but PM2.5 averages at near-road sites appeared higher than those at other nearby sites [10,11]. This indicates that on average, near-road areas are susceptible to relatively greater PM2.5 exposure. However, two studies conducted at near-road sites in the United States showed that traffic characteristics like traffic volume and traffic speed did not correlate well with PM2.5, but meteorological factors and their interaction with site characteristics had profound impacts on PM2.5 [10,11]. In both prior studies, a stable meteorological condition at night showed elevated PM2.5 concentration, whereas higher wind speed reduced the PM2.5 level. A source apportionment analysis at a roadside location in Maryland showed that on-road traffic contributed only 12.5–17% of PM2.5 [10]. High PM2.5 was associated with wind nearly perpendicular to the road, since the wind swept the maximum surface area of the road towards the monitor compared to when it blew precisely perpendicular to the road segment [11]. Moreover, the monitoring station with buildings on the windward side consistently displayed high PM2.5, suggesting that the presence of buildings trapped the pollutants coming from the roads [11].
Thus, besides emission sources, geographic features, spatial forms, and their interactions with wind played a role in determining PM2.5 concentration at locations. Mobile monitoring platforms allow for studying PM2.5 distributions in different urban microenvironments, such as the effect of trees and buildings near emission sources (e.g., roads) on PM2.5 concentration [3,15]. However, PM2.5 research in low-density and non-near-road sites appears to be lacking but can offer insights into the dynamics of PM2.5 concentration in microenvironments typical of pedestrian areas or shopping plazas. This study aimed to investigate PM2.5 variations in such a built environment.

2. Materials and Methods

2.1. Study Area

We chose the University of Texas at Dallas (UTD) campus, located in the Dallas–Fort Worth Metroplex, Texas, USA (Figure 1) for the study. The Campus spanned across Dallas and Collin counties. The figure depicts the metroplex boundary, associated counties, and point source PM2.5 emissions (year 2019) in the area obtained by the Texas Commission on Environmental Quality from State of Texas Air Reporting (STAR) system [16]. Around the campus, several manufacturing facilities contributed as point sources of PM2.5 emissions. However, emissions from most sources were minimal, falling within a range of 0–3 tons per year. Among non-point-emission sources, dust from construction and roads (paved or unpaved), emissions from waste disposal, commercial cooking, and residential fuel combustion (wood) were common sources of PM2.5 in Dallas and Collin counties. In Collin County, dust from agriculture activities related to crops and livestock was also one of the major sources of non-point PM2.5 emissions [17].
A part of the University of Texas at Dallas was chosen as an example of pedestrian-dominant urban settings where city dwellers are more likely to experience direct exposure to air pollution than when indoors or in a vehicle in other urban settings. Four roads surrounded the University: W Campbell Road, Waterview Parkway, N Floyd Road, and Synergy Park Blvd. Figure 2 displays an annual average daily and peak-hour traffic count for these roads. Specifically, the study focused on five paths (A, B, C, D, and E) in the interior of the University, away from the major roads mentioned above (Figure 2). Each path was surrounded by buildings with various characteristics in terms of size, spacing, and spatial configuration. PM2.5 data collection followed these paths when traffic was light.

2.2. Data Sources

2.2.1. PM2.5 Data

The study used DR1000 from Scentroid (Stouffville, ON, USA) [19], a portable PM2.5-monitoring instrument, to collect the PM2.5 data. Before data collection (in December 2019 and February 2020), we calibrated DR1000 using a reference instrument, Fidas® Frog (Palas GmbH, Karlsruhe, Germany), a fine dust measurement device using optical light scattering of a single particle technology. During its production, Fidas® Frog was validated for quality by comparing it with a certified Fidas® 200. Fidas® Frog measures particle number, size distributions, and mass (for PM1, PM2.5, PM4, PM10, and total suspended particles) [20] and has been used in other research studies for direct measurements or calibration [21,22,23].
During the calibration, two hours of co-located measurements (on 3 September 2019) from DR1000 and Fidas® Frog were collected. DR1000 recorded PM2.5 measurements at an interval of 3–4 s, whereas Fidas® Frog every 5 s. After the calibration, both instruments simultaneously collected data for about an hour (on 16 September 2019). For comparison, 1-min averages were calculated for each instrument. Figure 3 displays 1-min average measurements from both instruments before and after calibration. The goodness of fit (R2 value) improved from 0.55 to 0.81 after calibration.
Each collection of PM2.5 observations ran along all five paths with DR1000 mounted on a bicycle at a height of about 1.2 m above the ground. A total of nine data collection runs were completed: eight runs in December 2019 and one run in February 2020. Out of these nine runs, three runs collected data in the mornings, three in the afternoons, and three in the evenings. Each data-collection run took about an hour. Each run collected data from the north-to-south direction, was repeated from the south-to-north, and then repeated three times in each direction.

2.2.2. Weather Data

Variables such as wind direction and speed, temperature, and relative humidity could affect the PM2.5 observations. The weather station on the roof of Residential Hall West (Figure 2) on campus supplied weather data at the observation frequency of one minute. The weather station was about 320 m from path A in the northwest. The study retrieved weather data closest to the PM2.5 timestamps.

2.2.3. Building Data

The study used building footprints released by Microsoft (version 2.0) [24] and validated in the study. However, footprints for new construction on campus were missing in the data. Such footprints were digitized manually.

2.2.4. Digital Surface Model (DSM)

A high-resolution LiDAR point cloud for the study area was downloaded from USGS Geological Survey [25]. A DSM with the 50 cm resolution was derived from the LiDAR point cloud using ArcGIS Pro 2.5 (Esri, Redlands, CA, USA)

2.3. Data Preparation

Figure 4 depicts the flowchart for data preparation, while the following subsection describes the steps involved in preparing the final dataset for the analysis.

2.3.1. Segment-Wise PM2.5 Data

The PM2.5 data came from the data collection paths shown in Figure 2. Data processing included steps to divide each path into segments of 50 m and average PM2.5 observations along each 50 m segment as a representative measure. The divide-and-average process smoothened the measurement uncertainty from the instrument and variations in movement speed during the collection, while being sufficient to capture environmental variabilities within and between all paths.

2.3.2. Building Morphology

Building morphological parameters included building coverage ratio (ratio of the area occupied by buildings to total buffer area), mean building area, number of buildings, the average distance between nearest-neighbor buildings, and mean building height. All these parameters could affect pollution dispersion [12,13]. Additionally, Shi et al. [3] considered frontal area index and sky view factor to model PM2.5 and PM10. Both these factors were related to openness. They found that the frontal area index, representing building area in the direction of wind, played an important role in determining PM2.5 concentration. However, since this study already considered openness in the direction of the wind (see directional viewsheds in the next section), the frontal area index was excluded from the analysis.
This study assumed that the 100-m buffer around each 50-m road-segment, covering an area of 41,394 m2(equivalent to approximately eight American football fields) was sufficient to study the effect of buildings’ morphology on PM2.5 in 50-m road-segments. Therefore, buffers of 100 m around each 50-m segment delineated the proximity around roads to extract buildings’ characteristics for calculating morphological parameters.
Depending on the spatial relation of buildings with the 100 m buffer, we categorized buildings into three classes:
(1)
Buildings with their centroid falling within the buffer
(2)
Part of the building/buildings intersecting with the buffer
(3)
Buildings within or touching the buffer
Figure 5 depicts buildings in all three classes highlighted in orange used to compute various morphological measures of the buildings in relation to a 100-m buffer around each 50-m road-segment. Buildings with their centroid falling within the buffer (class 1) were used to compute mean building height, mean building area, and the number of buildings per 1000 m2. This was to avoid considering the buildings with only a small portion in the buffer. In order to compute the building coverage ratio, this study considered the portions of the building intersecting with the buffer (class 2). Furthermore, to determine the average nearest-neighbor distance between buildings in the buffer, this study considered all the buildings either within the buffer or touching the buffer (class 3), as a partial account of buildings would not reflect the true distance between them. Table 1 provides the relevance and computation of each of these parameters:
Table 2 summarizes the building morphological characteristics along the data collection paths. Among all the paths, the building coverage ratios of paths A and E were smaller. Many small buildings surrounded path A. In contrast, a few bigger and taller buildings surrounded path E. Path C had the largest site coverage ratio, followed by path D, and followed by path B. Still, on average, buildings around paths C and D were bigger in size, taller in height, and fewer in number compared to path B. In short, path A constituted smaller, shorter, and denser buildings with intermediate building coverage. Path B consisted of moderately dense, medium-sized buildings with median height and median building coverage. Path C was composed of relatively low-density, large-sized, and taller buildings but had the highest building coverage. Path D contained the largest-sized low-density and distantly placed buildings with medium height and relatively high building coverage. Finally, path E was surrounded by small-sized, sparsely placed and short buildings. Path E also had the lowest building coverage among all paths.

2.3.3. Directional Viewshed

The amount of open area around a location could impact the PM2.5 value observed at a location. Brown et al. [11] provided evidence for the effect of the interaction between wind direction and site characteristics on PM2.5. Thus, this study calculated the open area in the direction of the wind for each segment centroid (Figure 6a) at varying distances of 100 m, 200 m, 400 m, 800 m, and 1500 m by applying viewshed analysis to determine the amount of area visible from any given location in all directions. The viewshed analysis considered 12 different wind directions, starting with 0° in an increment of 30°, resulting in viewsheds at 15°, 45°, 75°, …, 345°. These directional viewsheds served as the basis to analyze the interactions of PM2.5 and openness based on the corresponding wind direction.

2.3.4. Data Integration

In each data collection run, the mean PM2.5 for each 50 m segment was obtained by averaging all the PM2.5 data points associated with the respective segment. Based on the timestamp of the segment-wise PM2.5, the weather data closest in time was combined with the segment-wise PM2.5 data. Further, based on the wind direction, directional viewshed data were integrated with PM2.5 and weather data. In combination, the study datasets consisted of the segment-wise PM2.5, weather-related variables, and directional viewshed. The resulting dataset was further combined with the corresponding building morphological data.

2.4. Modeling

The data collected in the study formed spatial panel data because it consisted of the same set of locations observed multiple times (n = 9). The study included time-varying variables to explain temporal variations in the PM2.5 observations. Location-specific variables, such as building morphological characteristics, would not change over time, and they would have fixed effects on the response variable. Therefore, the study built a fixed-effects panel model to assess the effect of time-varying variables on PM2.5 and extract the individual or location-specific fixed effects of time-invariant variables.
The equation for the fixed-effects panel model can be expressed as follows:
Yit = α𝑖 + β × Xit + eit
where i refers to a fixed location; t refers to the time at which an observation was collected; Yit is an observed PM2.5 at location i and time t; Xit is a k × 1 vector where k is the number of dependent variables specified in Equation (3); β is a 1 × k vector of parameters; e i t is the error term. Intercept α i represents a fixed effect that is time-invariant but varied across locations. Following fixed-effects models, the study estimated coefficient β by “de-meaning”, which removed the average over time from each observation:
Y i t Y i ¯ = α i α i ¯ + β × X i t X i ¯ + e i t e i ¯  
where Y i ¯ = 1 T t = 1 T Y i t , X i ¯ = 1 T t = 1 T X i t , and e i ¯ = 1 T t = 1 T e i t .
The variables with bars (i.e., ( Y i ¯ ,   α i ¯ , . .   ) represent temporal means. De-meaning removes the fixed effect α i since a constant equates its mean over time ( α i ¯ ). Thus, only time-varying effects remain in the equation. The study was subjected to the drawback of the fixed-effects model that the model could not include individual time-invariant covariates in the model because de-meaning the estimated α i (i.e., the sum of all fixed effects) effectively removed the time-invariant observations.
Specifically, this study used the fixed-effects model to assess the effect of weather-related variables and the openness in the direction of wind on PM2.5. Although openness in the wind direction was not a time-varying variable, this variable might not exhibit a full time-invariant nature since its value changed according to the wind direction. Wind direction was a circular variable; therefore, it was used in the form of sine and cosine components in the model to account for its cyclical nature. Another wind-related variable resulted from the interaction of wind direction and the inlet of the DR1000 instrument. When the wind blew into the instrument’s inlet, DR1000 recorded an elevated PM2.5 compared to the situation when the wind blew away from the inlet. Therefore, in addition to weather-related variables and directional viewshed, this fixed-effects model included the angle between the instrument travel direction and wind direction. When the wind blew into the inlet, this angle was 180°, whereas when it blew away from the inlet, it was 0°. In other cases, the angle varied between 0° to 180°. Consequently, PM2.5 maxima coincided with the angle at 180° and minima at 0°, confirming the cosine of this angle as the appropriate measure to include in the model.
Considering all the variables explained above, the fixed-effects model in Equation (1) was specified as Equation (3) to fit the panel data model using the R package ‘plm’ (version 2.6.4) [26]. Further, this study also tested spatial autocorrelation in the errors for each data collection run. The distance used to determine neighbors was set to 51.5 m, ensuring each data point had at least one neighbor. The weight matrix for computing spatial autocorrelation was derived using the ‘spdep’ package in R (version 1.3.4) [27], and spatial autocorrelation (Moran’s I) using the ‘ape’ package (version 5.8) [28].
P M 2.5 i t = α i + W i n d   S p e e d i t + C o s   W i n d   D i r i t + S i n   W i n d   D i r i t + T e m p e r a t u r e i t + R e l a t i v e   H u m i d i t y i t + C o s i n e   o f   t h e   a n g l e   b e t w e e n   t r a v e l   d i r e c t i o n   a n d   w i n d   d i r e c t i o n i t + D i r e c t i o n a l   V i e w s h e d   i t
The fixed effects in the model represented location-specific characteristics that were not included in the model. The study related these fixed effects with building morphological characteristics: building coverage ratio, mean building area, number of buildings, average distance between nearest-neighbor buildings, and mean building height. The study assumed that the key time-invariant factor was building morphology and correlated the estimated fixed effects with time-invariant building morphological parameters. In order to assess the contribution of building morphological characteristics, the study developed a regression model relating building morphological characteristics with the fixed effects.

3. Results

3.1. PM2.5 and Weather Data

Figure 7a shows the distribution of PM2.5 across all segments and all data collection runs. The PM2.5 data were positively skewed and adjusted with a Box–Cox transformation (Figure 7b) in order to use them in the proposed fixed-effects panel model. The Box–Cox transformation parameter, lambda, was −0.6057.
Table 3 provides the time, average PM2.5, and weather data for each data collection run. Across the runs, the maximum and minimum observed PM2.5 values were 13.10 μg/m3 and 4.21 μg/m3, respectively. There was a slight variation between the maximum and minimum PM2.5 recorded in each run, ranging from 1.2 μg/m3 to 4.5 μg/m3.

3.2. Variation in Runs across Paths

Figure 8 shows the PM2.5 distribution during each run. While the study focused on five data-collection paths, some runs covered additional nearby segments. All data collection runs encountered low variations in PM2.5 measurements. Considering the interquartile range, runs five and six exhibited comparably higher variations. Figure 9 shows PM2.5 distributions at individual data collection paths during all data collection runs. When the wind came from the south, the PM2.5 values along all paths were highly similar (runs three and four). All the paths oriented in the north–south direction with open passages in the same orientation could facilitate dispersion and explain the similar PM2.5 observations.
While variations in PM2.5 observations across all paths appeared low, there were a few exceptions. During northwest winds (runs five, six, and eight), path E collected higher PM2.5 values compared to the other paths. A possible explanation was that the open area to the north of path E allowed the efficient dispersion of pollutants from surrounding areas, and the pollutants were subsequently trapped by large buildings along this path. In other cases, when the wind came from different directions, these buildings shielded the incoming airflow. In general, path A had low PM2.5 values, except when the wind came from the northeast (run two) across a sizable parking lot. A higher wind speed led to better dispersion. At the wind speed of 7.8 mph (run seven), the highest wind speed observed among all the data collection runs, PM2.5 values were uniform across all paths.

3.3. Results from the Fixed-Effects Model

This study calculated the directional viewshed at multiple distances to account for the unknown effective distance over which openness in a viewshed would influence PM2.5 dispersion using Equation (3). An experiment using five fixed-effects models, each with an effective distance, examined the distance effect on model performance. All other weather-related variables remained the same in all five models. The results showed comparable model performances (R2) at varying viewshed distances (Table 4).
Nevertheless, small viewshed distances up to 200 m were not statistically significant in the model. The model R2 value increased up to 800 m, and then it started decreasing. Besides the directional viewshed, all other weather-related variables used in the model appeared to be statistically significant. Table 5 summarizes the detailed results of the model (Equation (3)) built with a viewshed of 800 m.
The small transformed PM2.5 values (0.84–1.32) in the dependent variable contributed to the small coefficients for all explanatory variables. Wind speed had a negative impact on PM2.5 as the model asserted a negative coefficient. Both the components of the wind direction, sine and cosine, had negative coefficients in the model, whereas both temperature and humidity had positive coefficients. This study used a digital surface model with a 50-cm resolution to determine the directional viewshed of a location, and 50 × 50 cm2 cells in a viewshed up to 800 m were used in the regression model. Depending on the location’s visibility, the number of cells in the viewshed varied from hundreds to thousands, which explained the small coefficient for directional viewsheds. The amount of open area in the wind direction had a small but positive impact on PM2.5.

3.3.1. Error Analysis

Moran’s I was calculated to evaluate the spatial autocorrelation in the errors for each data collection run. All but run three had a significant moderate to strong spatial autocorrelation in the errors (Table 6), suggesting the possible omission of variables responsible for the PM2.5 variations across the study area. To account for this spatial dependence, this study considered both spatial lag models and spatial error models. Spatial lag models examine the existing autocorrelation in the response variable [29]. As the study intended to mediate the model bias due to spatially autocorrelated errors and investigate the effect of predictors on PM2.5 estimates, the study opted for the spatial error modeling approach. Spatial autocorrelation in the errors increased with spatial specification in each panel. Moran’s I values and their significance to the errors observed in each data collection round are outlined in Table 6, evaluated within both the panel model (aspatial specification) and the spatial panel model. The results suggested that aspatial models suffered less than spatial models from spatial dependence in the errors, as indicated by Moran’s I. Therefore, the study selected the fixed-effects panel model over the fixed-effects spatial panel data model for further analysis.

3.3.2. Relation between Fixed Effects and Building Morphology

The correlation between the fixed effects and building morphological characteristics varied from 0.79 to −0.51 (Figure 10). The fixed effects were weak to moderately correlated with all building morphological characteristics, except for the average nearest-neighbor distance between buildings. Mean building area, mean building height, and building coverage ratio were strongly positively correlated with each other. On the other hand, while the average nearest-neighbor distance between buildings was only weakly correlated with all the other building morphological characteristics, the number of buildings was moderately negatively correlated with mean building height.
As such, the study excluded the average nearest-neighbor distance between buildings in regression modeling, as it was almost uncorrelated with the fixed effects. Table 7 presents the regression results. The study found only the building coverage ratio and the number of buildings to be significant variables that explained variations in the fixed effects. Overall, building morphological characteristics explained 33.22% variation in the fixed effects.

4. Discussion

Despite the moderate variation in building morphology, the study observed relatively uniform PM2.5 levels within an area of about 0.5 km2 across all data collection runs. The consistent PM2.5 levels observed in this study agree with the findings of a study by Harrison et al. [2] aimed at understanding spatial variability at a neighborhood scale. Harrison et al. [2] collected daily on-road PM2.5 measurements over a month within an area of 100 km2 (equivalent to MODIS aerosol product’s single pixel), encompassing the University of Texas at Dallas. Their study found that depending on the weather condition, the spatial scale of PM2.5 variation in the area varied from 0.8 km to 5.2 km, suggesting uniform PM2.5 levels over a distance of 0.8 km or smaller.
The building morphology accounted for 33.22% of the variations in the fixed effects, suggesting that the built environment could affect PM2.5. Regression results from the fixed effects showed that the building coverage ratio had a positive impact on PM2.5. Thus, the greater built-up area likely hindered the dispersion of the pollutants. Despite the negative impact of the number of buildings per 1000 m2 on PM2.5, a similar inference could be drawn. A negative correlation between the number of buildings per 1000 m2 and mean building area, mean building height, and building coverage ratio in the study area denoted the contrast between areas with many smaller buildings and areas with fewer larger buildings, suggesting areas with smaller and shorter buildings experienced better dispersion than other areas. Nevertheless, the study could not assess each variable in isolation, as an inference about it could change based on its relationships with other variables in the model. The insignificant impact of mean building height and mean building area could be attributed to limited variability among building morphological characteristics in the study area.
The negative effect of wind speed on PM2.5 in the model reaffirmed findings from the previous studies that increasing wind speed improved dispersion and hence reduced the PM2.5 concentration [10,11]. Increased temperature promoted air circulation and was negatively related to PM2.5 measures [10]. Therefore, the temperature was expected to have a negative effect on PM2.5 contrary to the positive effect found in this study. Several factors could contribute to this discrepancy. One possible explanation is that elevated temperature contributed to chemical reactions that might have formed new particles including PM2.5 [30]. Other reasons could be the blocking of the airflow by surrounding buildings or the desorption of particles from different surfaces due to increased temperature. As for the relative humidity, depending on its value, it could increase or decrease in PM2.5 mass. Increased humidity promoted particle size growth, but after a certain threshold, very high humidity could lead to particle deposition due to heavy particle growth [30]. Moreover, Lou et al. [31], in their 3-year study in the Yangtze River delta, described the relationship between relative humidity and PM2.5 as an inverted U-shape, where the hygroscopic growth of particles continued with an increase in relative humidity until it reached 70%, and then it started decreasing. The relative humidity observed in the study ranged from 39 to 70%—not high enough to cause particle deposition—and resulted in a positive correlation with PM2.5. The positive coefficient associated with directional viewshed suggested that a more open area in the direction of the wind was consistent with additional incoming pollution from the other areas, as confirmed by Brown et al. [11]. Of the different viewshed distances used in the model, model performance consistently increased up to 800 m distance, and then it decreased, suggesting openness up to 800 m distance affected PM2.5. However, increasing distance led to only marginal improvements in the model, and openness within 400 m also yielded satisfactory results.

5. Conclusions

The study collected the PM2.5 data using a mobile monitoring platform at a non-near-road site with different building morphological characteristics. Building morphological characteristics varied from high-density, small-sized, shorter buildings to low-density, medium to large-sized, taller buildings. The data collection strategy included multiple runs across paths representing varied building morphology. The fixed-effects panel model was used to investigate the effect of weather-related variables and building morphological characteristics on PM2.5. Unlike meteorological variables, building morphological characteristics did not change with time. A regression model was developed to find the contribution of building morphological characteristics to fixed effects extracted from the panel data model. While the weather-related variables explained variations in PM2.5, building morphological characteristics also showed positive effects on PM2.5. Furthermore, openness in the direction of the wind allowed pollutants from other areas and raised PM2.5 concentration in the area.
The model presented in the study carried spatially autocorrelated errors even with spatial specifications, possibly due to the small study area with low PM 2.5 variability. A larger study extent and coarser unit of spatial analysis might mediate the issues of spatially autocorrelated errors. Moreover, the fixed-effects model applied the same spatial weight matrix across all runs. Depending on the weather conditions and its interactions with the surrounding built environment, the nature of spatial dependence could change. Therefore, a dynamic spatial weight matrix would be more appropriate to account for this spatial dependence. Alternatively, space–time convolution could also address the complex space–time dependence in the data overlooked in the current model specification. Nevertheless, this study showed small spatial variation in PM2.5 over a small area (<1 km2), typical of non-near-road sites with moderate variation in building morphology. The study considered interactions between wind and openness in the wind direction and overall building morphological characteristics within a 100 m buffer. Future studies can incorporate interactions between wind and building heights and measures related to the spatial arrangement of buildings to further understand the effects of buildings on PM2.5. Furthermore, variation in PM2.5 depends on many factors like weather parameters and emission sources that vary across different parts of the year. All the scenarios cannot be explored at a single site with a limited number of runs. More sites with different urban settings need to be studied over an extended period to improve our understanding of the microscale dynamics of PM2.5.

Author Contributions

Conceptualization, M.Y. and Y.K.; methodology, Y.K.; software, Y.K.; validation, Y.K.; formal analysis, Y.K.; investigation, Y.K.; resources, M.Y.; data curation, Y.K.; writing—original draft preparation, Y.K.; writing—review and editing, Y.K. and M.Y.; visualization, Y.K.; supervision, M.Y.; funding acquisition, M.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the University of Texas STAR program through the University of Texas at Dallas. The work was finalized when M.Y. was serving at the US National Science Foundation. Any opinions, findings, and conclusions or recommendations expressed in this material are those of the author(s) and do not necessarily reflect the views of the University of Texas at Dallas or the US National Science Foundation.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors on request.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Environmental Protection Agency. 40 CFR Appendix E to Part 58—Probe and Monitoring Path Siting Criteria for Ambient Air Quality Monitoring. In Code of Federal Regulations; U.S. Government Publishing Office: Dallas, TX, USA, 2019; pp. 310–317. Available online: https://www.govinfo.gov/app/details/CFR-2019-title40-vol6/CFR-2019-title40-vol6-part58-appE/summary (accessed on 1 June 2024).
  2. Harrison, W.A.; Lary, D.; Nathan, B.; Moore, A.G. The Neighborhood Scale Variability of Airborne Particulates. J. Environ. Prot. 2015, 6, 464–476. [Google Scholar] [CrossRef]
  3. Shi, Y.; Lau, K.K.L.; Ng, E. Developing Street-Level PM2.5 and PM10 Land Use Regression Models in High-Density Hong Kong with Urban Morphological Factors. Environ. Sci. Technol. 2016, 50, 8178–8187. [Google Scholar] [CrossRef] [PubMed]
  4. Hankey, S.; Marshall, J.D. Land Use Regression Models of On-Road Particulate Air Pollution (Particle Number, Black Carbon, PM2.5, Particle Size) Using Mobile Monitoring. Environ. Sci. Technol. 2015, 49, 9194–9202. [Google Scholar] [CrossRef] [PubMed]
  5. Li, Z.; Fung, J.C.H.; Lau, A.K.H. High Spatiotemporal Characterization of On-Road PM2.5 Concentrations in High-Density Urban Areas Using Mobile Monitoring. Build. Environ. 2018, 143, 196–205. [Google Scholar] [CrossRef]
  6. Zhou, S.; Lin, R. Spatial-Temporal Heterogeneity of Air Pollution: The Relationship between Built Environment and on-Road PM2.5 at Micro Scale. Transp. Res. D Transp. Environ. 2019, 76, 305–322. [Google Scholar] [CrossRef]
  7. Harrison, W.A.; Lary, D.J.; Nathan, B.J.; Moore, A.G. Using Remote Control Aerial Vehicles to Study Variability of Airborne Particulates. Air Soil Water Res. 2015, 8, 43–51. [Google Scholar] [CrossRef]
  8. Peng, Z.R.; Wang, D.; Wang, Z.; Gao, Y.; Lu, S. A Study of Vertical Distribution Patterns of PM2.5 Concentrations Based on Ambient Monitoring with Unmanned Aerial Vehicles: A Case in Hangzhou, China. Atmos. Environ. 2015, 123, 357–369. [Google Scholar] [CrossRef]
  9. Vos, P.E.J.; Maiheu, B.; Vankerkom, J.; Janssen, S. Improving Local Air Quality in Cities: To Tree or Not to Tree? Environ. Pollut. 2013, 183, 113–122. [Google Scholar] [CrossRef] [PubMed]
  10. Ginzburg, H.; Liu, X.; Baker, M.; Shreeve, R.; Jayanty, R.K.M.; Campbell, D.; Zielinska, B. Monitoring Study of the Near-Road PM2.5 Concentrations in Maryland. J. Air Waste Manag. Assoc. 2015, 65, 1062–1071. [Google Scholar] [CrossRef] [PubMed]
  11. Brown, S.G.; Penfold, B.; Mukherjee, A.; Landsberg, K.; Eisinger, D.S. Conditions Leading to Elevated PM2.5 at near-Road Monitoring Sites: Case Studies in Denver and Indianapolis. Int. J. Environ. Res. Public. Health 2019, 16, 1634. [Google Scholar] [CrossRef] [PubMed]
  12. Edussuriya, P.; Chan, A.; Malvin, A. Urban Morphology and Air Quality in Dense Residential Environments: Correlations between Morphological Parameters and Air Pollution at Street-Level. J. Eng. Sci. Technol. 2014, 9, 64–80. [Google Scholar]
  13. Yang, J.; Shi, B.; Zheng, Y.; Shi, Y.; Xia, G. Urban Form and Air Pollution Disperse: Key Indexes and Mitigation Strategies. Sustain. Cities Soc. 2020, 57, 101955. [Google Scholar] [CrossRef]
  14. Wang, Z.; Zhong, S.; He, H.-d.; Peng, Z.R.; Cai, M. Fine-Scale Variations in PM2.5 and Black Carbon Concentrations and Corresponding Influential Factors at an Urban Road Intersection. Build. Environ. 2018, 141, 215–225. [Google Scholar] [CrossRef]
  15. Deshmukh, P.; Isakov, V.; Venkatram, A.; Yang, B.; Zhang, K.M.; Logan, R.; Baldauf, R. The Effects of Roadside Vegetation Characteristics on Local, near-Road Air Quality. Air Qual. Atmos. Health 2019, 12, 259–270. [Google Scholar] [CrossRef] [PubMed]
  16. Point Source Emission Inventory. Available online: https://www.tceq.texas.gov/airquality/point-source-ei/psei.html (accessed on 1 June 2024).
  17. U.S. Environmental Protection Agency. 2017 National Emissions Inventory (NEI) Data. Available online: https://www.epa.gov/air-emissions-inventories/2017-national-emissions-inventory-nei-data (accessed on 30 May 2024).
  18. Traffic Count Program Annual Report; 2019. Available online: https://www.cor.net/home/showdocument?id=28259 (accessed on 7 April 2021).
  19. Scentroid DR1000 Flying Lab. Available online: https://scentroid.com/products/analyzers/dr1000-flying-lab/ (accessed on 2 June 2024).
  20. Palas GmbH Fidas® Frog. Available online: https://www.palas.de/en/product/fidasfrog (accessed on 31 May 2024).
  21. Wijeratne, L. Coupling Physical Measurement with Machine Learning for Holistic Environmental Sensing. Ph.D. Thesis, University of Texas at Dallas, Richardson, TX, USA, 2021. [Google Scholar]
  22. Fernando, B.A.; Talebi, S.; Wijeratne, L.; Waczak, J.; Sooriyaarachchi, V.; Ruwali, S.; Hathurusinghe, P.; Lary, D.; Sadler, J.; Lary, T.; et al. Gauging Size Resolved Ambient Particulate Matter Concentration Solely Using Biometric Observations: A Machine Learning and Causal Approach. Med. Res. Arch. 2024, 12. [Google Scholar] [CrossRef]
  23. Talebi, S.; Lary, D.J.; Wijeratne, L.O.H.; Fernando, B.; Lary, T.; Lary, M.; Sadler, J.; Sridhar, A.; Waczak, J.; Aker, A.; et al. Decoding Physical and Cognitive Impacts of Particulate Matter Concentrations at Ultra-Fine Scales. Sensors 2022, 22, 4240. [Google Scholar] [CrossRef] [PubMed]
  24. Microsoft (Licensed under the Open Data Commons Open Database License (ODbL)) U.S. Building Footprints. Available online: https://github.com/Microsoft/USBuildingFootprints (accessed on 1 November 2020).
  25. U.S. Geological Survey. Lidar Point Cloud—USGS National Map 3DEP Downloadable Data Collection; U.S. Geological Survey: Reston, VA, USA, 2023.
  26. Croissant, Y.; Millo, G. Panel Data Econometrics in R: The Plm Package. J. Stat. Softw. 2008, 27, 1–43. [Google Scholar] [CrossRef]
  27. Bivand, R.; Pebesma, E.; Gómez-Rubio, V. Applied Spatial Data Analysis With R, 2nd ed.; Springer: New York, NY, USA, 2013. [Google Scholar]
  28. Paradis, E.; Schliep, K. Ape 5.0: An Environment for Modern Phylogenetics and Evolutionary Analyses in R. Bioinformatics 2019, 35, 526–528. [Google Scholar] [CrossRef] [PubMed]
  29. Anselin, L. Spatial Econometrics. In A Companion to Theoretical Econometrics; Blackwell Publishing Ltd.: Malden, MA, USA, 2003; pp. 310–330. [Google Scholar]
  30. Wang, J.; Ogawa, S. Effects of Meteorological Conditions on PM2.5 Concentrations in Nagasaki, Japan. Int. J. Environ. Res. Public. Health 2015, 12, 9089–9101. [Google Scholar] [CrossRef] [PubMed]
  31. Lou, C.; Liu, H.; Li, Y.; Peng, Y.; Wang, J.; Dai, L. Relationships of Relative Humidity with PM2.5 and PM10 in the Yangtze River Delta, China. Environ. Monit. Assess. 2017, 189, 582. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Study area context: UTD location (on the boundary of Dallas and Collin counties), metroplex and its counties, and PM2.5 emissions from point sources.
Figure 1. Study area context: UTD location (on the boundary of Dallas and Collin counties), metroplex and its counties, and PM2.5 emissions from point sources.
Applsci 14 05149 g001
Figure 2. PM2.5 data collection paths A, B, C, D, and E; pedestrian-friendly area within the UTD campus. The figure also depicts building footprints and major roads (Waterview Parkway, W Campbell, N Floyd, and Synergy Park Blvd) surrounding the campus. Brown road labels represent the average annual 24 h traffic count including peak-hour traffic data (in the brackets) from the City of Richardson’s 2019 Annual Report for Traffic Count Program [18].
Figure 2. PM2.5 data collection paths A, B, C, D, and E; pedestrian-friendly area within the UTD campus. The figure also depicts building footprints and major roads (Waterview Parkway, W Campbell, N Floyd, and Synergy Park Blvd) surrounding the campus. Brown road labels represent the average annual 24 h traffic count including peak-hour traffic data (in the brackets) from the City of Richardson’s 2019 Annual Report for Traffic Count Program [18].
Applsci 14 05149 g002
Figure 3. Linear regression between reference instrument (Fidas Frog) and DR1000 from Scentroid, a portable device measuring PM2.5 (a) before calibration and (b) after calibration.
Figure 3. Linear regression between reference instrument (Fidas Frog) and DR1000 from Scentroid, a portable device measuring PM2.5 (a) before calibration and (b) after calibration.
Applsci 14 05149 g003
Figure 4. The Flowchart for data preparation.
Figure 4. The Flowchart for data preparation.
Applsci 14 05149 g004
Figure 5. Buffer and buildings used in building morphology calculation. The orange color indicates that the buildings or their parts meet criteria for specified class (class 1, 2 or 3).
Figure 5. Buffer and buildings used in building morphology calculation. The orange color indicates that the buildings or their parts meet criteria for specified class (class 1, 2 or 3).
Applsci 14 05149 g005
Figure 6. Buildings’ morphological characteristics along data collection paths: (a) 50 m segments along data collection paths, segment centroids, and a 100 m buffer surrounding each segment with blue buffers of road segments along each collection path. As the buffers of road segments overlap, a single segment and its buffer boundary are highlighted in red for illustration. Buildings’ morphological characteristics for each segment are computed considering a surrounding 100 m buffer area and are displayed in (bf) using segment centroids: (b) mean building area (meters), (c) the number of buildings per 1000 m2 (d) mean building height (meters), (e) an average distance between nearest-neighbor building (meters), and (f) building coverage ratio.
Figure 6. Buildings’ morphological characteristics along data collection paths: (a) 50 m segments along data collection paths, segment centroids, and a 100 m buffer surrounding each segment with blue buffers of road segments along each collection path. As the buffers of road segments overlap, a single segment and its buffer boundary are highlighted in red for illustration. Buildings’ morphological characteristics for each segment are computed considering a surrounding 100 m buffer area and are displayed in (bf) using segment centroids: (b) mean building area (meters), (c) the number of buildings per 1000 m2 (d) mean building height (meters), (e) an average distance between nearest-neighbor building (meters), and (f) building coverage ratio.
Applsci 14 05149 g006
Figure 7. Distribution of PM2.5. (a) Raw PM2.5 observations. (b) Transformed PM2.5 data.
Figure 7. Distribution of PM2.5. (a) Raw PM2.5 observations. (b) Transformed PM2.5 data.
Applsci 14 05149 g007
Figure 8. Distribution of PM2.5 in the study area during each data collection run, along with corresponding wind direction indicated by arrows at the top.
Figure 8. Distribution of PM2.5 in the study area during each data collection run, along with corresponding wind direction indicated by arrows at the top.
Applsci 14 05149 g008
Figure 9. PM2.5 at each path during data collection runs, and the corresponding wind direction indicated by an arrow at the top.
Figure 9. PM2.5 at each path during data collection runs, and the corresponding wind direction indicated by an arrow at the top.
Applsci 14 05149 g009
Figure 10. Correlation between fixed effects and building morphological variables.
Figure 10. Correlation between fixed effects and building morphological variables.
Applsci 14 05149 g010
Table 1. Building morphology calculation and relevance.
Table 1. Building morphology calculation and relevance.
ParameterFormulaRelevance
Mean Building AreaThe sum of areas of all building footprints in a 100-m buffer/Number of buildings in 100-m bufferMeasures the average size of buildings. Larger buildings provide greater horizontal enclosure, hindering PM2.5 dispersion, and vice versa.
Number of Buildings per 1000 Square MetersNumber of buildings in a 100-m buffer × 1000/Area of buffer in square metersStandardized to number of buildings per 1000 m2 to account for buffer size variations. The number of buildings along with the mean building area indicates the amount of built-up area. A large number of large-sized buildings occupy more space and hence leave less room for dispersion, and vice versa.
Mean Building HeightThe sum of heights of all buildings in a 100-m buffer/Number of buildings in a 100-m bufferTaller buildings provide greater vertical enclosure.
Average Nearest-Neighbor Distance between BuildingsAverage nearest-neighbor distance between buildings. Computed using the generate near table tool from ArcGIS Pro 2.5Reflects building density. Tightly packed buildings hinder dispersion; greater distance between buildings allows more room for dispersion.
Building Coverage Ratio Area of a 100-m buffer occupied by buildings/Area of a 100-m bufferIndicates the proportion of buffer area covered by buildings. A greater coverage ratio limits room for dispersion; and vice versa.
Table 2. Summary of building morphological characteristics around data collection paths.
Table 2. Summary of building morphological characteristics around data collection paths.
PathBuilding Coverage RatioMean Building Area
(m2)
Number of Buildings per 1000 m2Average Distance
between Nearest-Neighbor Buildings (m)
Mean Building Height
(m)
Path A0.198860.23126.70
Path B0.3328370.1299.70
Path C0.3841320.071314.80
Path D0.3656610.061412.00
Path E0.1625200.04168.40
Table 3. Details of the date and time of data collection rounds.
Table 3. Details of the date and time of data collection rounds.
Run Date and TimePM2.5 (µg/m3)Wind Speed
(mph)
Wind
Direction
Temperature
(Degrees
Fahrenheit)
Relative
Humidity
(%)
MeanRangeInterquartile Range
110 December 2019 13:32:164.214.000.603.859044.5039
210 December 2019 17:42:525.202.000.904.425850.0042
312 December 2019 11:02:426.202.000.607.0017451.0045
412 December 2019 14:23:487.402.751.103.4018742.5053
513 December 2019 11:21:2613.103.001.002.0029675.6051
613 December 2019 17:22:5211.504.501.704.4027065.0058
714 December 2019 11:26:385.001.200.407.808749.0054
815 December 2019 16:35:228.403.400.805.6029740.0070
921 February 2020 15:24:035.002.000.302.6025020.0047
Table 4. Fixed-effects model results with varying viewshed distance.
Table 4. Fixed-effects model results with varying viewshed distance.
Directional Viewshed DistanceR2Is Directional Viewshed Significant?
100 m0.82496No
200 m0.82516No
400 m0.82658Yes
800 m0.82812Yes
1500 m0.82778Yes
Table 5. Fixed-effects model results.
Table 5. Fixed-effects model results.
VariableCoefficientp-value
Wind speed−0.00038942.358 × 10−6
Cosine of wind direction−0.0142960.0001239
Sine of wind direction−0.062830<2.2 × 10−16
Temperature0.006214<2.2 × 10−16
Relative humidity0.0040963<2.2 × 10−16
Cosine of angle between travel direction and wind direction−0.018950.0024263
Directional viewshed (800 m)3.3831 × 10−70.0013893
Table 6. Spatial autocorrelation in errors of the fixed-effects models.
Table 6. Spatial autocorrelation in errors of the fixed-effects models.
Data
Collection Run
Number of ObservationsMoran’s I
(Panel Data)
SignificanceMoran’s I
(Spatial Panel Data Model)
Significance
1570.47Significant0.68Significant
2570.87Significant0.90Significant
3570.14Insignificant0.30Significant
4570.31Significant0.26Significant
5570.31Significant0.55Significant
6570.48Significant0.80Significant
7570.43Significant0.75Significant
8570.40Significant0.68Significant
9570.48Significant0.59Significant
Table 7. Results of regression between fixed effects and building morphological characteristics.
Table 7. Results of regression between fixed effects and building morphological characteristics.
VariableCoefficientp-Value
Mean building area−2.33 × 10−60.1939
Building coverage ratio0.07310.0046 **
Mean building height−0.001620.1329
Number of buildings per 1000 m2−0.10990.0001 ***
Model R2: 0.3322
** and *** denote significance at the level of 0.01 and 0.0001, respectively.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Karale, Y.; Yuan, M. How May Building Morphology Influence Pedestrians’ Exposure to PM2.5? Appl. Sci. 2024, 14, 5149. https://doi.org/10.3390/app14125149

AMA Style

Karale Y, Yuan M. How May Building Morphology Influence Pedestrians’ Exposure to PM2.5? Applied Sciences. 2024; 14(12):5149. https://doi.org/10.3390/app14125149

Chicago/Turabian Style

Karale, Yogita, and May Yuan. 2024. "How May Building Morphology Influence Pedestrians’ Exposure to PM2.5?" Applied Sciences 14, no. 12: 5149. https://doi.org/10.3390/app14125149

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

Article Metrics

Back to TopTop