Evaluation of isoprene degradation in the detailed tropospheric chemical mechanism, MCM v3, using environmental chamber data

The isoprene degradation mechanism included in version 3 of the Master Chemical Mechanism (MCM v3) has been evaluated and refined, using the Statewide Air Pollution Research Center (SAPRC) environmental chamber datasets on the photo-oxidation of isoprene and its degradation products, methacrolein (MACR) and methylvinyl ketone (MVK). Prior to this, the MCM v3 butane degradation chemistry was also evaluated using chamber data on the photo-oxidation of butane, and its degradation products, methylethyl ketone (MEK), acetaldehyde (CH3CHO) and formaldehyde (HCHO), in conjunction with an initial evaluation of the chamber-dependent auxiliary mechanisms for the series of relevant chambers. The MCM v3 mechanisms for both isoprene and butane generally performed well and were found to provide an acceptable reaction framework for describing the NOx-photo-oxidation experiments on the above systems, although a number of parameter modifications and refinements were identified which resulted in an improved performance. All these relate to the magnitude of sources of free radicals from organic chemical process, such as carbonyl photolysis rates and the yields of radicals from the reactions of O3 with unsaturated oxygenates, and specific recommendations are made for refinements. In addition to this, it was necessary to include a representation of the reactions of O(3P) with isoprene, MACR and MVK (which were not previously treated in MCM v3), and conclusions are drawn concerning the required extent of free radical formation from these reactions. Throughout the study, the performance of MCM v3 was also compared with that of the SAPRC-99 mechanism, which was developed and optimized in conjunction with the chamber datasets.

Evaluation of alkene degradation in the detailed tropospheric chemistry mechanism, MCM v3, using environmental chamber data

The representation of alkene degradation in version 3 of the Master Chemical Mechanism (MCM v3) has been evaluated, using environmental chamber data on the photo-oxidation of ethene, propene, 1-butene and 1-hexene in the presence of NOx, from up to five chambers at the Statewide Air Pollution Research Center (SAPRC) at the University of California. As part of this evaluation, it was necessary to include a representation of the reactions of the alkenes with O(3P), which are significant under chamber conditions but generally insignificant under atmospheric conditions. The simulations for the ethene and propene systems, in particular, were found to be sensitive to the branching ratios assigned to molecular and free radical forming pathways of the O(3P) reactions, with the extent of radical formation required for proper fitting of the model to the chamber data being substantially lower than the reported consensus. With this constraint, the MCM v3 mechanisms for ethene and propene generally performed well. The sensitivity of the simulations to the parameters applied to a series of other radical sources and sink reactions (radical formation from the alkene ozonolysis reactions and product carbonyl photolysis; radical removal from the reaction of OH with NO2 and β-hydroxynitrate formation) were also considered, and the implications of these results are discussed. Evaluation of the MCM v3 1-butene and 1-hexene degradation mechanisms, using a more limited dataset from only one chamber, was found to be inconclusive. The results of sensitivity studies demonstrate that it is impossible to reconcile the simulated and observed formation of ozone in these systems for ranges of parameter values which can currently be justified on the basis of the literature. As a result of this work, gaps and uncertainties in the kinetic, mechanistic and chamber database are identified and discussed, in relation to both tropospheric chemistry and chemistry important under chamber conditions which may compromise the evaluation procedure, and recommendations are made for future experimental studies. Throughout the study, the performance of the MCM v3 chemistry was also simultaneously compared with that of the corresponding chemistry in the SAPRC-99 mechanism, which was developed and optimized in conjunction with the chamber datasets.

Modelling high particulate matter and ozone episodes in the Trans-Manche region

