INTRODUCTION
Around the world, PM2.5 (fine particles with aerodynamic diameters less than 2.5 µm) are significant pollutants that affect air quality in urban areas (WHO 2016). Moreover, high concentrations of PM2.5 in ambient air are associated with visibility, economic, and human health impacts (Pope et al. 2011, Grahame et al. 2014, Yang et al., 2019). Evidence suggests that considerable exposure to PM2.5 increases the risk of suffering from cardiovascular and respiratory diseases (Anake et al. 2020). Recently, the black carbon (BC) content of PM2.5 was significantly associated with cardiovascular disease mortality (Kirrane et al. 2019). However, very few comparable data are available for other locations around the world (Hopke et al. 2020).
There is continuous interest in accurately representing the chemical composition of PM2.5 because such knowledge will enable possible emission sources to be identified. This is important because different regions have varied chemical compositions; in addition, such chemical compositions are associated with PM2.5 toxicities (Anake et al. 2020). Chemical markers that are commonly used to trace source contributions include trace elements, soluble inorganic compounds, elemental carbon (EC), and compounds found in the organic carbon (OC) fraction of aerosols (Li et al. 2010). The specific composition of emission from different sources, such as vehicle exhaust, cooking operations, industries, biomass burning, vegetative detritus, and fugitive dusts should be determined to integrate a chemical source profile and accurately estimate the relative contribution of major emission sources to ambient PM2.5 concentrations (Pernigotti et al. 2016, Singh et al. 2017, Liu et al. 2018a, Miller et al. 2019, Hopke et al. 2020). These profiles can be used for different applications: 1) performing receptor models based on chemical profiles, such as chemical mass balance (CMB) models; 2) establishing a better speciated emissions inventory that could serve as input data to dispersion models; 3) interpreting the output of receptor models based exclusively on speciated ambient samples (such as positive matrix factorization [PMF]). CMB is a powerful and widely used receptor model that estimates the relative contribution of an emission source to ambient PM2.5 concentrations (USEPA 2004, Karagulian et al. 2015, Hopke et al. 2020). However, CMB’s results are often limited because the source profiles used might be site-specific; given that emission sources are quite different between sites, they cannot be easily generalized. Thus, the development and availability of source profiles relevant for different study areas are highly recommended (Xie et al. 2012, Belis et al. 2015, Karagulian et al. 2015). In spite of this limitation, several studies around the world used profiles from other sites to conduct source apportionment studies in order to identify tendencies in source contributions (Mohr et al. 2009, Karagulian et al. 2015, Mancilla et al. 2016, Kalaiarasan et al. 2018, Rehman et al. 2020).
Therefore, in the international context, the data of chemical source profiles must be increased. In Mexico, public resources are allocated to determine the causes of air quality deterioration using reference stations to monitor criteria pollutants in 83 cities (INECC 2019a). However, few studies have been conducted to chemically characterize fine particles in Mexican cities such as San Luis Potosí (INECC 2018a), Mexico (INECC 2019b) and Juárez (INECC 2018b). The Monterrey Metropolitan Area (MMA), which is the third most populated urban center in Mexico, hast to deal with air quality problems due to high concentrations of ambient PM2.5 (INECC 2019a). The chemical composition of ambient PM2.5 in this urban area has been widely studied (Badillo-Castañeda et al. 2015, Blanco-Jiménez et al. 2015, Martínez-Cinco et al. 2016, Mancilla et al. 2018). Moreover, source apportionment studies using data available for other locations around the world were conducted (Mancilla et al. 2016, Martínez-Cinco et al. 2016). These studies show that vehicle exhaust, cooking operations, and industries are major contributors to ambient PM2.5. The registered vehicle fleet of the MMA consists of 1.7 million passenger cars, 10 000 buses, 468 000 heavy trucks and 84 000 motorcycles (INEGI 2020). In relation to cooking operations, the average consumption of meat in the MMA is 29 kg/person, which is 71 % more than the national average. The culinary industry of the region is composed of more than 16 000 restaurants (INEGI 2015). Finally, the industrial activity in the MMA involves 135 000 economical units and is dominated by manufacturing (23 %), stores (17 %), construction (9 %), transport (9 %), and others (19 %) (INEGI 2017).
This work provides, for the first time, data of chemical source profiles based on the composition of trace elements, EC and OC for diesel-fueled vehicle exhausts, meat cooking operations, oil refinery and manufacturing industries, open-waste burning, and construction dust in the MMA. Moreover, for the second time, this study presents a chemical source profile for gasoline-fueled vehicle exhausts. These profiles provide necessary region-specific source profile data which are typically not measured in developing countries. Because these profiles were obtained based on a near-source measurement approach, correction by possible paved-road dust resuspension interference was applied to reflect a better composition of emissions and reduce errors in resulting source apportionment studies (Hopke et al. 2020). This is of particular relevance, given that the MMA is located in a semi-arid region where dry conditions foster the presence of dust in the ambient air. In addition, statistical tests were performed to establish a degree of similarity of the obtained profiles, thus ensuring source representativeness. The obtained source profiles can be used to conduct source apportionment studies in the MMA or other locations with similar emission sources.
MATERIALS AND METHODS
Sources and sampling locations
In the current investigation, the potential major source types that contribute to PM2.5 concentrations in the MMA were identified from previous studies (Mancilla et al. 2016, Martínez-Cinco et al. 2016) and official emission inventories (SEDESU 2016). Source samples were collected from 12 locations in the MMA (Mexico) from June 19, 2014 to July 22, 2014. The collected samples represent (1) residential and commercial meat cooking operations from natural gas and charcoal; (2) motor vehicle exhausts from gasoline- and diesel-fueled vehicles, including passenger cars and urban and freight transport; (3) open-waste burning; (4) oil refinery and manufacturing industries, and (5) urban background profiles including suburban area and urban construction sites.
All the samples were collected following the ground-based source-dominated method for source profile characterization (Chow et al. 2004, Zhao et al. 2019). Two filter-based instruments were used at the same time to collect the samples from each source emission. A low-volume sampler (MiniVol, Airmetrics) was used to collect samples on 47 mm Teflon filters with pore size of 0.2 µm (Whatman PTFE membrane circles, TE36) at an average flow rate of 0.34 m3/h (~5.7 Ll/min). The oother equipment was a high-volume sampler (TE-6070-2.5, Tisch Environmental) used to collect samples on 0.20 × 0.25 m pre-fired micro-quartz filters with pore size 2.2 µm (Whatman QM-A) at an average flow rate of 69 m3/h. The two samplers were calibrated before starting each sampling test by using certificated orifice calibrators (NIST traceable standard) MNF-1236 and TE-5040 for the MiniVol and Hi-Vol samplers, respectively. These calibrations were conducted to estimate a more precise flow rate sampling for each sampler. Quartz filters were collected by using Hi-Vol samplers instead of a MiniVol sampler due to further chemical analysis that will require higher amounts (> 400 µg) of organic carbon as a complement study of this work. The sampling locations, descriptions and times are shown in table I and illustrated in figure 1; preliminary sampling tests were conducted to define the most adequate sampling times for each chemical source profile. The same number of samples and sampling time per filter type were collected for each emission source profile. The average sampling volume can be calculated using the flow rates for each sampler and the average sampling time; for the Teflon filters these volumes range from 0.3-5.7 m3 while for Quartz filters they range from 78-1162 m3 with a correlation (r) between volumes of 0.88.
Source | Profile (samples used) | Code | Number of samples per filter type | Average sampling time (h) |
Meat cooking operations | Restaurant charbroiled and grilled meat: late afternoon samples. (25º39’18.9”N, 100º17’29.2”W) | CAR | 2 | 3.99 |
Residential charbroiled meat: late afternoon samples including start and end of burn and cooking process. (25º39’00.5”N, 100º17’05.5”W) | CAF | 2 | 2.90 | |
Supermarket charbroiled meat service. Afternoon samples including start and end of burn and cooking process. (25º38’21.7”N, 100º17’00.6”W) | CAS | 2 | 2.28 | |
Vehicle exhausts | Gasoline-fueled vehicles: heavy traffic in a bore tunnel with a +3.5% slope (Loma Larga Tunnel). (25º39’30.25”N, 100º20’11.85”W) | TLL | 3 | 3.60 |
Urban Transport: ground-based samples from intersection of Av. Juárez and Arteaga Street. (25º40’59.2”N, 100º18’48.6”W) | TP | 3 | 3.43 | |
Freight transport: ground-based samples from intersection of Av. Félix U. Gómez and Av. Ciudad Los Ángeles. (25º42’28.7”N, 100º16’55.9”W) | CC | 3 | 2.62 | |
Industry | Oil refinery: daytime samples in a hot spot near to the refinery. (25º36’10”N, 99º59’19.7”W) | PMX | 2 | 12.94 |
Manufacturing industry: daytime samples in a hot spot near to an industrial park. (25º39’30.25” N, 100º20’11.85” W) | PSN | 4 | 8.00 | |
Biomass and trash burning | Open-waste burning: field burn samples including start and end burn, active burn and heavy smoke. (25º41’24.1”N, 100º10’00.1”W) | QB | 2 | 0.89 |
Urban background | Suburban area in daytime: daytime samples during spring-summer time. (25º29’34.1”N, 100º10’52.8”W) | RTD | 3 | 6.80 |
Suburban area in nighttime: nighttime samples during spring-summer time. (25º29’34.1”N, 100º10’52.8”W) | RTN | 2 | 15.24 | |
Urban construction site: line 3 of the city subway construction site on the Av. Félix U. Gómez. (25º41’24.1”N, 100º17’47.5”W) | RP | 2 | 2.21 |
Sampling sites description
Motor vehicle exhaust samples were collected for gasoline- and diesel-fueled vehicles from a road tunnel and through side-road sampling, respectively. In all cases, the samplers were placed on a sidewalk within 2 m of the traffic lane, with the sampling inlet placed at 1.5 m above ground level. For the case of the road tunnel, general guidelines from a previous study in the same tunnel were followed (Araizaga et al. 2013). Three samples were collected in a two-bore urban tunnel (532 m long) in which the traffic is composed of 98 % gasoline-fueled vehicles, including motorcycles and passenger cars (Araizaga et al. 2013). The samplers were located 200 m downstream from the inlet of the bore; the traffic on the selected bore is subject to a 3.5 % positive slope. The vehicle exhausts for diesel trucks were divided into urban (public and private buses) and freight transport (at least four-axle single units and trailers). Information about the streets in the MMA with high urban and freight transport traffic was provided by the Secretaría de Seguridad Pública y Viabilidad de Monterrey (Monterrey Department of Public Safety and Roads). Based on this information, two locations were selected for sampling these source emissions. Samples were collected from these public road-side locations during the hours in which the fleet of gasoline-powered vehicles was low, which was determined by previously conducted vehicle count tests in these avenues. During the tests, the average vehicular capacities through each avenue were 400 buses and 200 freight trucks per hour for the urban and freight transport avenues, respectively.
Meat cooking emissions were collected based on the most frequent forms of cooking operations in the MMA. These samples included emissions from restaurants that use natural gas for combustion purposes and seed oils for heating (hot plate) and cooking food, supermarket charbroiled meat services, and residential charbroiled meat activities. In all cases, sampling was conducted under real conditions, which implied uncontrolled, non-conducted, open operations. In particular, the residential cooking test was comprised by a specific combination of meat types and other food complements (e.g., sausage, onions, cheese) using an outdoor charbroiling system. All tests involved the cooking of meat from beef, pork, sausage, and other typical food dishes associated with traditional charbroiled cooking operations in the MMA. The sampling equipment was positioned at 5 m downwind from each source; in the particular case of restaurant emissions, this implied positioning the sampling equipment in open, public spaces in the surroundings of the cooking operations.
The biomass and trash burning profile was obtained by ground-based source-dominated sampling in plumes of real open-waste burning emissions. The waste burning consisted of wood, grass, leaves, plastics, cardboard, paper, textiles, and household items. Two samples were collected at 2 m downwind from the open fire, where the plume directly impacted the sampler. This procedure was conducted in the suburbs of the MMA, where these burnings often occur. The site was located from an opportunity chase process: drive arounds and fire occurrence identification in open spaces.
The urban background profile was obtained by sampling in an open space (~60 000 m2), property of the Tecnológico de Monterrey, located in the suburban area of the MMA at ~30 km southeast from the downtown area. This type-of-villa space is far from the urban area, thus allowing the isolation of biogenic emissions. This site is a natural ecosystem for recreational purposes (uninhabited) and consists of shrubs, grass, and trees. Several activities, such as cooking operations and swimming, are frequently carried out in this site. Sampling was conducted during periods in which these activities were not occurring in order to minimize their influence on the profile. The photochemical oxidation of particles produced over natural areas leads to secondary organic aerosol (SOA) formation (Xu et al. 2020). Three samples were collected for daytime and two for nighttime to distinguish SOA formation. The samplers were placed on a small building’s rooftop, with the sampling inlet at 4.5 m above the ground.
Industrial profiles were obtained for oil refinery and manufacturing industry emissions. Sampling locations were set based on the locations of reference stations deployed to monitor these ambient hot spots. For the oil refinery, two samples were collected from a sampling site located at 3 km downwind from the source, positioned in a public open space. This oil refinery is composed of more than 30 facilities that process around 110 thousand barrels of petroleum per day to produce mainly gasoline and diesel fuels. For the industrial emissions, the selected site was an industrial park that houses several manufacturing companies in the field of automotive assembly, electrical appliances, cosmetics, and engineering products. Unlike the case of the oil refinery, four samples were collected from a sampling site located within 1 km downwind from these industrial sources, also positioned in a public open space.
The urban construction profile was obtained by positioning the samplers nearby an urban construction area. This area is influenced by particles coming from pave roads and the construction operations of a new subway line for local transportation. Significant resuspension of particles occurred due to construction activities that required the use of materials and machinery, such as concrete and off-road vehicles, respectively.
Chemical characterization
The Teflon filters were analyzed for mass by gravimetry using a Sartorius ME5 microbalance (1 mg readability) in a weighing chamber with controlled temperature and relative humidity to obtain the total collected particulate matter. The allowable weighing temperature was a 24-h average between 20 and 23 ºC (standard deviation less than 2 ºC), and the 24-h mean relative humidity was maintained between 30 and 40 % (standard deviation less than 5 %). Then, the same filters were analyzed for 38 trace elements (Na to Pb) by X-ray fluorescence (XRF) and using a KEVES EDX-771 X-ray spectrometer following protocol six of the non-destructive analytical method (IO-3.3) approved by the US-EPA (EPA/625/R-96/010a). In this method, a filter is bombarded with X-ray energy; then, the excitement of electrons release energy that can be measured. Each element has a characteristic energy discharge which is measured to determine the quantity of each element (Method IO-3.3).
The micro-quartz filters were analyzed for OC and EC by thermal-optical transmittance (TOT) using a Sunset Laboratory thermo-optical carbon analyzer (Model 4L); samples were analyzed using the normal parameters for the NIOSH 5040 method with modified parameters that increased and lengthened the heating steps in order to keep the instrument from going over scale. A 1 cm × 1.5 m standard punch (1.5 cm2) was extracted from a loaded quartz microfiber filter and placed in an oven using a quartz boat. The transmittance of a diode laser light through the filter was monitored during the analysis. The oven was initially purged with He. OC was then analyzed over varying time steps between 45 s and 300 s during the OC evolution at 310 ºC, 475 ºC, 615 ºC, and 870 ºC in an inert atmosphere (pure He). The heating desorbed OC thermally. The evolved OC was catalytically oxidized to CO2 in a manganese dioxide (MnO2) oxidizing oven; subsequently, the CO2 was swept out of the oxidizing oven with the He stream and reduced to CH4 in a (Ni/firebrick) methanator and quantified as CH4 by a flame ionization detector (FID). The EC analysis was conducted using temperature profiles of 550 ºC, 625 ºC, 700 ºC, 775 ºC, and 850 ºC with a withhold time of 45 s and a final holding time of 120 s at 870 ºC in an oxidizing atmosphere (He: O2 90:10 v/v). EC was oxidized from the filter into the oxidation oven, converted into CO2, reduced to CH4, and detected by FID as CH4. During this stage, a pyrolysis correction was made. A split point should be defined for OC and EC to quantify these components. The split point is defined as the point at which the light transmittance of the sample returns to the initial value. The carbon that evolved before or after the split point was considered OC or EC, respectively (Birch and Cary 1996). All laboratory analyses were conducted by Chester LabNet (Tigard, OR).
Correction of chemical source profile by dust resuspension
The source particle mass concentrations of trace metals and carbonaceous species can be impacted by dust resuspension given the ground-based source-dominated method utilized in this study. Hence, a correction for each profile sample based on mass percent was implemented as per equation 1.
where (m i,k )corrected is the corrected mass percent of element i in profile k, (m i,k )source is the source mass percent estimated from each sample of element i in profile k, (m i /m si )dust profile is the ratio of the mass percentages of element i to the concentration of silicon in the dust profile selected as the reference, and m si,k , is the mass percent of the silicon element in profile k. The dust profile used as the reference was the profile developed for a tunnel study in the MMA used in Mancilla et al. (2016).
In addition, the uncertainty corresponding to each (C i,k )corrected was identified. In cases where measurements were combined to form a “pooled” profile, an error propagation technique was followed to estimate the final uncertainties. The uncertainties of the individual profiles were considered independent and random in such a way that follows the sum in the quadrature technique (Taylor 1997). The estimated final uncertainties were calculated by adding the squares of each measured value and obtaining the square root.
Carbonaceous indicator (OC/EC)
In addition, a useful indicator for primary emission sources is the OC/EC ratio. Values between 1 and 3 for this ratio can be an indicator of primary emissions that come from fossil fuel combustion (Chow et al. 2010, Yan et al. 2019). However, high values of OC/EC ratios (> 3) can also be an indicator of contributions from other primary sources (Cabada et al. 2004, Chow et al. 2010, Hoyle et al. 2011). For example, some studies have reported high OC/EC ratios of 17 to 40 for wood combustion (Feng et al. 2009), approximately 3 to 7 for general biomass burning (Tao et al. 2013), and 2 to 54 for cooking operations (Wang et al. 2015a). Conversely, high values of OC/EC ratios (> 3) in ambient air samples can also be an indicator of SOA due to the photochemical oxidation of OC in the atmosphere (Hayes et al. 2013, Mancilla et al. 2015, Xu et al. 2020).
Similarity tests between chemical source profiles
A tool is needed to assess the homogeneity of chemical source profiles (Karagulian et al. 2015). The similarities between sources were tested with both the Pearson distance (PD) and the similarity identity distance (SID) in accordance with Belis et al. (2015). PD is equal to 1 - r 2 , where r 2 is the Pearson coefficient, and SID is defined by equation 2 (Weber et al. 2019):
where x j and y j are the relative masses of species j to the PM2.5 of two different source profiles and m is the number of common species between x and y. These two metrics (PD and SID) aim to compare two profiles based on their common chemical relative mass composition. PD is highly sensitive to variations in the major components of PM2.5, while SID is evenly sensitive to all components (Weber et al. 2019). PD < 0.4 and SID < 0.8 are acceptable criteria for profile similarity (Weber et al. 2019).
RESULTS AND DISCUSSION
PM2.5 composition for individual source profiles
Based on the chemical characterization of each source profile, the main chemical species and diagnostic ratios were identified and estimated, respectively, to assess the uniqueness of the profiles. The main chemical species for each profile are summarized in table II, and the average source profiles (mass percent) are shown in figures 2-5. The average source profiles were calculated by averaging the samples collected for each emission source. In addition, the Si was not reported for each profile due to correction by dust resuspension; the Si contribution before the correction ranged from 0.03 % to 9.07 % while for the dust profile used for the correction was 25.7%. Finally, the developed source profiles in this study are representative for spring-summer sampling periods.
Source Profile | Chemical abundances in mass percent | |||
< 0.1% | 0.1-1.0% | 1.0-10% | >10% | |
Residential charbroiled (CAF) | Zn, Pb, Ba, Cu, V, Mn, Co, Cr | Na, S | Cl, K, EC | OC |
Restaurant charbroiled (CAR) | Zn, Mg, Cu, V, Mn, Cr, Ni | EC | Na, S, Cl, K, Ca | OC |
Supermarket charbroiled service (CAS) | Zn, Pb, Mg, Cu, V, Co, Cr, Ni | Al, Sn | Na, S, Cl, K, EC | OC |
Gasoline-fueled vehicles (TLL) | P, Ti, V, Cr, Mn, Cu, Zn, Sb, Sr, Pb | Na, Mg, Al, Cl, Fe, Zn | S, EC | OC |
Urban transport (TP) | V, Cr, Mn, Zn, Se, As, Sr, Zr, Mo, Sn, Sb, Pb, Ti | Al, Cl, Fe, Ag, | Na, S | OC, EC |
Freight transport (CC) | V, Zn, Cr, Mn, Se, As, Sr, Zr, Mo, Hg, Hg | Na, Mg, Al, Cl, Fe, Sn | S | OC, EC |
Oil refinery (PMX) | Cr, Mn, Cu, Zr, Cd, Hg | Mg, Al, Cl, K, Ti, Fe, Sb | Na, S, EC | OC |
Manufacturing industry (PSN) | P, V, Mn, Cu, Zn, Ge, As, Se, Br, Sr, Zr, Mo, Pd, In, Sn, Pb | Na, Mg, Al, Cl, K, Cr, Fe, Pb | EC | S, OC |
Biomass and trash burning (QB) | Al, Ti, Cr, Mn, Fe, Sr, Cd, In, Ba, Hg, Pb | Na, Mg, S, Zn, Br, Sn | Cl, K | OC, EC |
Suburban area (RT) | Ti, Cr, Mn, Cu, Zn, Mo, In, Hg, Pb | Na, Mg, Al, Cl, K | S, EC | OC |
Urban construction site (RP) | P, Cr, As, Zr, Mo, Sn, Hg, Pb, Ni, Mn | Mg, Cl, K, Ti, Zn, Sb, Pb | Ca, S, EC | OC |
Figure 2 compares the average abundances of the profiles for gasoline- and diesel-fueled vehicles. The three gasoline-fueled vehicle profiles (TLL) are dominated by OC (6-31 %), EC (2-10 %), and S (3-8 %). Other appreciable elements (0.01-1.0 %) in the profiles from this study were Na, Mg, Al, Cl, Fe, Cu, Sb, Ti, and Zn. This dominance was different in Mexico City, where OC, EC, Br and Pb are the key components of vehicle emissions (INECC 2019b). Unlike in Hao et al. (2019), Hg was not detected in the gasoline profiles; in addition, the average abundances of Al, Pb, OC, and EC were 55-85, 70-95, 40-70, and 5-30 % lower than those reported in Kang et al. (2010), Zhang et al. (2015) and Hao et al. (2019), respectively.
Both diesel-fueled vehicle profiles (freight transport [CC] and urban transport profiles [TP]) were dominated by OC (CC: 17-18 %; TP: 8-24 %), EC (CC: 11-12 %; TP: 3-24 %) and S (CC: 2.2-3.6 %; TP: 2.0-3.7 %). Other noticeable trace elements (0.1-1.0 %) were Sn, Ag, Na, Cl, Fe, and Al. The average abundance of Al was 70 and 97 % lower than those reported by Liu et al. (2018 b) and Pulong et al. (2017), respectively. The average abundance of Fe was 2, 4, and 16 times higher than those reported by Hao et al. (2019), Zhang et al. (2016), and Cui et al. (2017), respectively, whereas it was 80 % lower than that reported by Pulong et al. (2017). Unlike in the urban transport profile, Ag was not detected in other studies, and the average abundance of EC was 35-40 % lower than those reported for other diesel profiles by Kang et al. (2010), Cui et al. (2017), Liu et al. (2018b), and Hao et al. (2019). In the course of this study the quality of gasoline and diesel followed Euro IV standards, however, differences in the abundances of some species can also be associated with vehicle maintenance and age (INECC 2019c).
Figure 3 compares the average abundances of the profiles for meat cooking operations. These profiles were divided into three types: residential charbroiled meat (CAF), restaurant charbroiled meat (CAR), and supermarket charbroiled meat service (CAS). Firstly, for the residential charbroiled meat, the profiles were dominated by OC (26-86 %), Cl (2-5 %), EC (0.4-4 %), K (3 %), Na (0.5 %), and S (0.3-0.4 %). Secondly, for the restaurant profiles, the most abundant species were OC (2-18 %), Cl (4-17 %), K (8-16 %), S (2-3 %), Na (1 %), and EC (0.1-0.3 %). Finally, the supermarket profiles were dominated by OC (73-99 %), Cl (1.8-4.2 %), K (1.9-4.8 %), S (1-2%), Na (1 %), and EC (1.4-4.3 %). S could be one of the impurities in charcoal and can be released to the atmosphere in the form of sulfur dioxide (SO2) (Susaya et al. 2010, Vicente et al. 2018). In addition, S is a major element emitted from cooking operations (Zhao et al. 2019); its values can range from 0.14 to 4.2 %, depending on the kind of biomass being burned (Sun et al. 2019a). Cr contributions typically come from stainless steel utensils (Zhao et al. 2019), while Mn and Pb could be emitted from cooking oil (Wang et al. 2011). Unlike in Zhao et al. (2019), the Ca in the current study was not representative (< 0.1 %). The Al, Mg, and Fe abundances (Al, Mg, and Fe < 0.1 %; Na < 1.0 %) were lower than those from Wang et al. (2015a) (Al: 0.5 %; Mg: 0.6 %; Fe: 1.0 %; Na: 1.4 %), Pulong et al. (2017) (Al: 10 %; Mg: 0.7 %; Fe:1.0 %; Na: 5.0 %), and Wang et al. (2020) (Al: 1.3 %; Mg: 11 %; Fe: 0.4 %; Na: 6.0 %). In contrast, the S, Cl, and K abundances (1-10 %) were higher than those from Pulong et al. (2017) (Cl: 0.45 %; K: 0.64 %) and Wang et al. (2020) (S: 0.5 %; Cl: 0.8 %; K: 0.4 %).
In figure 4, the average abundances of the profiles for the two categories of industries are compared. The oil refinery profiles (PMX) were dominated by OC (18-29 %), EC (1.2-3.3 %), S (3.7-5.5 %), and Na (1.0-2.5 %). Apart from the carbonaceous species, Mg, Al, Cl, K, Ti, Cr, Cu, Mn, Fe, Cd, Sb, and Pb were observable (< 1.0 %). The manufacturing industry profiles (PSN) were dominated by OC (21-22 %), S (12-14 %), and EC (1.1-2.2 %). In this case, the following elements were noticeable (< 1.0 %): Na, Mg, Al, P, Cl, K, Cr, Cu, Mn, Fe, Zn, As, Br, Mo, Pd, Ag, In, Sn, and Pb. This dominance for both was also observed for the industrial emissions in Mexico City (INECC 2019b). The S contribution in the manufacturing profiles suggests that the manufacturing industries might be using fossil fuels with S content (e.g., diesel in ancillary systems). In addition, the industries may be using petrochemical materials in their processes. For the oil refinery profiles, the average abundances of EC (2.23 %), Al (0.84 %), Cr (0.03 %), Fe (0.06 %), Pb (0.005 %), and Zn (0.002 %) were lower than those from Lou et al. (2018), which were 3.0, 0.97, 0.012, 1.18, 0.4, and 0.43 %, respectively; in contrast, the levels of OC and S were 20 and 40 % higher than those from Luo et al. (2018). Additionally, as a tracer of heavy oil burning emissions, V was not found in the oil refinery profiles, unlike in cases with similar source profiles (Taiwo et al. 2014, Zhang et al. 2016, Luo et al. 2018, Amoatey et al. 2019). These results are representative of weekdays when industrial activity is normal.
The average abundances for the open-waste burning and urban construction dust profiles are shown in figures 4 and 5, respectively. The open-waste burning profiles (QB) were dominated by OC (50-52 %) and EC (44-49 %). Apart from organic tracers, trace elements that are characteristics of organic matter (e.g., wood, grass, meals) and inorganic matter (e.g., plastics, textiles) were found in these emission profiles. Cl and K were present (1-10 %) as well as Na, Mg, and Br (0.1-1.0 %) for the organic matter (compostable materials) (Maasikmets et al. 2016), while for the inorganic matter (non-compostable materials), the noticeable trace elements (0.1-1.0 %) included Cr, As, Cd, Sn, Sb, Hg, and Pb. This is different than that observed in Mexico City, where K and Pb are the key components of biomass burning (INECC 2019b). The presence of Hg is associated with the burning or incineration processes of waste (Taiwo et al. 2014, Riffault et al. 2015) that could include discarded household items such as alkaline batteries, button cells, fluorescent lamps, thermometers, and sphygmomanometers, which may contain this element at largely varying levels (Chalkidis et al. 2020). In addition, the presence of OC, EC, Cl, Na, and K is an indicator of biomass burning processes (Tao et al. 2013, Khan et al. 2021). K/EC ratios were used to assess the biomass burning contribution (Khan et al. 2021). The average K/EC ratio in this study was 0.033; this suggests field grass and wood combustion (Tao et al. 2013, Khan et al. 2021), which agrees with the sampling site for this profile. Unlike in similar source profiles (Zhang et al. 2016, Sun et al. 2019a, Wang et al. 2020), Hg and other non-compostable tracers were found in open-waste burning profiles in the present study. Additionally, the average mass contribution of K (1.60 %) was lower than those from Sun et al. (2019a) and Wang et al. (2020), which were 7.10 and 2.39 %, respectively.
The urban construction profiles (RP) were dominated by Ca (5-7 %), OC (10-21 %), and EC (3-8 %). Other trace elements that were appreciable in these profiles were Al (2-3 %), Fe (1 %), and S (1-7 %). As can be seen in these profiles, the contributions of Ca were noticeable and can be an indicator of cement dust emissions, given that cement composition includes calcium silicates along with the presence of Fe, Al, Mg, Zn, As, and Pb (Taiwo et al. 2014, Zhang et al. 2014). A combination of Zn and Pb contributions indicates a possible paved road dust resuspension (Riffault et al. 2015). Diagnostic ratios between the trace elements that comprise fugitive dust profiles can be used as markers to track the origins of dust in urban areas (Arimoto et al. 1996, Alfaro et al. 2003, Sun et al. 2019b). In this study, the average ratios of Ca, Zn, Pb, and Fe to Al were estimated with the following values: Ca/Al = 11, Fe/Al = 0.4, Zn/Al = 0.22, and Pb/Al = 0.08. These values are different from those reported for construction dust by Sun et al. (2019b), which were 2.4-4.3, 0.7-1.2, 0.01-0.02, and 0.003-0.01, respectively. The high ratio of Ca/Al suggests a contribution from cement dust given the range values reported by Sun et al. (2019b), which were 6.3-29.5. Finally, different from other construction and fugitive dust profiles (Wang et al. 2015b, Zhang et al. 2016), the profiles in this study exhibited elemental S that can be associated with the chemical binding of sulfate with crustal elements in the road (Kim and Seinfeld 1995, Wang et al. 2015b).
The average abundances for the suburban background profiles (RT) are shown in figure 5. This profile, more than an emission profile, has the characteristic chemical abundance of air masses being transported from the eastern background region of the MMA. These profiles were dominated by OC (12-19 %), S (7-13 %), and EC (1-4 %). The contribution of elements was found to be in the order of S > Na > Mg > Al >Cl > K > Cu > Cr > Zn > Ti > Mn > Mo > Pb. Of these species, sulfur was the most abundant element. Its major chemical form in these profiles may be sulfate, which would be in line with the high OC/EC ratios (photochemically aged background air masses). Zn and Pb can be associated with emission coming from paved roads, open-waste burning, or mining industry (Riffault et al. 2015). Unlike the urban construction profile, Ca and Fe were not present in the suburban profiles because they represent outskirts areas of the MMA, while the urban construction profiles are associated to an anthropogenic activity inside the MMA.
OC/EC diagnostic ratios
The OC/EC ratios in table III are consistent within source categories but differ between categories. These consistencies and differences were statistically significant and were estimated by t-tests assuming a normal distribution and unknown variances at significance levels of 0.01, 0.05, and 0.10.
Source profile | OC/EC | Reference |
Residential charbroiled meat (CAF) | 23.74±3.67 | This studya |
Restaurant charbroiled and grilled meat (CAR) | 44.16±12.02 | |
Supermarket charbroiled meat service (CAS) | 29.88±4.80 | |
Gasoline-fueled vehicles (TLL) | 3.30±1.16 | |
Freight transport - diesel (CC) | 2.16±1.39 | |
Urban transport - diesel (TP) | 1.16±0.35 | |
Oil refinery (PMX) | 10.42±5.61 | |
Manufacturing industry (PSN) | 13.35±5.40 | |
Open-waste burning (QB) | 1.11±0.17 | |
Suburban area in daytime (RTD) | 12.85±5.0 | |
Suburban area in nighttime (RTN) | 5.23±1.15 | |
Urban construction site (RP) | 2.81±0.76 | |
Cooking smoke | 9.29 | Pulong et al. (2017) |
Gasoline-fueled vehicles | 13.11 | |
Garbage-fired power plant | 2.36 | |
Construction dust | 31.7, 5.5 | |
Barbeque (charcoal & meat) | 11.85 | Wang et al. (2020) |
Household garbage | 7.51 | |
Gasoline-fueled vehicles | 6.08 | Hao et al. (2019) |
Diesel-fueled vehicles | 1.13 | |
Biomass burning | 3, 6.5, 5.0 | Tao et al. (2013) |
Gasoline-fueled vehicles | 8.64 | Zhang et al. (2016) |
Diesel-fueled vehicles | 2.82 | |
Oil-fired power plants | 6.1 | |
Biomass burning | 13.3 | |
Construction dust | 19.8 | |
Gasoline-fueled vehicles | 2.65, 10.9 | Kang et al. (2010) |
Diesel-fueled vehicles | 0.33, 0.41, 0.67 | |
Diesel-fueled vehicles | 0.33 | Cui et al. (2017) |
Vehicles | 0.99, 2.7 | Pant et al. (2017) |
Diesel-fueled vehicles | 0.55 | Liu et al. (2018b) |
Petrochemical industry | 7.0 | Luo et al. (2018) |
Biomass burning (wood & grass) | 1.9 | Sun et al. (2019a) |
Construction dust | 3.05, 6.85 | Sun et al. (2019b) |
Barbeque (charcoal & meat) | 1.96 | Wang et al. (2015a) |
aAverage OC/EC ratios using individual values for each source sample.
The ratios for gasoline-fueled vehicles were 1.5-3.0 times higher than diesel-fueled vehicle ratios. This is related to the fact that EC levels for diesel vehicle exhausts are lower than those for gasoline vehicle exhausts because diesel engines typically have better efficiencies. The ratio for gasoline-fueled vehicles was 50-75 % lower than that reported by other studies. Meanwhile, the OC/EC ratios for diesel-fueled vehicles were both higher and lower than those reported from other studies. These differences between profiles can be associated with vehicle fleet composition, car models, fuel quality, traffic, and driving style, among others.
The average OC/EC ratios for the residential, restaurant, and supermarket cooking operations profiles were higher than those observed from table III; the high ratios in this study can be associated with the quality of the charcoal used and some other complements used to cook meat and burn charcoal. In their cooking operations, supermarkets use natural gas and charcoal, which substantially contribute OC. Additionally, the average EC fraction in the residential and supermarket cooking operations was 10 times higher than that in the restaurant cooking operations.
For the oil refinery and manufacturing industry profiles, the average OC/EC ratio was higher than that reported by other studies. Meanwhile, the open-waste burning profiles exhibited lower OC/EC ratios. For the latter profiles, the major difference lies in the presence of inorganic elements (non-compostable materials) which, along with the incomplete combustion of organic compounds, increase EC emissions. The average OC/EC ratio in the construction dust profile was much lower than that from similar profiles. Given that the construction site it could a possible influence of vehicle exhaust and fuel-powered machinery used in the construction site.
Finally, the daytime suburban OC/EC ratio (RTD profile) was three times higher than the nighttime one (RTN profile), which could be related to a secondary organic contribution due to transport emissions, local photochemical oxidation processes in the atmosphere, or presence of mixed vehicle exhaust in the background (Li et al. 2010, Mo et al. 2015, Lai et al. 2016). In addition, these ratios were 2-4 times higher than those for the urban construction profile, supporting the origin of these kind of emissions, anthropogenic for urban construction and natural-transport for suburban profiles.
Source profile similarity analysis
The main elemental composition fractions (as a percentage of the total mass) and results from the test of similarity between the source profiles in this study are shown in figures 6 and 7, respectively. Figure 6 indicates broadly that similar emission sources have a similar breakdown of main categories, or similar categories dominate the emissions. Metals dominated the elemental emissions of cooking operations, followed by non-metals. Heavy-duty diesel truck emissions (urban and freight transport) categories were similar, with a major difference in crustal presence with respect to gasoline emissions. The suburban (background) and urban construction site profiles based on the main elemental categories seemed unique and distinct from the other profiles.
Interestingly, the mass concentrations of Ti, V, Cr, Fe, and Zn for the diesel profiles were 1.5-2.0 times higher than those for the gasoline profiles. According to Alves et al. (2015), greater emissions of these elements were generated from heavy vehicles. Cu, Zn, and Sb could be emitted from brake and tire wear (Bukowiecki et al. 2009, Srimuruganandam and Nagendra 2011, Riffault et al. 2015, Pant et al. 2017). In the present study, for gasoline-fueled vehicles the Zn and Sb contributions were 6 and 1.5 times higher than those for diesel-fueled vehicles, rspectively. For the freight transport profiles (CC), one distinguishing element was Sn (0.09-0.16 %); for the urban transport profiles (TP), Ag was the distinguishing element (0.09-0.27 %). Meanwhile, for cooking operations, a remarkable result is that the S fraction in the restaurant and supermarket profiles was up to six times higher than that in the overall contribution to residential cooking operations. This suggests that the content of S in the charcoal used by the restaurants and supermarkets is relatively high. For the suburban profile, the average OC contribution during nighttime was 19 % higher than that during daytime. One reason for this could be the lower ambient temperatures during the night, which reduce the volatilization of the organic compounds emitted from vegetative detritus (Xu et al. 2020). In addition, the average EC contribution during nighttime was twice that during daytime, which could be associated with a local combustion source.
The similarities of all the chemical source profiles were investigated. The profiles were compared with each other using both the PD and SID similarity indicators. This assessment analysis helped to determine whether they were mixing with other emission sources. Figure 7 reveals that the average source profiles presented high PD and SID with values outside the acceptance box (PD < 0.4 and SID < 0.8), indicating that they were dissimilar (heterogeneous profiles) over the different source profiles of this study. This was expected, as such emissions were collected directly from the source. The variability in the PD-SID space, for those that reached the acceptance box, suggests that some profiles can be slightly mixed, or they may be different from most elements emitted from multiple types of sources, as is expected when the trace element approach is used to distinguish source profiles (Kalaiarasan et al. 2018). The most heterogeneous profiles were those for cooking operations, and the least heterogeneous were those for vehicle exhaust and suburban background. This last result indicates that, indeed, the background profile contains air masses that have received contributions from several sources, including vehicle emissions.
CONCLUSIONS
The chemical composition of PM2.5 samples collected from 12 anthropogenic sources in the MMA was obtained measuring 38 elements (Na to Pb) and two carbonaceous species (OC and EC). The source profiles of these sources were generally distinguishable, as determined by PD and SID for each pair of source profiles. Gasoline-fueled vehicles and construction dust were characterized by a particularly large contribution of crustal elements in comparison with the contribution of the other source profiles. Meanwhile, in the meat cooking operations, the transition metals were barely noticeable, unlike those in the other source profiles. The OC/EC diagnosis helped characterize primary emission sources; these ratios were found to be in the order of CAR > CAS > CAF > PSN > RTD > PMX > RTN > TLL > RP > CC > TP > QB, and they were statistically different between source emissions. Furthermore, the key components of each source profile seem to be slightly different from those in other Mexican cities in which source apportionment studies have been conducted. These profiles could further be used in receptor models, such as CMB models, to conduct source apportionment for ambient PM2.5 samples collected at different locations within the MMA or urban centers with similar emission source conditions.