Total organic carbon measurements reveal major gaps in petrochemical emissions reporting.

Anthropogenic organic carbon emissions reporting has been largely limited to subsets of chemically speciated volatile organic compounds. However, new aircraft-based measurements revealed total gas-phase organic carbon emissions that exceed oil sands industry-reported values by 1900% to over 6300%, the bulk of which was due to unaccounted-for intermediate-volatility and semivolatile organic compounds. Measured facility-wide emissions represented approximately 1% of extracted petroleum, resulting in total organic carbon emissions equivalent to that from all other sources across Canada combined. These real-world observations demonstrate total organic carbon measurements as a means of detecting unknown or underreported carbon emissions regardless of chemical features. Because reporting gaps may include hazardous, reactive, or secondary air pollutants, fully constraining the impact of anthropogenic emissions necessitates routine, comprehensive total organic carbon monitoring as an inherent check on mass closure.

Aircraft and satellite observations reveal historical gap between top-down and bottom-up CO2 emissions from Canadian oil sands.

Measurement-based estimates of greenhouse gas (GHG) emissions from complex industrial operations are challenging to obtain, but serve as an important, independent check on inventory-reported emissions. Such top-down estimates, while important for oil and gas (O&G) emissions globally, are particularly relevant for Canadian oil sands (OS) operations, which represent the largest O&G contributor to national GHG emissions. We present a multifaceted top-down approach for estimating CO2 emissions that combines aircraft-measured CO2/NOx emission ratios (ERs) with inventory and satellite-derived NOx emissions from Ozone Monitoring Instrument (OMI) and TROPOspheric Ozone Monitoring Instrument (TROPOMI) and apply it to the Athabasca Oil Sands Region (AOSR) in Alberta, Canada. Historical CO2 emissions were reconstructed for the surface mining region, and average top-down estimates were found to be >65% higher than facility-reported, bottom-up estimates from 2005 to 2020. Higher top-down vs. bottom-up emissions estimates were also consistently obtained for individual surface mining and in situ extraction facilities, which represent a growing category of energy-intensive OS operations. Although the magnitudes of the measured discrepancies vary between facilities, they combine such that the observed reporting gap for total AOSR emissions is ≥(31 ± 8) Mt for each of the last 3 years (2018-2020). This potential underestimation is large and broadly highlights the importance of continued review and refinement of bottom-up estimation methodologies and inventories. The ER method herein offers a powerful approach for upscaling measured facility-level or regional fossil fuel CO2 emissions by taking advantage of satellite remote sensing observations.

A newly developed Lagrangian chemical transport scheme: Part 1. Simulation of a boreal forest fire plume.

Forest fire research over the last several decades has improved the understanding of fire emissions and impacts. Nevertheless, the evolution of forest fire plumes remains poorly quantified and understood. Here, a Lagrangian chemical transport model, the Forward Atmospheric Stochastic Transport model coupled with the Master Chemical Mechanism (FAST-MCM), has been developed to simulate the transport and chemical transformations of plumes from a boreal forest fire over several hours since their emission. The model results for NOx (NO and NO2), O3, HONO, HNO3, pNO3 and 70 VOC species are compared with airborne in-situ measurements within plume centers and their surrounding portions during the transport. Comparisons between simulation results and measurements show that the FAST-MCM model can properly reproduce the physical and chemical evolution of forest fire plumes. The results indicate that the model can be an important tool used to aid the understanding of the downwind impacts of forest fire plumes.

Uncovering global-scale risks from commercial chemicals in air.