There is concern in the Cross-Channel region of Nord-Pas-de-Calais (France) and Kent (Great Britain), regarding the extent of atmospheric pollution detected in the area from emitted gaseous (VOC, NOx, S02)and particulate substances. In particular, the air quality of the Cross-Channel or "Trans-Manche" region is highly affected by the heavily industrial area of Dunkerque, in addition to transportation sources linked to cross-channel traffic in Kent and Calais, posing threats to the environment and human health. In the framework of the cross-border EU Interreg IIIA activity, the joint Anglo-French project, ATTMA, has been commissioned to study Aerosol Transport in the Trans-Manche Atmosphere. Using ground monitoring data from UK and French networks and with the assistance of satellite images the project aims to determine dispersion patterns. and identify sources responsible for the pollutants. The findings of this study will increase awareness and have a bearing on future air quality policy in the region. Public interest is evident by the presence of local authorities on both sides of the English Channel as collaborators. The research is based on pollution transport simulations using (a) Lagrangian Particle Dispersion (LPD) models, (b) an Eulerian Receptor Based model. This paper is concerned with part (a), the LPD Models. Lagrangian Particle Dispersion (LPD) models are often used to numerically simulate the dispersion of a passive tracer in the planetary boundary layer by calculating the Lagrangian trajectories of thousands of notional particles. In this contribution, the project investigated the use of two widely used particle dispersion models: the Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) model and the model FLEXPART. In both models forward tracking and inverse (or·. receptor-based) modes are possible. Certain distinct pollution episodes have been selected from the monitor database EXPER/PF and from UK monitoring stations, and their likely trajectory predicted using prevailing weather data. Global meteorological datasets were downloaded from the ECMWF MARS archive. Part of the difficulty in identifying pollution sources arises from the fact that much of the pollution outside the monitoring area. For example heightened particulate concentrations are to originate from sand storms in the Sahara, or volcanic activity in Iceland or the Caribbean work identifies such long range influences. The output of the simulations shows that there are notable differences between the formulations of and Hysplit, although both models used the same meteorological data and source input, suggesting that the identification of the primary emissions during air pollution episodes may be rather uncertain.

Biogenic nitric oxide emissions from soils: impact on NOx and ozone over West Africa during AMMA (African Monsoon Multidisciplinary Analysis) - Modelling study

Nitrogen oxide biogenic emissions from soils are driven by soil and environmental parameters. The relationship between these parameters and NO fluxes is highly non linear. A new algorithm, based on a neural network calculation, is used to reproduce the NO biogenic emissions linked to precipitations in the Sahel on the 6 August 2006 during the AMMA campaign. This algorithm has been coupled in the surface scheme of a coupled chemistry dynamics model (MesoNH Chemistry) to estimate the impact of the NO emissions on NOx and O3 formation in the lower troposphere for this particular episode. Four different simulations on the same domain and at the same period are compared: one with anthropogenic emissions only, one with soil NO emissions from a static inventory, at low time and space resolution, one with NO emissions from neural network, and one with NO from neural network plus lightning NOx. The influence of NOx from lightning is limited to the upper troposphere. The NO emission from soils calculated with neural network responds to changes in soil moisture giving enhanced emissions over the wetted soil, as observed by aircraft measurements after the passing of a convective system. The subsequent enhancement of NOx and ozone is limited to the lowest layers of the atmosphere in modelling, whereas measurements show higher concentrations above 1000 m. The neural network algorithm, applied in the Sahel region for one particular day of the wet season, allows an immediate response of fluxes to environmental parameters, unlike static emission inventories. Stewart et al (2008) is a companion paper to this one which looks at NOx and ozone concentrations in the boundary layer as measured on a research aircraft, examines how they vary with respect to the soil moisture, as indicated by surface temperature anomalies, and deduces NOx fluxes. In this current paper the model-derived results are compared to the observations and calculated fluxes presented by Stewart et al (2008).

Modelling study of the impact of deep convection on the utls air composition - Part I: Analysis of ozone precursors

