Evaluation and error apportionment of an ensemble of atmospheric chemistry transport modeling systems: multivariable temporal and spatial breakdown

What Impact Does Net Zero Action on Road Transport and Building Heating Have on Exposure to UK Air Pollution?

This study explores the cobenefits of reduced nitrogen dioxide (NO2), ozone (O3), and particulate matter (PM), through net zero (NZ) climate policy in the UK. Two alternative NZ scenarios, the balanced net zero (BNZP) and widespread innovation (WI) pathways, from the UK Climate Change Committee's Sixth Carbon Budget, were examined using a chemical transport model (CTM). Under the UK existing policy, Business as Usual (BAU), reductions in NO2 and PM were predicted by 2030 due to new vehicle technologies but plateau by 2040. The BNZP and WI scenarios show further reductions particularly by 2040, driven by accelerated electric vehicle (EV) uptake and low-carbon heating in buildings, with the building contribution to PM reduction being 2-3 times greater than road transport. The results demonstrate that the NZ transition to EVs (cars and vans) reduces both exhaust and nonexhaust emissions, as well as reducing traffic volumes. O3 trends are complex with a small overall increase by 2030 and a decrease by 2040. Although uncertain, 2050 predictions of BNZP showed important additional air pollution benefits. Our findings highlight the efficacy of NZ strategies, providing insights for UK and international policymakers interested in the air pollution cobenefits of climate policy.

Can the UK meet the World Health Organization PM2.5 interim target of 10 μg m-3 by 2030?

The recent United Kingdom (UK) Environment Act consultation had the intention of setting two targets for PM2.5 (particles with an aerodynamic diameter less than 2.5 μm), one related to meeting an annual average concentration and the second to reducing population exposure. As part of the consultation, predictions of PM2.5 concentrations in 2030 were made by combining European Union (EU) and UK government's emissions forecasts, with the Climate Change Committee's (CCC) Net Zero vehicle forecasts, and in London with the addition of local policies based on the London Environment Strategy (LES). Predictions in 2018 showed 6.4% of the UK's area and 82.6% of London's area had PM2.5 concentrations above the World Health Organization (WHO) interim target of 10 μg m-3, but by 2030, over 99% of the UK's area was predicted to be below it. However, kerbside concentrations in London and other major cities were still at risk of exceeding 10 μg m-3. With local action on PM2.5 in London, population weighted concentrations showed full compliance with the WHO interim target of 10 μg m-3 in 2030. However, predicting future PM2.5 concentrations and interpreting the results will always be difficult and uncertain for many reasons, such as imperfect models and the difficulty in estimating future emissions. To help understand the sensitivity of the model's PM2.5 predictions in 2030, current uncertainty was quantified using PM2.5 measurements and showed large areas in the UK that were still at risk of exceeding the WHO interim target despite the model predictions being below 10 μg m-3. Our results do however point to the benefits that policy at EU, UK and city level can have on achieving the WHO interim target of 10 μg m-3. These results were submitted to the UK Environment Act consultation. Nevertheless, the issues addressed here could be applicable to other European cities.

Linking gas, particulate, and toxic endpoints to air emissions in the Community Regional Atmospheric Chemistry Multiphase Mechanism (CRACMM)

