Skill and uncertainty of a regional air quality model ensemble

Recently several regional air quality projects were carried out to support the negotiation under the Clean Air For Europe (CAFE) programme by predicting the impact of emission control policies with an ensemble of models. Within these projects, CITYDELTA and EURODELTA, the fate of air quality at the scale of European cities or that of the European continent was studied using several models. In this article we focus on the results of EURODELTA. The predictive skill of the ensemble of models is described for ozone, nitrogen dioxide and secondary inorganic compounds, and the uncertainty in air quality modelling is examined through the model ensemble spread of concentrations.For ozone daily maxima the ensemble spread origin differs from one region to another. In the neighbourhood of cities or in mountainous areas the spread of predicted values does not span the range of observed data, due to poorly resolved emissions or complex-terrain meteorology. By contrast in Atlantic and North Sea coastal areas the spread of predicted values is found to be larger than the observations. This is attributed to large differences in the boundary conditions used in the different models. For NO₂ daily averages the ensemble spread is generally too small compared with observations. This is because models miss highest values occurring in stagnant meteorology in stable boundary layers near cities. For secondary particulate matter compounds the simulated concentration spread is more balanced, observations falling nearly equiprobably within the ensemble, and the spread originates both from meteorology and aerosol chemistry and thermodynamics.     

Two-stage fuzzy-stochastic robust programming: a hybrid model for regional air quality management

In this study, a hybrid two-stage fuzzy-stochastic robust programming (TFSRP) model is developed and applied to the planning of an air-quality management system. As an extension of existing fuzzy-robust programming and two-stage stochastic programming methods, the TFSRP can explicitly address complexities and uncertainties of the study system without unrealistic simplifications. Uncertain parameters can be expressed as probability density and/or fuzzy membership functions, such that robustness of the optimization efforts can be enhanced. Moreover, economic penalties as corrective measures against any infeasibilities arising from the uncertainties are taken into account. This method can, thus, provide a linkage to predefined policies determined by authorities that have to be respected when a modeling effort is undertaken. In its solution algorithm, the fuzzy decision space can be delimited through specification of the uncertainties using dimensional enlargement of the original fuzzy constraints. The developed model is applied to a case study of regional air quality management. The results indicate that reasonable solutions have been obtained. The solutions can be used for further generating pollution-mitigation alternatives with minimized system costs and for providing a more solid support for sound environmental decisions.     

Ground-level ozone in the Pearl River Delta region: Analysis of data from a recently established regional air quality monitoring network

The ambient air quality monitoring data of 2006 and 2007 from a recently established Pearl River Delta (PRD) regional air quality monitoring network are analyzed to investigate the characteristics of ground-level ozone in the region. Four sites covering urban, suburban, rural and coastal areas are selected as representatives for detailed analysis in this paper. The results show that there are distinct seasonal and diurnal cycles in ground-level ozone across the PRD region. Low ozone concentrations are generally observed in summer, while high O₃ levels are typically found in autumn. The O₃ diurnal variations in the urban areas are larger than those at the rural sites. The O₃ concentrations showed no statistically significant difference between weekend and weekdays in contrast to the findings in many other urban areas in the world. The average ozone concentrations are lower in urban areas compared to the sites outside urban centers. Back trajectories are used to show the major air-mass transport patterns and to examine the changes in ozone from the respective upwind sites to a site in the center of the PRD (Wanqingsha). The results show higher average ozone concentrations at the upwind sites in the continental and coastal air masses, but higher 1 h-max O₃ concentrations (by 8–16 ppbv) at the center PRD site under each of air-mass category, suggesting that the ozone pollution in the PRD region exhibits both regional and super-regional characteristics.     

Evaluation of near surface ozone in air quality simulations forced by a regional climate model over Europe for the period 1991–2000

