scholarly journals Inconsistency of ammonium–sulfate aerosol ratios with thermodynamic models in the eastern US: a possible role of organic aerosol

2017 ◽  
Vol 17 (8) ◽  
pp. 5107-5118 ◽  
Author(s):  
Rachel F. Silvern ◽  
Daniel J. Jacob ◽  
Patrick S. Kim ◽  
Eloise A. Marais ◽  
Jay R. Turner ◽  
...  

Abstract. Thermodynamic models predict that sulfate aerosol (S(VI)  ≡  H2SO4(aq) + HSO4−+ SO42−) should take up available ammonia (NH3) quantitatively as ammonium (NH4+) until the ammonium sulfate stoichiometry (NH4)2SO4 is close to being reached. This uptake of ammonia has important implications for aerosol mass, hygroscopicity, and acidity. When ammonia is in excess, the ammonium–sulfate aerosol ratio R =  [NH4+] ∕ [S(VI)] should approach 2, with excess ammonia remaining in the gas phase. When ammonia is in deficit, it should be fully taken up by the aerosol as ammonium and no significant ammonia should remain in the gas phase. Here we report that sulfate aerosol in the eastern US in summer has a low ammonium–sulfate ratio despite excess ammonia, and we show that this is at odds with thermodynamic models. The ammonium–sulfate ratio averages only 1.04 ± 0.21 mol mol−1 in the Southeast, even though ammonia is in large excess, as shown by the ammonium–sulfate ratio in wet deposition and by the presence of gas-phase ammonia. It further appears that the ammonium–sulfate aerosol ratio is insensitive to the supply of ammonia, remaining low even as the wet deposition ratio exceeds 6 mol mol−1. While the ammonium–sulfate ratio in wet deposition has increased by 5.8 % yr−1 from 2003 to 2013 in the Southeast, consistent with SO2 emission controls, the ammonium–sulfate aerosol ratio decreased by 1.4–3.0 % yr−1. Thus, the aerosol is becoming more acidic even as SO2 emissions decrease and ammonia emissions stay constant; this is incompatible with simple sulfate–ammonium thermodynamics. A tentative explanation is that sulfate particles are increasingly coated by organic material, retarding the uptake of ammonia. Indeed, the ratio of organic aerosol (OA) to sulfate in the Southeast increased from 1.1 to 2.4 g g−1 over the 2003–2013 period as sulfate decreased. We implement a simple kinetic mass transfer limitation for ammonia uptake to sulfate aerosols in the GEOS-Chem chemical transport model and find that we can reproduce both the observed ammonium–sulfate aerosol ratios and the concurrent presence of gas-phase ammonia. If sulfate aerosol becomes more acidic as OA ∕ sulfate ratios increase, then controlling SO2 emissions to decrease sulfate aerosol will not have the co-benefit of suppressing acid-catalyzed secondary organic aerosol (SOA) formation.

2016 ◽  
Author(s):  
Rachel F. Silvern ◽  
Daniel J. Jacob ◽  
Patrick S. Kim ◽  
Eloise A. Marais ◽  
Jay R. Turner