Chemical mechanisms describe the atmospheric transformations of organic and inorganic species and connect air emissions to secondary species such as ozone, fine particles, and hazardous air pollutants (HAPs) like formaldehyde. Recent advances in our understanding of several chemical systems and shifts in the drivers of atmospheric chemistry warrant updates to mechanisms used in chemical transport models such as the Community Multiscale Air Quality (CMAQ) modeling system. This work builds on the Regional Atmospheric Chemistry Mechanism version 2 (RACM2) and develops the Community Regional Atmospheric Chemistry Multiphase Mechanism (CRACMM) version 1.0, which demonstrates a fully coupled representation of chemistry leading to ozone and secondary organic aerosol (SOA) with consideration of HAPs. CRACMMv1.0 includes 178 gas-phase species, 51 particulate species, and 508 reactions spanning gas-phase and heterogeneous pathways. To support estimation of health risks associated with HAPs, nine species in CRACMM cover 50 % of the total cancer and 60 % of the total non-cancer emission-weighted toxicity estimated for primary HAPs from anthropogenic and biomass burning sources in the US, with the coverage of toxicity higher (>80 %) when secondary formaldehyde and acrolein are considered. In addition, new mechanism species were added based on the importance of their emissions for the ozone, organic aerosol, or atmospheric burden of total reactive organic carbon (ROC): sesquiterpenes, furans, propylene glycol, alkane-like low- to intermediate-volatility organic compounds (9 species), low- to intermediate-volatility oxygenated species (16 species), intermediate-volatility aromatic hydrocarbons (2 species), and slowly reacting organic carbon. Intermediate- and lower-volatility organic compounds were estimated to increase the coverage of anthropogenic and biomass burning ROC emissions by 40 % compared to current operational mechanisms. Autoxidation, a gas-phase reaction particularly effective in producing SOA, was added for C10 and larger alkanes, aromatic hydrocarbons, sesquiterpenes, and monoterpene systems including second-generation aldehydes. Integrating the radical and SOA chemistry put additional constraints on both systems and enabled the implementation of previously unconsidered SOA pathways from phenolic and furanone compounds, which were predicted to account for ~ 30 % of total aromatic hydrocarbon SOA under typical atmospheric conditions. CRACMM organic aerosol species were found to span the atmospherically relevant range of species carbon number, number of oxygens per carbon, and oxidation state with a slight high bias in the number of hydrogens per carbon. In total, 11 new emitted species were implemented as precursors to SOA compared to current CMAQv5.3.3 representations, resulting in a bottom-up prediction of SOA, which is required for accurate source attribution and the design of control strategies. CRACMMv1.0 is available in CMAQv5.4.

Evaluation and Projection of Surface PM2.5 and Its Exposure on Population in Asia Based on the CMIP6 GCMs

This paper evaluates the historical simulated surface concentrations of particulate matter small than 2.5 µm in diameter (PM2.5) and its components (black carbon (BC), dust, SO4, and organic aerosol (OA)) in Asia, which come from Coupled Model Intercomparison Project Phase 6 (CMIP6). In addition, future projected changes of surface PM2.5 and its components, as well as their exposure to population, under the different Shared Socioeconomic Pathway (SSP) scenarios are also provided. Results show that the simulated spatial distribution of surface PM2.5 concentrations is consistent with the Modern-Era Retrospective Analysis for Research and Applications version 2 (MERRA-2) and Socioeconomic Data and Applications Center (SEDAC). The model spreads are small/large over the regions with low/high climatic mean surface PM2.5 concentrations, i.e., Northern Asia/Saudi Arabia, Iran, and Xinjiang Province of China. The multi-model ensemble of CMIP6 reproduces the main features of annual cycles and seasonal variations in Asia and its sub-regions. Under the scenarios of SSP1-2.6, SSP2-4.5, and SSP5-8.5, compared to the present-day period of 1995-2014, annual mean surface PM2.5 concentrations are projected to decrease in Asia, with obvious differences among the scenarios. Meanwhile, the magnitudes and timings of changes at the regional scale are quite different, with the largest decreases in South Asia (SAS). Under SSP3-7.0, the increase of surface PM2.5 concentrations in SAS is the largest, with the increase value of 8 μg/m3 in 2050; while under SSP370-lowNTCF, which assumes stronger levels of air quality control measures relative to the SSP3-7.0, the decreases of surface PM2.5 concentrations in SAS, East Asia (EAS) and Southeast Asia (SEAS) are the largest. The characteristics of seasonal trends are consistent with that of the annual trend. The trends in the concentrations of surface PM2.5 and its components are similar. The population-weighted average values of surface PM2.5 concentrations are projected to decrease in Central Asia (CAS), EAS, North Asia (NAS), and SEAS, and it indicates that the surface PM2.5 concentrations over the most populated area of Asia will decrease. In SAS, because of its large population, the impact of air pollutants on human health is still disastrous in the future. In summary, the surface PM2.5 concentrations over the most area of Asia will decrease, which is beneficial to air quality and human health; under SSP370-lowNTCF, the reduction of short-lived climate forcers (SLCFs) will further improve air quality.

