scholarly journals In-cloud oxalate formation in the global troposphere: a 3-D modeling study

2011 ◽  
Vol 11 (1) ◽  
pp. 485-530 ◽  
Author(s):  
S. Myriokefalitakis ◽  
K. Tsigaridis ◽  
N. Mihalopoulos ◽  
J. Sciare ◽  
A. Nenes ◽  
...  

Abstract. Organic acids attract increasing attention as contributors to atmospheric acidity, secondary organic aerosol mass and aerosol hygroscopicity. Oxalic acid is globally the most abundant dicarboxylic acid, formed via chemical oxidation of gas-phase precursors in the aqueous phase of aerosols and droplets. Its lifecycle and atmospheric global distribution remain highly uncertain and are the focus of this study. The first global spatial and temporal distribution of oxalate, simulated using a state-of-the-art aqueous phase chemical scheme embedded within the global 3-dimensional chemistry/transport model TM4-ECPL, is here presented. The model accounts for comprehensive gas-phase chemistry and its coupling with major aerosol constituents (including secondary organic aerosol). Model results are consistent with ambient observations of oxalate at rural and remote locations (slope = 0.83 ± 0.06, r2 = 0.67, N = 106) and suggest that aqueous phase chemistry contributes significantly to the global atmospheric burden of secondary organic aerosol. In TM4-ECPL most oxalate is formed in-clouds and less than 10% is produced in aerosol water. About 61% of the oxalate is removed via wet deposition, 35% by in-cloud reaction with hydroxyl radical and 4% by dry deposition. The global oxalate net chemical production is calculated to be about 17–27 Tg yr−1 with almost 91% originating from biogenic hydrocarbons, mainly isoprene. This condensed phase net source of oxalate in conjunction with a global mean turnover time against deposition of about 5 days, maintain oxalate's global tropospheric burden of 0.24–0.39 Tg that is about 13–19% of calculated total organic aerosol burden.

2011 ◽  
Vol 11 (12) ◽  
pp. 5761-5782 ◽  
Author(s):  
S. Myriokefalitakis ◽  
K. Tsigaridis ◽  
N. Mihalopoulos ◽  
J. Sciare ◽  
A. Nenes ◽  
...  

Abstract. Organic acids attract increasing attention as contributors to atmospheric acidity, secondary organic aerosol mass and aerosol hygroscopicity. Oxalic acid is globally the most abundant dicarboxylic acid, formed via chemical oxidation of gas-phase precursors in the aqueous phase of aerosols and droplets. Its lifecycle and atmospheric global distribution remain highly uncertain and are the focus of this study. The first global spatial and temporal distribution of oxalate, simulated using a state-of-the-art aqueous-phase chemical scheme embedded within the global 3-dimensional chemistry/transport model TM4-ECPL, is here presented. The model accounts for comprehensive gas-phase chemistry and its coupling with major aerosol constituents (including secondary organic aerosol). Model results are consistent with ambient observations of oxalate at rural and remote locations (slope = 1.16 ± 0.14, r2 = 0.36, N = 114) and suggest that aqueous-phase chemistry contributes significantly to the global atmospheric burden of secondary organic aerosol. In TM4-ECPL most oxalate is formed in-cloud and less than 5 % is produced in aerosol water. About 62 % of the oxalate is removed via wet deposition, 30 % by in-cloud reaction with hydroxyl radical, 4 % by in-cloud reaction with nitrate radical and 4 % by dry deposition. The in-cloud global oxalate net chemical production is calculated to be about 21–37 Tg yr−1 with almost 79 % originating from biogenic hydrocarbons, mainly isoprene. This condensed phase net source of oxalate in conjunction with a global mean turnover time against deposition of about 5 days, maintain oxalate's global tropospheric burden of 0.2–0.3 Tg, i.e. 0.05–0.1 Tg-C that is about 5–9 % of model-calculated water soluble organic carbon burden.


