Restaurant Impacts on Outdoor Air Quality: Elevated Organic Aerosol

Aug 2, 2018 - Organic aerosol (OA) is a major component of fine particulate matter (PM2.5) in urban environments. We performed in-motion ambient ...
0 downloads 0 Views 2MB Size
Article Cite This: Environ. Sci. Technol. XXXX, XXX, XXX−XXX

pubs.acs.org/est

Restaurant Impacts on Outdoor Air Quality: Elevated Organic Aerosol Mass from Restaurant Cooking with Neighborhood-Scale Plume Extents Ellis Shipley Robinson,†,‡ Peishi Gu,†,‡ Qing Ye,§,∥,‡ Hugh Z. Li,†,‡ Rishabh Urvesh Shah,†,‡ Joshua Schulz Apte,⊥ Allen L. Robinson,†,‡,∥ and Albert A. Presto*,‡,† †

Department of Mechanical Engineering, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States Center for Atmospheric Particle Studies, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States § Department of Chemistry, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States ∥ Department of Engineering & Public Policy, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States ⊥ Department of Civil, Architectural & Environmental Engineering, University of Texas at Austin, Austin, Texas 78705, United States

Downloaded via UNIV OF CAMBRIDGE on August 2, 2018 at 14:09:19 (UTC). See https://pubs.acs.org/sharingguidelines for options on how to legitimately share published articles.



S Supporting Information *

ABSTRACT: Organic aerosol (OA) is a major component of fine particulate matter (PM2.5) in urban environments. We performed in-motion ambient sampling from a mobile platform with an aerosol mass spectrometer (AMS) to investigate the spatial variability and sources of OA concentrations in Pittsburgh, Pennsylvania, a midsize, largely postindustrial American city. To characterize the relative importance of cooking and traffic sources, we sampled in some of the most populated areas (∼18 km2) in and around Pittsburgh during afternoon rush hour and evening mealtime, including congested highways, major local roads, areas with high densities of restaurants, and urban background locations. We found greatly elevated OA concentrations (10s of μg m−3) in the vicinity of numerous individual restaurants and commercial districts containing multiple restaurants. The AMS mass spectral information indicates that majority of the high concentration plumes (71%) were from cooking sources. Areas containing both busy roads and restaurants had systematically higher OA concentrations than areas with only busy roads and urban background locations. Elevated OA concentrations were measured hundreds of meters downwind of some restaurants, indicating that these sources can influence air quality on neighborhood scales. Approximately 20% of the population (∼250 000 people) in the Pittsburgh area lives within 200 m of a restaurant; therefore, restaurant emissions are potentially an important source of outdoor PM exposures for this large population.



modeling7,8 and land-use regression (LUR) modeling9 are used to describe spatial variability at scales of tens to hundreds of meters in urban environments. LUR requires measurements for model construction and validation, and so measurement strategies that resolve pollutant variation at high spatial resolution have been employed. Sampling from mobile platforms (e.g., measurements on-board a motor vehicle10 or bicycle11) or distributed networks of sensors12 are both commonly used to characterize spatial variation of pollutant concentrations. Distributed sampling can be cost- and laborintensive, and therefore typically focuses on a small number of pollutants (e.g., PM2.5 and NO2), often using passive samplers13 or integrated filter samples.14 More recently, lowcost sensor networks have also been deployed to quantify spatial variations of these pollutants in real time.15 An

INTRODUCTION Exposure to particulate matter with diameters less than 2.5 μm (PM2.5) is a global risk factor for death and other negative health effects.1,2 Organic aerosol (OA) is a major component to PM2.5 mass, especially in urban locations. Important sources of primary OA are vehicles and food-cooking. Despite the consequences of PM2.5 exposure for human health, there remain important gaps in our understanding of its spatial distribution, especially in complex urban areas.3 Even in the United States, where there has been extensive air quality monitoring for the past half-century,4 estimates using data from fixed monitors can fail to capture the variability of pollutants in source-rich environments, and thus fail to accurately assess exposure.5 Recent work has shown that for some pollutants, for example various traffic-related pollutants (e.g., nitric oxide, ultrafine particles), intraurban variability in concentrations may be much greater than interurban differences measured at urban background monitors.6 There have been substantial efforts to better understand the spatial patterns of air pollution at fine scales. Dispersion © XXXX American Chemical Society

Received: May 21, 2018 Revised: July 13, 2018 Accepted: July 20, 2018

A

DOI: 10.1021/acs.est.8b02654 Environ. Sci. Technol. XXXX, XXX, XXX−XXX

Article

Environmental Science & Technology

Figure 1. Map of sampling domain showing all grid cells containing mobile data (map created using ArcMap software v.10.4, ESRI, Inc.). Grid cells are 100 m × 100 m, and are colored according to the land-use classifications in the Methods section. The inset box shows the sampling area for the time series discussed in Figure 2. The CMU campus background site is marked with a yellow star.

different sources, we conducted significant sampling in urban background locations, on busy roads, and near commercial districts with substantial numbers of restaurants. Source attribution was performed using both spatial and mass spectral information. We quantified the prevalence of OA plumes from different source categories, as well as the spatial scale for which primary emissions from restaurants are elevated downwind of the source. These measurements illustrate how restaurants affect local concentrations of OA, and thus human exposure to outdoor PM2.5.

advantage of mobile sampling is the potential use of sophisticated and high time resolution instrumentation, which has led to increased understanding of the spatial scales affecting pollutant concentrations3 and can offer insight into the sources driving observed spatial patterns.16 The spatial patterns of traffic-related air pollution have been extensively studied.6,17 Species such as CO, NO, and particulate black carbon (BC) are directly emitted from onroad vehicles. Concentrations of these pollutants are elevated at the road edge and fall to background levels within a few hundred meters downwind of the roadway. Little work has been done to establish the spatial patterns of PM related to food cooking or other urban sources. This is surprising given that in many locations the organic fraction of PM (organic aerosol, OA) attributable to cooking is comparable to or greater than the contribution from fresh traffic emissions.18−22 A few studies have examined the spatial distribution of cooking-related PM in the context of outdoor air quality. Elser et al.23 used mobile sampling to show that cooking OA is enhanced in the city center of Talinn, Estonia relative to the city’s outskirts. They attributed the enhancement to restaurants. Similarly, Vert et al.24 found that particle number concentrations are elevated in front of restaurants in Utrecht compared to nonrestaurant locations. Abernethy et al.25 found restaurant density to be a significant predictor of ultrafine particle concentrations in Vancouver, Canada. These studies show that restaurants contribute to particulate air pollution. Relative to traffic-related pollutants though, our understanding of how restaurant cooking-related PM concentrations and exposures are spatially distributed remains low. In this work, we systematically investigated the spatial distribution of OA from food cooking and traffic in an urban environment (Pittsburgh, PA) by combining a spatially stratified sampling strategy with chemically resolved measurements. Our measurements of OA, conducted on a mobile platform with an Aerosol Mass Spectrometer (AMS), target different sources and land-use types, such as commercial areas containing restaurants, areas with high traffic, and urban background sites. In order to understand the contribution of