Evaluating the performance of chemical transport models for PM2.5 source apportionment: An integrated application of spectral analysis and grey incidence analysis

Evaluating the performance of source apportionment (SA) models is difficult due to the non-observable nature of source contribution in reality. Here we propose a new approach to assess the performance of Chemical Transport Models (CTMs) for SA based on wavelet time-frequency spectral analysis and Grey Incidence Analysis (GIA). For each source category, certain species that better reflect the periodic characteristics of the emission sources were selected as the chemical tracers. The consistency of the time series between the simulated source contributions and the observed source-specific chemical tracers was then examined using a GIA model based on the perspective of similarity, and characterized by the GIA scores. By applying this method to six typical pollution episodes, we evaluated the performance of the Comprehensive Air Quality Model with Extensions-Particle Source Apportionment Technology (CAMx-PSAT) model for PM_2.5 SA from different temporal and spatial scales. The source- and episode-dependent optimal average time and main source regions were obtained. This approach is robust for facilitating a relatively meticulous evaluation of the performance of CTMs for PM_2.5 SA, and provides additional insight for decision-making for heavy pollution emergencies.

What Are the Sectors Contributing to the Exceedance of European Air Quality Standards over the Iberian Peninsula? A Source Contribution Analysis

The Iberian Peninsula, located in southwestern Europe, is exposed to frequent exceedances of different threshold and limit values of air pollution, mainly related to particulate matter, ozone, and nitrous oxide. Source apportionment modeling represents a useful modeling tool for evaluating the contribution of different emission sources or sectors and for designing useful mitigation strategies. In this sense, this work assesses the impact of various emission sectors on air pollution levels over the Iberian Peninsula using a source contribution analysis (zero-out method). The methodology includes the use of the regional WRF + CHIMERE modeling system (coupled to EMEP emissions). In order to represent the sensitivity of the chemistry and transport of gas-phase pollutants and aerosols, several emission sectors have been zeroed-out to quantify the influence of different sources in the area, such as on-road traffic or other mobile sources, combustion in energy generation, industrial emissions or agriculture, among others. The sensitivity analysis indicates that large reductions of precursor emissions (coming mainly from energy generation, road traffic, and maritime-harbor emissions) are needed for improving air quality and attaining the thresholds set in the European Directive 2008/50/EC over the Iberian Peninsula.

Establishing the Suitability of the Model for Prediction Across Scales for Global Retrospective Air Quality Modeling

The U.S. EPA is leveraging recent advances in meteorological modeling to construct an air quality modeling system to allow consistency from global to local scales. The Model for Prediction Across Scales-Atmosphere (MPAS-A or MPAS) has been developed by the National Center for Atmospheric Research (NCAR) as a global complement to the Weather Research and Forecasting model (WRF). Patterned after a regional coupled system with WRF, the Community Multiscale Air Quality (CMAQ) modeling system has been coupled within MPAS to explore global-to-local chemical transport modeling. Several options were implemented into MPAS for retrospective applications. Nudging-based data assimilation was added to support continuous simulations of past weather to minimize error growth that exists with a weather forecast configuration. The Pleim-Xiu land-surface model, the Asymmetric Convective Model 2 boundary layer scheme, and the Pleim surface layer scheme were added as the preferred options for retrospective air quality applications with WRF. Annual simulations were conducted using this EPA-enhanced MPAS configuration on two different mesh structures and compared against WRF. MPAS generally compares well with WRF over the conterminous United States. Errors in MPAS surface meteorology are comparable to WRF throughout the year. Precipitation statistics indicate MPAS performs slightly better than WRF. Solar radiation in MPAS is higher than WRF and measurements, suggesting fewer clouds in MPAS than WRF. Upper-air meteorology is well-simulated by MPAS, but errors are slightly higher than WRF. These comparisons lend confidence to use MPAS for retrospective air quality modeling and suggest ways it can be further improved in the future.