The aim of this work is to study the local impact on the upper troposphere/lower stratosphere air composition of an extreme deep convective system. For this purpose, we performed a simulation of a convective cluster composed of many individual deep convective cells that occurred near Bauru (Brazil). The simulation is performed using the 3-D mesoscale model RAMS coupled on-line with a chemistry model. The comparisons with meteorological measurements show that the model produces meteorological fields generally consistent with the observations. The present paper (part I) is devoted to the analysis of the ozone precursors (CO, NO x and non-methane volatile organic compounds) and HO x in the UTLS. The simulation results show that the distribution of CO with altitude is closely related to the upward convective motions and consecutive outflow at the top of the convective cells leading to a bulge of CO between 7 km altitude and the tropopause (around 17km altitude). The model results for CO are consistent with satellite-borne measurements at 700 hPa. The simulation also indicates enhanced amounts of NO x up to 2 ppbv in the 7-17 km altitude layer mainly produced by the lightning associated with the intense convective activity. For insoluble non-methane volatile organic compounds, the convective activity tends to significantly increase their amount in the 7-17km layer by dynamical effects. During daytime in the presence of lightning NO x, this bulge is largely reduced in the upper part of the layer for reactive species (e.g. isoprene, ethene) because of their reactions with OH that is increased on average during daytime. Lightning NO x also impacts on the oxydizing capacity of the upper troposphere by reducing on average HO x, HO 2, H 2O 2 and organic hydroperoxides. During the simulation time, the impact of convection on the air composition of the lower stratosphere is negligible for all ozone precursors although several of the simulated convective cells nearly reach the tropopause. There is no significant transport from the upper troposphere to the lower stratosphere, the isentropic barrier not being crossed by convection. The impact of the increase of ozone precursors and HO x in the upper troposphere on the ozone budget in the LS is discussed in part II of this series of papers.

Comparison of different stomatal conductance algorithms for ozone flux modelling

A multiplicative and a semi-mechanistic, BWB-type [Ball, J.T., Woodrow, I.E., Berry, J.A., 1987. A model predicting stomatalconductance and its contribution to the control of photosynthesis under different environmental conditions. In: Biggens, J. (Ed.), Progress in Photosynthesis Research, vol. IV. Martinus Nijhoff, Dordrecht, pp. 221–224.] algorithm for calculating stomatalconductance (gs) at the leaf level have been parameterised for two crop and two tree species to test their use in regional scale ozone deposition modelling. The algorithms were tested against measured, site-specific data for durum wheat, grapevine, beech and birch of different European provenances. A direct comparison of both algorithms showed a similar performance in predicting hourly means and daily time-courses of gs, whereas the multiplicative algorithm outperformed the BWB-type algorithm in modelling seasonal time-courses due to the inclusion of a phenology function. The re-parameterisation of the algorithms for local conditions in order to validate ozone deposition modelling on a European scale reveals the higher input requirements of the BWB-type algorithm as compared to the multiplicative algorithm because of the need of the former to model net photosynthesis (An)

Implementation of a module for risk of ozone impacts assessment to vegetation in the integrated assessment modelling system for the Iberian peninsula. Evaluation for wheat and holm oak

A module to estimate risks of ozone damage to vegetation has been implemented in the Integrated Assessment Modelling system for the Iberian Peninsula. It was applied to compute three different indexes for wheat and Holm oak; daylight AOT40 (cumulative ozone concentration over 40 ppb), cumulative ozone exposure index according to the Directive 2008/50/EC (AOT40-D) and PODY (Phytotoxic Ozone Dose over a given threshold of Y nmol m−2 s−1). The use of these indexes led to remarkable differences in spatial patterns of relative ozone risks on vegetation. Ozone critical levels were exceeded in most of the modelling domain and soil moisture content was found to have a significant impact on the results. According to the outputs of the model, daylight AOT40 constitutes a more conservative index than the AOT40-D. Additionally, flux-based estimations indicate high risk areas in Portugal for both wheat and Holm oak that are not identified by AOT-based methods.

Modelling ozone stomatal flux over wheat under Mediterranean conditions