Abstract. Acid-base neutralization of sulfate aerosol (S(VI) ≡ H2SO4(aq) + HSO4− + SO42−) by ammonia (NH3) has important implications for aerosol mass, hygroscopicity, and acidity. Surface network and aircraft observations across the eastern US show that sulfate aerosol is not fully neutralized even in the presence of excess ammonia, at odds with thermodynamic equilibrium models. The sulfate aerosol neutralization ratio (f = [NH4+]/2[S(VI)]) averages only 0.51 ± 0.11 mol mol−1 at sites in the Southeast and 0.78 ± 0.13 mol mol−1 in the Northeast in summer 2013, even though ammonia is in large excess as shown by the corresponding [NH4+]/2[S(VI)] ratio in wet deposition fluxes. There is in fact no site-to-site correlation between the two quantities; the aerosol neutralization ratio in the Southeast remains in a range 0.3–0.6 mol mol−1 even as the wet deposition neutralization ratio exceeds 3 mol mol−1. While the wet deposition neutralization ratio has increased by 4.6 % a−1 from 2003 to 2013 in the Southeast US, consistent with SO2 emission controls, the aerosol neutralization ratio has decreased by 1.0–3.2 % a−1. Thus the aerosol is becoming more acidic even as SO2 emissions decrease. One possible explanation is that sulfate particles are increasingly coated by organic material, retarding the uptake of ammonia. The ratio of organic aerosol (OA) to sulfate increases over the 2003–2013 period as sulfate decreases. We implement a kinetic mass transfer limitation for ammonia uptake to sulfate aerosols in the GEOS-Chem chemical transport model and find improved agreement with surface and aircraft observations of the aerosol neutralization ratio. If sulfate aerosol becomes more acidic as OA/sulfate ratios increase, then controlling SO2 emissions to decrease sulfate aerosol will not have the co-benefit of suppressing acid-catalyzed secondary organic aerosol (SOA) formation.


2010 ◽  
Vol 10 (9) ◽  
pp. 4377-4401 ◽  
Author(s):  
H. O. T. Pye ◽  
J. H. Seinfeld

Abstract. Global production of organic aerosol from primary emissions of semivolatile (SVOCs) and intermediate (IVOCs) volatility organic compounds is estimated using the global chemical transport model, GEOS-Chem. SVOC oxidation is predicted to be a larger global source of net aerosol production than oxidation of traditional parent hydrocarbons (terpenes, isoprene, and aromatics). Using a prescribed rate constant and reduction in volatility for atmospheric oxidation, the yield of aerosol from SVOCs is predicted to be about 75% on a global, annually-averaged basis. For IVOCs, the use of a naphthalene-like surrogate with different high-NOx and low-NOx parameterizations produces a global aerosol yield of about 30%, or roughly 5 Tg/yr of aerosol. Estimates of the total global organic aerosol source presented here range between 60 and 100 Tg/yr. This range reflects uncertainty in the parameters for SVOC volatility, SVOC oxidation, SVOC emissions, and IVOC emissions, as well as wet deposition. The highest estimates result if SVOC emissions are significantly underestimated (by more than a factor of 2) or if wet deposition of the gas-phase semivolatile species is less effective than previous estimates. A significant increase in SVOC emissions, a reduction of the volatility of the SVOC emissions, or an increase in the enthalpy of vaporization of the organic aerosol all lead to an appreciable reduction of prediction/measurement discrepancy. In addition, if current primary organic aerosol (POA) inventories capture only about one-half of the SVOC emission and the Henrys Law coefficient for oxidized semivolatiles is on the order of 103 M/atm, a global estimate of OA production is not inconsistent with the top-down estimate of 140 Tg/yr by (Goldstein and Galbally, 2007). Additional information is needed to constrain the emissions and treatment of SVOCs and IVOCs, which have traditionally not been included in models.


2019 ◽  
Vol 12 (7) ◽  
pp. 2983-3000 ◽  
Author(s):  
Duseong S. Jo ◽  
Alma Hodzic ◽  
Louisa K. Emmons ◽  
Eloise A. Marais ◽  
Zhe Peng ◽  
...  