Commercial chemicals are used extensively across urban centres worldwide1, posing a potential exposure risk to 4.2 billion people2. Harmful chemicals are often assessed on the basis of their environmental persistence, accumulation in biological organisms and toxic properties, under international and national initiatives such as the Stockholm Convention3. However, existing regulatory frameworks rely largely upon knowledge of the properties of the parent chemicals, with minimal consideration given to the products of their transformation in the atmosphere. This is mainly due to a dearth of experimental data, as identifying transformation products in complex mixtures of airborne chemicals is an immense analytical challenge4. Here we develop a new framework-combining laboratory and field experiments, advanced techniques for screening suspect chemicals, and in silico modelling-to assess the risks of airborne chemicals, while accounting for atmospheric chemical reactions. By applying this framework to organophosphate flame retardants, as representative chemicals of emerging concern5, we find that their transformation products are globally distributed across 18 megacities, representing a previously unrecognized exposure risk for the world's urban populations. More importantly, individual transformation products can be more toxic and up to an order-of-magnitude more persistent than the parent chemicals, such that the overall risks associated with the mixture of transformation products are also higher than those of the parent flame retardants. Together our results highlight the need to consider atmospheric transformations when assessing the risks of commercial chemicals.

Evolution of Atmospheric Total Organic Carbon from Petrochemical Mixtures.

Reactive organic compounds play a central role in the formation of ozone and secondary organic aerosols. The ability to accurately predict their fate, in part, relies upon quantitative knowledge of the chemical and physical parameters associated with the total organic carbon (TOC), which includes both precursors and oxidation products that evolve in the atmosphere over short to long time scales. However, such knowledge, obtained via limited carbon closure experiments, has not been attained for complex anthropogenic emissions. Here we present the first comprehensive characterization of TOC in the atmospheric oxidation of organic vapors from light and heavy oil mixtures associated with oil sand operations. Despite the complexity of the investigated oil mixtures, we are able to achieve carbon closure (83-116%) within the uncertainties (±20%), with the degree of the closure being dependent upon the vapor composition and NOx levels. In contrast to biogenic precursors (e.g., α-pinene), the photochemical time scale required for a largely complete oxidation and evolution of chemical parameters is very long for the petrochemical vapors (i.e., ∼7-10 days vs ∼1 day), likely due to the lower initial precursor reactivity. This suggests that petrochemical emissions and their impacts are likely to extend further spatially than biogenic emissions, and retain more of their complex composition and reactivity for many days. The results of this work provide key parameters to regional models for further improving the representation of the chemical evolution of petrochemical emissions.

Fugitive Emissions of Volatile Organic Compounds from a Tailings Pond in the Oil Sands Region of Alberta.

Tailings ponds in the oil sands (OS) region in Alberta, Canada, have been associated with fugitive emissions of volatile organic compounds (VOCs) and other pollutants to the atmosphere. However, the contribution of tailings ponds to the total fugitive emissions of VOCs from OS operations remains uncertain. To address this knowledge gap, a field study was conducted in the summer of 2017 at Suncor's Pond 2/3 to estimate emissions of a suite of pollutants including 68 VOCs using a combination of micrometeorological methods and measurements from a flux tower. The results indicate that in 2017, Pond 2/3 was an emission source of 3322 ± 727 tons of VOCs including alkanes, aromatics, and oxygenated and sulfur-containing organics. While the total VOC emissions were approximately a factor of 2 higher than those reported by Suncor, the individual VOC species emissions varied by up to a factor of 12. A chemical mass balance (CMB) receptor model was used to estimate the contribution of the tailings pond to VOC pollution events in a nearby First Nations and Metis community in Fort McKay. CMB results indicate that Suncor Pond 2/3 contributed up to 57% to the total mass of VOCs measured at Fort McKay, reinforcing the importance of accurate VOC emission estimation methods for tailings ponds.

Top-Down Determination of Black Carbon Emissions from Oil Sand Facilities in Alberta, Canada Using Aircraft Measurements.

Black carbon (BC) emissions from the Canadian oil sand (OS) surface mining facilities in Alberta were investigated using aircraft measurements. BC emission rates were derived with a top-down mass balance approach and were found to be linearly related to the volume of oil sand ore mined at each facility. Two emission factors were determined from the measurements; production-based BC emission factors were in the range of 0.6-1.7 g/tonne mined OS ore, whereas fuel-based BC emission factors were between 95 and 190 mg/kg-fuel, depending upon the facility. The annual BC emission, at 707 ± 117 tonnes/year for the facilities, was determined using the production-based emission factors and annual production data. Although this annual emission is in reasonable agreement with the BC annual emissions reported in the latest version of the Canadian national BC inventory (within 16%), the relative split between off-road diesel and stack sources is significantly different between the measurements and the inventory. This measurement evidence highlights the fact that the stack sources of BC may be overestimated and the off-road diesel sources may be underestimated in the inventory and points to the need for improved BC emission data from diesel sources within facilities.