This study aims to evaluate near surface ozone simulated with the modelling system RegCM3/CAMx against ozone measurements from the EMEP database for the recent decade 1991–2000. The RegCM3/CAMx simulations were performed on a 50 km × 50 km grid over Europe driven either by ERA-40 reanalysis (hereafter referred as ERA simulation) or the global circulation model (GCM) ECHAM5 (hereafter referred as ECHAM simulation). A set of statistical metrics is used for the model evaluation, including temporal correlation coefficient, the ratio of the standard deviations and the bias of simulated versus observed values. Overall, a good agreement is found for both ERA and ECHAM simulations at the majority of the selected EMEP stations in all metrics throughout the year based either on monthly or daily ozone values. Based on these results, it is assessed that the modelling system RegCM3/CAMx is suitable to be used for present and future regional climate-air quality simulations with emphasis on near surface ozone. The ERA simulations reproduce more accurately the observed ozone values in comparison to ECHAM simulations because the meteorology of the ERA experiment is closer to real atmospheric conditions than the GCM based experiment. On a seasonal basis, both ERA and ECHAM simulations exhibit a seasonally dependent bias, with winter and spring ozone values being generally under-estimated by the model and summer and autumn values being slightly overestimated. This seasonally dependent bias is also evident from median and peak midday ozone values. However, the highest observed midday ozone peaks in summer, with values higher than 80 ppbv, could not be captured either by ERA or ECHAM simulations. An analysis of day-time and night-time ERA and ECHAM modelled ozone values shows that CAMx performs better during the day-time.     

A regional air quality model for the Kwinana industrial area of Western Australia

The development and application of a three dimensional regional air dispersion model to the Kwinana industrial area is described. Spatial variation in the flow and variations in surface roughness are incorporated and the model results are shown to be in close agreement with visual observations.     

Use of an indoor air quality model (IAQM) to estimate indoor ozone levels

Currently, outdoor ozone levels in many U.S. cities exceed the primary health-based national ambient air quality standard. While outdoor ozone levels are an important measure of the severity of those exceedances, people typically spend more than 80 percent of their time indoors, where ozone levels are lower. Indoor ozone levels range from 10 to 80 percent of outdoor levels, with many people receiving a substantial portion of their ozone exposure while indoors. This paper uses an indoor air quality model (IAQM) to estimate indoor ozone levels by microenvironment type (home, office, and vehicle) and configuration (windows open, windows closed, older construction, weatherized, and air conditioned). The formulation of IAQM is discussed, along with specification of model parameters for ozone. The multicompartment version of IAQM is described, with a single-compartment version used for the analyses. IAQM-calculated ozone indoor-outdoor ratios compare well with research-reported values. Results indicate that ozone peak-concentration indoor-outdoor ratios range as follows: home--0.65 (windows open), 0.36 (air conditioned), 0.23 (typical construction, windows closed), and 0.05 (energy-efficient construction, windows closed); office--0.82 (heating, ventilation and air conditioning systems supplying 100 percent outdoor air), 0.60 (typical HVAC), and 0.32 (energy-efficient HVAC); and vehicle--0.41 (85 mph), 0.33 (55 mph), and 0.21 (10 mph). Analysis results are presented to characterize IAQM's sensitivity to assumed model parameters.     

Vegetation exposure to ozone over the continental United States: Assessment of exposure indices by the Eta-CMAQ air quality forecast model

The main use of air quality forecast (AQF) models is to predict ozone (O₃) exceedances of the primary O₃ standard for informing the public of potential health concerns. This study presents the first evaluation of the performance of the Eta-CMAQ air quality forecast model to predict a variety of widely used seasonal mean and cumulative O₃ exposure indices associated with vegetation using the U.S. AIRNow O₃ observations. These exposure indices include two concentration-based O₃ indices, M7 and M12 (the seasonal means of daytime 7-h and 12-h O₃ concentrations, respectively), and three cumulative exposure-based indices, SUM06 (the sum of all hourly O₃ concentrations ≥ 0.06 ppm), W126 (hourly concentrations weighed by a sigmoidal weighting function), and AOT40 (O₃ concentrations accumulated over a threshold of 40 ppb during daylight hours). During a three-month simulation (July–September 2005), the model over predicted the M7 and M12 values by 8–9 ppb, or a NMB value of 19% and a NME value of 21%. The model predicts a central belt of high O₃ extending from Southern California to Middle Atlantic where the seasonal means, M7 and M12 (the seasonal means of daytime 7-h and 12-h O₃ concentrations), are higher than 50 ppbv. In contrast, the model is less capable of reproducing the observed cumulative indices. For AOT40, SUM06 and W126, the NMB and NME values are two- to three-fold of that for M7, M12 or peak 8-h O₃ concentrations. The AOT40 values range from 2 to 33 ppm h by the model and from 1 to 40 ppm h by the monitors. There is a significantly higher AOT40 value experienced in the United States in comparison to Europe. The domain-wide mean SUM06 value is 14.4 ppm h, which is about 30% higher than W126, and 40% higher than AOT40 calculated from the same 3-month hourly O₃ data. This suggests that SUM06 and W126 represent a more stringent standard than AOT40 if either the SUM06 or the W126 was used as a secondary O₃ standard. Although CMAQ considerably over predicts SUM06 and W126 values at the low end, the model under predicts the extreme high exposure values (>50 ppm h). Most of these extreme high values are found at inland California sites. Based on our analysis, further improvement of the model is needed to better capture cumulative exposure indices.     

