Global modelling of secondary organic aerosol in the troposphere: A sensitivity analysis

International audienceA global 3-dimensional chemistry/transport model able to describe O3, NOx, Volatile Organic Compounds (VOC), sulphur and NH3 chemistry has been extended to simulate the temporal and spatial distribution of primary and secondary carbonaceous aerosols in the troposphere focusing on Secondary Organic Aerosol (SOA) formation. A number of global simulations have been performed to determine a possible range of annual global SOA production and investigate uncertainties associated with the model results. Uncertainties in the model calculations related to the enthalpy of vaporization, the solubility and the activity coefficient of the condensable species, the chemical fate of the first generation low volatility oxidation products, the ageing of particles with regard to their hydrophilic properties, the partitioning of SOA on various aerosol surfaces and the evaporation of semi-volatiles from aerosol surfaces have been evaluated. Although not all sources of uncertainties have been investigated, according to our calculations, the above factors within the experimental range of variations could result to an overall uncertainty of about a factor of 20 in the global SOA budget. The global annual SOA production from biogenic VOC might range from 2.5 to 44.5 Tg of organic matter per year, whereas that from anthropogenic VOC ranges from 0.05 to 2.62 Tg of organic matter per year. These estimates can be considered as a lower limit, since partitioning on coarse particles like nitrate, dust or sea-salt, together with the partitioning and the dissociation of the semi-volatile products in aerosol water has been neglected. The large uncertainties associated with the emissions of VOC have not been investigated in this study. Comparison of model results to observations, where available, shows a better agreement for the upper budget estimates than for the lower ones

Evaluating Ammonium, Nitrate and Sulfate Aerosols in 3-Dimensions

The effect aerosols have on climate and air quality is a func-on of their chemical composi-on, concentra-on and spa-al distribu-on. These parameters are controlled by emissions, heterogeneous and homogeneous chemistry, where thermodynamics plays a key role, transport, which includes stratospheric- tropospheric exchange, and deposi-onal sinks. In this work we demonstrate the effect of some of these processes on the SO4NH4NO3 system using the GISS ModelE2 Global Circula-on Model (GCM)

Impact of Improvements in Volcanic Implementation on Atmospheric Chemistry and Climate in the GISS-E2 Model

The representation of volcanic eruptions in climate models introduces some of the largest errors when evaluating historical simulations, partly due to the crude model parameterizations. We will show preliminary results from the Goddard Institute for Space Studies (GISS)-E2 model comparing traditional highly parameterized volcanic implementation (specified Aerosol Optical Depth, Effective Radius) to deploying the full aerosol microphysics module MATRIX and directly emitting SO2 allowing us the prognosically determine the chemistry and climate impact. We show a reasonable match in aerosol optical depth, effective radius, and forcing between the full aerosol implementation and reconstructions/observations of the Mt. Pinatubo 1991 eruption, with a few areas as targets for future improvement. This allows us to investigate not only the climate impact of the injection of volcanic aerosols, but also influences on regional water vapor, O3, and OH distributions. With the skill of the MATRIX volcano implementation established, we explore (1) how the height of the injection column of SO2 influence atmospheric chemistry and climate response, (2) how the initial condition of the atmosphere influences the climate and chemistry impact of the eruption with a particular focus on how ENSO and QBO and (3) how the coupled chemistry could mitigate the climate signal for much larger eruptions (i.e. the 1258 eruption, reconstructed to be approximately 10x Pinatubo). During each sensitivity experiment we assess the impact on profiles of water vapor, O3, and OH, and assess how the eruption impacts the budget of each

Formation of secondary organic aerosol from isoprene oxidation over Europe

