scholarly journals Improving the representation of secondary organic aerosol (SOA) in the MOZART-4 global chemical transport model

2012 ◽  
Vol 5 (4) ◽  
pp. 4187-4232 ◽  
Author(s):  
A. Mahmud ◽  
K. C. Barsanti

Abstract. The secondary organic aerosol (SOA) module in the Model for Ozone and Related chemical Tracers, version 4 (MOZART-4) has been updated by replacing existing two-product (2p) parameters with those obtained from two-product volatility basis set (2p-VBS) fits, and by treating SOA formation from the following volatile organic compounds (VOCs): isoprene, propene and lumped alkenes. Strong seasonal and spatial variations in global SOA distributions were demonstrated, with significant differences in the predicted concentrations between the base-case and updated model versions. The base-case MOZART-4 predicted annual average SOA of 0.36 ± 0.50 μg m−3 in South America, 0.31 ± 0.38 μg m−3 in Indonesia, 0.09 ± 0.05 μg m−3 in the USA, and 0.12 ± 0.07 μg m−3 in Europe. Concentrations from the updated versions of the model showed a~marked increase in annual average SOA. Using the updated set of parameters alone (MZ4-v1) increased annual average SOA by ~8%, ~16%, ~56%, and ~108% from the base-case in South America, Indonesia, USA, and Europe, respectively. Treatment of additional parent VOCs (MZ4-v2) resulted in an even more dramatic increase of ~178–406% in annual average SOA for these regions over the base-case. The increases in predicted SOA concentrations further resulted in increases in corresponding SOA contributions to annual average total aerosol optical depth (AOD) by <1% for MZ4-v1 and ~1–6% for MZ4-v2. Estimated global SOA production was ~6.6 Tg yr−1 and ~19.1 Tg yr−1 with corresponding burdens of ~0.24 Tg and ~0.59 Tg using MZ4-v1 and MZ4-v2, respectively. The SOA budgets predicted in the current study fall well within reported ranges for similar modeling studies, 6.7 to 96 Tg yr−1, but are lower than recently reported observationally-constrained values, 50 to 380 Tg yr−1. With MZ4-v2, simulated SOA concentrations at the surface were also in reasonable agreement with comparable modeling studies and observations. Concentrations of estimated organic aerosol (OA) at the surface, however, showed under-prediction in Europe and over-prediction in the Amazonian regions and Malaysian Borneo during certain months of the year. Overall, the updated version of MOZART-4, MZ4-v2, showed consistently better skill in predicting SOA and OA levels and spatial distributions as compared with unmodified MOZART-4. The MZ4-v2 updates may be particularly important when MOZART-4 output is used to generate boundary conditions for regional air quality simulations that require more accurate representation of SOA concentrations and distributions.

2013 ◽  
Vol 6 (4) ◽  
pp. 961-980 ◽  
Author(s):  
A. Mahmud ◽  
K. Barsanti

Abstract. The secondary organic aerosol (SOA) module in the Model for Ozone and Related Chemical Tracers, version 4 (MOZART-4) was updated by replacing existing two-product (2p) parameters with those obtained from two-product volatility basis set (2p-VBS) fits (MZ4-C1), and by treating SOA formation from the following additional volatile organic compounds (VOCs): isoprene, propene and lumped alkenes (MZ4-C2). Strong seasonal and spatial variations in global SOA distributions were demonstrated, with significant differences in the predicted concentrations between the base case and updated model simulations. Updates to the model resulted in significant increases in annual average SOA mass concentrations, particularly for the MZ4-C2 simulation in which the additional SOA precursor VOCs were treated. Annual average SOA concentrations predicted by the MZ4-C2 simulation were 1.00 ± 1.04 μg m−3 in South America, 1.57 ± 1.88 μg m−3 in Indonesia, 0.37 ± 0.27 μg m−3 in the USA, and 0.47 ± 0.29 μg m−3 in Europe with corresponding increases of 178, 406, 311 and 292% over the base-case simulation, respectively, primarily due to inclusion of isoprene. The increases in predicted SOA mass concentrations resulted in corresponding increases in SOA contributions to annual average total aerosol optical depth (AOD) by ~ 1–6%. Estimated global SOA production was 5.8, 6.6 and 19.1 Tg yr−1 with corresponding burdens of 0.22, 0.24 and 0.59 Tg for the base-case, MZ4-C1 and MZ4-C2 simulations, respectively. The predicted SOA budgets fell well within reported ranges for comparable modeling studies, 6.7 to 96 Tg yr−1, but were lower than recently reported observationally constrained values, 50 to 380 Tg yr−1. For MZ4-C2, simulated SOA concentrations at the surface also were in reasonable agreement with comparable modeling studies and observations. Total organic aerosol (OA) mass concentrations at the surface, however, were slightly over-predicted in Europe, Amazonian regions and Malaysian Borneo (Southeast Asia) during certain months of the year, and under-predicted in most sites in Asia; relative to those regions, the model performed better for sites in North America. Overall, with the inclusion of additional SOA precursors (MZ4-C2), namely isoprene, MOZART-4 showed consistently better skill (NMB (normalized mean bias) of −11 vs. −26%) in predicting total OA levels and spatial distributions of SOA as compared with unmodified MOZART-4. Treatment of SOA formation by these known precursors (isoprene, propene and lumped alkenes) may be particularly important when MOZART-4 output is used to generate boundary conditions for regional air quality simulations that require more accurate representation of SOA concentrations and distributions.