Coupling GEOS-CHEM with a regional air pollution model for Greece

The sensitivity of regional air quality modeling simulations to boundary conditions over Greece is investigated, for various synoptic conditions. For this purpose, a global to mesoscale model-chain is developed and applied, coupling the individual models' simulations. The global chemical transport model GEOS-CHEM, applied in a one-way nested procedure, is used to drive the regional UAM-V chemical dispersion model with time-varying lateral and top boundary conditions. The results of the coupling procedure are compared with the MINOS campaign measurements at Finokalia (Southern Greece) during the period from 1 to 16 August 2001 which is mainly characterized by an interchange of two synoptic types, High-Low and Long Wave trough.The comparison between the simulation results and the measurements reveals that the coupling procedure captures satisfactorily the range of observed CO concentrations at the southern part of Greece. The most severe deviations are observed under strongly variable atmospheric circulation, when no distinct synoptic circulation is allowed to be established in the area. Regarding O₃, the highest, though underestimated, surface concentrations are simulated under Long Wave trough conditions due to the influence of the ozone inflow predicted by GEOS-CHEM at the western boundary of the innermost domain and/or under enhanced NOy emissions arriving at Finokalia from urban and ships plumes.     

Variability of ozone air quality indicators in selected metropolitan statistical areas

This article analyzes numerical variability in ozone air quality data to understand how this variability affects the number of violations seen each year in metropolitan statistical areas (MSAs). Three commonly cited violation indices are used: 1) the annual number of expected exceedances averaged over 3 years is greater than 1; 2) the n+ 1th hourly value in n years of data is greater than 0.12 ppm; and 3) the annual number of expected exceedances is greater than 1. Only the first index is consistent with applicable regulations. The analyses indicate that about 23 percent of all MSAs with valid data had one or more change in their ozone violation status between 1979 and 1987. This change in status occurred for approximately 7 percent of all MSA-years of available data. This statistic was about one-third of the value usually obtained when the two incorrect, but commonly used, criteria of ozone violations are used.     

Economic damages of ozone air pollution to crops using combined air quality and GIS modelling

This study aims at presenting a combined air quality and GIS modelling methodological approach in order to estimate crop damages from photochemical air pollution, depict their spatial resolution and assess the order of magnitude regarding the corresponding economic damages. The analysis is conducted within the Greater Thessaloniki Area, Greece, a Mediterranean territory which is characterised by high levels of photochemical air pollution and considerable agricultural activity. Ozone concentration fields for 2002 and for specific emission reduction scenarios for the year 2010 were estimated with the Ozone Fine Structure model in the area under consideration. Total economic damage to crops turns out to be significant and estimated to be approximately 43 M€ for the reference year. Production of cotton presents the highest economic loss, which is over 16 M€, followed by table tomato (9 M€), rice (4.2 M€), wheat (4 M€) and oilseed rape (2.8 M€) cultivations. Losses are not spread uniformly among farmers and the major losses occur in areas with valuable ozone-sensitive crops. The results are very useful for highlighting the magnitude of the total economic impacts of photochemical air pollution to the area’s agricultural sector and can potentially be used for comparison with studies worldwide. Furthermore, spatial analysis of the economic damage could be of importance for governmental authorities and decision makers since it provides an indicative insight, especially if the economic instruments such as financial incentives or state subsidies to farmers are considered.     

Evaluation of air quality models for the simulation of a high ozone episode in the Houston metropolitan area

