Search

Biogenic volatile organic compounds (BVOCs) are emitted by vegetation and react with other compounds to form ozone and secondary organic matter (OM). In regional air quality models, biogenic emissions are often calculated using a plant functional type approach, which depends on the land use category. However, over cities, the land use is urban, so trees and their emissions are not represented. Here, we develop a bottom-up inventory of urban tree biogenic emissions in which the location of trees and their characteristics are derived from the tree database of the Paris city combined with allometric equations. Biogenic emissions are then computed for each tree based on their leaf dry biomass, tree-species-dependent emission factors, and activity factors representing the effects of light and temperature. Emissions are integrated in WRF-CHIMERE air quality simulations performed over June–July 2022. Over Paris city, the urban tree emissions have a significant impact on OM, inducing an average increase in the OM of about 5 %, reaching 14 % locally during the heatwaves. Ozone concentrations increase by 1.0 % on average and by 2.4 % during heatwaves, with a local increase of up to 6 %. The concentration increase remains spatially localized over Paris, extending to the Paris suburbs in the case of ozone during heatwaves. The inclusion of urban tree emissions improves the estimation of OM concentrations compared to in situ measurements, but they are still underestimated as trees are still missing from the inventory. OM concentrations are sensitive to terpene emissions, highlighting the importance of favoring urban tree species with low-terpene emissions.

Toluene represents a large fraction of anthropogenic emissions and significantly contributes to tropospheric ozone and secondary organic aerosol (SOA) formation. Despite the fact that toluene is one of the most studied aromatic compounds, detailed chemical mechanisms still fail to correctly reproduce the speciation of toluene gaseous and condensed oxidation products. This study aims to elucidate the role of initial experimental conditions in toluene SOA mass loadings and to investigate gas–particle partitioning of its reaction products at different relevant temperatures. Gaseous and particulate reaction products were identified and quantified using a proton transfer reaction time-of-flight mass spectrometer (PTR-ToF-MS) coupled to a CHemical Analysis of aeRosol ONline (CHARON) inlet. The chemical system exhibited a volatility distribution mostly in the semi-volatile regime. Temperature decrease caused a shift of saturation concentration towards lower values. The CHARON–PTR-ToF-MS instrument identified and quantified approximately 60 %–80 % of the total organic mass measured by an aerosol mass spectrometer. A detailed mechanism for toluene gaseous oxidation was developed based on the Master Chemical Mechanism (MCM) and Generator for Explicit Chemistry and Kinetics of Organics in the Atmosphere (GECKO-A) deterministic mechanisms, modified following the literature. The new mechanism showed improvements in modeling oxidation product speciation with more observed species represented and more representative concentrations compared to the MCM–GECKO-A reference. Tests on partitioning processes, nonideality, and wall losses highlighted the high dependency of SOA formation on the considered processes. Our results underline the fact that volatility is not sufficient to explain the gas–particle partitioning: the organic and the aqueous phases need to be considered as well as the interactions between compounds in the particle phase.

Air pollution in the central zone of Chile is not only a public health concern but also threatens water resources and climate, in connection with the transport and deposition of black carbon (BC) from urban centers onto the glaciers of the Andes. Chemistry-transport simulations reveal a seasonal dichotomy in the flux and latitudinal pattern of BC deposition on glaciers of the central Chilean Andes. The average deposition flux of BC on glaciers between 30 and 37∘ S is 4 times larger in winter, affecting mostly low-elevation glaciers, whereas the smaller summertime flux affects glaciers evenly, irrespective of their elevation. The contribution of emissions from the city of Santiago is dominant in summertime with more than 50 % along the Andes but minor in wintertime with less than 20 % even close to the capital city. Transport at larger scales and more local sources likely account for the remaining flux. The superimposition of synoptic-scale circulation and local mountain-valley circulation along the Andes drives the differences between summertime and wintertime deposition fluxes and generates a greater meteorological export potential during summer months. Future emissions and climate projections suggest that under the RCP8.5 scenario the gap between summertime and wintertime BC export and deposition flux could decrease, thereby pointing to summertime emission control gaining relevance. The chemistry-transport modeling approach for BC deposition on the Andes sheds light on the importance of the often disregarded summertime emissions on the radiative balance of its glaciers, particularly in the vicinity of Santiago.

With the growing evidence that high particle number concentrations may impact health, modelling their emissions and understanding formation processes is necessary, especially in cities where many people are exposed. As emission inventories of particle numbers and size distribution over cities are usually not available, a methodology is defined to estimate them from PM2.5 emissions and ratios of PM1 / PM2.5 and PM0.1 / PM2.5 by activity sector. In this methodology, a fitting parameter αem is used to redistribute the number concentrations in the lowest emission diameter range. This parameter is chosen by comparing measured and simulated number concentrations during non-nucleation days. The emission size distribution is then finely discretised by conserving both mass and number in each of the size ranges where emissions are specified. The methodology is applied over Greater Paris during the MEGAPOLI campaign (July 2009). Three-dimensional simulations are performed using the chemistry transport model Polair3D/Polyphemus coupled to the aerosol module SSH-aerosol to represent the evolution of particles by condensation, evaporation, coagulation, and nucleation, with a sectional approach for the size distribution. The model is first compared to measurements during non-nucleation days, and the influence over the month of July 2009 of three different nucleation parameterisations is assessed, i.e. binary (sulfuric acid, water), ternary (sulfuric acid, ammonia, water), and heteromolecular (extremely low-volatility organic compounds (ELVOCs) from monoterpenes and sulfuric acid). The modelled number concentrations compare very well to measurements, with an average normalised mean error of 42 % for the daily number concentrations of particles larger than 10 nm and 37 % for the number concentrations of particles larger than 100 nm. The influence of the binary nucleation is low, and the ternary nucleation scheme leads to better simulated number concentrations (in terms of bias and error) at only one site out of three, but it systematically reduces the model to measurement correlation, suggesting that ternary nucleation may not be the dominant process in new particle formation. However, the relative bias and error, as well as the correlation at suburban sites, are systematically improved using the heteromolecular nucleation scheme involving sulfuric acid and ELVOCs from monoterpenes. This suggests that heteromolecular nucleation may be important in cities, especially at suburban sites in summer, and that a better characterisation of the emissions of ELVOC precursors from traffic is needed.

