Search

Reliable quantification of the sources and sinks of atmospheric carbon dioxide (CO2), including that of their trends and uncertainties, is essential to monitoring the progress in mitigating anthropogenic emissions under the Kyoto Protocol and the Paris Agreement. This study provides a consolidated synthesis of estimates for all anthropogenic and natural sources and sinks of CO2 for the European Union and UK (EU27 + UK), derived from a combination of state-of-the-art bottom-up (BU) and top-down (TD) data sources and models. Given the wide scope of the work and the variety of datasets involved, this study focuses on identifying essential questions which need to be answered to properly understand the differences between various datasets, in particular with regards to the less-well-characterized fluxes from managed ecosystems. The work integrates recent emission inventory data, process-based ecosystem model results, data-driven sector model results and inverse modeling estimates over the period 1990–2018. BU and TD products are compared with European national greenhouse gas inventories (NGHGIs) reported under the UNFCCC in 2019, aiming to assess and understand the differences between approaches. For the uncertainties in NGHGIs, we used the standard deviation obtained by varying parameters of inventory calculations, reported by the member states following the IPCC Guidelines. Variation in estimates produced with other methods, like atmospheric inversion models (TD) or spatially disaggregated inventory datasets (BU), arises from diverse sources including within-model uncertainty related to parameterization as well as structural differences between models. In comparing NGHGIs with other approaches, a key source of uncertainty is that related to different system boundaries and emission categories (CO2 fossil) and the use of different land use definitions for reporting emissions from land use, land use change and forestry (LULUCF) activities (CO2 land). At the EU27 + UK level, the NGHGI (2019) fossil CO2 emissions (including cement production) account for 2624 Tg CO2 in 2014 while all the other seven bottom-up sources are consistent with the NGHGIs and report a mean of 2588 (± 463 Tg CO2). The inversion reports 2700 Tg CO2 (± 480 Tg CO2), which is well in line with the national inventories. Over 2011–2015, the CO2 land sources and sinks from NGHGI estimates report −90 Tg C yr−1 ± 30 Tg C yr−1 while all other BU approaches report a mean sink of −98 Tg C yr−1 (± 362 Tg of C from dynamic global vegetation models only). For the TD model ensemble results, we observe a much larger spread for regional inversions (i.e., mean of 253 Tg C yr−1 ± 400 Tg C yr−1). This concludes that (a) current independent approaches are consistent with NGHGIs and (b) their uncertainty is too large to allow a verification because of model differences and probably also because of the definition of “CO2 flux” obtained from different approaches. The referenced datasets related to figures are visualized at https://doi.org/10.5281/zenodo.4626578 (Petrescu et al., 2020a).

In 2015 and 2016, atmospheric transport modeling challenges were conducted in the context of the Comprehensive Nuclear-Test-Ban Treaty (CTBT) verification, however, with a more limited scope with respect to emission inventories, simulation period and number of relevant samples (i.e., those above the Minimum Detectable Concentration (MDC)) involved. Therefore, a more comprehensive atmospheric transport modeling challenge was organized in 2019. Stack release data of Xe-133 were provided by the Institut National des Radioéléments/IRE (Belgium) and the Canadian Nuclear Laboratories/CNL (Canada) and accounted for in the simulations over a three (mandatory) or six (optional) months period. Best estimate emissions of additional facilities (radiopharmaceutical production and nuclear research facilities, commercial reactors or relevant research reactors) of the Northern Hemisphere were included as well. Model results were compared with observed atmospheric activity concentrations at four International Monitoring System (IMS) stations located in Europe and North America with overall considerable influence of IRE and/or CNL emissions for evaluation of the participants' runs. Participants were prompted to work with controlled and harmonized model set-ups to make runs more comparable, but also to increase diversity. It was found that using the stack emissions of IRE and CNL with daily resolution does not lead to better results than disaggregating annual emissions of these two facilities taken from the literature if an overall score for all stations covering all valid observed samples is considered. A moderate benefit of roughly 10% is visible in statistical scores for samples influenced by IRE and/or CNL to at least 50% and there can be considerable benefit for individual samples. Effects of transport errors, not properly characterized remaining emitters and long IMS sampling times (12-24 h) undoubtedly are in contrast to and reduce the benefit of high-quality IRE and CNL stack data. Complementary best estimates for remaining emitters push the scores up by 18% compared to just considering IRE and CNL emissions alone. Despite the efforts undertaken the full multi-model ensemble built is highly redundant. An ensemble based on a few arbitrary runs is sufficient to model the Xe-133 background at the stations investigated. The effective ensemble size is below five. An optimized ensemble at each station has on average slightly higher skill compared to the full ensemble. However, the improvement (maximum of 20% and minimum of 3% in RMSE) in skill is likely being too small for being exploited for an independent period.