Experimental Study of OH-Initiated Heterogeneous Oxidation of Organophosphate Flame Retardants: Kinetics, Mechanism, and Toxicity.

The environmental risks and health impacts associated with particulate organophosphate flame retardants (OPFRs), which are ubiquitous in the global atmosphere, have not been adequately assessed due to the lack of data on the reaction kinetics, products, and toxicity associated with their atmospheric transformations. Here, the importance of such transformations for OPFRs are explored by investigating the reaction kinetics, degradation chemical mechanisms, and toxicological evolution of two OPFRs (2-ethylhexyl diphenyl phosphate (EHDP) and diphenyl phosphate (DPhP)) coated on (NH4)2SO4 particles upon heterogeneous OH oxidation. The derived reaction rate constants for the heterogeneous loss of EHDP and DPhP are (1.12 ± 0.22) × 10-12 and (2.33 ± 0.14) × 10-12 cm3 molecules-1 s-1, respectively. Using recently developed real-time particle chemical composition measurements, particulate products from heterogeneous photooxidation and the associated degradation mechanisms for particulate OPFRs are reported for the first time. Subsequent cytotoxicity analysis of the unreacted and oxidized OPFR particles indicated that the overall particle cytotoxicity was reduced by up to 94% with heterogeneous photooxidation, likely due to a significantly lower cytotoxicity associated with the oxidized OPFR products relative to the parent OPFRs. The present work not only provides guidance for future field sampling for the detection of transformation products of OPFRs, but also strongly supports the ongoing risk assessment of these emerging chemicals and most critically, their products.

Understanding the Impact of High-NOx Conditions on the Formation of Secondary Organic Aerosol in the Photooxidation of Oil Sand-Related Precursors.

Oil sands (OS) are an important type of heavy oil deposit, for which operations in Alberta, Canada, were recently found to be a large source of secondary organic aerosol (SOA). However, SOA formation from the OS mining, processing, and subsequent tailings, especially in the presence of NOx, remains unclear. Here, photooxidation experiments for OS-related precursors under high-NOx conditions were performed using an oxidation flow reactor, in which ∼95% of peroxy radicals (RO2) react with NO. The SOA yields under high-NOx conditions were found to be lower than yields under low-NOx conditions for all precursors, which is likely due to the higher volatilities of the products from the RO2 + NO pathway compared with RO2 + HO2. The SOA yields under high-NOx conditions show a strong dependence on pre-existing surface area (not observed in previous low-NOx experiments), again attributed to the higher product volatilities. Comparing the mass spectra of SOA formed from different precursors, we conclude that the fraction of m/z > 80 (F80) can be used as a parameter to separate different types of SOA in the region. In addition, particle-phase organic nitrate was found to be an important component (9-23%) of OS SOA formed under high-NOx conditions. These results have implications for better understanding the atmospheric processing of OS emissions.

A comprehensive risk assessment of human inhalation exposure to atmospheric halogenated flame retardants and organophosphate esters in an urban zone.

Inhalation exposure to flame retardants used as additives to minimize fire risk and plasticizers is ubiquitous in human daily activities, but has not been adequately assessed. To address this research gap, the present study conducted an assessment of human health risk for four age groups through inhalation exposure to size fractionated particle-bound and gaseous halogenated flame retardants (polybrominated diphenyl ethers (PBDEs) and alternative halogenated flame retardants (AHFRs)) and organophosphate esters (OPEs) at indoor and outdoor environments (school, office, and residence) in three districts of a megacity (Guangzhou, China). Results demonstrated that OPEs were the dominant components among all targets. Indoor daily intakes of PBDEs and OPEs were 13-16 times greater than outdoor levels for all age groups. Gaseous OPEs contributed significantly greater than particle-bound compounds to daily intakes of all target compounds. Based on the different life scenarios, hazard quotient (HQ) and incremental life cancer risk (ILCR) from adults exposure to PBDEs and OPEs in indoor and outdoor settings were the greatest, followed by adolescents, children, and seniors. The estimated HQ and ILCR for all age groups both indoors and outdoors were lower than the safe level (HQ = 1 and ILCR = 10-6), indicating that the potential health risk for local residents in Guangzhou via inhalation exposure to atmospheric halogenated flame retardants and OPEs was low.