The Eurodelta-Trends (EDT) multi-model experiment, aimed at assessing the efficiency of emission mitigation measures in improving air quality in Europe during 1990–2010, was designed to answer a series of questions regarding European pollution trends; i.e. were there significant trends detected by observations? Do the models manage to reproduce observed trends? How close is the agreement between the models and how large are the deviations from observations? In this paper, we address these issues with respect to particulate matter (PM) pollution. An in-depth trend analysis has been performed for PM10 and PM2.5 for the period of 2000–2010, based on results from six chemical transport models and observational data from the EMEP (Cooperative Programme for Monitoring and Evaluation of the Long-range Transmission of Air Pollutants in Europe) monitoring network. Given harmonization of set-up and main input data, the differences in model results should mainly result from differences in the process formulations within the models themselves, and the spread in the model-simulated trends could be regarded as an indicator for modelling uncertainty. The model ensemble simulations indicate overall decreasing trends in PM10 and PM2.5 from 2000 to 2010, with the total reductions of annual mean concentrations by between 2 and 5 (7 for PM10) µg m−3 (or between 10 % and 30 %) across most of Europe (by 0.5–2 µg m−3 in Fennoscandia, the north-west of Russia and eastern Europe) during the studied period. Compared to PM2.5, relative PM10 trends are weaker due to large inter-annual variability of natural coarse PM within the former. The changes in the concentrations of PM individual components are in general consistent with emission reductions. There is reasonable agreement in PM trends estimated by the individual models, with the inter-model variability below 30 %–40 % over most of Europe, increasing to 50 %–60 % in the northern and eastern parts of the EDT domain. Averaged over measurement sites (26 for PM10 and 13 for PM2.5), the mean ensemble-simulated trends are −0.24 and −0.22 µg m−3 yr−1 for PM10 and PM2.5, which are somewhat weaker than the observed trends of −0.35 and −0.40 µg m−3 yr−1 respectively, partly due to model underestimation of PM concentrations. The correspondence is better in relative PM10 and PM2.5 trends, which are −1.7 % yr−1 and −2.0 % yr−1 from the model ensemble and −2.1 % yr−1 and −2.9 % yr−1 from the observations respectively. The observations identify significant trends (at the 95 % confidence level) for PM10 at 56 % of the sites and for PM2.5 at 36 % of the sites, which is somewhat less that the fractions of significant modelled trends. Further, we find somewhat smaller spatial variability of modelled PM trends with respect to the observed ones across Europe and also within individual countries. The strongest decreasing PM trends and the largest number of sites with significant trends are found for the summer season, according to both the model ensemble and observations. The winter PM trends are very weak and mostly insignificant. Important reasons for that are the very modest reductions and even increases in the emissions of primary PM from residential heating in winter. It should be kept in mind that all findings regarding modelled versus observed PM trends are limited to the regions where the sites are located. The analysis reveals considerable variability of the role of the individual aerosols in PM10 trends across European countries. The multi-model simulations, supported by available observations, point to decreases in SO4-2 concentrations playing an overall dominant role. Also, we see relatively large contributions of the trends of NH4+ and NO3- to PM10 decreasing trends in Germany, Denmark, Poland and the Po Valley, while the reductions of primary PM emissions appear to be a dominant factor in bringing down PM10 in France, Norway, Portugal, Greece and parts of the UK and Russia. Further discussions are given with respect to emission uncertainties (including the implications of not accounting for forest fires and natural mineral dust by some of the models) and the effect of inter-annual meteorological variability on the trend analysis.