Abstract. Secondary organic aerosol derived from isoprene epoxydiols (IEPOX-SOA) is thought to contribute the dominant fraction of total isoprene SOA, but the current volatility-based lumped SOA parameterizations are not appropriate to represent the reactive uptake of IEPOX onto acidified aerosols. A full explicit modeling of this chemistry is however computationally expensive owing to the many species and reactions tracked, which makes it difficult to include it in chemistry–climate models for long-term studies. Here we present three simplified parameterizations (version 1.0) for IEPOX-SOA simulation, based on an approximate analytical/fitting solution of the IEPOX-SOA yield and formation timescale. The yield and timescale can then be directly calculated using the global model fields of oxidants, NO, aerosol pH and other key properties, and dry deposition rates. The advantage of the proposed parameterizations is that they do not require the simulation of the intermediates while retaining the key physicochemical dependencies. We have implemented the new parameterizations into the GEOS-Chem v11-02-rc chemical transport model, which has two empirical treatments for isoprene SOA (the volatility-basis-set, VBS, approach and a fixed 3 % yield parameterization), and compared all of them to the case with detailed fully explicit chemistry. The best parameterization (PAR3) captures the global tropospheric burden of IEPOX-SOA and its spatiotemporal distribution (R2=0.94) vs. those simulated by the full chemistry, while being more computationally efficient (∼5 times faster), and accurately captures the response to changes in NOx and SO2 emissions. On the other hand, the constant 3 % yield that is now the default in GEOS-Chem deviates strongly (R2=0.66), as does the VBS (R2=0.47, 49 % underestimation), with neither parameterization capturing the response to emission changes. With the advent of new mass spectrometry instrumentation, many detailed SOA mechanisms are being developed, which will challenge global and especially climate models with their computational cost. The methods developed in this study can be applied to other SOA pathways, which can allow including accurate SOA simulations in climate and global modeling studies in the future.


2010 ◽  
Vol 10 (2) ◽  
pp. 4079-4141 ◽  
Author(s):  
H. O. T. Pye ◽  
J. H. Seinfeld

Abstract. Organic aerosol from primary semivolatile and intermediate volatility compounds is estimated using a global chemical transport model. Semivolatile organic compound (SVOC, saturation concentrations between about 0.1 and 104 μg/m3) oxidation is predicted to be a much larger global source of net aerosol production than oxidation of traditional parent hydrocarbons (terpenes, isoprene, and aromatics). Using a prescribed rate constant and reduction in volatility, the yield of aerosol (defined as the net mass of aerosol formed divided by the total mass of the parent hydrocarbon emitted) from SVOCs is predicted to be about 75% on a global, annually averaged basis. Intermediate volatility compound (IVOC, saturation concentrations between about 104 and 106 μg/m3) emissions and oxidation are highly uncertain since they are not typically measured. The use of a naphthalene-like surrogate with different high-NOx and low-NOx parameterizations produces an aerosol yield of about 30% or roughly 5 Tg/yr of aerosol from IVOC oxidation on a global basis. Estimates of the total global organic aerosol source presented here range between 60 and 100 Tg/yr. This range reflects uncertainty in the parameters for SVOC volatility, SVOC oxidation, SVOC emissions, and IVOC emissions, as well as wet deposition. The highest estimates result if SVOC emissions are significantly underestimated (by more than a factor of 2) or if wet deposition of the gas-phase semivolatile species is less effective than previous estimates. Compared to a traditional non-volatile primary organic aerosol model without IVOCs, the global estimate of organic aerosol production is at most roughly 10% higher than previous studies. Additional information is needed to constrain the emissions and treatment of SVOCs and IVOCs, which have traditionally not been included in models. Comparisons to winter organic carbon observations over the US indicate that SVOC emissions are significantly underestimated by the traditional POA inventories. The degree to which IVOC emissions or other parameters are uncertain is unknown.


2021 ◽  
Vol 21 (22) ◽  
pp. 17115-17132
Author(s):  
Ksakousti Skyllakou ◽  
Pablo Garcia Rivera ◽  
Brian Dinkelacker ◽  
Eleni Karnezi ◽  
Ioannis Kioutsioukis ◽  
...  