Assessing the manageable portion of ground-level ozone in the contiguous United States

Regional air quality models are widely being used to understand the spatial extent and magnitude of the ozone non-attainment problem and to design emission control strategies needed to comply with the relevant ozone standard through direct emission perturbations. In this study, we examine the manageable portion of ground-level ozone using two simulations of the Community Multiscale Air Quality (CMAQ) model for the year 2010 and a probabilistic analysis approach involving 29 years (1990-2018) of historical ozone observations. The modeling results reveal that the reduction in the peak ozone levels from total elimination of anthropogenic emissions within the model domain is around 13-21 ppb for the 90th-100th percentile range of the daily maximum 8-hr ozone concentrations across the contiguous United States (CONUS). Large reductions in the 4th highest 8-hr ozone are seen in the regions of West (interquartile range (IQR) of 17-33%), South (IQR 22-34%), Central (IQR 19-31%), Southeast (IQR 25-34%), and Northeast (IQR 24-37%). However, sites in the western portion of the domain generally show smaller reductions even when all anthropogenic emissions are removed, possibly due to the strong influence of global background ozone, including sources such as intercontinental ozone transport, stratospheric ozone intrusions, wildfires, and biogenic precursor emissions. Probabilistic estimates of the exceedances for several hypothetical thresholds of the 4th highest 8-hr ozone indicate that, in some areas, exceedances of such hypothetical thresholds may occur even with no anthropogenic emissions due to the ever-present atmospheric stochasticity and the current global tropospheric ozone burden. Implications: Because air pollution is intricately linked to adverse health effects, National Ambient Air Quality Standards (NAAQS) have been established for criteria pollutants to safeguard human health and the environment. Areas not in compliance with the relevant standards are required to develop plans and policies to reduce their air pollution levels. Regional-scale air quality models are currently being used routinely to inform policies to identify the emissions reduction required to meet and maintain the NAAQS throughout the country. This paper examines the feasibility of the 4th highest ozone, which is used to derive the ozone design value for NAAQS, complying with various current and hypothetical 8-hr ozone thresholds over CONUS based on the information embedded in 29 years of historical ozone observations and two modeling scenarios with and without anthropogenic emissions loading.

Impact of the ‘13th Five-Year Plan’ Policy on Air Quality in Pearl River Delta, China: A Case Study of Haizhu District in Guangzhou City Using WRF-Chem

Due to increasingly stringent control policy, air quality has generally improved in major cities in China during the past decade. However, the standards of national regulation and the World Health Organization are yet to be fulfilled in certain areas (in some urban districts among the cities) and/or certain periods (during pollution episode event). A further control policy, hence, has been issued in the 13th Five-Year Plan (2016–2020, hereafter 13th FYP). It will be of interest to evaluate the air quality before the 13th FYP (2015) and to estimate the potential air quality by the end of the 13th FYP (2020) with a focus on the area of an urban district and the periods of severe pollution episodes. Based on observation data of major air pollutants, including SO2 (sulphur dioxide), NO2 (nitrogen dioxide), CO (carbon monoxide), PM10 (particulate matter with aerodynamic diameter equal to or less than 10 μm), PM2.5 (particulate matter with aerodynamic diameter equal to or less than 2.5 µm) and O3 (Ozone), the air quality of Haizhu district [an urban district in the Pearl River Delta (PRD), China] in 2015 suggested that typical heavy pollution occurred in winter and the hot season, with NO2 or PM2.5 as the key pollutants in winter and O3 as the key pollutant in the hot season. We also adopted a state-of-the-art chemical transport model, the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem), to predict the air quality in Haizhu District 2020 under different scenarios. The simulation results suggested that among the emission control scenarios, comprehensive measures taken in the whole of Guangzhou city would improve air quality more significantly than measures taken just in Haizhu, under all conditions. In the urban district, vehicle emission control would account more than half of the influence of all source emission control on air quality. Based on our simulation, by the end of the 13th FYP, it is noticeable that O3 pollution would increase, which indicates that the control ratio of volatile organic compounds (VOCs) and nitrogen oxides (NOx) may be unsuitable and therefore should be adjusted. Our study highlights the significance of evaluating the efficacy of current policy in reducing the air pollutants and recommends possible directions for further air pollution control for urban areas during the 13th FYP.