Agriculture is the main source of ammonia (NH3) in France, an important gaseous precursor of atmospheric particulate matter (PM). National and global emission inventories are known to have difficulty representing the large spatial and temporal variability inherent to atmospheric NH3. In this study, we compare NH3 emissions in France during spring 2011 from one reference inventory, the TNO inventory, and two alternative inventories that account in different manners for both the spatial and temporal variabilities of the emissions: (i) the NH3SAT satellite-derived inventory based on IASI NH3 columns and (ii) the CADASTRE-CIT inventory that combines NH3 emissions due to nitrogen fertilization calculated with the mechanistic model VOLT'AIR on the database of the CADASTRE_NH3 framework and other source emissions from the CITEPA. The total spring budgets, from March to May 2011, at the national level are higher when calculated with both alternative inventories than with the reference, the difference being more marked with CADASTRE-CIT. NH3SAT and CADASTRE-CIT inventories both yield to large NH3 spring emissions due to fertilization on soils with high pH in the northeastern part of France (65 and 135 kt NH3, respectively, vs. 48 kt NH3 for TNO-GEN), while soil properties are not accounted for by the TNO-GEN methodology. For the other parts of France, the differences are smaller. The timing of fertilization and associated ammonia emissions is closely related to the nitrogen requirements and hence the phenological stage of the crops, and therefore to the crop year's specific weather conditions. Maximum emissions are observed in March for 2011 for some regions for both alternative inventories, while April is the period with maximum emissions for the reference inventory regardless of the region or the year. Comparing the inventories at finer temporal resolutions, typically at daily scale, large differences are found. The convergence of alternative, independent and complementary methods on the spatiotemporal representation of the spring NH3 emissions, particularly over areas where the contribution of mineral fertilizer spreading to the spring budget is strong, encourages further developments in both prospective complementary directions, as this will help improve national NH3 emission inventories.

The VBS-GECKO (volatility basis set – Generator for Explicit Chemistry and Kinetics of Organics in the Atmosphere) parameterization for secondary organic aerosol (SOA) formation was integrated into the chemistry-transport model CHIMERE. Concentrations of organic aerosol (OA) and SOA were simulated over Europe for the July–August 2013 period. Simulated concentrations with VBS-GECKO were compared to results obtained with the former H2O parameterization implemented in CHIMERE and to observations from EMEP, ACTRIS and other observations available in the EBAS database. The model configuration using the VBS-GECKO parameterization slightly improves the performances compared to the model configuration using the former H2O parameterization. The VBS-GECKO model configuration performs well for stations showing a large SOA concentration from biogenic sources, especially in northern Europe, but underestimates OA concentrations over stations close to urban areas. Simulated OA was found to be mainly secondary (∼85 %) and from terpene oxidation. Simulations show negligible contribution of the oxidation of mono-aromatic compounds to SOA production. Tests performed to examine the sensitivity of simulated OA concentrations to hydro-solubility, volatility, aging rates and NOx regime have shown that the VBS-GECKO parameterization provides consistent results, with a weak sensitivity to changes in the parameters provided by the gas-phase mechanism included in CHIMERE (e.g., HOx or NOx concentrations). Different scenarios considering intermediate-volatility organic compound (IVOC) emissions were tested to examine the contribution of IVOC oxidation to SOA production. At the continental scale, these simulations show a weak sensitivity of OA concentrations to IVOC emission variations. At the local scale, accounting for IVOC emissions was found to lead to a substantial increase in OA concentrations in the plume from urban areas. This additional OA source remains too small to explain the gap between simulated and measured values at stations where anthropogenic sources are dominant.

This work aims at quantifying the relative contribution of secondary organic aerosol (SOA) precursors emitted by wildfires to organic aerosol (OA) formation during summer of 2007 over the Euro-Mediterranean region, where intense wildfires occurred. A new SOA formation mechanism, H2Oaro, including recently identified aromatic volatile organic compounds (VOCs) emitted from wildfires, is developed based on smog chamber experiment measurements under low- and high-NOx regimes. The aromatic VOCs included in the mechanism are toluene, xylene, benzene, phenol, cresol, catechol, furan, naphthalene, methylnaphthalene, syringol, guaiacol, and structurally assigned and unassigned compounds with at least six carbon atoms per molecule (USC>6). This mechanism H2Oaro is an extension of the H2O (hydrophilic–hydrophobic organic) aerosol mechanism: the oxidation of the precursor forms surrogate species with specific thermodynamic properties (volatility, oxidation degree and affinity to water). The SOA concentrations over the Euro-Mediterranean region in summer of 2007 are simulated using the chemistry transport model (CTM) Polair3D of the air-quality platform Polyphemus, where the mechanism H2Oaro was implemented. To estimate the relative contribution of the aromatic VOCs, intermediate volatility, semi-volatile and low-volatility organic compounds (I/S/L-VOCs), to wildfires OA concentrations, different estimations of the gaseous I/S/L-VOC emissions (from primary organic aerosol – POA – using a factor of 1.5 or from non-methanic organic gas – NMOG – using a factor of 0.36) and their ageing (one-step oxidation vs. multi-generational oxidation) are also tested in the CTM. Most of the particle OA concentrations are formed from I/S/L-VOCs. On average during the summer of 2007 and over the Euro-Mediterranean domain, they are about 10 times higher than the OA concentrations formed from VOCs. However, locally, the OA concentrations formed from VOCs can represent up to 30 % of the OA concentrations from biomass burning. Amongst the VOCs, the main contributors to SOA formation are phenol, benzene and catechol (CAT; 47 %); USC>6 compounds (23 %); and toluene and xylene (12 %). Sensitivity studies of the influence of the VOCs and the I/S/L-VOC emissions and chemical ageing mechanisms on PM2.5 concentrations show that surface PM2.5 concentrations are more sensitive to the parameterization used for gaseous I/S/L-VOC emissions than for ageing. Estimating the gaseous I/S/L-VOC emissions from POA or from NMOG has a high impact on local surface PM2.5 concentrations (reaching −30 % in the Balkans, −8 % to −16 % in the fire plume and +8 % to +16 % in Greece). Considering the VOC as SOA precursors results in a moderate increase in PM2.5 concentrations mainly in the Balkans (up to 24 %) and in the fire plume (+10 %).