Enhancement of aqueous sulfate formation by the coexistence of NO2/NH3 under high ionic strengths in aerosol water.

Current air quality models usually underestimate the concentration of ambient air sulfate, but the cause of this underestimation remains unclear. One reason for the underestimation is that the sulfate formation mechanism in the models is incomplete, and does not adequately consider the impact of the synergistic effects of high concentrations of multiple pollutants on sulfate formation. In this work, the roles of gaseous NO2, NH3 and solution ionic strength in the formation of sulfate in the aqueous phase were quantitatively investigated using a glass reactor and a 30 m3 smog chamber, separately. The results showed that sulfate formation was enhanced to different degrees in the presence of gas-phase NO2, NH3 and their coexistence as solutes in both liquid solution and aerosol water. NH3 enhances the aqueous oxidation of SO2 by NO2 mainly by accelerating the uptake of SO2 through increased solubility. More importantly, we found that high ionic strength in aerosol water could significantly accelerate the aqueous oxidation of SO2, resulting in unexpectedly high S(VI) formation rates. We estimate that under severe haze conditions, heterogeneous oxidation of SO2 by NO2 on aerosols may be much shorter than that through gas phase oxidation by OH, aided by high ionic strengths in aerosols. Considering the existence of complex air pollution conditions with high concentrations of NO2, NH3 and aerosol water, as expected in typical urban and suburban settings, the sulfate formation mechanisms revealed in the present work should be incorporated into air quality models to improve the prediction of sulfate concentrations.

Understanding the Impact of Relative Humidity and Coexisting Soluble Iron on the OH-Initiated Heterogeneous Oxidation of Organophosphate Flame Retardants.

The current uncertainties in the reactivity and atmospheric persistence of particle-associated chemicals present a challenge for the prediction of long-range transport and deposition of emerging chemicals such as organophosphate flame retardants, which are ubiquitous in the global environment. Here, the OH-initiated heterogeneous oxidation kinetics of organophosphate flame retardants (OPFRs) coated on inert (NH4)2SO4 and redox-active FeSO4 particles were systematically determined as a function of relative humidity (RH). The derived reaction rate constants for the heterogeneous loss of tricresyl phosphate (TCP; kTCP) and tris(2-butoxyethyl) phosphate (TBEP; kTBEP) were in the range of (2.69-3.57) × 10-12 and (3.06-5.55) × 10-12 cm3 molecules-1 s-1, respectively, depending on the RH and coexisting Fe(II) content. The kTCP (coated on (NH4)2SO4) was relatively constant over the investigated RH range while kTBEP was enhanced by up to 19% with increasing RH. For both OPFRs, the presence of Fe(II) enhanced their k by up to 53% over inert (NH4)2SO4. These enhancement effects (RH and Fe(II)) were attributed to fundamental changes in the organic phase state (higher RH lowered particle viscosity) and Fenton-type chemistry which resulted in the formation of reactive oxygen species, respectively. Such findings serve to emphasize the importance of ambient RH, the phase state of particle-bound organics in general, and the presence of coexisting metallic species for an accurate description of the degradation kinetics and aging of particulate OPFRs in models used to evaluate their atmospheric persistence.

Measured Canadian oil sands CO2 emissions are higher than estimates made using internationally recommended methods.