2012 ◽  
Vol 12 (4) ◽  
pp. 9857-9901 ◽  
Author(s):  
B. N. Murphy ◽  
N. M. Donahue ◽  
C. Fountoukis ◽  
M. Dall'Osto ◽  
C. O'Dowd ◽  
...  

Abstract. Multigenerational oxidation chemistry of atmospheric organic compounds and its effects on aerosol loadings and chemical composition is investigated by implementing the Two-Dimensional Volatility Basis Set (2-D-VBS) in a Lagrangian host chemical transport model. Three model formulations were chosen to explore the complex interactions between functionalization and fragmentation processes during gas-phase oxidation of organic compounds by the hydroxyl radical. The base case model employs a conservative transformation by assuming a reduction of one order of magnitude in effective saturation concentration and an increase of oxygen content by one or two oxygen atoms per oxidation generation. A second scheme simulates functionalization in more detail using group contribution theory to estimate the effects of oxygen addition to the carbon backbone on the compound volatility. Finally, a fragmentation scheme is added to the detailed functionalization scheme to create a functionalization-fragmentation parameterization. Two condensed-phase chemistry pathways are also implemented as additional sensitivity tests to simulate (1) heterogeneous oxidation via OH uptake to the particle-phase and (2) aqueous-phase chemistry of glyoxal and methylglyoxal. The model is applied to summer and winter periods at three sites where observations of organic aerosol (OA) mass and O:C were obtained during the European Integrated Project on Aerosol Cloud Climate and Air Quality Interactions (EUCAARI) campaigns. The base case model reproduces observed mass concentrations and O:C well, with fractional errors (FE) lower than 55% and 25%, respectively. The detailed functionalization scheme tends to overpredict OA concentrations, especially in the summertime, and also underpredicts O:C by approximately a factor of 2. The detailed functionalization model with fragmentation agrees well with the observations for OA concentration, but still underpredicts O:C. Both heterogeneous oxidation and aqueous-phase processing have small effects on OA levels but heterogeneous oxidation, as implemented here, does enhance O:C by about 0.1. The different schemes result in very different fractional attribution for OA between anthropogenic and biogenic sources.


2011 ◽  
Vol 45 (5) ◽  
pp. 1126-1137 ◽  
Author(s):  
Harshal M. Parikh ◽  
Annmarie G. Carlton ◽  
William Vizuete ◽  
Richard M. Kamens

2016 ◽  
Vol 16 (3) ◽  
pp. 1603-1618 ◽  
Author(s):  
E. A. Marais ◽  
D. J. Jacob ◽  
J. L. Jimenez ◽  
P. Campuzano-Jost ◽  
D. A. Day ◽  
...  

Abstract. Isoprene emitted by vegetation is an important precursor of secondary organic aerosol (SOA), but the mechanism and yields are uncertain. Aerosol is prevailingly aqueous under the humid conditions typical of isoprene-emitting regions. Here we develop an aqueous-phase mechanism for isoprene SOA formation coupled to a detailed gas-phase isoprene oxidation scheme. The mechanism is based on aerosol reactive uptake coefficients (γ) for water-soluble isoprene oxidation products, including sensitivity to aerosol acidity and nucleophile concentrations. We apply this mechanism to simulation of aircraft (SEAC4RS) and ground-based (SOAS) observations over the southeast US in summer 2013 using the GEOS-Chem chemical transport model. Emissions of nitrogen oxides (NOx  ≡  NO + NO2) over the southeast US are such that the peroxy radicals produced from isoprene oxidation (ISOPO2) react significantly with both NO (high-NOx pathway) and HO2 (low-NOx pathway), leading to different suites of isoprene SOA precursors. We find a mean SOA mass yield of 3.3 % from isoprene oxidation, consistent with the observed relationship of total fine organic aerosol (OA) and formaldehyde (a product of isoprene oxidation). Isoprene SOA production is mainly contributed by two immediate gas-phase precursors, isoprene epoxydiols (IEPOX, 58 % of isoprene SOA) from the low-NOx pathway and glyoxal (28 %) from both low- and high-NOx pathways. This speciation is consistent with observations of IEPOX SOA from SOAS and SEAC4RS. Observations show a strong relationship between IEPOX SOA and sulfate aerosol that we explain as due to the effect of sulfate on aerosol acidity and volume. Isoprene SOA concentrations increase as NOx emissions decrease (favoring the low-NOx pathway for isoprene oxidation), but decrease more strongly as SO2 emissions decrease (due to the effect of sulfate on aerosol acidity and volume). The US Environmental Protection Agency (EPA) projects 2013–2025 decreases in anthropogenic emissions of 34 % for NOx (leading to a 7 % increase in isoprene SOA) and 48 % for SO2 (35 % decrease in isoprene SOA). Reducing SO2 emissions decreases sulfate and isoprene SOA by a similar magnitude, representing a factor of 2 co-benefit for PM2.5 from SO2 emission controls.