The wet deposition of nitrogen and sulfur in Europe for the period 1990–2010 was estimated by six atmospheric chemistry transport models (CHIMERE, CMAQ, EMEP MSC-W, LOTOS-EUROS, MATCH and MINNI) within the framework of the EURODELTA-Trends model intercomparison. The simulated wet deposition and its trends for two 11-year periods (1990–2000 and 2000–2010) were evaluated using data from observations from the EMEP European monitoring network. For annual wet deposition of oxidised nitrogen (WNOx), model bias was within 30 % of the average of the observations for most models. There was a tendency for most models to underestimate annual wet deposition of reduced nitrogen (WNHx), although the model bias was within 40 % of the average of the observations. Model bias for WNHx was inversely correlated with model bias for atmospheric concentrations of NH3+NH4+, suggesting that an underestimation of wet deposition partially contributed to an overestimation of atmospheric concentrations. Model bias was also within about 40 % of the average of the observations for the annual wet deposition of sulfur (WSOx) for most models. Decreasing trends in WNOx were observed at most sites for both 11-year periods, with larger trends, on average, for the second period. The models also estimated predominantly decreasing trends at the monitoring sites and all but one of the models estimated larger trends, on average, for the second period. Decreasing trends were also observed at most sites for WNHx, although larger trends, on average, were observed for the first period. This pattern was not reproduced by the models, which estimated smaller decreasing trends, on average, than those observed or even small increasing trends. The largest observed trends were for WSOx, with decreasing trends at more than 80 % of the sites. On average, the observed trends were larger for the first period. All models were able to reproduce this pattern, although some models underestimated the trends (by up to a factor of 4) and others overestimated them (by up to 40 %), on average. These biases in modelled trends were directly related to the tendency of the models to under- or overestimate annual wet deposition and were smaller for the relative trends (expressed as % yr−1 relative to the deposition at the start of the period). The fact that model biases were fairly constant throughout the time series makes it possible to improve the predictions of wet deposition for future scenarios by adjusting the model estimates using a bias correction calculated from past observations. An analysis of the contributions of various factors to the modelled trends suggests that the predominantly decreasing trends in wet deposition are mostly due to reductions in emissions of the precursors NOx, NH3 and SOx. However, changes in meteorology (e.g. precipitation) and other (non-linear) interactions partially offset the decreasing trends due to emission reductions during the first period but not the second. This suggests that the emission reduction measures had a relatively larger effect on wet deposition during the second period, at least for the sites with observations.

In this study, assumptions (ideality and thermodynamic equilibrium) commonly made in three-dimensional (3-D) air quality models were reconsidered to evaluate their impacts on secondary organic aerosol (SOA) formation over Europe. To investigate the effects of non-ideality, dynamic mass transfer and aerosol viscosity on the SOA formation, the Secondary Organic Aerosol Processor (SOAP) model was implemented in the 3-D air quality model Polyphemus. This study presents the first 3-D modeling simulation which describes the impact of aerosol viscosity on the SOA formation. The model uses either the equilibrium approach or the dynamic approach with a method specially designed for 3-D air quality models to efficiently solve particle-phase diffusion when particles are viscous. Sensitivity simulations using two organic aerosol models implemented in Polyphemus to represent mass transfer between gas and particle phases show that the computation of the absorbing aerosol mass strongly influences the SOA formation. In particular, taking into account the concentrations of inorganic aerosols and hydrophilic organic aerosols in the absorbing mass of the aqueous phase increases the average SOA concentration by 5 % and 6 %, respectively. However, inorganic aerosols influence the SOA formation not only because they constitute an absorbing mass for hydrophilic SOA, but also because they interact with organic compounds. Non-ideality (short-, medium- and long-range interactions) was found to influence SOA concentrations by about 30 %. Concerning the dynamic mass transfer for the SOA formation, if the viscosity of SOA is not taken into account and if ideality of aerosols is assumed, the dynamic approach is found to give generally similar results to the equilibrium approach (indicating that equilibrium is an efficient hypothesis for inviscid and ideal aerosols). However, when a non-ideal aerosol is assumed, taking into account the dynamic mass transfer leads to a decrease of concentrations of the hydrophilic compounds (compared to equilibrium). This decrease is due to differences in the values of activity coefficients, which are different between values computed for bulk aerosols and those for each size section. This result indicates the importance of non-ideality on the dynamic evolution of SOA. For viscous aerosols, assuming a highly viscous organic phase leads to an increase in SOA concentrations during daytime (by preventing the evaporation of the most volatile organic compounds). The partitioning of nonvolatile compounds is not affected by viscosity, but the aging of more volatile compounds (that leads to the formation of the less volatile compounds) slows down as the evaporation of those compounds is stopped due to the viscosity of the particle. These results imply that aerosol concentrations may deviate significantly from equilibrium as the gas–particle partitioning could be higher than predicted by equilibrium. Furthermore, although a compound evaporates in the simulation using the equilibrium approach, the same compound can condense in the simulation using the dynamic approach if the particles are viscous. The results of this study emphasize the need for 3-D air quality models to take into account the effect of non-ideality on SOA formation and the effect of aerosol viscosity for the more volatile fraction of semi-volatile organic compounds.