The oil and gas (O&G) sector represents a large source of greenhouse gas (GHG) emissions globally. However, estimates of O&G emissions rely upon bottom-up approaches, and are rarely evaluated through atmospheric measurements. Here, we use aircraft measurements over the Canadian oil sands (OS) to derive the first top-down, measurement-based determination of the their annual CO2 emissions and intensities. The results indicate that CO2 emission intensities for OS facilities are 13-123% larger than those estimated using publically available data. This leads to 64% higher annual GHG emissions from surface mining operations, and 30% higher overall OS GHG emissions (17 Mt) compared to that reported by industry, despite emissions reporting which uses the most up to date and recommended bottom-up approaches. Given the similarity in bottom-up reporting methods across the entire O&G sector, these results suggest that O&G CO2 emissions inventory data may be more uncertain than previously considered.

Oxidative and Toxicological Evolution of Engineered Nanoparticles with Atmospherically Relevant Coatings.

The health impacts associated with engineered nanoparticles (ENPs) released into the atmosphere have not been adequately assessed. Such impacts could potentially arise from the toxicity associated with condensable atmospheric secondary organic material (SOM), or changes in the SOM composition induced by ENPs. Here, these possibilities are evaluated by investigating the oxidative and toxicological evolution of TiO2 and SiO2 nanoparticles which have been coated with SOM from the O3 or OH initiated oxidation of α-pinene. It was found that pristine SiO2 particles were significantly more cytotoxic compared to pristine TiO2 particles. TiO2 in the dark or under UV irradiation catalytically reacted with the SOM, increasing its O/C by up to 55% over photochemically inert SiO2 while having negligible effects on the overall cytotoxicity. Conversely, the cytotoxicity associated with SiO2 coated with SOM was markedly suppressed (by a factor of 9, at the highest exposure dose) with both increased SOM coating thickness and increased photochemical aging. These suppressing effects (organic coating and photo-oxidation of organics) were attributed to a physical hindrance of SiO2-cell interactions by the SOM and enhanced SOM viscosity and hydrophilicity with continued photo-oxidation, respectively. These findings highlight the importance of atmospheric processes in altering the cytotoxicity of ENPs.

Size Dependence of the Physical Characteristics of Particles Containing Refractory Black Carbon in Diesel Vehicle Exhaust.

The number and mass size distributions of refractory black carbon (rBC) cores in particles emitted from a diesel vehicle were investigated as a function of particle mobility diameter ( dmob) using a single particle soot photometer (SP2) and a differential mobility analyzer (DMA). The thickness and mass of coatings on the rBC cores were characterized. On the basis of the SP2 and DMA results, the physical properties of particles containing rBC, including effective density (ρeff), mass-mobility scaling exponent ( Dm), dynamic shape factor (χ), and mass absorption cross section (MAC), were derived as a function of dmob. At each dmob, the count median diameter (CMD) of the rBC cores was essentially the same as their mass median diameter (MMD), which increased linearly with dmob. The mass of the rBC cores was proportional to the cubic of their dmob. However, coating thickness on rBC cores remained unchanged with dmob, with an average thickness of 28.72 ± 4.81 nm. For particles containing rBC, ρeff decreased and χ increased with dmob. The Dm of particles containing rBC was calculated to be 2.09. At 355 and 532 nm wavelengths, the MAC of the diesel particles containing rBC was inversely dependent on dmob.

High resolution mapping of nitrogen dioxide with TROPOMI: First results and validation over the Canadian oil sands.