2020 ◽  
Vol 20 (24) ◽  
pp. 15513-15535
Author(s):  
Thomas Berkemeier ◽  
Masayuki Takeuchi ◽  
Gamze Eris ◽  
Nga L. Ng

Abstract. Organic aerosol constitutes a major fraction of the global aerosol burden and is predominantly formed as secondary organic aerosol (SOA). Environmental chambers have been used extensively to study aerosol formation and evolution under controlled conditions similar to the atmosphere, but quantitative prediction of the outcome of these experiments is generally not achieved, which signifies our lack in understanding of these results and limits their portability to large-scale models. In general, kinetic models employing state-of-the-art explicit chemical mechanisms fail to describe the mass concentration and composition of SOA obtained from chamber experiments. Specifically, chemical reactions including the nitrate radical (NO3) are a source of major uncertainty for assessing the chemical and physical properties of oxidation products. Here, we introduce a kinetic model that treats gas-phase chemistry, gas–particle partitioning, particle-phase oligomerization, and chamber vapor wall loss and use it to describe the oxidation of the monoterpenes α-pinene and limonene with NO3. The model can reproduce aerosol mass and nitration degrees in experiments using either pure precursors or their mixtures and infers volatility distributions of products, branching ratios of reactive intermediates and particle-phase reaction rates. The gas-phase chemistry in the model is based on the Master Chemical Mechanism (MCM) but trades speciation of single compounds for the overall ability of quantitatively describing SOA formation by using a lumped chemical mechanism. The complex branching into a multitude of individual products in MCM is replaced in this model with product volatility distributions and detailed peroxy (RO2) and alkoxy (RO) radical chemistry as well as amended by a particle-phase oligomerization scheme. The kinetic parameters obtained in this study are constrained by a set of SOA formation and evaporation experiments conducted in the Georgia Tech Environmental Chamber (GTEC) facility. For both precursors, we present volatility distributions of nitrated and non-nitrated reaction products that are obtained by fitting the kinetic model systematically to the experimental data using a global optimization method, the Monte Carlo genetic algorithm (MCGA). The results presented here provide new mechanistic insight into the processes leading to formation and evaporation of SOA. Most notably, the model suggests that the observed slow evaporation of SOA could be due to reversible oligomerization reactions in the particle phase. However, the observed non-linear behavior of precursor mixtures points towards a complex interplay of reversible oligomerization and kinetic limitations of mass transport in the particle phase, which is explored in a model sensitivity study. The methodologies described in this work provide a basis for quantitative analysis of multi-source data from environmental chamber experiments but also show that a large data pool is needed to fully resolve uncertainties in model parameters.


2014 ◽  
Vol 7 (1) ◽  
pp. 379-429 ◽  
Author(s):  
F. Couvidat ◽  
K. Sartelet