High resolution temporal profiles in the Emissions Database for Global Atmospheric Research

Emissions into the atmosphere from human activities show marked temporal variations, from inter-annual to hourly levels. The consolidated practice of calculating yearly emissions follows the same temporal allocation of the underlying annual statistics. However, yearly emissions might not reflect heavy pollution episodes, seasonal trends, or any time-dependant atmospheric process. This study develops high-time resolution profiles for air pollutants and greenhouse gases co- emitted by anthropogenic sources in support of atmospheric modelling, Earth observation communities and decision makers. The key novelties of the Emissions Database for Global Atmospheric Research (EDGAR) temporal profiles are the development of (i) country/region- and sector- specific yearly profiles for all sources, (ii) time dependent yearly profiles for sources with inter-annual variability of their seasonal pattern, (iii) country- specific weekly and daily profiles to represent hourly emissions, (iv) a flexible system to compute hourly emissions including input from different users. This work creates a harmonized emission temporal distribution to be applied to any emission database as input for atmospheric models, thus promoting homogeneity in inter-comparison exercises.

Urban pollution in the Danube and Western Balkans regions: The impact of major PM2.5 sources

The SHERPA tool was used to assess the major pollution sources and the geographical areas impacting on the PM_2.5 of the main cities in the Danube and Western Balkans regions. The activity sectors influencing most the PM_2.5 levels in the study area are energy production (22%), agriculture (19%), residential combustion (16%) and road transport (7%). The energy production in inefficient coal-fuelled power plants was identified as one of main source of PM_2.5 in the Western Balkans. As for the geographical origin of PM_2.5, the transboundary pollution is confirmed as the main origin of PM_2.5 (44%) in the investigated cities, while the city own emissions and the national sources outside the concerned city impact on average 22% and 15%, respectively. An association was observed between the long-range transport and the impact of agriculture and energy production, while both local urban emissions and long-range transport were associated with the residential sector. A special attention is given in this study to biomass, a renewable source, which use is often promoted in the frame of climate and energy policies. Nevertheless, the combustion of biomass in inefficient small appliances has considerable particulate matter emissions and therefore this type of practice impacts negatively on air quality. Considering that biomass is traditionally used in South-East Europe as fuel for residential heating, the interpretation of the model results was supported with the estimation of biomass burning contributions to PM_2.5 obtained with receptor models and data on biomass fuel consumption from the literature. The analysis of the contributions from biomass burning derived from receptor models suggests that biomass burning is the dominant source within the residential heating sector in the studied area and that the emissions from this source are likely underestimated. This study concludes that more effort is needed to improve the estimations of biomass burning emissions and that policies to improve air quality in the cities should involve a geographic context wider than the city level.

Modelling black carbon absorption of solar radiation: combining external and internal mixing assumptions

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 (E _abs) 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 E _abs to values more consistent with this limit. However, the spectral dependence of the absorption is not well reproduced, and the absorption Ångström exponent AAE 675 440 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.

A New Method for Assessing the Efficacy of Emission Control Strategies