TROPOMI, on-board the Sentinel-5 Precursor satellite is a nadir-viewing spectrometer measuring reflected sunlight in the ultraviolet, visible, near-infrared, and shortwave infrared spectral range. From these spectra several important air quality and climate-related atmospheric constituents are retrieved at an unprecedented high spatial resolution, including nitrogen dioxide (NO2). We present the first retrievals of TROPOMI NO2 over the Canadian Oil Sands, contrasting them with observations from the OMI satellite instrument, and demonstrate its ability to resolve individual plumes and highlight its potential for deriving emissions from individual mining facilities. Further, the first TROPOMI NO2 validation is presented, consisting of aircraft and surface in-situ NO2 observations, as well as ground-based remote-sensing measurements between March and May 2018. Our comparisons show that the TROPOMI NO2 vertical column densities are highly correlated with the aircraft and surface in-situ NO2 observations, and the ground-based remote-sensing measurements with a low bias (15-30 %) over the Canadian Oil Sands. PLAIN LANGUAGE SUMMARY: Nitrogen dioxide (NO2) is a pollutant that is linked to respiratory health issues and has negative environmental impacts such as soil and water acidification. Near the surface the most significant sources of NO2 are fossil fuel combustion and biomass burning. With a recently launched satellite instrument (TROPOspheric Monitoring Instrument; TROPOMI) NO2 can be measured with an unprecedented combination of accuracy, spatial coverage, and resolution. This work presents the first TROPOMI NO2 measurements near the Canadian Oil Sands and shows that these measurements have an outstanding ability to detect NO2 on a very high horizontal resolution that is unprecedented for satellite NO2 observations. Further, these satellite measurements are in excellent agreement with aircraft and ground-based measurements.

Diurnal Variability and Emission Pattern of Decamethylcyclopentasiloxane (D5) from the Application of Personal Care Products in Two North American Cities.

Decamethylcyclopentasiloxane (D5) is a cyclic volatile methyl siloxane (cVMS) that is widely used in consumer products and commonly observed in urban air. This study quantifies the ambient mixing ratios of D5 from ground sites in two North American cities (Boulder, CO, USA, and Toronto, ON, CA). From these data, we estimate the diurnal emission profile of D5 in Boulder, CO. Ambient mixing ratios were consistent with those measured at other urban locations; however, the diurnal pattern exhibited similarities with those of traffic-related compounds such as benzene. Mobile measurements and vehicle experiments demonstrate that emissions of D5 from personal care products are coincident in time and place with emissions of benzene from motor vehicles. During peak commuter times, the D5/benzene ratio (w/w) is in excess of 0.3, suggesting that the mass emission rate of D5 from personal care product usage is comparable to that of benzene due to traffic. The diurnal emission pattern of D5 is estimated using the measured D5/benzene ratio and inventory estimates of benzene emission rates in Boulder. The hourly D5 emission rate is observed to peak between 6:00 and 7:00 AM and subsequently follow an exponential decay with a time constant of 9.2 h. This profile could be used by models to constrain temporal emission patterns of personal care products.

Assessment of airborne polycyclic aromatic hydrocarbons in a megacity of South China: Spatiotemporal variability, indoor-outdoor interplay and potential human health risk.

Although a number of studies have assessed the occurrence of atmospheric polycyclic aromatic hydrocarbons (PAHs) in indoor environment, few studies have systemically examined the indoor-outdoor interplay of size-dependent particulate PAHs and potential health risk based on daily lifestyles. In the present study, size-dependent particle and gaseous samples were collected both indoors and outdoors within selected schools, offices and residences located in three districts of Guangzhou, China with different urbanization levels during the dry and wet weather seasons. Results from measurements of PAHs showed that higher total PAH concentrations occurred in residential areas than in other settings and in indoor than in outdoor environments. Compositional profiles and size distribution patterns of particle-bound PAHs were similar indoors and outdoors, predominated by 4-and 5-ring PAHs and the 0.56-1.0 µm particle fraction. Statistical analyses indicated that outdoor sources may have contributed to 38-99% and 62-100% of the variations for indoor particle-bound and gaseous PAH concentrations, respectively. Incremental life cancer risk (ILCR) from human exposure to indoor and outdoor PAHs based on different lifestyles followed the order of adults > children > adolescents > seniors. All average ILCR values for four age groups were below the lower limit of the Safe Acceptable Range (10-6). In addition, the ILCR value for adults (average: 7.2 × 10-7; 95% CI: 5.4 × 10-8‒2.5 × 10-6), estimated from outdoor air PAH levels with 24-h exposure time, was significantly higher than our assessment results (average: 5.9 × 10-7; 95% CI: 6.3 × 10-8‒1.9 × 10-6), suggesting the significance of assessing human inhalation exposure risks of indoor and outdoor PAHs in urban air based on daily lifestyles.