Abstract. The Secondary Organic Aerosol Processor (SOAP v1.0) model is presented. This model is designed to be modular with different user options depending on the computing time and the complexity required by the user. This model is based on the molecular surrogate approach, in which each surrogate compound is associated with a molecular structure to estimate some properties and parameters (hygroscopicity, absorption on the aqueous phase of particles, activity coefficients, phase separation). Each surrogate can be hydrophilic (condenses only on the aqueous phase of particles), hydrophobic (condenses only on the organic phase of particles) or both (condenses on both the aqueous and the organic phases of particles). Activity coefficients are computed with the UNIFAC thermodynamic model for short-range interactions and with the AIOMFAC parameterization for medium and long-range interactions between electrolytes and organic compounds. Phase separation is determined by Gibbs energy minimization. The user can choose between an equilibrium and a dynamic representation of the organic aerosol. In the equilibrium representation, compounds in the particle phase are assumed to be at equilibrium with the gas phase. However, recent studies show that the organic aerosol (OA) is not at equilibrium with the gas phase because the organic phase could be semi-solid (very viscous liquid phase). The condensation or evaporation of organic compounds could then be limited by the diffusion in the organic phase due to the high viscosity. A dynamic representation of secondary organic aerosols (SOA) is used with OA divided into layers, the first layer at the center of the particle (slowly reaches equilibrium) and the final layer near the interface with the gas phase (quickly reaches equilibrium).


2015 ◽  
Vol 8 (4) ◽  
pp. 1111-1138 ◽  
Author(s):  
F. Couvidat ◽  
K. Sartelet

Abstract. In this paper the Secondary Organic Aerosol Processor (SOAP v1.0) model is presented. This model determines the partitioning of organic compounds between the gas and particle phases. It is designed to be modular with different user options depending on the computation time and the complexity required by the user. This model is based on the molecular surrogate approach, in which each surrogate compound is associated with a molecular structure to estimate some properties and parameters (hygroscopicity, absorption into the aqueous phase of particles, activity coefficients and phase separation). Each surrogate can be hydrophilic (condenses only into the aqueous phase of particles), hydrophobic (condenses only into the organic phases of particles) or both (condenses into both the aqueous and the organic phases of particles). Activity coefficients are computed with the UNIFAC (UNIversal Functional group Activity Coefficient; Fredenslund et al., 1975) thermodynamic model for short-range interactions and with the Aerosol Inorganic–Organic Mixtures Functional groups Activity Coefficients (AIOMFAC) parameterization for medium- and long-range interactions between electrolytes and organic compounds. Phase separation is determined by Gibbs energy minimization. The user can choose between an equilibrium representation and a dynamic representation of organic aerosols (OAs). In the equilibrium representation, compounds in the particle phase are assumed to be at equilibrium with the gas phase. However, recent studies show that the organic aerosol is not at equilibrium with the gas phase because the organic phases could be semi-solid (very viscous liquid phase). The condensation–evaporation of organic compounds could then be limited by the diffusion in the organic phases due to the high viscosity. An implicit dynamic representation of secondary organic aerosols (SOAs) is available in SOAP with OAs divided into layers, the first layer being at the center of the particle (slowly reaches equilibrium) and the final layer being near the interface with the gas phase (quickly reaches equilibrium). Although this dynamic implicit representation is a simplified approach to model condensation–evaporation with a low number of layers and short CPU (central processing unit) time, it shows good agreements with an explicit representation of condensation–evaporation (no significant differences after a few hours of condensation).


2012 ◽  
Vol 12 (22) ◽  
pp. 10797-10816 ◽  
Author(s):  
B. N. Murphy ◽  
N. M. Donahue ◽  
C. Fountoukis ◽  
M. Dall'Osto ◽  
C. O'Dowd ◽  
...  