Regional-scale air quality models and observations at routine air quality monitoring sites are used to determine attainment/non-attainment of the ozone air quality standard in the United States. In current regulatory applications, a regional-scale air quality model is applied for a base year and a future year with reduced emissions using the same meteorological conditions as those in the base year. Because of the stochastic nature of the atmosphere, the same meteorological conditions would not prevail in the future year. Therefore, we use multi-decadal observations to develop a new method for estimating the confidence bounds for the future ozone design value (based on the 4^th highest value in the daily maximum 8-hr ozone concentration time series, DM8HR) for each emission loading scenario along with the probability of the design value exceeding a given ozone threshold concentration at all monitoring sites in the contiguous United States. To this end, we spectrally decompose the observed DM8HR ozone time series covering the period from 1981 to 2014 using the Kolmogorov-Zurbenko (KZ) filter and examine the variability in the relative strengths of the short-term variations (induced by synoptic-scale weather fluctuations; referred to as synoptic component, SY) and the long-term component (dictated by changes in emissions, seasonality and other slow-changing processes such as climate change; referred to as baseline component, BL). Results indicate that combining the projected change in the ozone baseline level with the adjusted synoptic forcing in historical ozone observations enables us to provide a probabilistic assessment of the efficacy of a selected emissions control strategy in complying with the ozone standard in future years. In addition, attainment demonstration is illustrated with a real-world application of the proposed methodology by using air quality model simulations, thereby helping build confidence in the use of regional-scale air quality models for supporting regulatory policies.

Attributing differences in the fate of lateral boundary ozone in AQMEII3 models to physical process representations

Increasing emphasis has been placed on characterizing the contributions and the uncertainties of ozone imported from outside the US. In chemical transport models (CTMs), the ozone transported through lateral boundaries (referred to as LB ozone hereafter) undergoes a series of physical and chemical processes in CTMs, which are important sources of the uncertainty in estimating the impact of LB ozone on ozone levels at the surface. By implementing inert tracers for LB ozone, the study seeks to better understand how differing representations of physical processes in regional CTMs may lead to differences in the simulated LB ozone that eventually reaches the surface across the US. For all the simulations in this study (including WRF/CMAQ, WRF/CAMx, COSMO-CLM/CMAQ, and WRF/DEHM), three chemically inert tracers that generally represent the altitude ranges of the planetary boundary layer (BC1), free troposphere (BC2), and upper troposphere-lower stratosphere (BC3) are tracked to assess the simulated impact of LB specification. Comparing WRF/CAMx with WRF/CMAQ, their differences in vertical grid structure explain 10 %-60 % of their seasonally averaged differences in inert tracers at the surface. Vertical turbulent mixing is the primary contributor to the remaining differences in inert tracers across the US in all seasons. Stronger vertical mixing in WRF/CAMx brings more BC2 downward, leading to higher BCT (BCT = BC1+BC2+BC3) and BC2/BCT at the surface in WRF/CAMx. Meanwhile, the differences in inert tracers due to vertical mixing are partially counteracted by their difference in sub-grid cloud mixing over the southeastern US and the Gulf Coast region during summer. The process of dry deposition adds extra gradients to the spatial distribution of the differences in DM8A BCT by 5-10 ppb during winter and summer. COSMO-CLM/CMAQ and WRF/CMAQ show similar performance in inert tracers both at the surface and aloft through most seasons, which suggests similarity between the two models at process level. The largest difference is found in summer. Sub-grid cloud mixing plays a primary role in their differences in inert tracers over the southeastern US and the oceans in summer. Our analysis of the vertical profiles of inert tracers also suggests that the model differences in dry deposition over certain regions are offset by the model differences in vertical turbulent mixing, leading to small differences in inert tracers at the surface in these regions.

Seasonal ozone vertical profiles over North America using the AQMEII3 group of air quality models: model inter-comparison and stratospheric intrusions

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.

Modeled deposition of nitrogen and sulfur in Europe estimated by 14 air quality model systems: evaluation, effects of changes in emissions and implications for habitat protection

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 m² 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.

Assessment and economic valuation of air pollution impacts on human health over Europe and the United States as calculated by a multi-model ensemble in the framework of AQMEII3

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.

Two-scale multi-model ensemble: is a hybrid ensemble of opportunity telling us more?

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.

Impacts of different characterizations of large-scale background on simulated regional-scale ozone over the continental United States