To track progress towards keeping global warming well below 2 ∘C or even 1.5 ∘C, as agreed in the Paris Agreement, comprehensive up-to-date and reliable information on anthropogenic emissions and removals of greenhouse gas (GHG) emissions is required. Here we compile a new synthetic dataset on anthropogenic GHG emissions for 1970–2018 with a fast-track extension to 2019. Our dataset is global in coverage and includes CO2 emissions, CH4 emissions, N2O emissions, as well as those from fluorinated gases (F-gases: HFCs, PFCs, SF6, NF3) and provides country and sector details. We build this dataset from the version 6 release of the Emissions Database for Global Atmospheric Research (EDGAR v6) and three bookkeeping models for CO2 emissions from land use, land-use change, and forestry (LULUCF). We assess the uncertainties of global greenhouse gases at the 90 % confidence interval (5th–95th percentile range) by combining statistical analysis and comparisons of global emissions inventories and top-down atmospheric measurements with an expert judgement informed by the relevant scientific literature. We identify important data gaps for F-gas emissions. The agreement between our bottom-up inventory estimates and top-down atmospheric-based emissions estimates is relatively close for some F-gas species (∼ 10 % or less), but estimates can differ by an order of magnitude or more for others. Our aggregated F-gas estimate is about 10 % lower than top-down estimates in recent years. However, emissions from excluded F-gas species such as chlorofluorocarbons (CFCs) or hydrochlorofluorocarbons (HCFCs) are cumulatively larger than the sum of the reported species. Using global warming potential values with a 100-year time horizon from the Sixth Assessment Report by the Intergovernmental Panel on Climate Change (IPCC), global GHG emissions in 2018 amounted to 58 ± 6.1 GtCO2 eq. consisting of CO2 from fossil fuel combustion and industry (FFI) 38 ± 3.0 GtCO2, CO2-LULUCF 5.7 ± 4.0 GtCO2, CH4 10 ± 3.1 GtCO2 eq., N2O 2.6 ± 1.6 GtCO2 eq., and F-gases 1.3 ± 0.40 GtCO2 eq. Initial estimates suggest further growth of 1.3 GtCO2 eq. in GHG emissions to reach 59 ± 6.6 GtCO2 eq. by 2019. Our analysis of global trends in anthropogenic GHG emissions over the past 5 decades (1970–2018) highlights a pattern of varied but sustained emissions growth. There is high confidence that global anthropogenic GHG emissions have increased every decade, and emissions growth has been persistent across the different (groups of) gases. There is also high confidence that global anthropogenic GHG emissions levels were higher in 2009–2018 than in any previous decade and that GHG emissions levels grew throughout the most recent decade. While the average annual GHG emissions growth rate slowed between 2009 and 2018 (1.2 % yr−1) compared to 2000–2009 (2.4 % yr−1), the absolute increase in average annual GHG emissions by decade was never larger than between 2000–2009 and 2009–2018. Our analysis further reveals that there are no global sectors that show sustained reductions in GHG emissions. There are a number of countries that have reduced GHG emissions over the past decade, but these reductions are comparatively modest and outgrown by much larger emissions growth in some developing countries such as China, India, and Indonesia. There is a need to further develop independent, robust, and timely emissions estimates across all gases. As such, tracking progress in climate policy requires substantial investments in independent GHG emissions accounting and monitoring as well as in national and international statistical infrastructures. The data associated with this article (Minx et al., 2021) can be found at https://doi.org/10.5281/zenodo.5566761.