New parameterizations for the formation of organic aerosols were developed. These parameterizations cover secondary organic aerosol (SOA) formation from biogenic and anthropogenic precursors, NOx dependency, oligomerization and the reactive uptake of pinonaldehyde. These parameterizations were implemented in a box model in which the condensation and/or evaporation of semi-volatile organic compounds was simulated by the Secondary Organic Aerosol Processor (SOAP) model to take the dynamic evolution of concentrations into account.The parameterizations were tested against several experiments carried out in previous studies in the EUPHORE outdoor chamber. Two datasets of experiments were used: the anthropogenic experiments (in which SOA is formed mainly from a mixture of toluene, 1,3,5-trimethylbenzene and o-xylene) and the biogenic experiments (in which SOA is formed mainly from α-pinene and limonene).When assuming no wall deposition of organic vapors, satisfactory results (bias lower than 20 %) were obtained for the biogenic experiments and for most of the anthropogenic experiments. However, a decrease of SOA concentrations (up to 30 %) was found when taking wall deposition of organic vapors into account (with the parameters of Zhang et al., 2014). The anthropogenic experiments seem to indicate a complex NOx dependency that could not be reproduced by the model. Oligomerization was found to have a strong effect on SOA composition (oligomers were estimated to account for up to 78 % of the SOA mass) and could therefore have a strong effect on the formation of SOA. The uptake of pinonaldehyde (which is a high-volatility semi-volatile organic compound, SVOC) onto acidic aerosol was found to be too slow to be significant under atmospheric conditions (no significant amount of SOA formed after 3 days of evolution), indicating that the parameterization of Pun and Seigneur (2007) used in some air quality models may lead to an overestimation of SOA concentrations. The uptake of aldehydes could nevertheless be an important SOA formation pathway for less volatile or more reactive aldehydes than pinonaldehyde.Regarding viscosity, a low effect of viscosity on SOA concentrations was estimated by the model, although a decrease of SVOC evaporation was found when taking it into account, as well as a lower sensitivity of concentrations to changes of temperature during the experiments.

Atmospheric chambers have been widely used to study secondary organic aerosol (SOA) properties and formation from various precursors under different controlled environmental conditions and to develop parameterization to represent SOA formation in chemical transport models (CTMs). Chamber experiments are however limited in number, performed under conditions that differ from the atmosphere and can be subject to potential artefacts from chamber walls. Here, the Generator for Explicit Chemistry and Kinetics of Organics in the Atmosphere (GECKO-A) modelling tool has been used in a box model under various environmental conditions to (i) explore the sensitivity of SOA formation and properties to changes on physical and chemical conditions and (ii) develop a volatility basis set (VBS)-type parameterization. The set of parent hydrocarbons includes n-alkanes and 1-alkenes with 10, 14, 18, 22 and 26 carbon atoms, α-pinene, β-pinene and limonene, benzene, toluene, o-xylene, m-xylene and p-xylene. Simulated SOA yields and their dependences on the precursor structure, organic aerosol load, temperature and NOx levels are consistent with the literature. GECKO-A was used to explore the distribution of molar mass, vaporization enthalpy, OH reaction rate and Henry's law coefficient of the millions of secondary organic compounds formed during the oxidation of the different precursors and under various conditions. From these explicit simulations, a VBS-GECKO parameterization designed to be implemented in 3-D air quality models has been tuned to represent SOA formation from the 18 precursors using GECKO-A as a reference. In evaluating the ability of VBS-GECKO to capture the temporal evolution of SOA mass, the mean relative error is less than 20 % compared to GECKO-A. The optimization procedure has been automated to facilitate the update of the VBS-GECKO on the basis of the future GECKO-A versions, its extension to other precursors and/or its modification to carry additional information.

The implementation of European emission abatement strategies has led to a significant reduction in the emissions of ozone precursors during the last decade. Ground-level ozone is also influenced by meteorological factors such as temperature, which exhibit interannual variability and are expected to change in the future. The impacts of climate change on air quality are usually investigated through air-quality models that simulate interactions between emissions, meteorology and chemistry. Within a multi-model assessment, this study aims to better understand how air-quality models represent the relationship between meteorological variables and surface ozone concentrations over Europe. A multiple linear regression (MLR) approach is applied to observed and modelled time series across 10 European regions in springtime and summertime for the period of 2000–2010 for both models and observations. Overall, the air-quality models are in better agreement with observations in summertime than in springtime and particularly in certain regions, such as France, central Europe or eastern Europe, where local meteorological variables show a strong influence on surface ozone concentrations. Larger discrepancies are found for the southern regions, such as the Balkans, the Iberian Peninsula and the Mediterranean basin, especially in springtime. We show that the air-quality models do not properly reproduce the sensitivity of surface ozone to some of the main meteorological drivers, such as maximum temperature, relative humidity and surface solar radiation. Specifically, all air-quality models show more limitations in capturing the strength of the ozone–relative-humidity relationship detected in the observed time series in most of the regions, for both seasons. Here, we speculate that dry-deposition schemes in the air-quality models might play an essential role in capturing this relationship. We further quantify the relationship between ozone and maximum temperature (mo3 − T, climate penalty) in observations and air-quality models. In summertime, most of the air-quality models are able to reproduce the observed climate penalty reasonably well in certain regions such as France, central Europe and northern Italy. However, larger discrepancies are found in springtime, where air-quality models tend to overestimate the magnitude of the observed climate penalty.