2012 ◽  
Vol 12 (18) ◽  
pp. 8499-8527 ◽  
Author(s):  
R. Bergström ◽  
H. A. C. Denier van der Gon ◽  
A. S. H. Prévôt ◽  
K. E. Yttri ◽  
D. Simpson

Abstract. A new organic aerosol module has been implemented into the EMEP chemical transport model. Four different volatility basis set (VBS) schemes have been tested in long-term simulations for Europe, covering the six years 2002–2007. Different assumptions regarding partitioning of primary organic aerosol and aging of primary semi-volatile and intermediate volatility organic carbon (S/IVOC) species and secondary organic aerosol (SOA) have been explored. Model results are compared to filter measurements, aerosol mass spectrometry (AMS) data and source apportionment studies, as well as to other model studies. The present study indicates that many different sources contribute significantly to organic aerosol in Europe. Biogenic and anthropogenic SOA, residential wood combustion and vegetation fire emissions may all contribute more than 10% each over substantial parts of Europe. This study shows smaller contributions from biogenic SOA to organic aerosol in Europe than earlier work, but relatively greater anthropogenic SOA. Simple VBS based organic aerosol models can give reasonably good results for summer conditions but more observational studies are needed to constrain the VBS parameterisations and to help improve emission inventories. The volatility distribution of primary emissions is one important issue for further work. Emissions of volatile organic compounds from biogenic sources are also highly uncertain and need further validation. We can not reproduce winter levels of organic aerosol in Europe, and there are many indications that the present emission inventories substantially underestimate emissions from residential wood combustion in large parts of Europe.


2010 ◽  
Vol 10 (2) ◽  
pp. 525-546 ◽  
Author(s):  
A. P. Tsimpidi ◽  
V. A. Karydis ◽  
M. Zavala ◽  
W. Lei ◽  
L. Molina ◽  
...  

Abstract. New primary and secondary organic aerosol modules have been added to PMCAMx, a three dimensional chemical transport model (CTM), for use with the SAPRC99 chemistry mechanism based on recent smog chamber studies. The new modelling framework is based on the volatility basis-set approach: both primary and secondary organic components are assumed to be semivolatile and photochemically reactive and are distributed in logarithmically spaced volatility bins. This new framework with the use of the new volatility basis parameters for low-NOx and high-NOx conditions tends to predict 4–6 times higher anthropogenic SOA concentrations than those predicted with the older generation of models. The resulting PMCAMx-2008 was applied in Mexico City Metropolitan Area (MCMA) for approximately a week during April 2003 during a period of very low regional biomass burning impact. The emission inventory, which uses as a starting point the MCMA 2004 official inventory, is modified and the primary organic aerosol (POA) emissions are distributed by volatility based on dilution experiments. The predicted organic aerosol (OA) concentrations peak in the center of Mexico City, reaching values above 40 μg m−3. The model predictions are compared with the results of the Positive Matrix Factorization (PMF) analysis of the Aerosol Mass Spectrometry (AMS) observations. The model reproduces both Hydrocarbon-like Organic Aerosol (HOA) and Oxygenated Organic Aerosol (OOA) concentrations and diurnal profiles. The small OA underprediction during the rush-hour periods and overprediction in the afternoon suggest potential improvements to the description of fresh primary organic emissions and the formation of the oxygenated organic aerosols, respectively, although they may also be due to errors in the simulation of dispersion and vertical mixing. However, the AMS OOA data are not specific enough to prove that the model reproduces the organic aerosol observations for the right reasons. Other combinations of contributions of primary and secondary organic aerosol production rates may lead to similar results. The model results strongly suggest that, during the simulated period, transport of OA from outside the city was a significant contributor to the observed OA levels. Future simulations should use a larger domain in order to test whether the regional OA can be predicted with current SOA parameterizations. Sensitivity tests indicate that the predicted OA concentration is especially sensitive to the volatility distribution of the emissions in the lower volatility bins.


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.


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.


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.


