Examination of the Community Multiscale Air Quality (CMAQ) model performance over the North American and European domains

Copyright 2011 Elsevier B.V., All rights reserved.The CMAQ modeling system has been used to simulate the air quality for North America and Europe for the entire year of 2006 as part of the Air Quality Model Evaluation International Initiative (AQMEII). The operational model performance of tropospheric ozone (O), fine particulate matter (PM) and total particulate matter (PM) for the two continents has been assessed. The model underestimates daytime (8am-8pm LST) O mixing ratios by 13% in the winter for North America, primarily due to an underestimation of daytime O mixing ratios in the middle and lower troposphere from the lateral boundary conditions. The model overestimates winter daytime O mixing ratios in Europe by an average of 8.4%. The model underestimates daytime O by 4-5% in the spring for both continents, while in the summer daytime O is overestimated by 9.8% for North America and slightly underestimated by 1.6% for Europe. The model overestimates daytime O in the fall for both continents, grossly overestimating daytime O by over 30% for Europe. The performance for PM varies both seasonally and geographically for the two continents. For North American, PM is overestimated in the winter and fall, with an average Normalized Mean Bias (NMB) greater than -30%, while performance in the summer is relatively good, with an average NMB of -4.6%. For Europe, PM is underestimated throughout the entire year, with the NMB ranging from -24% in the fall to -55% in the winter. PM is underestimated throughout the year for both North America and Europe, with remarkably similar performance for both continents. The domain average NMB for PM ranges between -45% and -65% for the two continents, with the largest underestimation occurring in the summer for North American and the winter for Europe.Peer reviewedSubmitted Versio

Evaluation of the performance of four chemical transport models in predicting the aerosol chemical composition in Europe in 2005

© Author(s) 2016.Four regional chemistry transport models were applied to simulate the concentration and composition of particulate matter (PM) in Europe for 2005 with horizontal resolution 20 km. The modelled concentrations were compared with the measurements of PM chemical composition by the European Monitoring and Evaluation Programme (EMEP) monitoring network. All models systematically underestimated PM10 and PM2:5 by 10–60 %, depending on the model and the season of the year, when the calculated dry PM mass was compared with the measurements. The average water content at laboratory conditions was estimated between 5 and 20% for PM2:5 and between 10 and 25% for PM10. For majority of the PM chemical components, the relative underestimation was smaller than it was for total PM, exceptions being the carbonaceous particles and mineral dust. Some species, such as sea salt and NO3, were overpredicted by the models. There were notable differences between the models’ predictions of the seasonal variations of PM, mainly attributable to different treatments or omission of some source categories and aerosol processes. Benzo(a)pyrene concentrations were overestimated by all the models over the whole year. The study stresses the importance of improving the models’ skill in simulating mineral dust and carbonaceous compounds, necessity for high-quality emissions from wildland fires, as well as the need for an explicit consideration of aerosol water content in model–measurement comparison.Peer reviewedFinal Published versio

The impact of assimilation of satellite derived wind observations for the prediction of a monsoon depression over India using a mesoscale model

Copyright 2008 Elsevier B.V., All rights reserved.The Penn State/NCAR mesoscale model (MM5) has been used in this study to ingest and assimilate the INSAT-CMV (Indian National Satellite System-Cloud Motion Vector) wind observations using analysis nudging (four-dimensional data assimilation, FDDA) to improve the prediction of a monsoon depression which occurred over the Bay of Bengal, India during 28 July 2005 to 31 July 2005. To determine the impact of assimilation of INSAT-CMV winds on the prediction of a monsoon depression, three sets of numerical experiments (NOFDDA, FDDA and FDDA CMV) were designed. While the FDDA CMV run assimilated satellite derived winds only, the FDDA run assimilated both satellite and conventional observations. The NOFDDA run used neither satellite nor conventional observations. The results of the study indicate that the simulated sea level pressure field from the FDDA run is more consistent with the sea level pressure field from NCEP-FNL compared to the FDDA CMV and NOFDDA runs. The highest correlation and lowest rms error of the sea level pressure field are associated with the FDDA run, and this provides a quantitative verification of the improvement due to the assimilation of satellite derived winds and the conventional upper air observations for the prediction of monsoon depression. All the three model simulated winds are in good agreement with the analysis winds at 850 hPa, 500 hPa and 200 hPa levels. The simulated structure of the spatial precipitation pattern for the assimilation experiments (FDDA and FDDA CMV) are closer to the TRMM observations with more rainfall simulated over the east coast regions in the assimilation experiments. The rms errors of the wind speed for the FDDA run show lower values at 500 hPa for all the three model runs, with a reduction in all three levels of up to 0.8-1.4 m s for the FDDA run and 0.5-1.9 m s for the FDDA CMV run with respect to the NOFDDA run. The statistical significance of the sea level pressure and the precipitation differences between the FDDA and the NOFDDA as well as the differences between the FDDA CMV and the NOFDDA have been calculated using the two-tailed Student's t-test and were found to be statistically significant. The influence of varying the nudging coefficients in the FDDA experiment has been studied.Peer reviewe

