1. Introduction
The emission of particles derived from anthropogenic activities into the atmosphere affects the hydrological cycle, air quality, and global radiative balance (Ramanathan and Carmichael, 2008). Since atmospheric particles absorb and scatter light, they lead to changes in energy and temperature in the environment (Horvath, 1993); hence, the terrestrial radiation budget and global warming are intimately connected with atmospheric aerosols. Furthermore, in polluted urban atmospheres, particulate matter with an aerodynamic diameter of 2.5 µm or less (PM2.5) contains black carbon (BC), a critical short-life climate forcer. That is why some countries have included BC within their mitigation agendas.
The incomplete combustion of coal, trash, biomass, and fossil fuels produces gases and particles containing BC and other light-absorbing substances (Paredes-Miranda et al., 2009; Aiken et al., 2010). BC absorbs light efficiently at all visible wavelengths, but other carbonaceous substances (i.e., brown carbon, BrC) prefer UV and short wavelengths (Li et al., 2020). Although particles with BC are slowly reacting in the atmosphere, they interact with many other species leading to changes in optical properties (Peng et al., 2016). For instance, some authors have found fractal soot particles with non-absorbing material that modifies the particle’s physical, hygroscopic, and optical properties (Weingartner et al., 2003).
BC and elemental carbon (EC) are usually interchangeable terms, despite different analytical methods supporting the definitions of each species. Instruments measure BC based on the absorption coefficient of particles, and EC is measured with thermochemical analyses (Pöschl, 2005). The link between light absorption (b abs ) and BC/EC abundance is the mass absorption cross-section or mass absorption efficiency (MAE) (Cheng et al., 2011; Tao et al., 2018; Cappa et al., 2019); however, it depends on particle size, aging, mixing state, coating distribution of species, and other particle factors (Knox et al., 2009), so it is hardly a constant value due to the diversity and complexity of atmospheric particles in nature. We have defined the MAE as the optical interaction of a chemical compound in the particle (in this case, elemental carbon) with solar radiation.
Some air quality networks determine BC concentrations with monitors that measure the optical properties of the particle because it is faster and cheaper than chemical techniques and there is sometimes a bias, since instruments assume that BC is the only absorber and the particle has a constant MAE. Aerosols from Mexico City contain around 50% by mass of carbonaceous compounds (Torres-Jardón et al., 2009; Peralta et al., 2019) and this quantity may alter the optical properties of the particle.
The atmosphere of Mexico City is under a volatile organic compound regime. It is a highly reactive system due to ozone and other oxidants in the air that probably accelerate aging processes and modify the optical properties of the particles (García-Reynoso et al., 2009; Peralta et al., 2021). We measured the light absorption (babs) and scattering (bsca) coefficients of PM2.5 at two wavelengths, sampled PM2.5, quantified EC and organic carbon (OC) content, and described aerosol MAE changes as a function of the CO/CE ratio. The different values of MAE will reduce the uncertainty of BC measurements in cities of developing countries.
2. Method and instruments
The measurement campaign took place at the Instituto de Ciencias de la Atmósfera y Cambio Climático (Institute for Atmospheric Sciences and Climatic Change, ICAyCC) building on the Universidad Nacional Autónoma de México (National Autonomous University of Mexico, UNAM) campus, from January 19 to March 20, 2015, as a part of the Atmospheric Aerosols Campaign at University Campus (ACU15) within the Air Quality Project Study in Central Mexico (ECAIM) (Salcedo et al., 2018).
2.1 Meteorology
The Red Universitaria de Observatorios Atmosféricos (University Network of Atmospheric Observatories, RUOA) and the Programa de Estaciones Metereológicas del Bachillerato Universitario (University Program of High School Meteorological Stations, PEMBU) provided the meteorological data (wind speed and direction, temperature, solar radiation, pressure, and relative humidity) for the campaign. PEMBU has Vantage Pro 2 weather monitoring stations (Davis Instruments). Data can be consulted at https://www.ruoa.unam.mx/pembu/index.php?page=historical_facts
2.2 Particle optical properties
The instruments were located on the roof of the ICAyCC building, and the sampling probe was 5 m above the roof and 10 m above the ground to avoid resuspended soil particles. We used a photoacoustic spectrometer (PAS) and a photoacoustic extinctiometer (PAX) to measure the optical properties of PM2.5. Both instruments operated continuously from January 21 to February 6 and from February 17 to March 21, 2015. We connected both instruments to a TSI aerosol splitter and then to a 2.5 µm cyclone (URG). PAS and PAX sampled aerosols with a flow rate of 1.0 l min-1.
PAS and PAX have the same operation principle, where the quantity of absorption measured is proportional to the sound pressure produced in an acoustic resonator caused by the absorption of light. The acoustic resonator operates in plane wave mode with a resonant frequency of 500 Hz and a photoacoustic coefficient of 12.8 Pa (W m-1). The coherent acoustic noise is suppressed with filters at the pressure nodes of the resonator. The laser wavelengths selected for the photoacoustic spectrometers avoid absorption of gaseous atmospheric light, with a gaseous absorption coefficient of 0.1 Mm-1 (Liñán-Abanto et al., 2021). The instruments have also a wide-angle integrating reciprocal nephelometer to measure the scattering coefficient of light.
The lower detection limit for light absorption is 0.4 Mm-1, which corresponds to a mass density of elemental carbon of 40 ng m-3, assuming a light absorption efficiency of 10 m2 g-1. Also, the absorption and scattering coefficients of atmospheric gases like O2, NO2, O3, and H2O are two orders of magnitude lower compared to those of aerosols (Horvath, 1993).
The University of Nevada built the PAS, which has a laser operating at 532 nm (Arnott et al., 1999). The PAX is a commercial model (Droplet Measurement Technologies, USA) with a laser operating at 870 nm. Both spectrometers recorded the particle’s optical properties (absorption, scattering, and extinction coefficients) every minute. PAX also quantified BC using a default mass absorption efficiency (MAE) of 4.74 m2 g-1 (Prieto et al., 2023). The intercomparison of PAS and PAX data had a good determination coefficient (R2 = 0.97).
2.3 PM 2.5 sampling
Before sampling, 47 mm quartz filters (Tissue quartz 2500QAT-UO, 47 mm, PALL) were preheated at 600 ºC for 6 h to eliminate the organic material that could be present. Then, we conditioned the substrates for 24 h at constant relative humidity (less than 30%) and 25 ºC, weighed them on an analytical balance (Sartorius CPA225D-OCE), and boxed them individually in plastic Petri dishes. The plastic dishes were wrapped with parafilm tape and opened until the sampling time. Filters after sampling were stored at 4 ºC, stabilized for 24 h at the same conditions as before and weighed.
We calibrated the Minivol sampler (Airmetrics) flow rate with a Gillibrator at a constant 5.0 L min-1. The impactor collected PM2.5 on filters every 24 h. The Red Automática de Monitoreo Atmosférico (Automatic Atmospheric Monitoring Network, RAMA) of Mexico City provided continuous PM2.5 concentrations, with a Thermo Andresen model FH62C14 continuous particulate monitor based on beta attenuation.
2.4 Carbon analysis
The carbon analyzer CM5014 (UIC) measures OC and total carbon (TC) content. The instrument quantifies CO2 produced from the complete combustion of carbonaceous material in a coulombimetric cell and converts it into carbon content. We calibrated the instrument with NIST SRM 1649a Urban Dust (Álvarez-Ospina et al., 2016). OC evolved at 450 ºC and TC at 700 ºC. The difference between TC and OC corresponds to elemental carbon (EC), following Eq. (1):
2.5 Optical properties and statistical analysis
We measured the particle’s absorption and scattering coefficients every minute and averaged the values for 24 h to match the optical property databases with those of OC and EC analyses. With the daily averages, we calculated the MAE532 (532 nm) and MAE870 (870 nm) according to Eq. (2), which represents the ratio of particle absorption coefficient and the concentration of light-absorbing chemical species, basically EC. We also calculated the single scattering albedo (SSA) with Eq. (3):
We applied a linear regression of b abs as a function of EC, where the slope corresponds to the MAE (m2 g-1). To contrast the hypothesis of normality, we used the Kolmogorov-Smirnov statistical test (α = 0.05). The data had a non-parametric distribution; hence, we applied the Mann-Whitney statistical test to compare b abs and EC.
3. Results
3.1 Meteorology
January, February, and March are part of the cold-dry season in Mexico City, and cold air masses often blow from the S-SW to the city. Southerly winds are associated with air masses containing low temperatures and high relative humidity. These air masses are not polluted, unlike those coming from the north, where the industrial area of Mexico City is located.
From January to mid-March the mean wind speed was 1.8 m s-1. Wind speed was slow, so we believe local activities produced most of the airborne particles. By March 15, the average wind speed increased to 3.0 m s-1 (Fig. 1). In January and February, the W-NW wind direction prevailed and in the week of March 8 to 15, the wind changed its direction and blew from the south (Fig. 1). From January 19 to February 15, the average temperature was 13.9 ± 1.5 ºC and from March 8 it increased to 17.8 ± 1.6 ºC (Fig. 1). From January 19 to March 8 relative humidity was 50% and after March 8 it increased to 60% with southerly winds. Figure 1 shows the basic meteorological parameters (wind speed and direction, temperature, and relative humidity) recorded in the measurement campaign.
3.2 Particle optical properties
At 532 nm, the average babs was 13.68 ± 6.89 Mm-1, and at 870 nm it was 10.40 ± 5.33 Mm-1. From January 18 to 31, it decreased as the wind intensity increased. In February, babs varied from one week to another, ranging from 7.25 to 21.92 Mm-1 (532 nm) and from 5.10 to 17.44 Mm-1 (870 nm). From March 8 to 15, it decreased as the wind speed went up (Fig. 2).
The average bsca was 71.00 ± 40.56 Mm-1 at 532 nm and 21.37 ± 12.31 Mm-1 at 870 nm. Between January 21 and February 6, the scattering decreased to half of its average value, with the wind blowing from N-NW, probably transporting pollutants from the industrial part of the city. From March 8-15, bsca decreased associated with a cold front entering Mexico City (Fig. 2).
At both wavelengths, the SSA reached its minimum value at 07:00, when light absorption was at its maximum produced by the presence of dark particles. The SSA532 maximum was 0.89, and the SSA870 maximum was 0.77, around 13:00 LT, when the absorption coefficient dropped to its minimum value (Prieto et al., 2023). The average SSA532 was 0.83 ± 0with.35 Mm-1 and the average SSA870 was 0.66 ± 0.28 Mm-1. The difference in values is probably due to more chemical species in the particle being sensitive to lower wavelengths.
Particles are covered with different materials, oxidated species, water, and other coating substances changing the particle optical properties and promoting the formation of secondary aerosols. Furthermore, SSA532 values were close to 1.0 (Fig. 2), probably linked to secondary aerosols and other materials that scatter light (Eck et al., 1999; Cheng et al., 2011; Schumann, 2012; Tao et al., 2019).
3.3 Carbon content
The average mass concentration of PM2.5 was 18.45 ± 8.00 µg m-3, while the average concentration of OC was 7.96 ± 3.61 µg m-3, and for EC it was 2.44 ± 1.77 µg m-3. So, EC corresponded to 12% and OC to 39% of the total PM2.5 mass, which means that the carbonaceous material is 51% of the particle mass. Vehicles and commercial activities were probably the main sources of emission of carbonaceous material near the sampling site. The carbon composition had significant variations, perhaps due to the contribution of various sources and secondary atmospheric processes. The correlation coefficient between OC and EC was 0.10. The average OC/EC ratio ranged from 1.08 to 22.89, with a mean of 5.58 ± 1.08 (Turpin et al., 1990; Chow et al., 2011).
The was 1.08, and primary organic carbon (POC) was calculated following Eq. (4). We used Eq. (5) to calculate secondary organic carbon (SOC) (Castro et al., 1999). SOC average concentration was 5.32 ± 2.96 µg m-3, meaning that the secondary species contribution was more significant than POC. SOC is linked to temperature, but the statistical analysis concludes that there is no significant difference between SOC concentrations and monthly temperature (Zhang et al., 2013). In January, February, and March, the average temperature was 14.54, 14.94, and 15.88 ºC, respectively.
Figure 3 shows the daily variations of organic carbon (POC and SOC) and EC. Variations in POC and SOC concentrations are likely due to changes in emission sources, transport of airborne material, or secondary processes in the atmosphere.
3.4 Mass absorption efficiency (MAE)
MAE efficiency corresponds to the ratio of daily average absorption coefficient and EC concentration, following Eq. (2). The range of MAE532 was 0.11-2.44 m2 g-1 with an average value of 0.59 ± 0.57 m2 g-1. The range of MAE870 was 1.39-25.60 m2 g-1, with an average value of 6.78 ± 5.85 m2 g-1. The PAX uses a default MAE of 4.74 m2 g-1 that falls into the average MAE870.
The EC concentration and the daily average babs changed from one day to another, probably caused by changes in the emission rates of local sources, transport, or atmospheric aging processes. Hence, MAE at both wavelengths had considerable variations, and the correlation coefficients of EC versus babs were 0.41 (at 532 nm) and 0.39 (at 870 nm).
We sorted the data by OC/EC ratios in four groups assuming an external mixing of the particle. We assumed EC was the only absorbing material in the particle at both wavelengths and different concentrations of scattering material modified the absorbing properties of EC, like a shadowing effect (Weingartner et al., 2003). Then we calculated the linear regression for each group of data. Groups 1 and 5 corresponded to OC/EC ratios < 2.0, where OC and EC concentrations were similar. Groups 2 and 6 were for OC/EC ratios between 2.0 and 3.1. Groups 3 and 7 comprised OC/EC ratios between 3.1 and 8.0. Groups 4 and 8 were for particles with a notable predominance of organic compounds, with OC/EC ratios greater than 8.0. Table I shows the OC/EC ratio ranges, the population, the linear regression equation, and the determination coefficient for each group. Figure 4 shows the data and linear adjustment for each group. All regressions intersect the babs axis at EC = 0 because we assumed EC was the unique contribution to light-absorbing in the particle and OC modified the light-absorbing properties of EC but did not absorb light directly.
Group | OC/EC ratio range | n | y = mx | MAE (m2 g-1) | R2 | SOC/POC range |
Wavelength: 532 nm | ||||||
1 | OC/EC < 2.0 | 6 | y = 19.75x | 19.75 | 0.88 | 0.00 - 0.64 |
2 | 2.0 < OC/EC < 3.1 | 14 | y = 7.75x | 7.75 | 0.94 | 1.02 - 1.83 |
3 | 3.1 < OC/EC < 8. 0 | 19 | y = 4.46x | 4.46 | 0.97 | 2.06 - 5.67 |
4 | OC/EC > 8.0 | 8 | y = 2.67x | 2.67 | 0.97 | 7.88 - 20.16 |
Wavelength: 870 nm | ||||||
5 | OC/EC < 2.0 | 6 | y = 15.26x | 15.26 | 0.89 | 0.00 - 0.64 |
6 | 2.0 < OC/EC < 3.1 | 14 | y = 5.95x | 5.95 | 0.94 | 1.02 - 1.83 |
7 | 3.1 < OC/EC < 8.0 | 19 | y = 3.37x | 3.37 | 0.97 | 2.06 - 5.67 |
8 | OC/EC > 8.0 | 8 | y = 2.03x | 2.03 | 0.97 | 7.88 - 20.16 |
OC: organic carbon; EC: elemental carbon; MAE: mass absorption efficiency; SOC: secondary organic carbon; POC: primary organic carbon.
The carbonaceous material, with a high contribution of OC and other species, changed the optical properties of aerosols, and the amount of primary and secondary organic compounds modified the aerosols’ absorbing light capacity. When the amount of OC was like EC in the particle (OC/EC ratio < 2.0), we had a MAE532 of 19.75 m2 g-1, and MAE870 was 15.26 m2 g-1. As the amount of OC increased MAE diminished, and if the OC concentration was eight times the EC mass (OC/EC ratio > 8.0), then the MAE532 was equal to 2.67 m2 g-1, and the MAE870 was 2.03 m2 g-1.
4. Discussion
Activities at the UNAM campus, as well as nearby commercial areas and vehicles, impact the measuring site by emitting particles. The SSA suggests the presence of particulate matter with a significant coating of scattering substances since values are 0.82 at 532 and 0.66 at 880 nm at noon. SSA is linked to the radiative energy budget, and those values show the variety of absorbing and scattering material in urban aerosols.
It is also possible to find aerosols with a large and thick coat of scattering material like water and organic or inorganic compounds. Authors mention that the aging processes of aerosols modify the BC light absorption capacity, as observed in Mexico City (Salcedo et al., 2018) and in China (Wang et al., 2018, 2020). We found a large amount of OC in particles, probably resulting from natural and human atmospheric sources.
Linear regressions of the different groups show high determination coefficients (R2 > 0.85), meaning that organic compounds are important in the particle’s interaction with light. Table II shows MAE from other studies for comparison with the results shown in Table I. MAE values shown in Table II do not have an increasing sequence as the OC/EC ratio decreases. The samples were collected in different environments, and the particle optical properties were measured at different wavelengths.
OC/EC reported by other authors | MAE reported by other authors | Wavelength | Source | Reference |
8.17 | 8.11 (EC) | 880 nm | Biomass burning | (Tao et al., 2020) |
0.24 (OC) | ||||
8.50 | 11.3 (EC) | 660 nm | Urban | (Ram and Sarin, 2009) |
2.94 | 1.00 (WSOC) | 700 nm | Urban | (Li et al., 2019) |
17.62 | 6.10 (EC) | 632 nm | Biomass Burning | (Hu et al., 2017) |
1.19 | 8.10 (EC) | Vehicle emissions | ||
2.8 | 0.3 (WSOC) | 365 nm | Urban | (Srinivas and Sarin, 2013) |
1.5 | 0.3 (WSOC) | |||
0.7 | 0.2 (WSOC) |
OC: organic carbon; EC: elemental carbon; MAE: mass absorption efficiency; SOC: secondary organic carbon; POC: primary organic carbon; WSOC: water-soluble organic carbon. Abbreviations in parentheses indicate that compounds were used to calculate MAE.
The OC average concentration was three times higher than EC, and there were days when the OC fraction reached eight times the EC. OC is primary and secondary, and the OC/EC ratios change MAE532 and MAE870. The absorbing material of the particle probably has both BC and BrC; the first species is associated with vehicular activities and is sensitive to detection at 870 nm. The 532 nm wavelength is commonly used to detect BrC from biomass sources (i.e., natural emissions and biomass burning, among others). However, in our case, we adjusted both wavelengths to detect only EC, which is a chemical species that we can analyze with thermochemical techniques.
Since the OC/EC ratio modifies the particle’s optical properties, it affects BC measurements underestimating or overestimating its concentration during the day. Based on a constant MAE870 = 4.74 m2 g-1 (like PAX), the underestimation reaches 75% if the OC/EC ratio is below 2.0, and the overestimation is close to 100% if the OC/EC ratio is above 8.0. The determination of BC based on the particle optical properties and OC/EC ratios is probably influenced by emissions from different sources of combustion (gasoline, diesel, biomass, among others), meteorological changes (i.e., relative humidity, temperature), and other atmospheric reactants.
If the emission rates of pollutants and the photo-oxidation processes in the atmosphere change on time modifying the particle optical properties, it is not easy to estimate a general MAE without chemical analyses of both OC and EC. Some authors find that O3 is a proxy of SOC presence (Turpin and Huntzicker, 1995). Figure 5 (upper panel) shows BC readings by PAX, BC corrected by MAE870, and EC, and the same figure (lower panel) shows the OC/EC ratio and the daily median of O3. There is a tendency to have high OC/EC ratios with O3 concentrations above 30 ppb, but this does not always happen.
5. Conclusions
Based on the SOC/POC ratio and the the measuring site is a receptor point with aging aerosols prevailing over newly emitted particles, where OC is the predominant species, capable of modifying the light absorption properties of absorbing species (assuming EC is the only absorbing species). The particle undergoes daily transformations, modifying its optical properties. The MAEs obtained by other authors show a variety of values since they depend on the measuring wavelength, the type of carbonaceous compounds, and the predominant particle source. However, other reactive species, like O3, can operate as a proxy for the presence of SOC.