The evaluation and intercomparison of air quality models is key to reducing model errors and uncertainty. The projects AQMEII3 and EURODELTA-Trends, in the framework of the Task Force on Hemispheric Transport of Air Pollutants and the Task Force on Measurements and Modelling, respectively (both task forces under the UNECE Convention on the Long Range Transport of Air Pollution, LTRAP), have brought together various regional air quality models to analyze their performance in terms of air concentrations and wet deposition, as well as to address other specific objectives.This paper jointly examines the results from both project communities by intercomparing and evaluating the deposition estimates of reduced and oxidized nitrogen (N) and sulfur (S) in Europe simulated by 14 air quality model systems for the year 2010. An accurate estimate of deposition is key to an accurate simulation of atmospheric concentrations. In addition, deposition fluxes are increasingly being used to estimate ecological impacts. It is therefore important to know by how much model results differ and how well they agree with observed values, at least when comparison with observations is possible, such as in the case of wet deposition.This study reveals a large variability between the wet deposition estimates of the models, with some performing acceptably (according to previously defined criteria) and others underestimating wet deposition rates. For dry deposition, there are also considerable differences between the model estimates. An ensemble of the models with the best performance for N wet deposition was made and used to explore the implications of N deposition in the conservation of protected European habitats. Exceedances of empirical critical loads were calculated for the most common habitats at a resolution of 100  ×  100 m2 within the Natura 2000 network, and the habitats with the largest areas showing exceedances are determined.Moreover, simulations with reduced emissions in selected source areas indicated a fairly linear relationship between reductions in emissions and changes in the deposition rates of N and S. An approximate 20 % reduction in N and S deposition in Europe is found when emissions at a global scale are reduced by the same amount. European emissions are by far the main contributor to deposition in Europe, whereas the reduction in deposition due to a decrease in emissions in North America is very small and confined to the western part of the domain. Reductions in European emissions led to substantial decreases in the protected habitat areas with critical load exceedances (halving the exceeded area for certain habitats), whereas no change was found, on average, when reducing North American emissions in terms of average values per habitat.

The simulation of fine organic aerosols with CTMs (chemistry–transport models) in the western Mediterranean basin has not been studied until recently. The ChArMEx (the Chemistry-Aerosol Mediterranean Experiment) SOP 1b (Special Observation Period 1b) intensive field campaign in summer of 2013 gathered a large and comprehensive data set of observations, allowing the study of different aspects of the Mediterranean atmosphere including the formation of organic aerosols (OAs) in 3-D models. In this study, we used the CHIMERE CTM to perform simulations for the duration of the SAFMED (Secondary Aerosol Formation in the MEDiterranean) period (July to August 2013) of this campaign. In particular, we evaluated four schemes for the simulation of OA, including the CHIMERE standard scheme, the VBS (volatility basis set) standard scheme with two parameterizations including aging of biogenic secondary OA, and a modified version of the VBS scheme which includes fragmentation and formation of nonvolatile OA. The results from these four schemes are compared to observations at two stations in the western Mediterranean basin, located on Ersa, Cap Corse (Corsica, France), and at Cap Es Pinar (Mallorca, Spain). These observations include OA mass concentration, PMF (positive matrix factorization) results of different OA fractions, and 14C observations showing the fossil or nonfossil origins of carbonaceous particles. Because of the complex orography of the Ersa site, an original method for calculating an orographic representativeness error (ORE) has been developed. It is concluded that the modified VBS scheme is close to observations in all three aspects mentioned above; the standard VBS scheme without BSOA (biogenic secondary organic aerosol) aging also has a satisfactory performance in simulating the mass concentration of OA, but not for the source origin analysis comparisons. In addition, the OA sources over the western Mediterranean basin are explored. OA shows a major biogenic origin, especially at several hundred meters height from the surface; however over the Gulf of Genoa near the surface, the anthropogenic origin is of similar importance. A general assessment of other species was performed to evaluate the robustness of the simulations for this particular domain before evaluating OA simulation schemes. It is also shown that the Cap Corse site presents important orographic complexity, which makes comparison between model simulations and observations difficult. A method was designed to estimate an orographic representativeness error for species measured at Ersa and yields an uncertainty of between 50 and 85 % for primary pollutants, and around 2–10 % for secondary species.

A PM2.5-capable aerosol chemical speciation monitor (Q-ACSM) was deployed in urban Nanjing, China, for the first time to measure in situ non-refractory fine particle (NR-PM2.5) composition from 20 October to 19 November 2015, along with parallel measurements of submicron aerosol (PM1) species by a standard Q-ACSM. Our results show that the NR-PM2.5 species (organics, sulfate, nitrate, and ammonium) measured by the PM2.5-Q-ACSM are highly correlated (r2 > 0.9) with those measured by a Sunset Lab OC  /  EC analyzer and a Monitor for AeRosols and GAses (MARGA). The comparisons between the two Q-ACSMs illustrated similar temporal variations in all NR species between PM1 and PM2.5, yet substantial mass fractions of aerosol species were observed in the size range of 1–2.5 µm. On average, NR-PM1−2.5 contributed 53 % of the total NR-PM2.5, with sulfate and secondary organic aerosols (SOAs) being the two largest contributors (26 and 27 %, respectively). Positive matrix factorization of organic aerosol showed similar temporal variations in both primary and secondary OAs between PM1 and PM2.5, although the mass spectra were slightly different due to more thermal decomposition on the capture vaporizer of the PM2.5-Q-ACSM. We observed an enhancement of SOA under high relative humidity conditions, which is associated with simultaneous increases in aerosol pH, gas-phase species (NO2, SO2, and NH3) concentrations and aerosol water content driven by secondary inorganic aerosols. These results likely indicate an enhanced reactive uptake of SOA precursors upon aqueous particles. Therefore, reducing anthropogenic NOx, SO2, and NH3 emissions might not only reduce secondary inorganic aerosols but also the SOA burden during haze episodes in China.