A high ozone event in the Houston–Galveston–Brazoria area was utilized to study the shortcomings of the current air quality models. To improve the baseline simulations with the Comprehensive Air quality Model with Extensions (CAMx) for developing the state implementation plan, the Texas Commission on Environmental Quality (TCEQ) imputed emissions of highly reactive volatile organic compounds (HRVOCs) by scaling the amount of fugitive emissions of olefins to co-emitted NO_x from selected point sources, effectively multiplying by 3–12 times over the regular inventory values. In this paper, CAMx and the Community Multiscale Air Quality (CMAQ) model were used to determine if the imputed HRVOC emissions were consistent with the observed atmospheric conditions. With the base emissions, CMAQ and CAMx both with the Carbon-Bond 4 (CB-4) mechanism simulated similar ozone concentrations. But with the imputed HRVOC emissions, CMAQ predicted lower ozone peaks than CAMx in the vicinity and downwind of the Ship Channel and other highly HRVOC-rich areas. Based on analyses of sensitivity simulations of CMAQ with different emission inputs and vertical diffusion algorithms in the model, we found that the modeled atmosphere lacked reactivity to produce the observed high ozone event. Although the imputed HRVOC emissions improved ozone prediction at the surface sites, but the ethylene concentrations were not consistent with the measurements at the super sites (La Porte and Clinton) and by NOAA aircraft. Several sensitivity tests designed to provide additional radicals into the system and other research results suggested that the lack of reactivity may need to be corrected by targeted, and probably of episodic, increase of HRVOC emissions, from the sources in the Houston Ship Channel. Additional investigation of the ozone production efficiency for different chemical mechanisms is necessary to pinpoint the emissions uncertainty issues.     

Statistical and diagnostic evaluation of the ADMS-Urban model compared with an urban air quality monitoring network

The present study examines the behaviour of the ADMS-Urban air quality forecasting model in predicting dispersion of traffic-related pollutants in urban areas. The study has been carried out in Ravenna (NE Italy), a medium-sized town where pollution produced by vehicle traffic accounts for most of the emissions. ADMS-Urban performances have been assessed through statistical analysis, by comparing carbon monoxide concentrations (vehicle traffic tracing pollutant) estimated by the model with concentrations measured by stations of the air quality monitoring network. Although the correspondence of values estimated by ADMS-Urban with measured values turns out to be satisfactory, the study shows that the model tends to produce an underestimated value compared with the actual situation, and identifies a corrective method that makes it possible to improve the relevant performances. Furthermore, the diagnostic analysis highlights that the model performances depend upon some meteorological parameters.     

区域低碳交通评价模型研究

在中国机动车快速增长的前提下,低碳交通运输是实现节能减排、发展低碳经济的重要组成部分。借鉴戴利的消费和福利绩效公式,基于导向性、数据可获性与可比性等基本原则,构建区域低碳交通评价指标体系,该体系包括低碳交通能效指标和低碳交通运营指标。基于灰关联分析法和模糊数学方法,建立区域低碳交通模糊评价模型。最后,对江苏2010年的低碳交通发展水平进行评价,结果表明,该年江苏低碳交通发展水平为"良",今后发展的重点在于大力发展城市轨道交通、积极开展低碳公路建设、开展碳税制度和优化公交体系等。     

On the influence of chemical initial and boundary conditions on annual regional air quality model simulations for North America

The influence of chemical initial conditions and chemical lateral boundary conditions (CLBCs) on long-term regional air quality model simulations was investigated using outputs from an annual simulation of the year 2002 on a North American domain. This simulation was carried out using the AURAMS regional air quality model. It was subdivided into three multi-month segments with two overlap periods (May 15–30 and September 1–30) to allow the segments to be run in parallel. For this approach to work, model predictions had to match very closely by the end of the two-week and four-week overlap periods. The time required for the values of daily domain-average surface PM_2.5 concentration to match for the two simulation segments associated with each of the two overlap periods was four and six days, respectively. For individual locations within the model domain, however, the required spin-up period was as much as nine days, considerably longer than the 2–4-day spin-up period usually assumed in the literature. For ozone, on the other hand, the daily domain-average surface ozone concentration values did not converge for either overlap period. A zero-gradient CLBC had been used for all run segments and species. When a time-invariant CLBC for ozone was used instead, the daily domain-average surface ozone concentration values behaved more realistically and did converge after fewer than three days. A similar improvement was also obtained for individual locations, but with spin-up periods of up to nine days. Model chemical spin-up time thus seems to be dependent on the species considered, the time required for the influence of the inflow boundaries to reach all locations within the domain, and the impact of local emissions sources. These results suggest the use of a spin-up period of longer than one week for a large (continental) domain and long-term simulation of PM_2.5 and O₃ rather than the 2–4 days commonly assumed in the literature.     