This study analyzes simulated regional-scale ozone burdens both near the surface and aloft, estimates process contributions to these burdens, and calculates the sensitivity of the simulated regional-scale ozone burden to several key model inputs with a particular emphasis on boundary conditions derived from hemispheric or global-scale models. The Community Multiscale Air Quality (CMAQ) model simulations supporting this analysis were performed over the continental US for the year 2010 within the context of the Air Quality Model Evaluation International Initiative (AQMEII) and Task Force on Hemispheric Transport of Air Pollution (TF-HTAP) activities. CMAQ process analysis (PA) results highlight the dominant role of horizontal and vertical advection on the ozone burden in the mid-to-upper troposphere and lower stratosphere. Vertical mixing, including mixing by convective clouds, couples fluctuations in free-tropospheric ozone to ozone in lower layers. Hypothetical bounding scenarios were performed to quantify the effects of emissions, boundary conditions, and ozone dry deposition on the simulated ozone burden. Analysis of these simulations confirms that the characterization of ozone outside the regional-scale modeling domain can have a profound impact on simulated regional-scale ozone. This was further investigated by using data from four hemispheric or global modeling systems (Chemistry - Integrated Forecasting Model (C-IFS), CMAQ extended for hemispheric applications (H-CMAQ), the Goddard Earth Observing System model coupled to chemistry (GEOS-Chem), and AM3) to derive alternate boundary conditions for the regional-scale CMAQ simulations. The regional-scale CMAQ simulations using these four different boundary conditions showed that the largest ozone abundance in the upper layers was simulated when using boundary conditions from GEOS-Chem, followed by the simulations using C-IFS, AM3, and H-CMAQ boundary conditions, consistent with the analysis of the ozone fields from the global models along the CMAQ boundaries. Using boundary conditions from AM3 yielded higher springtime ozone columns burdens in the middle and lower troposphere compared to boundary conditions from the other models. For surface ozone, the differences between the AM3-driven CMAQ simulations and the CMAQ simulations driven by other large-scale models are especially pronounced during spring and winter where they can reach more than 10 ppb for seasonal mean ozone mixing ratios and as much as 15 ppb for domain-averaged daily maximum 8 h average ozone on individual days. In contrast, the differences between the C-IFS-, GEOS-Chem-, and H-CMAQ-driven regional-scale CMAQ simulations are typically smaller. Comparing simulated sur face ozone mixing ratios to observations and computing seasonal and regional model performance statistics revealed that boundary conditions can have a substantial impact on model performance. Further analysis showed that boundary conditions can affect model performance across the entire range of the observed distribution, although the impacts tend to be lower during summer and for the very highest observed percentiles. The results are discussed in the context of future model development and analysis opportunities.

Influence of anthropogenic emissions and boundary conditions on multi-model simulations of major air pollutants over Europe and North America in the framework of AQMEII3

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 nonlocal 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.

Source apportionment of PM2.5 in North India using source-oriented air quality models

In recent years, severe pollution events were observed frequently in India especially at its capital, New Delhi. However, limited studies have been conducted to understand the sources to high pollutant concentrations for designing effective control strategies. In this work, source-oriented versions of the Community Multi-scale Air Quality (CMAQ) model with Emissions Database for Global Atmospheric Research (EDGAR) were applied to quantify the contributions of eight source types (energy, industry, residential, on-road, off-road, agriculture, open burning and dust) to fine particulate matter (PM_2.5) and its components including primary PM (PPM) and secondary inorganic aerosol (SIA) i.e. sulfate, nitrate and ammonium ions, in Delhi and three surrounding cities, Chandigarh, Lucknow and Jaipur in 2015. PPM mass is dominated by industry and residential activities (>60%). Energy (∼39%) and industry (∼45%) sectors contribute significantly to PPM at south of Delhi, which reach a maximum of 200 μg/m³ during winter. Unlike PPM, SIA concentrations from different sources are more heterogeneous. High SIA concentrations (∼25 μg/m³) at south Delhi and central Uttar Pradesh were mainly attributed to energy, industry and residential sectors. Agriculture is more important for SIA than PPM and contributions of on-road and open burning to SIA are also higher than to PPM. Residential sector contributes highest to total PM_2.5 (∼80 μg/m³), followed by industry (∼70 μg/m³) in North India. Energy and agriculture contribute ∼25 μg/m³ and ∼16 μg/m³ to total PM_2.5, while SOA contributes <5 μg/m³. In Delhi, industry and residential activities contribute to 80% of total PM_2.5.