The role of isoprene as a precursor to secondary organic aerosol (SOA) over Europe is studied with the two-way nested global chemistry transport model TM5. The inclusion of the formation of SOA from isoprene oxidation in our model almost doubles the atmospheric burden of SOA over Europe compared to SOA formation from terpenes and aromatics. The reference simulation, which considers SOA formation from isoprene, terpenes and aromatics, predicts a yearly European production rate of 1.0 Tg SOA yr<sup>−1</sup> and an annual averaged atmospheric burden of about 50 Gg SOA over Europe. A fraction of 35% of the SOA produced in the boundary layer over Europe is transported to higher altitudes or to other world regions. Summertime measurements of organic matter (OM) during the extensive EMEP OC/EC campaign 2002/2003 are better reproduced when SOA formation from isoprene is taken into account, reflecting also the strong seasonality of isoprene and other biogenic volatile organic compounds (BVOC) emissions from vegetation. However, during winter, our model strongly underestimates OM, likely caused by missing wood burning in the emission inventories. Uncertainties in the parameterisation of isoprene SOA formation have been investigated. Maximum SOA production is found for irreversible sticking (non-equilibrium partitioning) of condensable vapours on particles, with tropospheric SOA production over Europe increased by a factor of 4 in summer compared to the reference case. Completely neglecting SOA formation from isoprene results in the lowest estimate (0.51 Tg SOA yr<sup>−1</sup>). The amount and the nature of the absorbing matter are shown to be another key uncertainty when predicting SOA levels. Consequently, smog chamber experiments on SOA formation should be performed with different types of seed aerosols and without seed aerosols in order to derive an improved treatment of the absorption of SOA in the models. Consideration of a number of recent insights in isoprene SOA formation mechanisms reduces the tropospheric production of isoprene derived SOA over Europe from 0.4 Tg yr<sup>−1</sup> in our reference simulation to 0.1 Tg yr<sup>−1</sup>

Formation of Secondary Organic Aerosol from Isoprene Oxidation over Europe

The role of isoprene as a precursor to secondary organic aerosol (SOA) over Europe is studied with the two-way nested global chemistry transport model TM5. The inclusion of the formation of SOA from isoprene oxidation in our model almost doubles the atmospheric burden of SOA over Europe compared to SOA formation from terpenes and aromatics. The reference simulation, which considers SOA formation from isoprene, terpenes and aromatics, predicts a yearly European production rate of 1.0 Tg SOA yr-1 and an annual averaged atmospheric burden of about 50 Gg SOA over Europe. A fraction of 35% of the SOA produced in the boundary layer over Europe is transported to higher altitudes or to other world regions. Summertime measurements of particulate organic matter (POM) during the extensive EMEP OC/EC campaign 2002/2003 are better reproduced when SOA formation from isoprene is taken into account, reflecting also the strong seasonality of isoprene and other biogenic volatile organic compounds (BVOC) emissions from vegetation. However, during winter, our model strongly underestimates POM, likely caused by missing wood burning in the emission inventories. Uncertainties in the parameterisation of isoprene SOA formation have been investigated. Maximum SOA production is found for irreversible sticking (non-equilibrium partitioning) of condensable vapours on particles, with tropospheric SOA production over Europe increased by a factor of 4 in summer compared to the reference case. Completely neglecting SOA formation from isoprene results in the lowest estimate (0.51 Tg SOA yr-1). The amount and the nature of the absorbing matter are shown to be another key uncertainty when predicting SOA levels. Tropospheric isoprene SOA production over Europe in summer more than doubles when, in addition to preexisting carbonaceous aerosols, condensation of semi volatile vapours on ammonium and sulphate aerosols is considered. Consequently, smog chamber experiments on SOA formation should be performed with different types of seed aerosols and without seed aerosols in order to derive an improved treatment of the absorption of SOA in the models. Consideration of a number of recent insights in isoprene SOA formation mechanisms reduces the tropospheric production of isoprene derived SOA over Europe from 0.4 Tg yr-1 in our reference simulation to 0.1 Tg yr-1JRC.DDG.H.2-Climate chang

Regional and global temperature response to anthropogenic SO2 emissions from China in three climate models

We use the HadGEM3-GA4, CESM1, and GISS ModelE2 climate models to investigate the global and regional aerosol burden, radiative flux, and surface temperature responses to removing anthropogenic sulfur dioxide (SO2) emissions from China. We find that the models differ by up to a factor of six in the simulated change in aerosol optical depth (AOD) and shortwave radiative flux over China that results from reduced sulfate aerosol, leading to a large range of magnitudes in the regional and global temperature responses. Two of the three models simulate a near-ubiquitous hemispheric warming due to the regional SO2 removal, with similarities in the local and remote pattern of response, but overall with a substantially different magnitude. The third model simulates almost no significant temperature response. We attribute the discrepancies in the response to a combination of substantial differences in the chemical conversion of SO2 to sulfate, translation of sulfate mass into AOD, cloud radiative interactions, and differences in the radiative forcing efficiency of sulfate aerosol in the models. The model with the strongest response (HadGEM3-GA4) compares best with observations of AOD regionally, however the other two models compare similarly (albeit poorly) and still disagree substantially in their simulated climate response, indicating th at total AOD observations are far from sufficient to determine which model response is more plausible. Our results highlight that there remains a large uncertainty in the representation of both aerosol chemistry as well as direct and indirect aerosol radiative effects in current climate models, and reinforces that caution must be applied when interpreting the results of modelling studies of aerosol influences on climate. Model studies that implicate aerosols in climate responses should ideally explore a range of radiative forcing strengths representative of this uncertainty, in addition to thoroughly evaluating the models used against observations