The impact of assimilation of MODIS data for the prediction of a tropical low-pressure system over India using a mesoscale model

Copyright 2009 Elsevier B.V., All rights reserved.A low-pressure system formed over the Bay of Bengal, India on 15 October 2003 and crossed the east coast of India during the late hours of 17 October 2003. The system, which provided copious rainfall over the Bay of Bengal and stations on the east coast, is investigated in this study using the fifth Generation Mesoscale Model (MM5). Three sets of numerical experiments are designed in this study. While the first set utilizes National Center for Environmental Prediction - Aviation (NCEP-AVN) analysis (for the initial conditions and lateral boundary conditions) only in the MM5 simulation, the second set utilizes the vertical profiles of temperature and humidity obtained from Moderate Resolution Imaging Spectroradiometer (MODIS) (as well as a few radiosonde station data) to provide an improved analysis. The third set used the vertical profiles of temperature and humidity from MODIS alone to provide an improved analysis. The results of the three sets of simulation are compared with one another as well as with the analysis and observations. It is found that the predicted sea level pressure of the MM5 simulation which utilized the improved analysis: reproduces the large-scale structure of the low-pressure system as manifested in the NCEP-AVN analysis; provides a stronger and deeper low-pressure system as seen from the sea level pressure field; and shows a larger northward extent of the associated precipitation pattern as compared with the simulation with just the analysis. The results of the third experiment (impact of vertical profiles of temperature and humidity using MODIS alone) compare well with the results of the second experiment except that in the former, the associated cyclonic circulation in the lower troposphere appears weaker. The results of this study, although restricted to a single case study, demonstrate that inclusion of MODIS derived vertical temperature and humidity profiles together with radiosonde data caused a favourable impact on the simulated structure of the low-pressure system.Peer reviewe

The effect of satellite and conventional meteorological data assimilation on the mesoscale modeling of monsoon depressions over India

Copyright 2008 Elsevier B.V., All rights reserved.The Fifth Generation Mesoscale Model (MM5) is used to study the effect of assimilated satellite and conventional data on the prediction of three monsoon depressions over India using analysis nudging. The satellite data comprised the vertical profiles of temperature and humidity (NOAA-TOVS: - National Oceanic and Atmospheric Administration-TIROS Operational Vertical Sounder; MODIS: - MODerate resolution Imaging Spectroradiometer) and the surface wind vector over the sea (QuikSCAT: - Quick Scatterometer); the conventional meteorological data included the upper-air and surface data from the India Meteorological Department (IMD). Two sets of numerical experiments are performed for each case: the first set, NOFDDA (no nudging), utilizes NCEP reanalysis (for the initial conditions and lateral boundary conditions) in the simulation, the second set, FDDA, employs the satellite and conventional meteorological data for an improved analysis through analysis nudging. Two additional experiments are performed to study the effect of increased vertical and horizontal resolution as well as convective parameterization for one of the depressions for which special fields observations were available. The results from the simulation are compared with each other and with the analysis and observations. The results show that the predicted sea level pressure (SLP), the lower tropospheric cyclonic circulation, and the precipitation of the FDDA simulation reproduced the large-scale structure of the depression as manifested in the NCEP reanalysis. The simulation of SLP using no assimilation high-resolution runs (HRSKF10KM, HRSKF3.3KM) with the Kain-Fritsch cumulus parameterization scheme appeared poor in comparison with the FDDA run, while the no assimilation high-resolution runs (HRSGR10KM, HRSGR3.3KM) with the Grell cumulus scheme provided better results. However, the space correlation and the root mean square (rms) error of SLP for the HRSKF10KM was better than the FDDA; the largest and smallest space correlation for HRSKF10KM, FDDA, and HRSGR10KM were 0.894 and 0.623, 0.663 and 0.195, and 0.733 and 0.338 respectively; the smallest and largest rms error for HRSKF10KM, FDDA and HRSGR10KM were 1.879 and 5.245, 2.308 and 4.242, and 2.055 and 4.909 respectively. The precipitation simulations with the 3.3 km high-resolution, no assimilation runs performed no better than the precipitation simulation with the FDDA run. Thus, a significant finding of this study is that over the Indian monsoon region, the improvements in the simulation using nudging in the FDDA run are of similar magnitude (or better) than the improvements in the simulation due to high-resolution and to cumulus parameterization sensitivity. The improvements in the FDDA run due to analysis nudging were also verified in two more depression cases. The current operational regional models in India do not incorporate the data assimilation of NOAA-TOVS/MODIS and QuikSCAT satellite data, and hence the results of this study are relevant to operational forecasters in India and their ongoing efforts to improve weather forecasting using satellite data sets.Peer reviewe