Abstract. Significant reductions in emissions of SO2, NOx, volatile organic compounds (VOCs), and primary particulate matter (PM) took place in the US from 1990 to 2010. We evaluate here our understanding of the links between these emissions changes and corresponding changes in concentrations and health outcomes using a chemical transport model, the Particulate Matter Comprehensive Air Quality Model with Extensions (PMCAMx), for 1990, 2001, and 2010. The use of the Particle Source Apportionment Algorithm (PSAT) allows us to link the concentration reductions to the sources of the corresponding primary and secondary PM. The reductions in SO2 emissions (64 %, mainly from electric-generating units) during these 20 years have dominated the reductions in PM2.5, leading to a 45 % reduction in sulfate levels. The predicted sulfate reductions are in excellent agreement with the available measurements. Also, the reductions in elemental carbon (EC) emissions (mainly from transportation) have led to a 30 % reduction in EC concentrations. The most important source of organic aerosol (OA) through the years according to PMCAMx is biomass burning, followed by biogenic secondary organic aerosol (SOA). OA from on-road transport has been reduced by more than a factor of 3. On the other hand, changes in biomass burning OA and biogenic SOA have been modest. In 1990, about half of the US population was exposed to annual average PM2.5 concentrations above 20 µg m−3, but by 2010 this fraction had dropped to practically zero. The predicted changes in concentrations are evaluated against the observed changes for 1990, 2001, and 2010 in order to understand whether the model represents reasonably well the corresponding processes caused by the changes in emissions.


2019 ◽  
Author(s):  
Duseong S. Jo ◽  
Alma Hodzic ◽  
Louisa K. Emmons ◽  
Eloise A. Marais ◽  
Zhe Peng ◽  
...  

Abstract. Secondary organic aerosol derived from isoprene epoxydiols (IEPOX-SOA) is thought to contribute the dominant fraction of total isoprene SOA, but the current volatility-based lumped SOA parameterizations are not appropriate to represent the reactive uptake of IEPOX onto acidified aerosols. A full explicit modelling of this chemistry is however computationally expensive owing to the many species and reactions tracked, which makes it difficult to include it in chemistry climate models for long-term studies. Here we present three simplified parameterizations for IEPOX-SOA simulation, based on an approximate analytical/fitting solution of the IEPOX-SOA yield and formation timescale. The yield and timescale can then be directly calculated using the global model fields of oxidants, NO, aerosol pH and other key properties, and dry deposition rates. The advantage of the proposed parameterizations is that they do not require the simulation of the intermediates while retaining the key physico-chemical dependencies. We have implemented the new parameterizations into the GEOS-Chem v11-02-rc chemical transport model, which has two empirical treatments for isoprene SOA (the volatility basis set (VBS) approach and a fixed 3 % yield parameterization) and compared all of them to the case with detailed full chemistry. The best parameterization (PAR3) captures the global tropospheric burden of IEPOX-SOA and its spatio-temporal distribution (R2 = 0.93) vs. those simulated by the full chemistry, while being more computationally efficient (~ 5 times faster), and accurately captures the response to changes on NOx and SO2 emissions. On the other hand, the constant 3 % yield that is now default in GEOS-Chem deviates strongly (R2 = 0.65, 63 % underestimation), as does the VBS (R2 = 0.45, 78 % underestimation), with neither parameterization capturing the response to emission changes. With the advent of new mass spectrometry instrumentation, many detailed SOA mechanisms are being developed, which will challenge global and especially climate models with their computational cost. The methods developed in this study can be applied to other SOA pathways, which can allow including accurate SOA simulations in climate and global modeling studies in the future.


2021 ◽  
Author(s):  
Ksakousti Skyllakou ◽  
Pablo Garcia Rivera ◽  
Brian Dinkelacker ◽  
Eleni Karnezi ◽  
Ioannis Kioutsioukis ◽  
...  