2020 ◽  
Author(s):  
Georgia N. Theodoritsi ◽  
Giancarlo Ciarelli ◽  
Spyros N. Pandis

Abstract. A source-resolved three-dimensional chemical transport model, PMCAMx-SR, was applied in the continental U.S. to investigate the contribution of the various components (primary and secondary) of biomass burning organic aerosol (bbOA) to organic aerosol levels. Two different schemes based on the volatility basis set were used for the simulation of the bbOA during different seasons. The first is the default scheme of PMCAMx-SR and the second is a recently developed scheme based on laboratory experiments of the bbOA evolution. The simulations with the alternative bbOA scheme predict much higher total bbOA concentrations when compared with the base case ones. This is mainly due to the high emissions of intermediate volatility organic compounds (IVOCs) assumed in the alternative scheme. The oxidation of these compounds is predicted to be a significant source of secondary organic aerosol. The impact of the other parameters that differ in the two schemes is low to negligible. The monthly average maximum predicted concentrations of the alternative bbOA scheme were approximately an order of magnitude higher than those of the default scheme during all seasons. The performance of the two schemes was evaluated against observed total organic aerosol concentrations from several measurement sites across the US. The results were mixed. The default scheme performed better during July and September while the alternative scheme performed a little better during April. These results illustrate the uncertainty of the corresponding predictions, the need to quantify the emissions and reactions of IVOCs from specific biomass sources, and to better constrain the total (primary and secondary) bbOA levels.


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.


2009 ◽  
Vol 9 (3) ◽  
pp. 13693-13737 ◽  
Author(s):  
A. P. Tsimpidi ◽  
V. A. Karydis ◽  
M. Zavala ◽  
W. Lei ◽  
L. Molina ◽  
...  

Abstract. New primary and secondary organic aerosol modules have been added to PMCAMx, a three dimensional chemical transport model (CTM), for use with the SAPRC99 chemistry mechanism based on recent smog chamber studies. The new modeling framework is based on the volatility basis-set approach: both primary and secondary organic components are assumed to be semivolatile and photochemically reactive and are distributed in logarithmically spaced volatility bins. This new framework with the use of the new volatility basis parameters for low-NOx and high-NOx conditions tends to predict 4–6 times higher anthropogenic SOA concentrations than those predicted with older generation of models. The resulting PMCAMx-2008 was applied in Mexico City Metropolitan Area (MCMA) for approximately a week during April of 2003. The emission inventory, which uses as starting point the MCMA 2004 official inventory, is modified and the primary organic aerosol (POA) emissions are distributed by volatility based on dilution experiments. The predicted organic aerosol (OA) concentrations peak in the center of Mexico City reaching values above 40 μg m−3. The model predictions are compared with Aerosol Mass Spectrometry (AMS) observations and their Positive Matrix Factorization (PMF) analysis. The model reproduces both Hydrocarbon-like Organic Aerosol (HOA) and Oxygenated Organic Aerosol (OOA) concentrations and diurnal profiles. The small OA underprediction during the rush hour periods and overprediction in the afternoon suggest potential improvements to the description of fresh primary organic emissions and the formation of the oxygenated organic aerosols respectively, although they may also be due to errors in the simulation of dispersion and vertical mixing. However, the AMS OOA data are not specific enough to prove that the model reproduces the organic aerosol observations for the right reasons. Other combinations of contributions of primary, aged primary, and secondary organic aerosol production rates may lead to similar results. The model results suggest strongly that during the simulated period transport of OA from outside the city was a significant contributor to the observed OA levels. Future simulations should use a larger domain in order to test whether the regional OA can be predicted with current SOA parameterizations. Sensitivity tests indicate that the predicted OA concentration is especially sensitive to the volatility distribution of the emissions in the lower volatility bins.


Sign in / Sign up

Export Citation Format

Share Document