Correct estimation of leaf-level stomatal conductance (gsto) is central for current ozone (O3) risk assessment of wheat yield loss based on the absorbed O3 phytotoxic dose (POD). The gsto model parameterizations developed in Europe must be checked in the different climatic regions where they are going to be applied in order to reduce the uncertainties associated with the POD approach. This work proposes a new gsto model parameterization for estimating POD of Triticum aestivum and Triticum durum under Mediterranean conditions, based on phenological observations over 25 years and gsto field measurements during 5 growing seasons. Results show that POD in the Mediterranean area might be higher than previously estimated. However, caution must be paid when assessing the risk of yield loss for wheat in this area since field validation of O3 impacts is still limited.

Modelling photochemical production of fine particulates in a toluene/NOx/water vapour system

This work investigates the computer modelling of the photochemical formation of smog products such as ozone and aerosol, in a system containing toluene, NOx and water vapour. In particular, the problem of modelling this process in the Commonwealth Scientific and Industrial Research Organization (CSIRO) smog chambers, which utilize outdoor exposure, is addressed. The primary requirement for such modelling is a knowledge of the photolytic rate coefficients. Photolytic rate coefficients of species other than N02 are often related to JNo2 (rate coefficient for the photolysis ofN02) by a simple factor, but for outdoor chambers, this method is prone to error as the diurnal profiles may not be similar in shape. Three methods for the calculation of diurnal JNo2 are investigated. The most suitable method for incorporation into a general model, is found to be one which determines the photolytic rate coefficients for N02, as well as several other species, from actinic flux, absorption cross section and quantum yields. A computer model was developed, based on this method, to calculate in-chamber photolysis rate coefficients for the CSIRO smog chambers, in which ex-chamber rate coefficients are adjusted by accounting for variation in light intensity by transmittance through the Teflon walls, albedo from the chamber floor and radiation attenuation due to clouds. The photochemical formation of secondary aerosol is investigated in a series of toluene-NOx experiments, which were performed in the CSIRO smog chambers. Three stages of aerosol formation, in plots of total particulate volume versus time, are identified: a delay period in which no significant mass of aerosol is formed, a regime of rapid aerosol formation (regime 1) and a second regime of slowed aerosol formation (regime 2). Two models are presented which were developed from the experimental data. One model is empirically based on observations of discrete stages of aerosol formation and readily allows aerosol growth profiles to be calculated. The second model is based on an adaptation of published toluene photooxidation mechanisms and provides some chemical information about the oxidation products. Both models compare favorably against the experimental data. The gross effects of precursor concentrations (toluene, NOx and H20) and ambient conditions (temperature, photolysis rate) on the formation of secondary aerosol are also investigated, primarily using the mechanism model. An increase in [NOx]o results in increased delay time, rate of aerosol formation in regime 1 and volume of aerosol formed in regime 1. This is due to increased formation of dinitrocresol and furanone products. An increase in toluene results in a decrease in the delay time and an increase in the rate of aerosol formation in regime 1, due to enhanced reactivity from the toluene products, such as the radicals from the photolysis of benzaldehyde. Water vapor has very little effect on the formation of aerosol volume, except that rates are slightly increased due to more OH radicals from reaction with 0(1D) from ozone photolysis. Increased temperature results in increased volume of aerosol formed in regime 1 (increased dinitrocresol formation), while increased photolysis rate results in increased rate of aerosol formation in regime 1. Both the rate and volume of aerosol formed in regime 2 are increased by increased temperature or photolysis rate. Both models indicate that the yield of secondary particulates from hydrocarbons (mass concentration aerosol formed/mass concentration hydrocarbon precursor) is proportional to the ratio [NOx]0/[hydrocarbon]0

Regional and urban air quality modelling under a climate change scenario

