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

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

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

A Global Modeling Study on Carbonaceous Aerosol Microphysical Characteristics and Radiative Effects

Recently, attention has been drawn towards black carbon aerosols as a short-term climate warming mitigation candidate. However the global and regional impacts of the direct, indirect and semi-direct aerosol effects are highly uncertain, due to the complex nature of aerosol evolution and the way that mixed, aged aerosols interact with clouds and radiation. A detailed aerosol microphysical scheme, MATRIX, embedded within the GISS climate model is used in this study to present a quantitative assessment of the impact of microphysical processes involving black carbon, such as emission size distributions and optical properties on aerosol cloud activation and radiative effects. Our best estimate for net direct and indirect aerosol radiative flux change between 1750 and 2000 is -0.56 W/m2. However, the direct and indirect aerosol effects are quite sensitive to the black and organic carbon size distribution and consequential mixing state. The net radiative flux change can vary between -0.32 to -0.75 W/m2 depending on these carbonaceous particle properties at emission. Taking into account internally mixed black carbon particles let us simulate correct aerosol absorption. Absorption of black carbon aerosols is amplified by sulfate and nitrate coatings and, even more strongly, by organic coatings. Black carbon mitigation scenarios generally showed reduced radiative fluxeswhen sources with a large proportion of black carbon, such as diesel, are reduced; however reducing sources with a larger organic carbon component as well, such as bio-fuels, does not necessarily lead to a reduction in positive radiative flux

Can semi-volatile organic aerosols lead to fewer cloud particles?

The impact of condensing organic aerosols on activated cloud number concentration is examined in a new aerosol microphysics box model, MATRIX-VBS. The model includes the volatility basis set (VBS) framework coupled with the aerosol microphysical scheme MATRIX (Multiconfiguration Aerosol TRacker of mIXing state) that resolves aerosol mass and number concentrations and aerosol mixing state. By including the condensation of organic aerosols, the new model produces fewer activated particles compared to the original model, which treats organic aerosols as nonvolatile. Parameters such as aerosol chemical composition, mass and number concentrations, and particle sizes that affect activated cloud number concentration are thoroughly tested via a suite of Monte Carlo simulations. Results show that by considering semi-volatile organics in MATRIX-VBS, there is a lower activated particle number concentration, except in cases with low cloud updrafts, in clean environments at above-freezing temperatures, and in polluted environments at high temperatures (310&thinsp;K) and extremely low-humidity conditions.</p

Investigating hydroclimatic impacts of the 168–158&thinsp;BCE volcanic quartet and their relevance to the Nile River basin and Egyptian history

The Ptolemaic era (305–30 BCE) is an important period of Ancient Egyptian history known for its material and scientific advances, but also intermittent political and social unrest in the form of (sometimes widespread) revolts against the Ptolemaic elites. While the role of environmental pressures has long been overlooked in this period of Egyptian history, ice-core-based volcanic histories have identified the period as experiencing multiple notable eruptions, and a repeated temporal association between explosive volcanism and revolt has recently been noted. Here we analyze the global and regional (Nile River basin) hydroclimatic response to a unique historical sequence of four large and closely timed volcanic eruptions (first a tropical one, followed by three extratropical northern hemispheric events) between 168 and 158 BCE, a particularly troubled period in Ptolemaic history for which we now provide a more detailed hydroclimatic context. The NASA (National Aeronautics and Space Administration) GISS (Goddard Institute for Space Studies) ModelE2.1 Earth system model simulates a strong radiative response with a radiative forcing (top of atmosphere) of −7.5 W m−2 (following the first eruption) and −2.5 W m−2 (after each of the three remaining eruptions) at a global scale. Associated with this, we observe a global surface cooling of the order of 1.5 ∘C following the first (tropical) eruption, with the following three extratropical eruptions extending the cooling period for more than 15 years. Consequently, this series of eruptions is observed to constrain the northward migration of the inter-tropical convergence zone (ITCZ) during the Northern Hemisphere summer monsoon season, and major monsoon zones (African, South Asian, and East Asian) were seen to experience a suppression of rainfall of &gt;1 mm d−1 during the monsoon (JJAS) season averaged for 2 years after each eruption. A substantial suppression of the Indian and North African summer monsoon (over the Nile River headwater region) was seen to strongly affect the modeled river flow in the catchment and discharge at river mouth. River mass flow over the basin was observed to decrease by 29 % and 38 % relative to an unperturbed (non-volcanic) annual mean flow in the first and second year, respectively, after the first (i.e., tropical) eruption. A moderate decrease ranging between 5 % and 18 % was observed after the third and fourth (extratropical) eruptions. These results indicate, in sum, that the first eruption likely produced a strong hydroclimate response, with the following extratropical eruptions prolonging this. These results also support the recently hypothesized association between ice-core-based signals of explosive volcanism and hydroclimatic variability during the Ptolemaic era, including the suppression of the agriculturally critical Nile summer flooding.</p

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

Downscaling a Global Climate Model to Simulate Climate Change Impacts on U.S. Regional and Urban Air Quality

Climate change can exacerbate future regional air pollution events by making conditions more favorable to form high levels of ozone. In this study, we use spectral nudging with WRF to downscale NASA earth system GISS modelE2 results during the years 2006 to 2010 and 2048 to 2052 over the continental United States in order to compare the resulting meteorological fields from the air quality perspective during the four seasons of five-year historic and future climatological periods. GISS results are used as initial and boundary conditions by the WRF RCM to produce hourly meteorological fields. The downscaling technique and choice of physics parameterizations used are evaluated by comparing them with in situ observations. This study investigates changes of similar regional climate conditions down to a 12km by 12km resolution, as well as the effect of evolving climate conditions on the air quality at major U.S. cities. The high resolution simulations produce somewhat different results than the coarse resolution simulations in some regions. Also, through the analysis of the meteorological variables that most strongly influence air quality, we find consistent changes in regional climate that would enhance ozone levels in four regions of the U.S. during fall (Western U.S., Texas, Northeastern, and Southeastern U.S), one region during summer (Texas), and one region where changes potentially would lead to better air quality during spring (Northeast). We also find that daily peak temperatures tend to increase in most major cities in the U.S. which would increase the risk of health problems associated with heat stress. Future work will address a more comprehensive assessment of emissions and chemistry involved in the formation and removal of air pollutants