Reconciling Warming Trends

Climate models projected stronger warming over the past 15 years than has been seen in observations. Conspiring factors of errors in volcanic and solar inputs, representations of aerosols, and El NiNo evolution, may explain most of the discrepancy

MATRIX-VBS Condensing Organic Aerosols in an Aerosol Microphysics Model

The condensation of organic aerosols is represented in a newly developed box-model scheme, where its effect on the growth and composition of particles are examined. We implemented the volatility-basis set (VBS) framework into the aerosol mixing state resolving microphysical scheme Multiconfiguration Aerosol TRacker of mIXing state (MATRIX). This new scheme is unique and advances the representation of organic aerosols in models in that, contrary to the traditional treatment of organic aerosols as non-volatile in most climate models and in the original version of MATRIX, this new scheme treats them as semi-volatile. Such treatment is important because low-volatility organics contribute significantly to the growth of particles. The new scheme includes several classes of semi-volatile organic compounds from the VBS framework that can partition among aerosol populations in MATRIX, thus representing the growth of particles via condensation of low volatility organic vapors. Results from test cases representing Mexico City and a Finish forrest condistions show good representation of the time evolutions of concentration for VBS species in the gas phase and in the condensed particulate phase. Emitted semi-volatile primary organic aerosols evaporate almost completely in the high volatile range, and they condense more efficiently in the low volatility range

Naturally driven variability in the global secondary organic aerosol over a decade

International audienceIn order to investigate the variability of the secondary organic aerosol (SOA) distributions and budget and provide a measure for the robustness of the conclusions on human induced changes of SOA, a global 3-dimensional chemistry transport model describing both the gas and the particulate phase chemistry of the troposphere has been applied. The response of the global budget of SOA to temperature and moisture changes as well as to biogenic emission changes over a decade (1984-1993) has been evaluated. The considered emissions of biogenic non-methane volatile organic compounds (VOC) are driven by temperature, light and vegetation. They vary between 756 and 810 Tg Cy-1 and are therefore about 5.5 times higher than the anthropogenic VOC emissions. All secondary aerosols (sulphuric, nitrates and organics) are computed on-line together with the aerosol associated water. Over the studied decade, the computed natural variations (8%) in the chemical SOA production from biogenic VOC oxidation equal the chemical SOA production from anthropogenic VOC oxidation. Maximum values are calculated for 1990 (warmer and drier) and minimum values for 1986 (colder and wetter). The SOA computed variability results from a 7% increase in biogenic VOC emissions from 1986 to 1990 combined with 8.5% and 6% increases in the wet and dry deposition of SOA and leads to about 11.5% increase in the SOA burden of biogenic origin. The present study also demonstrates the importance of the hydrological cycle in determining the built up and fate of SOA in the atmosphere. It also reveals the existence of significant positive and negative feedback mechanisms in the atmosphere responsible for the non linear relationship between emissions of biogenic VOC and SOA burden

Past, present and future atmospheric nitrogen deposition

Reactive nitrogen emissions into the atmosphere are increasing due to human activities, affecting nitrogen deposition to the surface and impacting the productivity of terrestrial and marine ecosystems. An atmospheric chemistry-transport model (TM4-ECPL) is here used to calculate the global distribution of total nitrogen deposition, accounting for the first time for both its inorganic and organic fractions in gaseous and particulate phases, and past and projected changes due to anthropogenic activities. The anthropogenic and biomass burning ACCMIP historical and RCP6.0 and RCP8.5 emissions scenarios are used. Accounting for organic nitrogen (ON) primary emissions, the present-day global nitrogen atmospheric source is about 60% anthropogenic, while total N deposition increases by about 20% relative to simulations without ON primary emissions. About 20-25% of total deposited N is ON. About 10% of the emitted nitrogen oxides are deposited as ON instead of inorganic nitrogen (IN) as is considered in most global models. Almost a 3-fold increase over land (2-fold over the ocean) has been calculated for soluble N deposition due to human activities from 1850 to present. The investigated projections indicate significant changes in the regional distribution of N deposition and chemical composition, with reduced compounds gaining importance relative to oxidized ones, but very small changes in the global total flux. Sensitivity simulations quantify uncertainties due to the investigated model parameterizations of IN partitioning onto aerosols and of N chemically fixed on organics to be within 10% for the total soluble N deposition and between 25-35% for the dissolved ON deposition. Larger uncertainties are associated with N emissions