The better understanding of the interactions between climate change and air quality is an emerging priority for research and policy. Climate change will bring changes in the climate system, which will affect the concentration and dispersion of air pollutants. The main objective of the current study is to assess the impacts of climate change on air quality in 2050 over Portugal and Porto urban area. First, an evaluation and characterization of the air quality over mainland Portugal was performed for the period between 2002 and 2012. The results show that NO2, PM10 and O3 are the critical pollutants in Portugal. Also, the influence of meteorology on O3, NO2 and PM10 levels was investigate in the national main urban areas (Porto and Lisboa) and was verified that O3 has a statistically significant relationship with temperature in most of the components. The results also indicate that emission control strategies are primary regulators for NO2 and PM10 levels. After, understanding the national air quality problems and the influence that meteorology had in the historical air quality levels, the air quality modelling system WRF-CAMx was tested and the required inputs for the simulations were prepared to fulfil the main goal of this work. For the required air quality modelling inputs, an Emission Projections under RCP scenarios (EmiPro-RCP) model was developed to assist the estimation of future emission inventories for GHG and common air pollutants. Also, the current emissions were estimated for Portugal with a higher detailed disaggregation to improve the performance of the air quality simulations. The air quality modelling system WRF/CAMx was tested and evaluated over Portugal and Porto urban area and the results point out that is an adequate tool for the analysis of air quality under climate change. For this purpose, regional simulations of air quality during historical period and future (2045-2050) were conducted with CAMx version 6.0 to evaluate the impacts of simulated future climate and anthropogenic emission projections on air quality over the study area. The climate and the emission projections were produced under the RCP8.5 scenario. The results from the simulations point out, that if the anthropogenic emissions keep the same in 2050, the concentrations of NO2, PM10 and O3 will increase in Portugal. When, besides the climate change effects, is consider the projected anthropogenic emissions the annual mean concentrations of NO2 decrease significantly in Portugal and Porto urban area, and on the contrary the annual mean PM10 concentrations increases in Portugal and decrease in Porto urban area. The O3 results are mainly caused by the reduction of ozone precursors, getting the higher reductions in urban areas and increases in the surrounding areas. All the analysis performed for both simulations for Porto urban area support that, for PM10 and O3, there will be an increase in the occurrence of extreme values, surpassing the annual legislated parameters and having more daily exceedances. This study constitutes an innovative scientific tool to help in future air quality management in order to mitigate future climate change impacts on air quality.

Application of the WRF-Chem Model for a High Ozone Episode in SW Poland – Preliminary Results

The Weather Research and Forecasting model, integrated online with chemistry module, is a multi-scale model suitable for both research and operational forecasts of meteorology and air quality. It is used by many institutions for a variety of applications. In this study, the WRF v3.5 with chemistry (WRF-Chem) is applied to the area of Poland, for a period of 3-20 July 2006, when high concentrations of ground level ozone were observed. The meteorological and chemistry simulations were initiated with ERA-Interim reanalysis and TNO MACC II emissions database, respectively. The model physical parameterization includes RRTM shortwave radiation, Kain-Fritsch cumulus scheme, Purdue Lin microphysics and ACM2 PBL, established previously as the optimal configuration. Chemical mechanism used for the study was RADM2 with MADE/SORGAM aerosols. Simulations were performed for three one-way nested domains covering Europe (36 km x 36 km), Central Europe (12 km x 12 km) and Poland (4 km x 4 km). The results from the innermost domain were analyzed and compared to measurements of ozone concentration at three stations in different environments. The results show underestimation of observed values and daily amplitude of ozone concentrations.

The potential to narrow uncertainty in projections of stratospheric ozone over the 21st century

Future stratospheric ozone concentrations will be determined both by changes in the concentration of ozone depleting substances (ODSs) and by changes in stratospheric and tropospheric climate, including those caused by changes in anthropogenic greenhouse gases (GHGs). Since future economic development pathways and resultant emissions of GHGs are uncertain, anthropogenic climate change could be a significant source of uncertainty for future projections of stratospheric ozone. In this pilot study, using an "ensemble of opportunity" of chemistry-climate model (CCM) simulations, the contribution of scenario uncertainty from different plausible emissions pathways for ODSs and GHGs to future ozone projections is quantified relative to the contribution from model uncertainty and internal variability of the chemistry-climate system. For both the global, annual mean ozone concentration and for ozone in specific geographical regions, differences between CCMs are the dominant source of uncertainty for the first two-thirds of the 21st century, up-to and after the time when ozone concentrations return to 1980 values. In the last third of the 21st century, dependent upon the set of greenhouse gas scenarios used, scenario uncertainty can be the dominant contributor. This result suggests that investment in chemistry-climate modelling is likely to continue to refine projections of stratospheric ozone and estimates of the return of stratospheric ozone concentrations to pre-1980 levels.