In the framework of the Chemistry-Aerosol Mediterranean Experiment, a measurement site was set up at a remote site (Ersa) on Corsica Island in the northwestern Mediterranean Sea. Measurement campaigns performed during the summers of 2012 and 2013 showed high organic aerosol concentrations, mostly from biogenic origin. This work aims to represent the organic aerosol concentrations and properties (oxidation state and hydrophilicity) using the air-quality model Polyphemus with a surrogate approach for secondary organic aerosol (SOA) formation. Biogenic precursors are isoprene, monoterpenes and sesquiterpenes. In this work, the following model oxidation products of monoterpenes are added: (i) a carboxylic acid (MBTCA) to represent multi-generation oxidation products in the low-NOx regime, (ii) organic nitrate chemistry and (iii) extremely low-volatility organic compounds (ELVOCs) formed by ozonolysis. The model shows good agreement of measurements of organic concentrations for both 2012 and 2013 summer campaigns. The modelled oxidation property and hydrophilic organic carbon properties of the organic aerosols also agree reasonably well with the measurements. The influence of the different chemical processes added to the model on the oxidation level of organics is studied. Measured and simulated water-soluble organic carbon (WSOC) concentrations show that even at a remote site next to the sea, about 64 % of the organic carbon is soluble. The concentrations of WSOC vary with the origins of the air masses and the composition of organic aerosols. The marine organic emissions only contribute to a few percent of the organic mass in PM1, with maxima above the sea.

The EURODELTA III exercise has facilitated a comprehensive intercomparison and evaluation of chemistry transport model performances. Participating models performed calculations for four 1-month periods in different seasons in the years 2006 to 2009, allowing the influence of different meteorological conditions on model performances to be evaluated. The exercise was performed with strict requirements for the input data, with few exceptions. As a consequence, most of differences in the outputs will be attributed to the differences in model formulations of chemical and physical processes. The models were evaluated mainly for background rural stations in Europe. The performance was assessed in terms of bias, root mean square error and correlation with respect to the concentrations of air pollutants (NO2, O3, SO2, PM10 and PM2.5), as well as key meteorological variables. Though most of meteorological parameters were prescribed, some variables like the planetary boundary layer (PBL) height and the vertical diffusion coefficient were derived in the model preprocessors and can partly explain the spread in model results. In general, the daytime PBL height is underestimated by all models. The largest variability of predicted PBL is observed over the ocean and seas. For ozone, this study shows the importance of proper boundary conditions for accurate model calculations and then on the regime of the gas and particle chemistry. The models show similar and quite good performance for nitrogen dioxide, whereas they struggle to accurately reproduce measured sulfur dioxide concentrations (for which the agreement with observations is the poorest). In general, the models provide a close-to-observations map of particulate matter (PM2.5 and PM10) concentrations over Europe rather with correlations in the range 0.4–0.7 and a systematic underestimation reaching −10 µg m−3 for PM10. The highest concentrations are much more underestimated, particularly in wintertime. Further evaluation of the mean diurnal cycles of PM reveals a general model tendency to overestimate the effect of the PBL height rise on PM levels in the morning, while the intensity of afternoon chemistry leads formation of secondary species to be underestimated. This results in larger modelled PM diurnal variations than the observations for all seasons. The models tend to be too sensitive to the daily variation of the PBL. All in all, in most cases model performances are more influenced by the model setup than the season. The good representation of temporal evolution of wind speed is the most responsible for models' skillfulness in reproducing the daily variability of pollutant concentrations (e.g. the development of peak episodes), while the reconstruction of the PBL diurnal cycle seems to play a larger role in driving the corresponding pollutant diurnal cycle and hence determines the presence of systematic positive and negative biases detectable on daily basis.

The aerosol speciation and size distribution is modeled during the summer 2013 and over a large area encompassing Africa, Mediterranean and western Europe. The modeled aerosol is compared to available measurements such as the AERONET aerosol optical depth (AOD) and aerosol size distribution (ASD) and the EMEP network for surface concentrations of particulate matter PM2.5, PM10 and inorganic species (nitrate, sulfate and ammonium). The main goal of this study is to quantify the model ability to realistically model the speciation and size distribution of the aerosol. Results first showed that the long-range transport pathways are well reproduced and mainly constituted by mineral dust: spatial correlation is ≈ 0.9 for AOD and Ångström exponent, when temporal correlations show that the day-to-day variability is more difficult to reproduce. Over Europe, PM2.5 and PM10 have a mean temporal correlation of ≈ 0.4 but the lowest spatial correlation ( ≈ 0.25 and 0.62, respectively), showing that the fine particles are not well localized or transported. Being short-lived species, the uncertainties on meteorology and emissions induce these lowest scores. However, time series of PM2.5 with the speciation show a good agreement between model and measurements and are useful for discriminating the aerosol composition. Using a classification from the south (Africa) to the north (northern Europe), it is shown that mineral dust relative mass contribution decreases from 50 to 10 % when nitrate increases from 0 to 20 % and all other species, sulfate, sea salt, ammonium, elemental carbon, primary organic matter, are constant. The secondary organic aerosol contribution is between 10 and 20 % with a maximum at the latitude of the Mediterranean Sea (Spanish stations). For inorganic species, it is shown that nitrate, sulfate and ammonium have a mean temporal correlation of 0.25, 0.37 and 0.17, respectively. The spatial correlation is better (0.25, 0.5 and 0.87), showing that the mean values may be biased but the spatial localization of sulfate and ammonium is well reproduced. The size distribution is compared to the AERONET product and it is shown that the model fairly reproduces the main values for the fine and coarse mode. In particular, for the fine mode, the model overestimates the aerosol mass in Africa and underestimates it in Europe.