Application of chemical transport model CMAQ to policy decisions regarding PM2.5 in the UK

This document is the Accepted Manuscript, made available under the terms of the Creative Commons Attribution-NonCommercial-NoDerivatives License CC BY NC-ND 4.0 (http://creativecommons.org/licenses/by-nc-nd/4.0/), which permits non-commercial re-use, distribution, and reproduction in any medium, provided the original work is properly cited, and is not altered, transformed, or built upon in any way.This paper shows how the advanced chemical transport model CMAQ can be used to estimate future levels of PM2.5 in the UK, the key air pollutant in terms of human health effects, but one which is largely made up from the formation of secondary particulate in the atmosphere. By adding the primary particulate contribution from typical urban roads and including a margin for error, it is concluded that the current indicative limit value for PM2.5 will largely be met in 2020 assuming 2006 meteorological conditions. Contributions to annual average regional PM2.5 concentration from wild fires in Europe in 2006 and from possible climate change between 2006 and 2020 are shown to be small compared with the change in PM2.5 concentration arising from changes in emissions between 2006 and 2020. The contribution from emissions from major industrial sources regulated in the UK is estimated from additional CMAQ calculations. The potential source strength of these emissions is a useful indicator of the linearity of the response of the atmosphere to changes in emissions. Uncertainties in the modelling of regional and local sources are taken into account based on previous evaluations of the models. Future actual trends in emissions mean that exceedences of limit values may arise, and these and further research into PM2.5 health effects will need to be part of the future strategy to manage PM2.5 concentrationsPeer reviewedFinal Accepted Versio

Operational model evaluation for particulate matter in Europe and North America in the context of AQMEII

Ten state-of-the-science regional air quality (AQ) modeling systems have been applied to continental-scale domains in North America and Europe for full-year simulations of 2006 in the context of Air Quality Model Evaluation International Initiative (AQMEII), whose main goals are model intercomparison and evaluation. Standardised modeling outputs from each group have been shared on the web-distributed ENSEMBLE system, which allows statistical and ensemble analyses to be performed. In this study, the one-year model simulations are inter-compared and evaluated with a large set of observations for ground-level particulate matter (PK10 and PM2.5) and its chemical components. Modeled concentrations of gaseous PM precursors, SO2 and NO2, have also been evaluated against observational data for both continents. Furthermore, modeled deposition (dry and wet) and emissions of several species relevant to PM are also inter-compared. The unprecedented scale of the exercise (two continents, one full year, fifteen modeling groups) allows for a detailed description of AQ model skill and uncertainty with respect to PM. Analyses of PM10 yearly time series and mean diurnal cycle show a large underestimation throughout the year for the AQ models included in AQMEII. The possible causes of PM bias, including errors in the emissions and meteorological inputs (e.g., wind speed and precipitation), and the calculated deposition are investigated. Further analysis of the coarse PM components, PM2.5 and its major components (SO4, NH4, NO3, elemental carbon), have also been performed, and the model performance for each component evaluated against measurements. Finally, the ability of the models to capture high PM concentrations has been evaluated by examining two separate PM2.5 episodes in Europe and North America. A large variability among models in predicting emissions, deposition, and concentration of PM and its precursors during the episodes has been found. Major challenges still remain with regards to identifying and eliminating the sources of PM bias in the models. Although PM2.5 was found to be much better estimated by the models than PM10, no model was found to consistently match the observations for all locations throughout the entire yearPeer reviewedFinal Accepted Versio

Model evaluation and ensemble modelling of surface-level ozone in Europe and North America in the context of AQMEII

More than ten state-of-the-art regional air quality models have been applied as part of the Air Quality Model Evaluation International Initiative (AQMEII). These models were run by twenty independent groups in Europe and North America. Standardised modelling outputs over a full year (2006) from each group have been shared on the web-distributed ENSEMBLE system, which allows for statistical and ensemble analyses to be performed by each group. The estimated ground-level ozone mixing ratios from the models are collectively examined in an ensemble fashion and evaluated against a large set of observations from both continents. The scale of the exercise is unprecedented and offers a unique opportunity to investigate methodologies for generating skilful ensembles of regional air quality models outputs. Despite the remarkable progress of ensemble air quality modelling over the past decade, there are still outstanding questions regarding this technique. Among them, what is the best and most beneficial way to build an ensemble of members? And how should the optimum size of the ensemble be determined in order to capture data variability as well as keeping the error low? These questions are addressed here by looking at optimal ensemble size and quality of the members. The analysis carried out is based on systematic minimization of the model error and is important for performing diagnostic/probabilistic model evaluation. It is shown that the most commonly used multi-model approach, namely the average over all available members, can be outperformed by subsets of members optimally selected in terms of bias, error, and correlation. More importantly, this result does not strictly depend on the skill of the individual members, but may require the inclusion of low-ranking skill-score members. A clustering methodology is applied to discern among members and to build a skilful ensemble based on model association and data clustering, which makes no use of priori knowledge of model skill. Results show that, while the methodology needs further refinement, by optimally selecting the cluster distance and association criteria, this approach can be useful for model applications beyond those strictly related to model evaluation, such as air quality forecasting. (C) 2012 Elsevier Ltd. All rights reserved.Peer reviewedSubmitted Versio

Evaluating the capability of regional-scale air quality models to capture the vertical distribution of pollutants

© Author(s) 2013. This work is distributed under the Creative Commons Attribution 3.0 LicenseThis study is conducted in the framework of the Air Quality Modelling Evaluation International Initiative (AQMEII) and aims at the operational evaluation of an ensemble of 12 regional-scale chemical transport models used to predict air quality over the North American (NA) and European (EU) continents for 2006. The modelled concentrations of ozone and CO, along with the meteorological fields of wind speed (WS) and direction (WD), temperature (T), and relative humidity (RH), are compared against high-quality in-flight measurements collected by instrumented commercial aircraft as part of the Measurements of OZone, water vapour, carbon monoxide and nitrogen oxides by Airbus In-service airCraft (MOZAIC) programme. The evaluation is carried out for five model domains positioned around four major airports in NA (Portland, Philadelphia, Atlanta, and Dallas) and one in Europe (Frankfurt), from the surface to 8.5 km. We compare mean vertical profiles of modelled and measured variables for all airports to compute error and variability statistics, perform analysis of altitudinal error correlation, and examine the seasonal error distribution for ozone, including an estimation of the bias introduced by the lateral boundary conditions (BCs). The results indicate that model performance is highly dependent on the variable, location, season, and height (e.g. surface, planetary boundary layer (PBL) or free troposphere) being analysed. While model performance for T is satisfactory at all sites (correlation coefficient in excess of 0.90 and fractional bias a parts per thousand currency sign 0.01 K), WS is not replicated as well within the PBL (exhibiting a positive bias in the first 100 m and also underestimating observed variability), while above 1000 m, the model performance improves (correlation coefficient often above 0.9). The WD at NA airports is found to be biased in the PBL, primarily due to an overestimation of westerly winds. RH is modelled well within the PBL, but in the free troposphere large discrepancies among models are observed, especially in EU. CO mixing ratios show the largest range of modelled-to-observed standard deviations of all the examined species at all heights and for all airports. Correlation coefficients for CO are typically below 0.6 for all sites and heights, and large errors are present at all heights, particularly in the first 250 m. Model performance for ozone in the PBL is generally good, with both bias and error within 20%. Profiles of ozone mixing ratios depend strongly on surface processes, revealed by the sharp gradient in the first 2 km (10 to 20 ppb km(-1)). Modelled ozone in winter is biased low at all locations in the NA, primarily due to an underestimation of ozone from the BCs. Most of the model error in the PBL is due to surface processes (emissions, transport, photochemistry), while errors originating aloft appear to have relatively limited impact on model performance at the surface. Suggestions for future work include interpretation of the model-to-model variability and common sources of model bias, and linking CO and ozone bias to the bias in the meteorological fields. Based on the results from this study, we suggest possible in-depth, process-oriented and diagnostic investigations to be carried out next.Peer reviewe