The Emissions Database for Global Atmospheric Research (EDGAR) compiles anthropogenic emissions data for greenhouse gases (GHGs), and for multiple air pollutants, based on international statistics and emission factors. EDGAR data provide quantitative support for atmospheric modelling and for mitigation scenario and impact assessment analyses as well as for policy evaluation. The new version (v4.3.2) of the EDGAR emission inventory provides global estimates, broken down to IPCC-relevant source-sector levels, from 1970 (the year of the European Union's first Air Quality Directive) to 2012 (the end year of the first commitment period of the Kyoto Protocol, KP). Strengths of EDGAR v4.3.2 include global geo-coverage (226 countries), continuity in time, and comprehensiveness in activities. Emissions of multiple chemical compounds, GHGs as well as air pollutants, from relevant sources (fossil fuel activities but also, for example, fermentation processes in agricultural activities) are compiled following a bottom-up (BU), transparent and IPCC-compliant methodology. This paper describes EDGAR v4.3.2 developments with respect to three major long-lived GHGs (CO2, CH4, and N2O) derived from a wide range of human activities apart from the land-use, land-use change and forestry (LULUCF) sector and apart from savannah burning; a companion paper quantifies and discusses emissions of air pollutants. Detailed information is included for each of the IPCC-relevant source sectors, leading to global totals for 2010 (in the middle of the first KP commitment period) (with a 95 % confidence interval in parentheses): 33.6(±5.9) Pg CO2 yr−1, 0.34(±0.16) Pg CH4 yr−1, and 7.2(±3.7) Tg N2O yr−1. We provide uncertainty factors in emissions data for the different GHGs and for three different groups of countries: OECD countries of 1990, countries with economies in transition in 1990, and the remaining countries in development (the UNFCCC non-Annex I parties). We document trends for the major emitting countries together with the European Union in more detail, demonstrating that effects of fuel markets and financial instability have had greater impacts on GHG trends than effects of income or population. These data (https://doi.org/10.5281/zenodo.2658138, Janssens-Maenhout et al., 2019) are visualised with annual and monthly global emissions grid maps of 0.1∘×0.1∘ for each source sector.

This paper describes, documents and validates the TM5-Fast Scenario Screening Tool (TM5-FASST), a global reduced-form air quality source-receptor model that has been designed to compute ambient pollutant concentrations as well as broad range of pollutant-related impacts, related to human health, agricultural crop production, and short-lived pollutant climate metrics, taking as input annual pollutant emission data aggregated at the national or regional level. The TM5-FASST tool, providing a trade-off between accuracy and applicability, is based on linearized emission-concentration sensitivities derived with the full chemistry-transport model TM5. The tool has been extensively applied in various recent critical studies. Although informal and fragmented validation has already been performed in various publications, this paper provides a comprehensive documentation of all components of the model and a validation against the full TM5 model. We find that the simplifications introduced in order to generate immediate results from emission scenarios are not compromising the validity of the output and as such TM5-FASST is proven to be a useful tool in science-policy analysis. Furthermore, it constitutes a suitable architecture for implementing the ensemble of source-receptor relations obtained in the frame of the HTAP modelling exercises, thus creating a link between the scientific community and policy-oriented users.

In this study we introduce a hybrid ensemble consisting of air quality models operating at both the global and regional scale. The work is motivated by the fact that these different types of models treat specific portions of the atmospheric spectrum with different levels of detail, and it is hypothesized that their combination can generate an ensemble that performs better than mono-scale ensembles. A detailed analysis of the hybrid ensemble is carried out in the attempt to investigate this hypothesis and determine the real benefit it produces compared to ensembles constructed from only global-scale or only regional-scale models. The study utilizes 13 regional and 7 global models participating in the Hemispheric Transport of Air Pollutants phase 2 (HTAP2)–Air Quality Model Evaluation International Initiative phase 3 (AQMEII3) activity and focuses on surface ozone concentrations over Europe for the year 2010. Observations from 405 monitoring rural stations are used for the evaluation of the ensemble performance. The analysis first compares the modelled and measured power spectra of all models and then assesses the properties of the mono-scale ensembles, particularly their level of redundancy, in order to inform the process of constructing the hybrid ensemble. This study has been conducted in the attempt to identify that the improvements obtained by the hybrid ensemble relative to the mono-scale ensembles can be attributed to its hybrid nature. The improvements are visible in a slight increase of the diversity (4 % for the hourly time series, 10 % for the daily maximum time series) and a smaller improvement of the accuracy compared to diversity. Root mean square error (RMSE) improved by 13–16 % compared to G and by 2–3 % compared to R. Probability of detection (POD) and false-alarm rate (FAR) show a remarkable improvement, with a steep increase in the largest POD values and smallest values of FAR across the concentration ranges. The results show that the optimal set is constructed from an equal number of global and regional models at only 15 % of the stations. This implies that for the majority of the cases the regional-scale set of models governs the ensemble. However given the high degree of redundancy that characterizes the regional-scale models, no further improvement could be expected in the ensemble performance by adding yet more regional models to it. Therefore the improvement obtained with the hybrid set can confidently be attributed to the different nature of the global models. The study strongly reaffirms the importance of an in-depth inspection of any ensemble of opportunity in order to extract the maximum amount of information and to have full control over the data used in the construction of the ensemble.

An accurate simulation of the absorption properties is key for assessing the radiative effects of aerosol on meteorology and climate. The representation of how chemical species are mixed inside the particles (the mixing state) is one of the major uncertainty factors in the assessment of these effects. Here we compare aerosol optical properties simulations over Europe and North America, coordinated in the framework of the third phase of the Air Quality Model Evaluation International Initiative (AQMEII), to 1 year of AERONET sunphotometer retrievals, in an attempt to identify a mixing state representation that better reproduces the observed single scattering albedo and its spectral variation. We use a single post-processing tool (FlexAOD) to derive aerosol optical properties from simulated aerosol speciation profiles, and focus on the absorption enhancement of black carbon when it is internally mixed with more scattering material, discarding from the analysis scenes dominated by dust. We found that the single scattering albedo at 440 nm (ω0,440) is on average overestimated (underestimated) by 3–5 % when external (core-shell internal) mixing of particles is assumed, a bias comparable in magnitude with the typical variability of the quantity. The (unphysical) homogeneous internal mixing assumption underestimates ω0,440 by ∼14 %. The combination of external and core-shell configurations (partial internal mixing), parameterized using a simplified function of air mass aging, reduces the ω0,440 bias to -1/-3 %. The black carbon absorption enhancement (Eabs) in core-shell with respect to the externally mixed state is in the range 1.8–2.5, which is above the currently most accepted upper limit of ∼1.5. The partial internal mixing reduces Eabs to values more consistent with this limit. However, the spectral dependence of the absorption is not well reproduced, and the absorption Ångström exponent AAE675440 is overestimated by 70–120 %. Further testing against more comprehensive campaign data, including a full characterization of the aerosol profile in terms of chemical speciation, mixing state, and related optical properties, would help in putting a better constraint on these calculations.

In this study, methods are proposed to diagnose the causes of errors in air quality (AQ) modelling systems. We investigate the deviation between modelled and observed time series of surface ozone through a revised formulation for breaking down the mean square error (MSE) into bias, variance, and the minimum achievable MSE (mMSE). The bias measures the accuracy and implies the existence of systematic errors and poor representation of data complexity, the variance measures the precision and provides an estimate of the variability of the modelling results in relation to the observed data, and the mMSE reflects unsystematic errors and provides a measure of the associativity between the modelled and the observed fields through the correlation coefficient. Each of the error components is analysed independently and apportioned to resolved process based on the corresponding timescale (long scale, synoptic, diurnal, and intra-day) and as a function of model complexity. The apportionment of the error is applied to the AQMEII (Air Quality Model Evaluation International Initiative) group of models, which embrace the majority of regional AQ modelling systems currently used in Europe and North America. The proposed technique has proven to be a compact estimator of the operational metrics commonly used for model evaluation (bias, variance, and correlation coefficient), and has the further benefit of apportioning the error to the originating timescale, thus allowing for a clearer diagnosis of the process that caused the error.

This study evaluates simulated vertical ozone profiles produced in the framework of the third phase of the Air Quality Model Evaluation International Initiative (AQMEII3) against ozonesonde observations in North America for the year 2010. Four research groups from the United States (US) and Europe have provided modeled ozone vertical profiles to conduct this analysis. Because some of the modeling systems differ in their meteorological drivers, wind speed and temperature are also included in the analysis. In addition to the seasonal ozone profile evaluation for 2010, we also analyze chemically inert tracers designed to track the influence of lateral boundary conditions on simulated ozone profiles within the modeling domain. Finally, cases of stratospheric ozone intrusions during May–June 2010 are investigated by analyzing ozonesonde measurements and the corresponding model simulations at Intercontinental Chemical Transport Experiment Ozonesonde Network Study (IONS) experiment sites in the western United States. The evaluation of the seasonal ozone profiles reveals that, at a majority of the stations, ozone mixing ratios are underestimated in the 1–6 km range. The seasonal change noted in the errors follows the one seen in the variance of ozone mixing ratios, with the majority of the models exhibiting less variability than the observations. The analysis of chemically inert tracers highlights the importance of lateral boundary conditions up to 250 hPa for the lower-tropospheric ozone mixing ratios (0–2 km). Finally, for the stratospheric intrusions, the models are generally able to reproduce the location and timing of most intrusions but underestimate the magnitude of the maximum mixing ratios in the 2–6 km range and overestimate ozone up to the first kilometer possibly due to marine air influences that are not accurately described by the models. The choice of meteorological driver appears to be a greater predictor of model skill in this altitude range than the choice of air quality model.

In the framework of the third phase of the Air Quality Model Evaluation International Initiative (AQMEII3), and as contribution to the second phase of the Hemispheric Transport of Air Pollution (HTAP2) activities for Europe and North America, the impacts of a 20 % decrease of global and regional anthropogenic emissions on surface air pollutant levels in 2010 are simulated by an international community of regional-scale air quality modeling groups, using different state-of-the-art chemistry and transport models (CTMs). The emission perturbations at the global level, as well as over the HTAP2-defined regions of Europe, North America and East Asia, are first simulated by the global Composition Integrated Forecasting System (C-IFS) model from European Centre for Medium-Range Weather Forecasts (ECMWF), which provides boundary conditions to the various regional CTMs participating in AQMEII3. On top of the perturbed boundary conditions, the regional CTMs used the same set of perturbed emissions within the regional domain for the different perturbation scenarios that introduce a 20 % reduction of anthropogenic emissions globally as well as over the HTAP2-defined regions of Europe, North America and East Asia. Results show that the largest impacts over both domains are simulated in response to the global emission perturbation, mainly due to the impact of domestic emission reductions. The responses of NO2, SO2 and PM concentrations to a 20 % anthropogenic emission reduction are almost linear (∼ 20 % decrease) within the global perturbation scenario with, however, large differences in the geographical distribution of the effect. NO2, CO and SO2 levels are strongly affected over the emission hot spots. O3 levels generally decrease in all scenarios by up to ∼ 1 % over Europe, with increases over the hot spot regions, in particular in the Benelux region, by an increase up to ∼ 6 % due to the reduced effect of NOx titration. O3 daily maximum of 8 h running average decreases in all scenarios over Europe, by up to ∼ 1 %. Over the North American domain, the central-to-eastern part and the western coast of the US experience the largest response to emission perturbations. Similar but slightly smaller responses are found when domestic emissions are reduced. The impact of intercontinental transport is relatively small over both domains, however, still noticeable particularly close to the boundaries. The impact is noticeable up to a few percent, for the western parts of the North American domain in response to the emission reductions over East Asia. O3 daily maximum of 8 h running average decreases in all scenarios over north Europe by up to ∼ 5 %. Much larger reductions are calculated over North America compared to Europe. In addition, values of the Response to Extra-Regional Emission Reductions (RERER) metric have been calculated in order to quantify the differences in the strengths of non-local source contributions to different species among the different models. We found large RERER values for O3 (∼ 0.8) over both Europe and North America, indicating a large contribution from non-local sources, while for other pollutants including particles, low RERER values reflect a predominant control by local sources. A distinct seasonal variation in the local vs. non-local contributions has been found for both O3 and PM2.5, particularly reflecting the springtime long-range transport to both continents.

In this study we present a novel approach for improving the air quality predictions using an ensemble of air quality models generated in the context of AQMEII (Air Quality Model Evaluation International Initiative). The development of the forecasting method makes use of modelled and observed time series (either spatially aggregated or relative to single monitoring stations) of ozone concentrations over different areas of Europe and North America. The technique considers the underlying forcing mechanisms on ozone by means of spectrally decomposed previsions. With the use of diverse applications, we demonstrate how the approach screens the ensemble members, extracts the best components and generates bias-free forecasts with improved accuracy over the candidate models. Compared to more traditional forecasting methods such as the ensemble median, the approach reduces the forecast error and at the same time it clearly improves the modelled variance. Furthermore, the result is not a mere statistical outcome depended on the quality of the selected members. The few individual cases with degraded performance are also identified and analysed. Finally, we show the extensions of the approach to other pollutants, specifically particulate matter and nitrogen dioxide, and provide a framework for its operational implementation. *One out of many

The ongoing trend of urbanisation worldwide is leading to a growing requirement for detailed flow and transport parameterisations to be included within numerical weather prediction (NWP) models. Such models often employ a simple roughness parameterisation for urban areas, which is not particularly accurate in predicting or assessing the flow and dispersion at street scale. Moreover, this kind of parameterisation offers too poor a representation of the mechanical and thermal forcing exerted by urban areas on the larger scale flow. At present, high computational costs and long simulation running times are among the constraints for the implementation of more detailed urban sub-models within NWP models. To overcome such limitations, a downscaling procedure from the atmospheric flow at the synoptic scale to the neighbourhood scale and below, is presented in this study. This is achieved by means of a simple urban model based on a parameterised formulation of the drag exerted by the building on the airflow. Application of the urban model for estimating spatially-averaged mean wind speed and the urban heat island over a selected neighbourhood area in Lisbon, Portugal, is presented. The results show the capability of the urban model to provide more accurate mean wind and temperature profiles. Moreover, the urban model has the advantage of being cost effective, as it requires small computational resources, and thus is suitable to be adopted in an operational context. The model is simple enough to be also used to assess how the resolving of urban surface processes may affect those at the larger scales.

Simulations from chemical weather models are subject to uncertainties in the input data (e.g. emission inventory, initial and boundary conditions) as well as those intrinsic to the model (e.g. physical parameterization, chemical mechanism). Multi-model ensembles can improve the forecast skill, provided that certain mathematical conditions are fulfilled. In this work, four ensemble methods were applied to two different datasets, and their performance was compared for ozone (O3), nitrogen dioxide (NO2) and particulate matter (PM10). Apart from the unconditional ensemble average, the approach behind the other three methods relies on adding optimum weights to members or constraining the ensemble to those members that meet certain conditions in time or frequency domain. The two different datasets were created for the first and second phase of the Air Quality Model Evaluation International Initiative (AQMEII). The methods are evaluated against ground level observations collected from the EMEP (European Monitoring and Evaluation Programme) and AirBase databases. The goal of the study is to quantify to what extent we can extract predictable signals from an ensemble with superior skill over the single models and the ensemble mean. Verification statistics show that the deterministic models simulate better O3 than NO2 and PM10, linked to different levels of complexity in the represented processes. The unconditional ensemble mean achieves higher skill compared to each station's best deterministic model at no more than 60 % of the sites, indicating a combination of members with unbalanced skill difference and error dependence for the rest. The promotion of the right amount of accuracy and diversity within the ensemble results in an average additional skill of up to 31 % compared to using the full ensemble in an unconditional way. The skill improvements were higher for O3 and lower for PM10, associated with the extent of potential changes in the joint distribution of accuracy and diversity in the ensembles. The skill enhancement was superior using the weighting scheme, but the training period required to acquire representative weights was longer compared to the sub-selecting schemes. Further development of the method is discussed in the conclusion.

Urban heat island (UHI) is one of the most well known forms of localised anthropogenic climate modification. It causes a local alteration of atmospheric stability. According to a top-down methodology, mesoscale meteorological modelling is a commonly used approach to study the impact of UHI on atmospheric stability at the urban scale. Our work is an effort to investigate UHI using a bottom-up approach by looking at the UHI through a computational fluid dynamics (CFD) model applied to the street canyons of a neighbourhood area. The CFD code is set up to model the thermal response (structure surface temperature and ambient air temperature) of an urban system to the outside climate. The determination of the air temperature in an urban unit allows the calculation of the ΔTu-r factor representing the difference between the air temperature in the urban system (u) and the air temperature recorded at the closest meteorological station (r), generally situated in the countryside. This factor, introduced by Oke in Boundary Layer Climate, (1987), enables the analysis of the heat island generated by an urban system. The simulation results obtained from the CFD model allows the estimation of the ΔTu-r factor in relation to physical aspects and geometrical configurations. We apply this technique to study UHI of a Mediterranean city of which some urban temperature measurements and morphometry from a digital elevation model (DEM) are available.

This paper deals with the modelling of the flow in the urban canopy layer. It critically reviews a well-known formula for the spatially-averaged wind profile, originally proposed by Cionco in 1965, and provides a new interpretation for it. This opens up a number of new applications for modelling mean wind flow over the neighbourhood scale. The model is based on a balance equation between the obstacle drag force and the local shear stress as proposed by Cionco for a vegetative canopy. The buildings within the canopy are represented as a canopy element drag formulated in terms of morphological parameters such as λf and λp (the ratios of plan area and frontal area of buildings to the lot area). These parameters can be obtained from the analysis of urban digital elevation models. The shear stress is parameterised using a mixing length approach. Spatially-averaged velocity profiles for different values of building packing density corresponding to different flow regimes are obtained and analysed. The computed solutions are compared with published data from windtunnel and water-tunnel experiments over arrays of cubes. The model is used to estimate the spatially-averaged velocity profile within and above neighbourhood areas of real cities by using vertical profiles of λf .

The simulation of the coupled evolution of atmospheric dynamics, pollutant transport, chemical reactions and atmospheric composition is one of the most challenging tasks in environmental modelling, climate change studies, and weather forecasting for the next decades as they all involve strongly integrated processes. Weather strongly influences air quality (AQ) and atmospheric transport of hazardous materials, while atmospheric composition can influence both weather and climate by directly modifying the atmospheric radiation budget or indirectly affecting cloud formation. Until recently, however, due to the scientific complexities and lack of computational power, atmospheric chemistry and weather forecasting have developed as separate disciplines, leading to the development of separate modelling systems that are only loosely coupled. The continuous increase in computer power has now reached a stage that enables us to perform online coupling of regional meteorological models with atmospheric chemical transport models. The focus on integrated systems is timely, since recent research has shown that meteorology and chemistry feedbacks are important in the context of many research areas and applications, including numerical weather prediction (NWP), AQ forecasting as well as climate and Earth system modelling. However, the relative importance of online integration and its priorities, requirements and levels of detail necessary for representing different processes and feedbacks can greatly vary for these related communities: (i) NWP, (ii) AQ forecasting and assessments, (iii) climate and earth system modelling. Additional applications are likely to benefit from online modelling, e.g.: simulation of volcanic ash or forest fire plumes, pollen warnings, dust storms, oil/gas fires, geo-engineering tests involving changes in the radiation balance. The COST Action ES1004 – European framework for online integrated air quality and meteorology modelling (EuMetChem) – aims at paving the way towards a new generation of online integrated atmospheric chemical transport and meteorology modelling with two-way interactions between different atmospheric processes including dynamics, chemistry, clouds, radiation, boundary layer and emissions. As its first task, we summarise the current status of European modelling practices and experience with online coupled modelling of meteorology with atmospheric chemistry including feedback mechanisms and attempt reviewing the various issues connected to the different modules of such online coupled models but also providing recommendations for coping with them for the benefit of the modelling community at large.

We explicitly address the fundamental issue of member diversity in multi-model ensembles. To date no attempts in this direction are documented within the air quality (AQ) community, although the extensive use of ensembles in this field. Common biases and redundancy are the two issues directly deriving from lack of independence, undermining the significance of a multi-model ensemble, and are the subject of this study. Shared biases among models will determine a biased ensemble, making therefore essential the errors of the ensemble members to be independent so that bias can cancel out. Redundancy derives from having too large a portion of common variance among the members of the ensemble, producing overconfidence in the predictions and underestimation of the uncertainty. The two issues of common biases and redundancy are analysed in detail using the AQMEII ensemble of AQ model results for four air pollutants in two European regions. We show that models share large portions of bias and variance, extending well beyond those induced by common inputs. We make use of several techniques to further show that subsets of models can explain the same amount of variance as the full ensemble with the advantage of being poorly correlated. Selecting the members for generating skilful, non-redundant ensembles from such subsets proved, however, non-trivial. We propose and discuss various methods of member selection and rate the ensemble performance they produce. In most cases, the full ensemble is outscored by the reduced ones. We conclude that, although independence of outputs may not always guarantee enhancement of scores (but this depends upon the skill being investigated) we discourage selecting the members of the ensemble simply on the basis of scores, that is, independence and skills need to be considered disjointly.

The impact of air pollution on human health and the associated external costs in Europe and the United States (US) for the year 2010 are modeled by a multi-model ensemble of regional models in the frame of the third phase of the Air Quality Modelling Evaluation International Initiative (AQMEII3). The modeled surface concentrations of O3, CO, SO2 and PM2.5 are used as input to the Economic Valuation of Air Pollution (EVA) system to calculate the resulting health impacts and the associated external costs from each individual model. Along with a base case simulation, additional runs were performed introducing 20 % anthropogenic emission reductions both globally and regionally in Europe, North America and east Asia, as defined by the second phase of the Task Force on Hemispheric Transport of Air Pollution (TF-HTAP2). Health impacts estimated by using concentration inputs from different chemistry–transport models (CTMs) to the EVA system can vary up to a factor of 3 in Europe (12 models) and the United States (3 models). In Europe, the multi-model mean total number of premature deaths (acute and chronic) is calculated to be 414 000, while in the US, it is estimated to be 160 000, in agreement with previous global and regional studies. The economic valuation of these health impacts is calculated to be EUR 300 billion and 145 billion in Europe and the US, respectively. A subset of models that produce the smallest error compared to the surface observations at each time step against an all-model mean ensemble results in increase of health impacts by up to 30 % in Europe, while in the US, the optimal ensemble mean led to a decrease in the calculated health impacts by ∼ 11 %. A total of 54 000 and 27 500 premature deaths can be avoided by a 20 % reduction of global anthropogenic emissions in Europe and the US, respectively. A 20 % reduction of North American anthropogenic emissions avoids a total of ∼ 1000 premature deaths in Europe and 25 000 total premature deaths in the US. A 20 % decrease of anthropogenic emissions within the European source region avoids a total of 47 000 premature deaths in Europe. Reducing the east Asian anthropogenic emissions by 20 % avoids ∼ 2000 total premature deaths in the US. These results show that the domestic anthropogenic emissions make the largest impacts on premature deaths on a continental scale, while foreign sources make a minor contribution to adverse impacts of air pollution.

Air quality modelling near airports has received attention due to the impact of emissions from aircrafts near ground level. This work is part of the model inter-comparison study undertaken for the Department of Transport in connection with air quality near Heathrow Airport. Results formed part of a submission to the UK government in July 2006. The Emissions and Dispersion Modelling System (EDMS) was used. Our contribution required the setting up and running of EDMS for Heathrow Airport and its surroundings to simulate the year 2002. NO(X), NO(2) and PM(10) were chosen for the study; these being of particular concern.

This paper deals with the modelling of the flow in the urban canopy layer. It critically reviews a well-known formula for the spatially-averaged wind profile, originally proposed by Cionco in 1965, and provides a new interpretation for it. This opens up a number of new applications for modelling mean wind flow over the neighbourhood scale. The model is based on a balance equation between the obstacle drag force and the local shear stress as proposed by Cionco for a vegetative canopy. The buildings within the canopy are represented as a canopy element drag formulated in terms of morphological parameters such as λ f and λ p (the ratios of plan area and frontal area of buildings to the lot area). These parameters can be obtained from the analysis of urban digital elevation models. The shear stress is parameterised using a mixing length approach. Spatially-averaged velocity profiles for different values of building packing density corresponding to different flow regimes are obtained and analysed. The computed solutions are compared with published data from wind-tunnel and water-tunnel experiments over arrays of cubes. The model is used to estimate the spatially-averaged velocity profile within and above neighbourhood areas of real cities by using vertical profiles of λ f .

Deep street canyons are a common feature of densely populated urban areas in Southern Europe. In these locations, multi-storey residential and/or office buildings surround relatively narrow streets reducing natural air ventilation. That may result in high concentrations of traffic-related pollutants, such as CO, in heavily congested streets and/or streets with high density of two-stroke motorcycles. Well calibrated and user-friendly dispersion models are needed for assessing air quality under these conditions. In the present study, two different modelling techniques are used to simulate CO concentrations at different heights within a deep street canyon (Via Nardones) in Naples, Italy.