Beirut, the capital of Lebanon, which is located on the eastern shore of the Mediterranean basin, experiences high air pollution episodes. Annual average concentrations of coarse and fine particulate matter (PM2.5) as well as nitrogen oxides (NOx) often exceed the World Health Organization (WHO) guidelines. Therefore, improving air quality in this region is essential. The Polyphemus/Polair3D modeling system is used here to investigate air pollution episodes in Beirut during 2 to 18 July 2011. The modeling domain covers two nested grids of 1 and 5 km horizontal resolution over greater Beirut and Lebanon, respectively. The anthropogenic emission inventory was developed earlier (Waked et al., 2012). The Weather and Research Forecasting (WRF) model is used to generate the meteorological fields and the Model of Emissions of Gases and Aerosols from Nature (MEGAN) is used for biogenic emissions. The results of the study are compared to measurements from a field campaign conducted in the suburb of Beirut during 2–18 July 2011. The model reproduces satisfactorily the concentrations of most gaseous pollutants, the total mass of PM2.5 as well as PM2.5 elemental carbon (EC), organic carbon (OC), and sulfate. Ozone concentrations are overestimated and it appears that this overestimation results mainly from the boundary conditions.

A secondary organic aerosol (SOA) model, H2O (Hydrophilic/Hydrophobic Organic), is evaluated over the Paris area. This model treats the formation of SOA with two kinds of surrogate species: hydrophilic species (which condense preferentially on an aqueous phase) and hydrophobic species (which condense only on an organic phase). These surrogates species are formed from the oxidation in the atmosphere of volatile organic compounds (VOC) by radicals (HO and NO3) and ozone. These VOC are either biogenic (isoprene, monoterpenes and sesquiterpenes) or anthropogenic (mainly aromatic compounds). This model includes the formation of aerosols from different precursors (biogenic precursors, aromatics), and semi-volatile organic compounds (SVOC) from traffic. The H2O aerosol model was incorporated into the Polyphemus air quality modeling platform and applied to the Paris area and evaluated by comparison to measurements performed during the Megapoli campaign in July 2009. The comparison to measurements in the suburbs and in the city center of Paris shows that the model gives satisfactory results for both elemental carbon (EC) and organic carbon (OC). However, the model gives a peak of OC concentrations in the morning due to high emissions from traffic, which does not appear in measurements. Uncertainties in the modeled temperature, which can affect the gas-particle partitioning, in the partitioning of primary SVOC or underestimation of primary organic aerosol (POA) evaporation by the model could explain the differences between model and measurements. Moreover, using a theoretical mechanism for the oxidation of primary SVOC and intermediate volatility organic compounds (IVOC), POA concentrations were found to be likely overestimated by models due to the use of simple partitioning constants (which do not take into account the affinity of a compound with the liquid aerosol solution) or due to the assumption that the organic aerosol solution is a one-phase ideal solution. The organic aerosol in the city center of Paris was found to be originating mostly from distant sources with only 30 to 38% due to local sources.

A new model for the formation of secondary organic aerosol (SOA) from isoprene was developed. This model uses surrogate molecular species (hydroxy-hydroperoxides, tetrols, methylglyceric acid, organic nitrates) to represent SOA formation. The development of this model used available experimental data on yields and molecular composition of SOA from isoprene and methacrolein oxidation. This model reproduces the amount of particles measured in smog chambers under both low-NOx and high-NOx conditions. Under low-NOx conditions, the model reproduces the transitional formation of hydroxy-hydroperoxides particles, which are photolyzed and lead to SOA mass decrease after reaching a maximum. Under high-NOx conditions, particles are assumed to be formed mostly from the photo-oxidation of a PAN-type molecule derived from methacrolein (MPAN). This model successfully reproduces the complex NOx-dependence of isoprene oxidation and suggests a possible yield increase under some high-NOx conditions. Experimental data correspond to dry conditions (RH < 10%). However, particles formed from isoprene are expected to be highly hydrophilic, and isoprene oxidation products would likely partition between an aqueous phase and the gas phase at high humidity in the atmosphere. The model was extended to take into account the hydrophilic properties of SOA, which are relevant under atmospheric conditions, and investigate the effect of particulate liquid water on SOA formation. An important increase in SOA mass was estimated for humid conditions due to the hydrophilic properties. Experiments under high relative humidity conditions should be conducted to confirm the results of this study, which have implications for SOA modeling.