Abstract. Significant reductions of emissions of SO2, NOx, volatile organic compounds (VOCs) and primary particulate matter (PM) took place in the US from 1990 to 2010. We evaluate here our understanding of the links between these emissions changes and corresponding changes in concentrations and health outcomes using a chemical transport model, the Particulate Matter Comprehensive Air Quality Model with Extensions (PMCAMx) with the Particle Source Apportionment Algorithm (PSAT). Results for 1990, 2001 and 2010 are presented. The reductions in SO2 emissions (64 %, mainly from electric generating units) during these 20 years have dominated the reductions in PM2.5 leading to a 45 % reduction in sulfate levels. The predicted sulfate reductions are in excellent agreement with the available measurements. Also, the reductions in elemental carbon (EC) emissions (mainly from transportation) have led to a 30 % reduction of EC concentrations. The most important source of organic aerosol (OA) through the years according to PMCAMx is biomass burning, followed by biogenic secondary organic aerosol (SOA). OA from on-road transport has been reduced by more than a factor of three. On the other hand, changes in biomass burning OA and biogenic SOA have been modest. In 1990, about half of the US population was exposed to annual-average PM2.5 concentrations above 20 μg m−3, but by 2010 this fraction had dropped to practically zero. The predicted changes in concentrations are evaluated against the observed changes for 1990, 2001, and 2010, in order to understand if the model represents reasonably well the corresponding processes caused by the changes in emissions.


2012 ◽  
Vol 12 (2) ◽  
pp. 5939-6018
Author(s):  
C. A. Stroud ◽  
M. D. Moran ◽  
P. A. Makar ◽  
S. Gong ◽  
W. Gong ◽  
...  

Abstract. Observations from the 2007 Border Air Quality and Meteorology Study (BAQS-Met 2007) in southern Ontario (ON), Canada, were used to evaluate Environment Canada's regional chemical transport model predictions of primary organic aerosol (POA). Environment Canada's operational numerical weather prediction model and the 2006 Canadian and 2005 US national emissions inventories were used as input to the chemical transport model (named AURAMS). Particle-component-based factor analysis was applied to aerosol mass spectrometer measurements made at one urban site (Windsor, ON) and two rural sites (Harrow and Bear Creek, ON) to derive hydrocarbon-like organic aerosol (HOA) factors. Co-located carbon monoxide (CO), PM2.5 black carbon (BC), and PM1 SO4 measurements were also used for evaluation and interpretation, permitting a detailed diagnostic model evaluation. At the urban site, good agreement was observed for the comparison of daytime campaign PM1 POA and HOA mean values: 1.1 μg m−3 vs. 1.2 μg m−3, respectively. However, a POA overprediction was evident on calm nights due to an overly-stable model surface layer. Biases in model POA predictions trended from positive to negative with increasing HOA values. This trend has several possible explanations, including (1) underweighting of urban locations in particulate matter (PM) spatial surrogate fields, (2) overly-coarse model grid spacing for resolving urban-scale sources, and (3) lack of a model particle POA evaporation process during dilution of vehicular POA tail-pipe emissions to urban scales. Furthermore, a trend in POA bias was observed at the urban site as a function of the BC/HOA ratio, suggesting a possible association of POA underprediction for diesel combustion sources. For several time periods, POA overprediction was also observed for sulphate-rich plumes, suggesting that our model POA fractions for the PM2.5 chemical speciation profiles may be too high for these point sources. At the rural Harrow site, significant underpredictions in PM1 POA concentration were found compared to observed HOA concentration and were associated, based on back-trajectory analysis, with (1) transport from the Detroit/Windsor urban complex, (2) longer-range transport from the US Midwest, and (3) biomass burning. Daytime CO concentrations were significantly overpredicted at Windsor but were unbiased at Harrow. Collectively, these biases provide support for a hypothesis that combines a current underweighting of PM spatial surrogate fields for urban locations with insufficient model vertical mixing for sources close to the urban measurement sites. The magnitude of the area POA emissions sources in the US and Canadian inventories (e.g., food cooking, road and soil dust, waste disposal burning) suggests that more effort should be placed at reducing uncertainties in these sectors, especially spatial and temporal surrogates.


2007 ◽  
Vol 7 (3) ◽  
pp. 9053-9092 ◽  
Author(s):  
C. R. Hoyle ◽  
T. Berntsen ◽  
G. Myhre ◽  
I. S. A. Isaksen