developing a practical dispersion model for an air quality region

A series of computer models have been developed to predict air quality in the New York/New Jersey/Connecticut Air Quality Region. Efforts have been directed at models which have a shorter time scale than climatological models, and which are capable of providing better recommendations for effective abatement and planning, but use input data presently available.

The basic dispersion model for these investigations is a steady-state,nondivergent Gaussian-type model. A modified inventory of SO₂ sources,based on published data for the New York/New Jersey/Connecticut Air Quality Region, was prepared for use with the model. The basic model has been subjected to various internal sensitivity analyses, in which was isolated the variation produced in the pollutant concentration by a given change in each of the factors that contribute, e.g., wind speed, wind direction,mixing depth, stability conditions, source strengths, and grid size for the area sources.

To date, validation tests of the model have been made against the July and August 1969 data for the ten telemetering stations of the New York City Aerometric Network. Hourly as well as averaged concentrations were considered. Various sets of meteorological data from the network stations and the three area airports, were compared and tested. Additional tests, particularly for the winter season, are needed to substantiate the preliminary conclusions suggested by the results to date.

Considerable insight into the relative importance of model components has been acquired from the sensitivity studies. Furthermore the validation results lend support to the belief that a reasonably simple, practical dispersion model can be developed for the region.     

Integrated biomonitoring of air quality with plants and lichens: a case study on ambient ozone from central Italy

A biennial integrated survey, based on the use of vascular plants for the bioindication of the effects of tropospheric ozone together with the use of automatic analysers of ozone, as well as the mapping of lichen biodiversity was performed in the area of Castelfiorentino (Tuscany, central Italy). Photochemically produced ozone proved to be a fundamental presence during the warm season, with maximum hourly means reaching 114 ppb, exceeding the information threshold as fixed by EU: the use of supersensitive tobacco Bel-W3 confirmed the opportunity of carrying out detailed cost-effective monitoring surveys. The potential for didactical and educational implications of this methodology are appealing. Critical levels set up for the protection of vegetation have exceeded considerably. The comparison of biomass productivity in sensitive and resistant individuals (NC-S and NC-R white clover clones, in the framework of an European network) provided evidence that ambient ozone levels are associated with relevant reduction (up to 30%) in the performance of sensitive material; effects on flowering were also pronounced. The economic assessment of such an impact deserves attention. Mapping of epiphytic lichen biodiversity – which has been used to monitor air quality worldwide – was not related to ozone geographical distribution as depicted by tobacco response.     

Spatial inhomogeneity in pollutant concentrations, and their implications for air quality model evaluation

The Southern California Air Quality Study database provides a valuable resource with which to test urban-scale photochemical models and to achieve a better understanding of the atmospheric dynamics of pollutant formation. The CIT model was evaluated using the SCAQS database according to traditional model performance guidelines. A first application, reported previously, focused on model enhancement and application of the model to the 27–29 August 1987 episode. This study evaluates the CIT model using the 24–25 June SCAQS episode, providing further evaluation of the model. Results show that the CIT airshed model can follow the diurnal variations of reactive species and the transport for relatively unreactive species. The normalized gross error for ozone was 31 % in June compared to 38% in August. However, to fully judge model performance in proper perspective, a question arises: “How well do the measurements reflect the air quality surrounding the monitoring station, not just in that location?” This is an important but seldom quantitatively considered factor, not only in model evaluation but in the study of health effects as well. Analyses indicate that individual concentration measurements only approximately represent the true volume-averaged concentrations within a computational grid cell and that significant spatial variations exist. Thus any evaluation of models using these data sets should take these local variations into consideration. A series of tests found that the local inhomogeneities had a normalized gross error in the range of 25–45% depending on the pollutant. In this context, the performance of the CIT model is consistent with known modeling limitations such as emissions inventories and sub-grid scale variation of observations.