Persistence of initial conditions in continental scale air quality simulations

This study investigates the effect of initial conditions (IC) for pollutant concentrations in the atmosphere and soil on simulated air quality for two continental-scale Community Multiscale Air Quality (CMAQ) model applications. One of these applications was performed for springtime and the second for summertime. Results show that a spin-up period of ten days commonly used in regional-scale applications may not be sufficient to reduce the effects of initial conditions to less than 1% of seasonally-averaged surface ozone concentrations everywhere while 20 days were found to be sufficient for the entire domain for the spring case and almost the entire domain for the summer case. For the summer case, differences were found to persist longer aloft due to circulation of air masses and even a spin-up period of 30 days was not sufficient to reduce the effects of ICs to less than 1% of seasonally-averaged layer 34 ozone concentrations over the southwestern portion of the modeling domain. Analysis of the effect of soil initial conditions for the CMAQ bidirectional NH₃ exchange model shows that during springtime they can have an important effect on simulated inorganic aerosols concentrations for time periods of one month or longer. The effects are less pronounced during other seasons. The results, while specific to the modeling domain and time periods simulated here, suggest that modeling protocols need to be scrutinized for a given application and that it cannot be assumed that commonly-used spin-up periods are necessarily sufficient to reduce the effects of initial conditions on model results to an acceptable level. What constitutes an acceptable level of difference cannot be generalized and will depend on the particular application, time period and species of interest. Moreover, as the application of air quality models is being expanded to cover larger geographical domains and as these models are increasingly being coupled with other modeling systems to better represent air-surface-water exchanges, the effects of model initialization in such applications needs to be studied in future work.

Advanced error diagnostics of the CMAQ and Chimere modelling systems within the AQMEII3 model evaluation framework

The work here complements the overview analysis of the modelling systems participating in the third phase of the Air Quality Model Evaluation International Initiative (AQMEII3) by focusing on the performance for hourly surface ozone by two modelling systems, Chimere for Europe and CMAQ for North America. The evaluation strategy outlined in the course of the three phases of the AQMEII activity, aimed to build up a diagnostic methodology for model evaluation, is pursued here and novel diagnostic methods are proposed. In addition to evaluating the "base case" simulation in which all model components are configured in their standard mode, the analysis also makes use of sensitivity simulations in which the models have been applied by altering and/or zeroing lateral boundary conditions, emissions of anthropogenic precursors, and ozone dry deposition. To help understand of the causes of model deficiencies, the error components (bias, variance, and covariance) of the base case and of the sensitivity runs are analysed in conjunction with timescale considerations and error modelling using the available error fields of temperature, wind speed, and NO x concentration. The results reveal the effectiveness and diagnostic power of the methods devised (which remains the main scope of this study), allowing the detection of the timescale and the fields that the two models are most sensitive to. The representation of planetary boundary layer (PBL) dynamics is pivotal to both models. In particular, (i) the fluctuations slower than ~ 1.5 days account for 70-85 % of the mean square error of the full (undecomposed) ozone time series; (ii) a recursive, systematic error with daily periodicity is detected, responsible for 10-20 % of the quadratic total error; (iii) errors in representing the timing of the daily transition between stability regimes in the PBL are responsible for a covariance error as large as 9 ppb (as much as the standard deviation of the network-average ozone observations in summer in both Europe and North America); (iv) the CMAQ ozone error has a weak/negligible dependence on the errors in NO2, while the error in NO2 significantly impacts the ozone error produced by Chimere; (v) the response of the models to variations of anthropogenic emissions and boundary conditions show a pronounced spatial heterogeneity, while the seasonal variability of the response is found to be less marked. Only during the winter season does the zeroing of boundary values for North America produce a spatially uniform deterioration of the model accuracy across the majority of the continent.