Abstract. The global chemical transport model Oslo CTM2 has been extended to include the formation, transport and deposition of secondary organic aerosol (SOA). Precursor hydrocarbons which are oxidised to form condensible species include both biogenic species such as terpenes and isoprene, as well as species emitted predominantly by anthropogenic activities (toluene, m-xylene, methylbenzene and other aromatics). A model simulation for 2004 gives an annual global SOA production of approximately 55 Tg. Of this total, 2.5 Tg is found to consist of the oxidation products of anthropogenically emitted hydrocarbons, and about 15 Tg is formed by the oxidation products of isoprene. The global production of SOA is increased to about 76 Tg yr−1 by allowing semi-volatile species to condense on ammonium sulphate aerosol. This brings modelled organic aerosol values closer to those observed, however observations in Europe remain significantly underestimated, raising the possibility of an unaccounted for SOA source. Allowing SOA to form on ammonium sulphate aerosol increases the contribution of anthropogenic SOA from about 4.5% to almost 9% of the total production. The importance of NO3 as an oxidant of SOA precursors is found to vary regionally, causing up to 50%–60% of the total amount of SOA near the surface in polluted regions and less than 25% in more remote areas. This study underscores the need for SOA to be represented in a more realistic way in global aerosol models in order to better reproduce observations of organic aerosol burdens in industrialised and biomass burning regions.


2021 ◽  
Author(s):  
Ruochong Xu ◽  
Joel A. Thornton ◽  
Ben H. Lee ◽  
Yanxu Zhang ◽  
Lyatt Jaeglé ◽  
...  

Abstract. We evaluate monoterpene-derived peroxy radical (MT-RO2) unimolecular autoxidation and self and cross reactions with other RO2 in the GEOS-Chem global chemical transport model. Formation of associated highly oxygenated organic molecule (HOM) and accretion products are tracked in competition with other bimolecular reactions. Autoxidation is the dominant fate up to 6–8 km for first-generation MT-RO2 which can undergo unimolecular H-shifts. Reaction with NO can be a more common fate for H-shift rate constants < 0.1 s−1 or at altitudes higher than 8 km due to the imposed Arrhenius temperature dependence of unimolecular H-shifts. For MT-derived HOM-RO2, generated by multi-step autoxidation of first-generation MT-RO2, reaction with other RO2 is predicted to be the major fate throughout most of the boreal and tropical forested regions, while reaction with NO dominates in temperate and subtropical forests of the Northern Hemisphere. The newly added reactions result in ~4 % global average decrease of HO2 and RO2 mainly due to faster self-/cross-reactions of MT-RO2, but the impact upon HO2/OH/NOx abundances is only important in the planetary boundary layer (PBL) over portions of tropical forests. Within the bounds of formation kinetics and HOM photochemical lifetime constraints from laboratory studies, predicted HOM concentrations in MT-rich regions and seasons reach 10 % or even exceed total organic aerosol as predicted by the standard GEOS-Chem model. Comparisons to observations reveal large uncertainties remain for key reaction parameters and processes, especially the photochemical lifetime of HOM and associated accretion products. Using the highest reported yields and H-shift rate constants of MT-RO2 that undergo autoxidation, HOM concentrations tend to exceed the limited set of observations. Similarly, we infer that RO2 cross reactions rate constants near the gas-kinetic limit with accretion product branching greater than ~0.25 are inconsistent with total organic aerosol unless there is rapid decomposition of accretion products, the accretion products have saturation vapor concentrations > > 1 μg m−3, or modeled MT emission rates are overestimated. This work suggests further observations and laboratory studies related to MT-RO2 derived HOM and gas-phase accretion product formation kinetics, and especially their atmospheric fate, such as gas-particle partitioning, multi-phase chemistry, and net SOA formation, are needed.


Sign in / Sign up

Export Citation Format

Share Document