METHODS Mobile Sampling Strategy. We conducted mobile sampling to quantify the spatial variation in submicron mass concentration and composition of OA as part of the Center for Air, Climate, and Energy Solutions (CACES). Our sampling was stratified across the following land-use categories in Pittsburgh, PA and adjacent boroughs in Allegheny County: urban background locations, areas containing restaurants but not busy roads, areas containing busy roads but not restaurants, and areas containing both restaurants and busy roads. The sampling domain (18 km2, total land area of Pittsburgh ∼151 km2) and the distribution of these land-use categories is shown in Figure 1, and how we performed this categorization is discussed below. By sampling across these land-use types, we aimed to measure OA concentrations in locations influenced by both traffic and cooking activities, and quantify how cooking OA is spatially distributed in Pittsburgh. Because our focus is the spatial distribution of cooking-related OA, we sampled largely during evening mealtimes when cooking emissions from restaurants are likely at their peak. Roads were defined as “busy” if they were in the upper two quintiles of annual average daily traffic (AADT) counts of state-managed roads (AADT > 7058 vehicles/day). The source of road data and traffic volumes came from Pennsylvania Spatial Data Access (PASDA) database (http://www.pasda. psu.edu/uci/SearchResults.aspx?Shortcut=transportation). Restaurant locations were determined using a list of restaurants in Allegheny County from the Allegheny County Health B

DOI: 10.1021/acs.est.8b02654 Environ. Sci. Technol. XXXX, XXX, XXX−XXX

Article

Environmental Science & Technology

source attribution can be decoupled from location because of the chemical information provided by the AMS. Semistationary Sampling to Quantify the Spatial Extent of Plumes. For a limited number of periods, we used low-speed saturation sampling to quantify the spatial extent of elevated OA concentrations around small clusters of restaurants (containing 3−7 restaurants). These results are presented as “fall-off curves,” similar to those shown in Karner et al. for pollutant concentrations near roadways.17 In each of these four cases, the restaurant cluster was isolated from other potentially confounding OA sources (e.g., away from busy roads, bus routes, etc.). Downwind distance was calculated using the centroid of the restaurant cluster. We parked our mobile laboratory for 1−5 min periods, or drove at very slow speeds (0.82) and low 55:57 values ( 1.6 and similarity values between 0.6−0.9 were apportioned to cooking sources. Plumes that did not fit either of the above criteria are labeled “other.” AMS cooking spectra across a wide range of cooking methods and ingredients appear similar when grouped according to the above criteria. This idea is supported by recent work from Reyes-Villegas et al.37 who compared AMS spectra from a variety of cooked foods and cooking methods. They found similar mass spectra for fresh emissions from a variety of cooking sources. While AMS spectra may not be useful for distinguishing different cooking types of foods, the spectral differences between the broad POA categories (cooking, traffic, biomass-burning) are significant enough to accurately differentiate plumes into these classes. We tested our algorithm using reference spectra from the major sources of primary OA emissions from the AMS Spectral Database (http://cires1.colorado.edu/jimenez-group/AMSsd/ ): traffic, cooking, and biomass-burning. The reference spectra are from source-categorization experiments, smog chamber experiments, and ambient measurements. Ambient measurements included source-specific factors determined via positive matrix factorization (PMF; e.g. HOA, cooking OA, biomassburning OA). Near-source spectra from our mobile data set e.g., vehicle chase plumes or measurements collected directly outside of a restaurant exhaust ventwere used as well. Further details of all reference spectra (Table S2), and a visual presentation of the categorization method (Figure S1) are presented in the SI. As shown in Figure S1, all reference spectra (39 total), both from the AMS Spectral Database and our own source-profile sampling, were sorted correctly according to this empirical classification scheme, with the possible exception of one spectrum. The one exception was sampled downwind of a restaurant serving pizza and hamburgers from a wood-fired grill and oven. This plume was not sorted into the “cooking” category, but rather is located near biomass-burning plumes in this similarity vs 55:57 space. Two other data points from other wood-fired restaurants fall within the cooking category but are closer to

tions. All AMS data were corrected for residence time in the sampling inlet. AMS data points were then matched to the GPS location (sampled at 1 Hz) closest to the midpoint of the AMS collection cycle. We used the Squirrel 1.57I toolkit for analyzing AMS data at unit mass resolution (written by Donna Sueper at University of Colorado for Igor PRO software, v6.27, Wavemetrics, Inc.). We applied a collection efficiency (CE) of 0.5 to the measured organic mass, which is standard practice with ambient AMS data sets for OA quantification.32 CE is an estimate of particle losses within the AMS, which are largely driven by the mechanism of particle bounce upon impact with the vaporizer prior to electron impact ionization.32 Different OA types can have substantially different CE values,33,34 and the results of Ye et al.35 show clear external mixtures of OA in near-source environments in Pittsburgh. However, we were unable to measure the CE for individual components during this study and so apply a uniform CE of 0.5 in the absence of better information. On the basis of the results of Docherty et al.,36 who show CE as a function of oxidation state for lab-generated OA, CE = 0.5 seems like an appropriate choice for cookingrelated OA, though likely is too low for fresh vehicle emissions, possibly resulting in overestimating some OA mass attributable to traffic sources. Similarly, we use the standard AMS relative ionization efficiency (RIE) of 1.4 for all organic mass. Recent work from Reyes-Villegas et al.37 has suggested that this RIE is appropriate for fresh cooking-related aerosol emissions, but may overestimate OA mass for aged cooking-related aerosols. Given the focus on fresh cooking-related plumes in this work, we deemed RIE = 1.4 to be appropriate for this analysis. To quantify the contribution of local sources, we corrected the data to account for temporal changes in background concentrations within and between drives. First, we identified areas with no known sources within 200 m; these were either on our driving route (e.g., large city parks, residential streets with low vehicular activity, etc.) or at the CMU campus urban background site. OA data collected from these locations were averaged together, defining the background concentration for that time. For all drives, data collected at the CMU campus site before and after the drive were used to define the background. For most drives, we also collected background data during the middle of the drives at the aforementioned areas (parks, quiet streets). A linear interpolation in time was then applied between background averaging periods for each drive. This was subtracted from the total OA concentration to provide the above-background OA concentration. This above-background OA concentration should be indicative of local sources and processes affecting OA. We refer to this above-background OA concentration as “ΔOA” hereafter. We refer to “plumes” throughout the rest of the paper as those data where OA concentrations were greater than 2 μg m−3 above background (ΔOA > 2 μg m−3). The mass spectral information provided by the AMS allows us to attribute plumes to vehicle sources, cooking sources, and other sources. OA from vehicle sources, often referred to as HOA (hydrocarbon-like OA), is relatively unoxidized, with an oxygen-to-carbon (O:C) of ∼0.03.38 HOA characteristically has high signal at m/z 57 (largely C4H9+) and other unoxidized alkyl (CxHy+) fragments. OA from cooking activities is typically more oxidized (4−7×) than HOA38,39 (O:C ≈ 0.11−0.2), and has a high relative signal at m/z 55 (significant contributions from both C4H7+ and C3H3O+). We use the ratio of ions at m/ z 55 and 57 (denoted “55:57” hereafter) as one metric to D

DOI: 10.1021/acs.est.8b02654 Environ. Sci. Technol. XXXX, XXX, XXX−XXX

Article

Environmental Science & Technology

Figure 2. (a) Time series from a typical drive, where we sampled in three commercial districts containing restaurants and in afternoon rush hour on a highway. The plot shows OA mass concentrations (green line) in the three different districts (left axis) and similarity to HOA (circle markers, right axis) colored by the ratio of 55:57. Upper panels show two mass spectra of different plume event types: (b) mass spectra collected in the highway tunnel and (c) near the kitchen exhaust of a burger restaurant (c). Panel d shows a difference spectra illustrating the mass spectral differences between (c) and (b).

urban park (∼2 km2) provides a sampling area devoid of local OA sources, and a location where we collected AMS mass spectra for background averaging. The outbound highway (AADT of 20 685 vehicles/day) is one of the major commuter highways in the Pittsburgh area and includes a two-bore 1.3 km-long traffic tunnel within the stretch that we sampled. Each of the three neighborhoods contains one or multiple clusters of restaurants (ranging from 7 to 30 restaurants), which are located on streets with AADT between 4300 and 10 000 vehicles per day. The commercial districts of Neighborhoods A and B are situated within residential areas containing both apartment buildings and single-family houses, while the landuse of Neighborhood C is a mix of neighborhood commercial and neighborhood industrial. Our driving within these neighborhoods encircled all sides of their respective commercial district(s), in an effort to minimize any biases from dispersion conditions during a given drive (e.g., wind direction). In some cases during this drive, which are labeled in Figure 2a, we were able to pinpoint plumes to individual restaurants. The time series in Figure 2a illustrates multiple important findings from this study. First, concentrations of OA are highly variable in space and this variability is driven significantly by restaurants. Second, OA concentrations in commercial districts with restaurants are frequently much higher than on highways. For example, during the drive shown in Figure 2a, the highest OA concentration we measured during the 20 min of rushhour driving on the highway, on-ramps, and through a tunnel was 6 μg m−3, while there were many instances of OA concentrations exceeding 10 μg m−3 in the commercial business districts. Third, as discussed below, the AMS mass spectra had chemical signatures of cooking emissions. For

the biomass-burning plumes than the other cooking plumes in our data set. It seems reasonable that these plumes contain a mix of organic molecules from the incomplete combustion of wood as well as aerosolized constituents from meats or other oils. Thus, there is some potential for classifying restaurant plumes in the “other” category, if wood is used as a cooking fuel for some restaurants in this sampling domain. Similarly, there appears to be a continuous spectrum of plumes between our cooking and vehicle categories, and we determined the empirical boundary based on the reference spectra. Thus, there is some potential for misclassifying vehicle plumes as cooking and vice versa, but this was not an issue for the 39 reference spectra with which we tested our classification scheme. Our classification separates HOA from cooking plumes at 55:57 = 1.6, where there is a clear boundary between ambient HOA spectra and laboratory cooking spectra from the reference data in Table S2. Below, we present the sensitivity of our plume classification to the specific choice of 55:57 as the boundary between HOA and cooking plumes.



RESULTS Figure 2a presents measurements collected on a typical drive, both in terms of sampling design and the important data types collected. In this time series, we show OA concentration and similarity to HOA for plumes (ΔOA > 2 μg m−3). Points are colored by 55:57. For this sample drive, our route passed through the following areas after departing CMU’s campus: a large urban park, the congested outbound lanes of an interstate highway during rush-hour, and three neighborhoods with restaurantcontaining commercial districts and roads spanning a wide range of traffic volumes, before returning to campus. The large E

DOI: 10.1021/acs.est.8b02654 Environ. Sci. Technol. XXXX, XXX, XXX−XXX

Article

Environmental Science & Technology

Figure 3. Composition analysis of all OA plumes greater than 2 μg m−3 above the background (ΔOA > 2 μg m−3). (a) Similarity to HOA values for all plume mass spectra as a function of ΔOA mass concentration. Each point is colored by the ratio of 55:57, which serves as an indicator of whether a plume is likely from traffic or cooking sources. Letters correspond to plumes identified with certainty as being emitted from vehicles or restaurant sources. Red and purple colored bands indicate the range of Similarity values used to classify spectra. (b) Histogram showing the frequency of plumes as a function of similarity, which shares a left-axis with panel (a). (c) Histogram summarizing the frequency of measured plumes as a function of ΔOA for each plume category.

between 0.7 and 0.9 illustrated by the histogram in panel b of the figure. The summary of the plume categorization is shown in Figure 3c as a function of ΔOA concentration. The majority of OA plumes in Pittsburgh neighborhoods during the sampling time frame have cooking-like mass spectra (71%). We do see a significant fraction of plumes with HOA-like mass spectra (20%), but this is more than a factor of 3 less than the frequency of cooking-like spectra. Nine percent of plumes were attributed to other, unidentified sources. The fraction of plumes we attribute to cooking increases with the measured OA concentration, as shown in Figure 3c. Therefore, plumes with the highest concentrations were most frequently due to cooking. During evening mealtime (5 pm to 9 pm local time), cooking appears to be the dominant source of high concentration OA plumes in the Pittsburgh area. As described in the Methods section and Figure S1, trafficand cooking-related plumes seem to form a continuum, with apportionment between these sources relying in large part on the choice of 55:57 to use as the boundary between plume classes. We chose 55:57 = 1.6 to define the boundary between traffic- and cooking-related mass spectra based on laboratory and ambient spectra from the literature. To assess how sensitive our results were to this cutoff, we adjusted the boundary to 55:57 = 2.05, which corresponds to the lowest 55:57 value of direct source-sampling from this data set for cooking plumes. Using 55:57 = 2.05 results in 53% of plumes attributed to cooking, 26% to vehicles/traffic, and 21% to “other.” Figure S1 shows this sensitivity test in similarity-55:57 space. Even this boundary (55:57 = 2.05), which does not seem reasonable based on the available reference spectra, does not change our conclusion that we see many more cooking-like OA plumes compared to HOA-like plumes.

example, almost all of the plumes in Figure 2a have similarity values less than 0.8 and high 55:57 values, which are characteristic of cooking OA. Almost 90% of the plumes in Figure 2 are attributed to cooking emissions versus 5% to traffic. Examples of plume mass spectra (MS) collected within the highway tunnel (Figure 2b) and next to the kitchen exhaust of a hamburger restaurant (Figure 2c) are shown in the upper panels of Figure 2. Both MS are background-subtracted, as discussed in the Methods section. Figure 2d shows a difference MS between background-subtracted tunnel and restaurant spectra to further illustrate differences between these two plumes. There are distinct mass spectral signatures for trafficand cooking-related OA. While both MS have aliphatic fragments, the series of fragments enhanced in the traffic tunnel spectra (m/z 57, 69, 71, 83, 85) is distinct from the series of fragments that are more abundant in the restaurant emissions (m/z 41, 55, 67), as shown in this difference plot (Figure 2d). These spectral differences form the basis of our plume categorization, and are also visualized from the colored similarity points in Figure 2a. Plume Categorization. The patterns observed in Figure 2a, that cooking plumes are more numerous and larger in magnitude than traffic-related plumes, apply to the rest of our mobile sampling data set as well. The composition of all OA plumes is shown in Figure 3. Figure 3a shows similarity to HOA as a function of ΔOA concentration, colored by 55:57 values. Letters in Figure 3a are from unambiguous (known a priori) restaurant and traffic sources sampled in this data set and are labeled above the plot. We sort all plume data by their background-subtracted mass spectra into three categories, as described in the Methods section: HOA-like, cooking-like, and other. The majority (93%) of the plumes have similarity values F

DOI: 10.1021/acs.est.8b02654 Environ. Sci. Technol. XXXX, XXX, XXX−XXX

Article

Environmental Science & Technology

Figure 4. ΔOA concentration as a function of distance from centroid of restaurant cluster. The clusters are named and displayed as follows: (a) Restaurant cluster 1, (b) Restaurant cluster 2, (c and d) Restaurant cluster 3 on two different dates. Each example shows the raw data (gray markers), binned-averages of the raw data (colored circles), and an exponential decay fitted to the binned averages. The right axis shows the enhancement above the background (all data are normalized by the upwind average concentration). The distance at which the fit fell to half of the binned-average max (D50) is displayed for each example. Vertical bars indicate the minimum and maximum concentration measured within each bin. Horizontal bars indicate the width of the bin. Two visits were made to Restaurant cluster 3 (panel (c) and (d)), which had the highest and lowest D50 values. This likely reflects variability in source strength and dispersion conditions.

Spatial Scale of Elevated Concentrations around Restaurants. Figure 4 shows OA concentrations as a function of upwind and downwind distance from three restaurant clusters. In each example, OA is elevated downwind of restaurants, but not upwind, and we measure OA concentrations significantly above the background level hundreds of meters downwind of the restaurant clusters. The mass spectra of each downwind plume was categorized as cooking-like. Measured ΔOA concentrations from these plumes are highly variable across these three restaurant clusters, but still greatly elevated above upwind concentrations. For example, ΔOA concentrations in the first 50 m downwind of Restaurant Cluster 3 were as high as 200 μg m−3 (mean 54 μg m−3) on 2016/12/14, but only as high as 40 μg m−3 (mean 7 μg m−3) on 2016/11/16. Restaurant Clusters 1 and 2 had maximum concentrations just above 30 μg m−3. The differences in plume magnitude and spatial extent could be caused by the volume or style of food cooking and/or dispersion conditions for each measurement set. Despite the variability in the plume concentrations, the binned-averages are significantly enhanced above the background, by tens of μg m−3 in the first 50 m downwind of the sources. The right axis of Figure 4 shows the downwind enhancement, defined by normalizing all downwind data to the upwind average, which is between 3−10× greater than the upwind concentrations for these four cases. This enhancement is much higher than that shown for pollutants downwind of roadways in Karner et al.17 The falloff curves for near-road environments vary by pollutant but are usually less than 3× the background value for particulate pollutants immediately downwind of road sources.

To quantify the spatial extent of the elevated OA concentrations downwind of the restaurant cluster, we fit an exponential decay to spatially binned average concentrations in Figure 4. Measurements were grouped into fifty-meter bins closest to the restaurant cluster, but larger bin widths were used when data was sparse such that there were always three data points per bin at minimum. The distance where the fit fell to half of the max value (D50) varied from 52−368 m for these three examples. There was both day-to-day and locationspecific variations in D50. The lowest and highest values of the range of D50 were from the same restaurant cluster (#3), but on different days. The falloff curves in Figure 4 demonstrate that restaurants substantially impact neighborhood-scale OA concentrations. The length scale of restaurant plume impacts in Figure 4 (52− 368m) is similar to the fall off for traffic-related pollutants downwind of roadways. Concentrations of CO, ultrafine particles, and particulate black carbon all have near-roadway D50 between 50 and 150 m.17 Massoli et al. observed a similar D50, between 100 and 200 m, for total PM emissions from traffic (sum of HOA and BC) downwind of a highway.42 Figure 4 shows that there is variability in the OA distribution downwind of restaurant clusters. This scatter is likely due to multiple factors. Changing wind direction and turbulence caused by complex urban topography affect pollutant dispersion and prevents a precise definition of the plume centerline. The distance from the restaurant cluster is calculated from the center of the cluster, which may not be the location of the exhaust for those restaurant(s) emitting the OA. Lastly, restaurant cooking is done in batches on the order of minutes to match customer demand. Thus, some scatter in G

DOI: 10.1021/acs.est.8b02654 Environ. Sci. Technol. XXXX, XXX, XXX−XXX

Article

Environmental Science & Technology

and restaurants having the highest ΔOA concentrations. These results, combined with the source attribution using chemical composition, present a strong case that OA concentrations in areas with both roads and busy restaurants are likely to be meaningfully higher than concentrations in urban background locations, and that restaurant emissions are a major contributor to elevated OA concentrations in these areas. The results also show that restaurants can be important contributors to local OA hotspots during mealtimes.

the data is inherent to restaurants as sources of OA. The source strength can be highly dynamic and collecting the data to define up- and down-wind concentrations took 30 min to an hour for each case. However, even with this variability, it is clear restaurants contribute significantly to outdoor OA concentrations at neighborhood scales. Distribution of ΔOA by Land-Use Category. Figures 2 and 4 show the presence and spatial extent of plumes associated with individual restaurants or clusters of restaurants. Here, we compare OA concentrations across the four land-use classes defined above: Background, Roads-only, Restaurantsonly, and Both restaurants and busy roads. These concentration differences are summarized in the normalized cumulative distributions of grid-cell average ΔOA concentrations for each land-use type shown in Figure 5. There are clear variations in ΔOA across the land use classes, with the lowest median concentrations in Background grid cells and highest in the Both category.



IMPLICATIONS Restaurant emissions can create elevated neighborhood scale submicron organic PM concentrations. In Allegheny County, 260 000 people (21% of total population) live within 200 m of a restaurant (63 mi.2 total buffer area). This is comparable to the 432 000 people (35% of population) living within 200 m of a busy (>7058 AADT) road (164 mi.2 total buffer area). However, organic aerosol concentrations downwind of restaurant clusters often greatly exceed those near roadways. Busy roads are more spatially distributed than restaurants in Pittsburgh, but restaurants are often located in more-densely populated areas relative to busy roads. Despite busy roads being the larger areal source in Pittsburgh, we measured a much larger number of OA plumes from food cooking (71%) compared to vehicle sources (20%). A significant portion of people living near restaurants also live near busy roads (14.8% of people in Allegheny County live within 200 m of both). This represents the potential for compound exposure to primary pollutants from these two source classes. Restaurants may be important for driving acute ambient exposures to OA. Acute exposure to high PM2.5 concentrations can have strong associations with negative health outcomes43 as well as total mortality.44 Peak exposures to PM2.5 may have greater associations with some negative health outcomes than average PM2.5 exposures.45 The numerous, but localized, hotspots of elevated OA concentrations around many restaurants create acute exposures to high OA concentrations, up to ∼100 μg m−3, during times of heavy cooking. During evening meal-times, the largest plumes of OA (and thus PM) are attributable to cooking activities. The data presented here may not be representative of longterm, spatially resolved OA concentrations as the sampling was conducted only during winter and with mostly evening driving. However, we consistently observed OA plumes in the vicinity of many restaurants and in areas with high restaurant density. We confirmed that the origin of many of these plumes to be cooking through the chemical information provided by the AMS, which can identify the broad source categories of primary OA. Assuming that other urban areas have a similar source mix of vehicles and restaurants, our data suggest that restaurant emissions are likely a major contributor to large spatial gradients and high exposures to PM in many cities.

Figure 5. Normalized cumulative distributions of average ΔOA values for grid cells of each land-use category: Background, Roads-only, Restaurants-only, and Both.

For each of the four land-use types, the majority of grid cells have average ΔOA concentrations less than 2 μg m−3. Almost all (96%) of the grid-cells for the Background land-use type have ΔOA less than 2 μg m−3. Busy-roads only, Restaurantsonly, and Both land-use types have 86%, 79%, and 68% of their grid cells with ΔOA less than 2 μg m−3, respectively. While there are differences between the land-use types, the majority of cells within each land-use type do not meet the “plume” criteria (ΔOA > 2 μg m−3). Nonetheless, any grid-cell average ΔOA concentrations that are above zero still represent local variations in OA concentrations. For example, the median ΔOA grid cell concentration for the Both land-use type is 1.2 μg m−3, 39% higher than the mean OA concentration for all background-averaging periods in the data set. There were grid cells with significantly elevated average ΔOA concentrations, and the distribution of these highconcentration grid cells also appears to be correlated with differences in land use. For the top decile of ΔOA, for example, the difference between the Both land-use type and the other land-use types becomes more substantial, though the ordering remains the same. Both is the highest (ΔOA = 5.1 μg m−3), followed by the Restaurants-only (3.2 μg m−3), Busy roads only (2.5 μg m−3), and Background (1.3 μg m−3) land-use types. Land-use clearly has an effect on OA concentration patterns in Pittsburgh, with areas containing both busy roads



ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.est.8b02654. Table with details of individual drives, table of known sources shown in Figure 3, and plot further detailing the similarity/55:57 space and plume categorization method (PDF) H

DOI: 10.1021/acs.est.8b02654 Environ. Sci. Technol. XXXX, XXX, XXX−XXX

Article

Environmental Science & Technology



T.; Marcon, A.; Eriksen, K. T.; Raaschou-Nielsen, O.; Stephanou, E.; Patelarou, E.; Lanki, T.; Yli-Tuomi, T.; Declercq, C.; Falq, G.; Stempfelet, M.; Birk, M.; Cyrys, J.; von Klot, S.; Nádor, G.; Varró, M. J.; Dėdelė, A.; Gražulevičienė, R.; Mölter, A.; Lindley, S.; Madsen, C.; Cesaroni, G.; Ranzi, A.; Badaloni, C.; Hoffmann, B.; Nonnemacher, M.; Krämer, U.; Kuhlbusch, T.; Cirach, M.; de Nazelle, A.; Nieuwenhuijsen, M.; Bellander, T.; Korek, M.; Olsson, D.; Strömgren, M.; Dons, E.; Jerrett, M.; Fischer, P.; Wang, M.; Brunekreef, B.; de Hoogh, K. Development of NO2 and NOx land use regression models for estimating air pollution exposure in 36 study areas in Europe − The ESCAPE project. Atmos. Environ. 2013, 72, 10−23. (13) Cyrys, J.; Eeftens, M.; Heinrich, J.; Ampe, C.; Armengaud, A.; Beelen, R.; Bellander, T.; Beregszaszi, T.; Birk, M.; Cesaroni, G.; Cirach, M.; de Hoogh, K.; de Nazelle, A.; de Vocht, F.; Declercq, C.; Dėdelė, A.; Dimakopoulou, K.; Eriksen, K.; Galassi, C.; Gra̧ulevičiené, R.; Grivas, G.; Gruzieva, O.; Gustafsson, A. H.; Hoffmann, B.; Iakovides, M.; Ineichen, A.; Krämer, U.; Lanki, T.; Lozano, P.; Madsen, C.; Meliefste, K.; Modig, L.; Mölter, A.; Mosler, G.; Nieuwenhuijsen, M.; Nonnemacher, M.; Oldenwening, M.; Peters, A.; Pontet, S.; Probst-Hensch, N.; Quass, U.; Raaschou-Nielsen, O.; Ranzi, A.; Sugiri, D.; Stephanou, E. G.; Taimisto, P.; Tsai, M.-Y.; Vaskövi, É .; Villani, S.; Wang, M.; Brunekreef, B.; Hoek, G. Variation of NO2 and NOx concentrations between and within 36 European study areas: Results from the ESCAPE study. Atmos. Environ. 2012, 62, 374−390. (14) Eeftens, M.; Tsai, M.-Y.; Ampe, C.; Anwander, B.; Beelen, R.; Bellander, T.; Cesaroni, G.; Cirach, M.; Cyrys, J.; de Hoogh, K.; de Nazelle, A.; de Vocht, F.; Declercq, C.; Dėdelė, A.; Eriksen, K.; Galassi, C.; Gražulevičienė, R.; Grivas, G.; Heinrich, J.; Hoffmann, B.; Iakovides, M.; Ineichen, A.; Katsouyanni, K.; Korek, M.; Krämer, U.; Kuhlbusch, T.; Lanki, T.; Madsen, C.; Meliefste, K.; Mölter, A.; Mosler, G.; Nieuwenhuijsen, M.; Oldenwening, M.; Pennanen, A.; Probst-Hensch, N.; Quass, U.; Raaschou-Nielsen, O.; Ranzi, A.; Stephanou, E.; Sugiri, D.; Udvardy, O.; Vaskövi, É .; Weinmayr, G.; Brunekreef, B.; Hoek, G. Spatial variation of PM2.5, PM10, PM2.5 absorbance and PMcoarse concentrations between and within 20 European study areas and the relationship with NO2 - Results of the ESCAPE project. Atmos. Environ. 2012, 62, 303−317. (15) Zimmerman, N.; Presto, A. A.; Kumar, S. P. N.; Gu, J.; Hauryliuk, A.; Robinson, E. S.; Robinson, A. L.; Subramanian, R. A machine learning calibration model using random forests to improve sensor performance for lower-cost air quality monitoring. Atmos. Meas. Tech. 2018, 11, 291−313. (16) Mohr, C.; DeCarlo, P. F.; Heringa, M. F.; Chirico, R.; Richter, R.; Crippa, M.; Querol, X.; Baltensperger, U.; Prévôt, A. S. H. Spatial Variation of Aerosol Chemical Composition and Organic Components Identified by Positive Matrix Factorization in the Barcelona Region. Environ. Sci. Technol. 2015, 49, 10421−10430. (17) Karner, A. A.; Eisinger, D. S.; Niemeier, D. A. Near-Roadway Air Quality: Synthesizing the Findings from Real-World Data. Environ. Sci. Technol. 2010, 44, 5334−5344. (18) Hayes, P. L.; Ortega, A. M.; Cubison, M. J.; Froyd, K. D.; Zhao, Y.; Cliff, S. S.; Hu, W. W.; Toohey, D. W.; Flynn, J. H.; Lefer, B. L.; Grossberg, N.; Alvarez, S.; Rappenglück, B.; Taylor, J. W.; Allan, J. D.; Holloway, J. S.; Gilman, J. B.; Kuster, W. C.; de Gouw, J. A.; Massoli, P.; Zhang, X.; Liu, J.; Weber, R. J.; Corrigan, A. L.; Russell, L. M.; Isaacman, G.; Worton, D. R.; Kreisberg, N. M.; Goldstein, A. H.; Thalman, R.; Waxman, E. M.; Volkamer, R.; Lin, Y. H.; Surratt, J. D.; Kleindienst, T. E.; Offenberg, J. H.; Dusanter, S.; Griffith, S.; Stevens, P. S.; Brioude, J.; Angevine, W. M.; Jimenez, J. L. Organic aerosol composition and sources in Pasadena, California, during the 2010 CalNex campaign. Journal of Geophysical Research: Atmospheres 2013, 118, 9233−9257. (19) Sun, Y. L.; Zhang, Q.; Schwab, J. J.; Demerjian, K. L.; Chen, W. N.; Bae, M. S.; Hung, H. M.; Hogrefe, O.; Frank, B.; Rattigan, O. V.; Lin, Y. C. Characterization of the sources and processes of organic and inorganic aerosols in New York city with a high-resolution time-

AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected] (A.A.P.). ORCID

Ellis Shipley Robinson: 0000-0003-1695-6392 Rishabh Urvesh Shah: 0000-0002-4608-1972 Joshua Schulz Apte: 0000-0002-2796-3478 Allen L. Robinson: 0000-0002-1819-083X Albert A. Presto: 0000-0002-9156-1094 Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS This article was developed under Assistance Agreement No. RD83587301 awarded by the U.S. Environmental Protection Agency. It has not been formally reviewed by EPA. The views expressed in this document are solely those of authors and do not necessarily reflect those of the Agency. EPA does not endorse any products or commercial services mentioned in this publication. Additional funding comes from NSF grant number AGS1543786.



REFERENCES

(1) Pope, C. A., III; Dockery, D. W. Health Effects of Fine Particulate Air Pollution: Lines that Connect. J. Air Waste Manage. Assoc. 2006, 56, 709−742. (2) Lepeule, J.; Laden, F.; Dockery, D.; Schwartz, J. Chronic Exposure to Fine Particles and Mortality: An Extended Follow-up of the Harvard Six Cities Study from 1974 to 2009. Environ. Health Perspect. 2012, 120, 965−970. (3) Apte, J. S.; Messier, K. P.; Gani, S.; Brauer, M.; Kirchstetter, T. W.; Lunden, M. M.; Marshall, J. D.; Portier, C. J.; Vermeulen, R. C. H.; Hamburg, S. P. High-Resolution Air Pollution Mapping with Google Street View Cars: Exploiting Big Data. Environ. Sci. Technol. 2017, 51, 6999−7008. (4) US EPA. National Ambient Air Quality Standards for Particulate Matter: Final Rule; Federal Register, 1997. (5) Jerrett, M.; Arain, A.; Kanaroglou, P.; Beckerman, B.; Potoglou, D.; Sahsuvaroglu, T.; Morrison, J.; Giovis, C. A review and evaluation of intraurban air pollution exposure models. J. Exposure Sci. Environ. Epidemiol. 2005, 15, 185−204. (6) Zhu, Y. F.; Hinds, W. C.; Kim, S.; Shen, S.; Sioutas, C. Study of ultrafine particles near a major highway with heavy-duty diesel traffic. Atmos. Environ. 2002, 36, 4323−4335. (7) Holmes, N. S.; Morawska, L. A review of dispersion modelling and its application to the dispersion of particles: An overview of different dispersion models available. Atmos. Environ. 2006, 40, 5902− 5928. (8) Fallah-Shorshani, M.; Shekarrizfard, M.; Hatzopoulou, M. Integrating a street-canyon model with a regional Gaussian dispersion model for improved characterisation of near-road air pollution. Atmos. Environ. 2017, 153, 21−31. (9) Li, H. Z.; Dallmann, T. R.; Li, X.; Gu, P.; Presto, A. A. Urban Organic Aerosol Exposure: Spatial Variations in Composition and Source Impacts. Environ. Sci. Technol. 2018, 52, 415−426. (10) Hatzopoulou, M.; Valois, M. F.; Levy, I.; Mihele, C.; Lu, G.; Bagg, S.; Minet, L.; Brook, J. Robustness of Land-Use Regression Models Developed from Mobile Air Pollutant Measurements. Environ. Sci. Technol. 2017, 51, 3938−3947. (11) Hankey, S.; Marshall, J. D. Land Use Regression Models of OnRoad Particulate Air Pollution (Particle Number, Black Carbon, PM 2.5, Particle Size) Using Mobile Monitoring. Environ. Sci. Technol. 2015, 49, 9194−9202. (12) Beelen, R.; Hoek, G.; Vienneau, D.; Eeftens, M.; Dimakopoulou, K.; Pedeli, X.; Tsai, M.-Y.; Künzli, N.; Schikowski, I

DOI: 10.1021/acs.est.8b02654 Environ. Sci. Technol. XXXX, XXX, XXX−XXX

Article

Environmental Science & Technology of-flight aerosol mass apectrometer. Atmos. Chem. Phys. 2011, 11, 1581−1602. (20) Crippa, M.; El-Haddad, I.; Slowik, J. G.; DeCarlo, P. F.; Mohr, C.; Heringa, M. F.; Chirico, R.; Marchand, N.; Sciare, J.; Baltensperger, U.; Prévôt, A. S. H. Identification of marine and continental aerosol sources in Paris using high resolution aerosol mass spectrometry. Journal of Geophysical Research: Atmospheres 2013, 118, 1950−1963. (21) Ge, X.; Setyan, A.; Sun, Y.; Zhang, Q. Primary and secondary organic aerosols in Fresno, California during wintertime: Results from high resolution aerosol mass spectrometry. Journal of Geophysical Research: Atmospheres 2012, 117, D19301. (22) Mohr, C.; DeCarlo, P. F.; Heringa, M. F.; Chirico, R.; Slowik, J. G.; Richter, R.; Reche, C.; Alastuey, A.; Querol, X.; Seco, R.; Peñ uelas, J.; Jimenez, J. L.; Crippa, M.; Zimmermann, R.; Baltensperger, U.; Prévôt, A. S. H. Identification and quantification of organic aerosol from cooking and other sources in Barcelona using aerosol mass spectrometer data. Atmos. Chem. Phys. 2012, 12, 1649− 1665. (23) Elser, M.; Bozzetti, C.; El-Haddad, I.; Maasikmets, M.; Teinemaa, E.; Richter, R.; Wolf, R.; Slowik, J. G.; Baltensperger, U.; Prevot, A. S. H. Urban increments of gaseous and aerosol pollutants and their sources using mobile aerosol mass spectrometry measurements. Atmos. Chem. Phys. 2016, 16, 7117−7134. (24) Vert, C.; Meliefste, K.; Hoek, G. Outdoor ultrafine particle concentrations in front of fast food restaurants. J. Exposure Sci. Environ. Epidemiol. 2016, 26, 35−41. (25) Abernethy, R. C.; Allen, R. W.; McKendry, I. G.; Brauer, M. A Land Use Regression Model for Ultrafine Particles in Vancouver, Canada. Environ. Sci. Technol. 2013, 47, 5217−5225. (26) Ryan, P. H.; LeMasters, G. K. A review of land-use regression models for characterizing intraurban air pollution exposure. Inhalation Toxicol. 2007, 19, 127−133. (27) Li, H. Z.; Dallmann, T. R.; Gu, P.; Presto, A. A. Application of mobile sampling to investigate spatial variation in fine particle composition. Atmos. Environ. 2016, 142, 71−82. (28) Tan, Y.; Lipsky, E. M.; Saleh, R.; Robinson, A. L.; Presto, A. A. Characterizing the spatial variation of air pollutants and the contributions of high emitting vehicles in Pittsburgh, PA. Environ. Sci. Technol. 2014, 48, 14186−14194. (29) Tan, Y.; Dallmann, T. R.; Robinson, A. L.; Environment, A. P. A. Application of Plume Analysis to Build Land Use Regression Models from Mobile Sampling to Improve Model Transferability; Elsevier, 2016. (30) Jayne, J. T.; Leard, D. C.; Zhang, X. F.; Davidovits, P.; Smith, K. A.; Kolb, C. E.; Worsnop, D. R. Development of an aerosol mass spectrometer for size and composition analysis of submicron particles. Aerosol Sci. Technol. 2000, 33, 49−70. (31) DeCarlo, P. F.; Kimmel, J. R.; Trimborn, A.; Northway, M. J.; Jayne, J. T.; Aiken, A. C.; Gonin, M.; Fuhrer, K.; Horvath, T.; Docherty, K. S.; Worsnop, D. R.; Jimenez, J. L. Field-Deployable, High-Resolution, Time-of-Flight Aerosol Mass Spectrometer. Anal. Chem. 2006, 78, 8281−8289. (32) Huffman, J. A.; Ziemann, P. J.; Jayne, J. T.; Worsnop, D. R.; Jimenez, J. L. Development and Characterization of a Fast-Stepping/ Scanning Thermodenuder for Chemically-Resolved Aerosol Volatility Measurements. Aerosol Sci. Technol. 2008, 42, 395−407. (33) Middlebrook, A. M.; Bahreini, R.; Jimenez, J. L.; Canagaratna, M. R. Evaluation of Composition-Dependent Collection Efficiencies for the Aerodyne Aerosol Mass Spectrometer using Field Data. Aerosol Sci. Technol. 2012, 46, 258−271. (34) Robinson, E. S.; Onasch, T. B.; Worsnop, D.; Donahue, N. M. Collection efficiency of a-pinene secondary organic aerosol particles explored via light-scattering single-particle aerosol mass spectrometry. Atmos. Meas. Tech. 2017, 10, 1139−1154. (35) Ye, Q.; Gu, P.; Li, H. Z.; Robinson, E. S.; Lipsky, E.; Kaltsonoudis, C.; Lee, A. K. Y.; Apte, J. S.; Robinson, A. L.; Sullivan, R. C.; Presto, A. A.; Donahue, N. M. Spatial Variability of Sources and Mixing State of Atmospheric Particles in a Metropolitan Area. Environ. Sci. Technol. 2018, 52, 6807−6815.

(36) Docherty, K. S.; Jaoui, M.; Corse, E.; Jimenez, J. L.; Offenberg, J. H.; Lewandowski, M.; Kleindienst, T. E. Collection Efficiency of the Aerosol Mass Spectrometer for Chamber-Generated Secondary Organic Aerosols. Aerosol Sci. Technol. 2013, 47, 294−309. (37) Reyes-Villegas, E.; Bannan, T.; Le Breton, M.; Mehra, A.; Priestley, M.; Percival, C.; Coe, H.; Allan, J. D. Online Chemical Characterization of Food-Cooking Organic Aerosols: Implications for Source Apportionment. Environ. Sci. Technol. 2018, 52, 5308−5318. (38) Mohr, C.; Huffman, J. A.; Cubison, M. J.; Aiken, A. C.; Docherty, K. S.; Kimmel, J. R.; Ulbrich, I. M.; Hannigan, M.; Jimenez, J. L. Characterization of Primary Organic Aerosol Emissions from Meat Cooking, Trash Burning, and Motor Vehicles with HighResolution Aerosol Mass Spectrometry and Comparison with Ambient and Chamber Observations. Environ. Sci. Technol. 2009, 43, 2443−2449. (39) Ng, N. L.; Canagaratna, M. R.; Jimenez, J. L.; Chhabra, P. S.; Seinfeld, J. H.; Worsnop, D. R. Changes in organic aerosol composition with aging inferred from aerosol mass spectra. Atmos. Chem. Phys. 2011, 11, 6465−6474. (40) Zhang, Q.; Jimenez, J. L.; Canagaratna, M. R.; Ulbrich, I. M.; Ng, N. L.; Worsnop, D. R.; Sun, Y. Understanding atmospheric organic aerosols via factor analysis of aerosol mass spectrometry: a review. Anal. Bioanal. Chem. 2011, 401, 3045−3067. (41) Liu, S.; Russell, L. M.; Sueper, D. T.; Onasch, T. B. Organic particle types by single-particle measurements using a time-of-flight aerosol mass spectrometer coupled with a light scattering module. Atmos. Meas. Tech. 2013, 6, 187−197. (42) Massoli, P.; Fortner, E. C.; Canagaratna, M. R. Pollution gradients and chemical characterization of particulate matter from vehicular traffic near major roadways: Results from the 2009 Queens College Air Quality. Aerosol Sci. Technol. 2012, 46, 1201−1218. (43) Zanobetti, A.; Dominici, F.; Wang, Y.; Schwartz, J. D. A national case-crossover analysis of the short-term effect of PM2.5 on hospitalizations and mortality in subjects with diabetes and neurological disorders. Environ. Health 2014, 13, 898−12. (44) Son, J.-Y.; Bell, M. L. The relationships between short-term exposure to particulate matter and mortality in Korea: impact of particulate matter exposure metrics for sub-daily exposures. Environ. Res. Lett. 2013, 8, 014015−9. (45) Lin, H.; Liu, T.; Xiao, J.; Zeng, W.; Guo, L.; Li, X.; Xu, Y.; Zhang, Y.; Chang, J. J.; Vaughn, M. G.; Qian, Z. M.; Ma, W. Hourly peak PM2.5 concentration associated with increased cardiovascular mortality in Guangzhou, China. J. Exposure Sci. Environ. Epidemiol. 2017, 27, 333−338.

J

DOI: 10.1021/acs.est.8b02654 Environ. Sci. Technol. XXXX, XXX, XXX−XXX