Abstract. Multigenerational oxidation chemistry of atmospheric organic compounds and its effects on aerosol loadings and chemical composition is investigated by implementing the Two-Dimensional Volatility Basis Set (2-D-VBS) in a Lagrangian host chemical transport model. Three model formulations were chosen to explore the complex interactions between functionalization and fragmentation processes during gas-phase oxidation of organic compounds by the hydroxyl radical. The base case model employs a conservative transformation by assuming a reduction of one order of magnitude in effective saturation concentration and an increase of oxygen content by one or two oxygen atoms per oxidation generation. A second scheme simulates functionalization in more detail using group contribution theory to estimate the effects of oxygen addition to the carbon backbone on the compound volatility. Finally, a fragmentation scheme is added to the detailed functionalization scheme to create a functionalization-fragmentation parameterization. Two condensed-phase chemistry pathways are also implemented as additional sensitivity tests to simulate (1) heterogeneous oxidation via OH uptake to the particle-phase and (2) aqueous-phase chemistry of glyoxal and methylglyoxal. The model is applied to summer and winter periods at three sites where observations of organic aerosol (OA) mass and O:C were obtained during the European Integrated Project on Aerosol Cloud Climate and Air Quality Interactions (EUCAARI) campaigns. The base case model reproduces observed mass concentrations and O:C well, with fractional errors (FE) lower than 55% and 25%, respectively. The detailed functionalization scheme tends to overpredict OA concentrations, especially in the summertime, and also underpredicts O:C by approximately a factor of 2. The detailed functionalization model with fragmentation agrees well with the observations for OA concentration, but still underpredicts O:C. Both heterogeneous oxidation and aqueous-phase processing have small effects on OA levels but heterogeneous oxidation, as implemented here, does enhance O:C by about 0.1. The different schemes result in very different fractional attribution for OA between anthropogenic and biogenic sources.


2015 ◽  
Vol 15 (21) ◽  
pp. 32005-32047 ◽  
Author(s):  
E. A. Marais ◽  
D. J. Jacob ◽  
J. L. Jimenez ◽  
P. Campuzano-Jost ◽  
D. A. Day ◽  
...  

Abstract. Isoprene emitted by vegetation is an important precursor of secondary organic aerosol (SOA), but the mechanism and yields are uncertain. Aerosol is prevailingly aqueous under the humid conditions typical of isoprene-emitting regions. Here we develop an aqueous-phase mechanism for isoprene SOA formation coupled to a detailed gas-phase isoprene oxidation scheme. The mechanism is based on aerosol reactive uptake probabilities (γ) for water-soluble isoprene oxidation products, including sensitivity to aerosol acidity and nucleophile concentrations. We apply this mechanism to simulation of aircraft (SEAC4RS) and ground-based (SOAS) observations over the Southeast US in summer 2013 using the GEOS-Chem chemical transport model. Emissions of nitrogen oxides (NOx ≡ NO + NO2) over the Southeast US are such that the peroxy radicals produced from isoprene oxidation (ISOPO2) react significantly with both NO (high-NOx pathway) and HO2 (low-NOx pathway), leading to different suites of isoprene SOA precursors. We find a mean SOA mass yield of 3.3 % from isoprene oxidation, consistent with the observed relationship of OA and formaldehyde (a product of isoprene oxidation). The yield is mainly contributed by two immediate gas-phase precursors, isoprene epoxydiols (IEPOX, 58 % of isoprene SOA) from the low-NOx pathway and glyoxal (28 %) from both low- and high-NOx pathways. This speciation is consistent with observations of IEPOX SOA from SOAS and SEAC4RS. Observations show a strong relationship between IEPOX SOA and sulfate aerosol that we explain as due to the indirect effect of sulfate on aerosol acidity and volume, rather than a direct mechanistic role for sulfate. Isoprene SOA concentrations increase as NOx emissions decrease (favoring the low-NOx pathway for isoprene oxidation), but decrease as SO2 emissions decrease (due to the effect of sulfate on aerosol acidity and volume). The US EPA projects 2013–2025 decreases in anthropogenic emissions of 34 % for NOx (leading to 7 % increase in isoprene SOA) and 48 % for SO2 (35 % decrease in isoprene SOA). The combined projected decreases in NOx and SO2 emissions reduce isoprene SOA yields from 3.3 to 2.3 %. Reducing SO2 emissions decreases sulfate and isoprene SOA by a similar magnitude, representing a factor of 2 co-benefit for PM2.5 from SO2 emission controls.


Sign in / Sign up

Export Citation Format

Share Document