Coupling aerosol surface and bulk chemistry with a kinetic double layer model (K2-SUB): oxidation of oleic acid by ozone

We present a kinetic double layer model coupling aerosol surface and bulk chemistry (K2-SUB) based on the PRA framework of gas-particle interactions (Poschl-Rudich-Ammann, 2007). K2-SUB is applied to a popular model system of atmospheric heterogeneous chemistry: the interaction of ozone with oleic acid. We show that our modelling approach allows de-convoluting surface and bulk processes, which has been a controversial topic and remains an important challenge for the understanding and description of atmospheric aerosol transformation. In particular, we demonstrate how a detailed treatment of adsorption and reaction at the surface can be coupled to a description of bulk reaction and transport that is consistent with traditional resistor model formulations. From literature data we have derived a consistent set of kinetic parameters that characterise mass transport and chemical reaction of ozone at the surface and in the bulk of oleic acid droplets. Due to the wide range of rate coefficients reported from different experimental studies, the exact proportions between surface and bulk reaction rates remain uncertain. Nevertheless, the model results suggest an important role of chemical reaction in the bulk and an approximate upper limit of similar to 10(-11) cm(2) s(-1) for the surface reaction rate coefficient. Sensitivity studies show that the surface accommodation coefficient of the gas-phase reactant has a strong non-linear influence on both surface and bulk chemical reactions. We suggest that K2-SUB may be used to design, interpret and analyse future experiments for better discrimination between surface and bulk processes in the oleic acid-ozone system as well as in other heterogeneous reaction systems of atmospheric relevance.

Aircraft routing with minimal climate impact: the REACT4C climate cost function modelling approach (V1.0)

In addition to CO2, the climate impact of aviation is strongly influenced by non-CO2 emissions, such as nitrogen oxides, influencing ozone and methane, and water vapour, which can lead to the formation of persistent contrails in ice-supersaturated regions. Because these non-CO2 emission effects are characterised by a short lifetime, their climate impact largely depends on emission location and time; that is to say, emissions in certain locations (or times) can lead to a greater climate impact (even on the global average) than the same emission in other locations (or times). Avoiding these climate-sensitive regions might thus be beneficial to climate. Here, we describe a modelling chain for investigating this climate impact mitigation option. This modelling chain forms a multi-step modelling approach, starting with the simulation of the fate of emissions released at a certain location and time (time-region grid points). This is performed with the chemistry–climate model EMAC, extended via the two submodels AIRTRAC (V1.0) and CONTRAIL (V1.0), which describe the contribution of emissions to the composition of the atmosphere and to contrail formation, respectively. The impact of emissions from the large number of time-region grid points is efficiently calculated by applying a Lagrangian scheme. EMAC also includes the calculation of radiative impacts, which are, in a second step, the input to climate metric formulas describing the global climate impact of the emission at each time-region grid point. The result of the modelling chain comprises a four-dimensional data set in space and time, which we call climate cost functions and which describes the global climate impact of an emission at each grid point and each point in time. In a third step, these climate cost functions are used in an air traffic simulator (SAAM) coupled to an emission tool (AEM) to optimise aircraft trajectories for the North Atlantic region. Here, we describe the details of this new modelling approach and show some example results. A number of sensitivity analyses are performed to motivate the settings of individual parameters. A stepwise sanity check of the results of the modelling chain is undertaken to demonstrate the plausibility of the climate cost functions.

Local scale finite element modelling of stack pollutant emissions

Congresos y conferencias