An Assessment of the Emissions Inventory Processing Systems EMS-2001 and SMOKE in Grid-Based Air Quality Models

Abstract

In the United States, emission processing models such as Emissions Modeling System-2001 (EMS-2001), Emissions Preprocessor System-Version 2.5 (EPS2.5), and the Sparse Matrix Operator Kernel Emissions (SMOKE) model are currently being used to generate gridded, hourly, speciated emission inputs for urban and regional-scale photochemical models from aggregated pollutant inventories. In this study, two models, EMS-2001 and SMOKE, were applied with their default internal data sets to process a common inventory database for a high ozone (O₃) episode over the eastern United States using the Carbon Bond IV (CB4) chemical speciation mechanism. A comparison of the emissions processed by these systems shows differences in all three of the major processing steps performed by the two models (i.e., in temporal allocation, spatial allocation, and chemical speciation). Results from a simulation with a photochemical model using these two sets of emissions indicate differences on the order of ±20 ppb in the predicted 1-hr daily maximum O₃ concentrations. It is therefore critical to develop and implement more common and synchronized temporal, spatial, and speciation cross-reference systems such that the processes within each emissions model converge toward reasonably similar results. This would also help to increase confidence in the validity of photochemical grid model results by reducing one aspect of modeling uncertainty.     

Space-time analysis of the Air Quality Model Evaluation International Initiative (AQMEII) Phase 1 air quality simulations

This study presents an evaluation of summertime ozone concentrations over North America (NA) and Europe (EU) using the database generated from Phase 1 of the Air Quality Model Evaluation International Initiative (AQMEII). The analysis focuses on identifying temporal and spatial features that can be used to stratify operational model evaluation metrics and to test the extent to which the various modeling systems can replicate the features seen in the observations. Using a synoptic map typing approach, it is demonstrated that model performance varies with meteorological conditions associated with specific synoptic-scale flow patterns over both eastern NA and EU. For example, the root mean square error of simulated daily maximum 8-hr ozone was twice as high when cloud fractions were high compared with when cloud fractions were low over eastern NA. Furthermore, results show that over both NA and EU the regional models participating in AQMEII were able to better reproduce the observed variance in ambient ozone levels than the global model used to specify chemical boundary conditions, although the variance simulated by almost all regional models is still less that the observed variance on all spatiotemporal scales. In addition, all modeling systems showed poor correlations with observed fluctuations on the intraday time scale over both NA and EU. Furthermore, a methodology is introduced to distinguish between locally influenced and regionally representative sites for the purpose of model evaluation. Results reveal that all models have worse model performance at locally influenced sites. Overall, the analyses presented in this paper show how observed temporal and spatial information can be used to stratify operational model performance statistics and to test the modeling systems’ ability to replicate observed temporal and spatial features, especially at scales the modeling systems are designed to capture.

Implications: The analyses presented in this paper demonstrate how observed temporal and spatial information can be used to stratify operational model performance and to test the modeling systems’ ability to replicate observed temporal and spatial features. Decisions for the improvement of regional air quality models should be based on the information derived from only regionally representative sites.     

Demonstrating Attainment of the Air Quality Standards: Integration of Observations and Model Predictions into the Probabilistic Framework

ABSTRACT

This paper introduces an integrated observational-modeling approach to transform the deterministic nature of attainment demonstrations of the National Ambient Air Quality Standard (NAAQS) into the probabilistic framework. While the methods presented here can be used to address any air quality standard that is based on extreme values, this paper focuses on the application to the 1-hr and 8-hr NAAQS for ozone. Extreme value statistics and resampling techniques are applied to estimate the probability of exceeding the NAAQS for both 1-hr and 8-hr ozone concentrations. Within the integrated observation-modeling analysis approach, we show that the model-to-model differences in the predicted responses to emission reductions are smaller than the model-to-model differences in predicted absolute ozone concentrations. We illustrate that the emission reductions stemming from a real-world emission control strategy would substantially reduce the probability of exceeding the NAAQS over a large portion of the eastern United States, especially for the 8-hr average ozone concentrations.     

A retrospective comparison of model-based forecasted PM2.5 concentrations with measurements

This study presents an assessment of the performance of the Community Multiscale Air Quality (CMAQ) photochemical model in forecasting daily PM2.5 (particulate matter < or = 2.5 microm in aerodynamic diameter) mass concentrations over most of the eastern United States for a 2-yr period from June 14, 2006 to June 13, 2008. Model predictions were compared with filter-based and continuous measurements of PM2.5 mass and species on a seasonal and regional basis. Results indicate an underprediction of PM2.5 mass in spring and summer, resulting from under-predictions in sulfate and total carbon concentrations. During winter, the model overpredicted mass concentrations, mostly at the urban sites in the northeastern United States because of overpredictions in unspeciated PM2.5 (suggesting possible overestimation of primary emissions) and sulfate. A comparison of observed and predicted diurnal profiles of PM2.5 mass at five sites in the domain showed significant discrepancies. Sulfate diurnal profiles agreed in shape across three sites in the southern portion of the domain but differed at two sites in the northern portion of the domain. Predicted organic carbon (OC) profiles were similar in shape to mass, suggesting that discrepancies in mass profiles probably resulted from the underprediction in OC. The diurnal profiles at a highly urbanized site in New York City suggested that the overpredictions at that site might be resulting from overpredictions during the morning and evening hours, displayed as sharp peaks in predicted profiles. An examination of the predicted planetary boundary layer (PBL) heights also showed possible issues in the modeling of PBL.     

Parameterization of Albedo over Heterogeneous Surfaces in Coupled Land-Atmosphere Schemes for Environmental Modeling. Part I: Theoretical Background

In grid-based environmental models, the underlying surface consists of patches of solid and liquid parts and different plant communities, creating a very heterogeneous picture in the grid cell. In these cases, numerical modelers usually use a simple arithmetic average to determine the grid-cell albedo, a key variable in the parameterization of the land-surface radiative transfer over the grid cell. The object of this paper is to consider the assumptions for aggregating the albedo over a very heterogeneous surface where various surfaces occur at different heights, and, then propose a method for deriving a general expression for it. The suggested expression for the albedo is compared with the conventional approach, for the two-patches grid-cell with a simple geometrical distribution and different heights of its components. A numerical test is performed to compare the two approaches by numerical simulation of the evolution of the surface temperature over the particular grid-cell. Specifically, a one-dimensional land-surface model was applied to an isolated rocky grid-cell with a hole in the center; the model was forced with meteorological observations taken on July 17, 1999 in Philadelphia, PA.     

Examination of model predictions at different horizontal grid resolutions

While fluctuations in meteorological and air quality variables occur on a continuum of spatial scales, the horizontal grid spacing of coupled meteorological and photochemical models sets a lower limit on the spatial scales that they can resolve. However, both computational costs and data requirements increase significantly with increasing grid resolution. Therefore, it is important to examine, for any given application, whether the expected benefit of increased grid resolution justifies the extra costs. In this study, we examine temperature and ozone observations and model predictions for three high ozone episodes that occurred over the northeastern United States during the summer of 1995. In the first set of simulations, the meteorological model RAMS4a was run with three two-way nested grids of 108/36/12 km grid spacing covering the United States and the photochemical model UAM-V was run with two grids of 36/12 km grid spacing covering the eastern United States. In the second set of simulations, RAMS4a was run with four two-way nested grids of 108/36/12/4 km grid spacing and UAM-V was run with three grids of 36/12/4 km grid spacing with the finest resolution covering the northeastern United States. Our analysis focuses on the comparison of model predictions for the finest grid domain of the simulations, namely, the region overlapping the 12 km and 4 km domains. A comparison of 12 km versus 4 km fields shows that the increased grid resolution leads to finer texture in the model predictions; however, comparisons of model predictions with observations do not reveal the expected improvement in the predictions. While high-resolution modeling has scientific merit and potential uses, the currently available monitoring networks, in conjunction with the scarceness of highly resolved spatial input data and the limitations of model formulation, do not allow confirmation of the expected superiority of the high-resolution model predictions.The U.S. Governments right to retain a non-exclusive royalty-free licence in and to any copyright is acknowledged.     

Semicontinuous PM2.5 sulfate and nitrate measurements at an urban and a rural location in New York: PMTACS-NY summer 2001 and 2002 campaigns

Several collocated semicontinuous instruments measuring particulate matter with particle sizes < or =2.5 microm (PM2.5) sulfate (SO4(2-)) and nitrate (NO3-) were intercompared during two intensive field campaigns as part of the PM2.5 Technology Assessment and Characterization Study. The summer 2001 urban campaign in Queens, NY, and the summer 2002 rural campaign in upstate New York (Whiteface Mountain) hosted an operation of an Aerosol Mass Spectrometer, Ambient Particulate Sulfate and Nitrate Monitors, a Continuous Ambient Sulfate Monitor, and a Particle-Into-Liquid Sampler with Ion Chromatographs (PILS-IC). These instruments provided near real-time particulate SO4(2-) and NO3- mass concentration data, allowing the study of particulate SO4(2-)/NO3- diurnal patterns and detection of short-term events. Typical particulate SO4(2-) concentrations were comparable at both sites (ranging from 0 to 20 microg/m3), while ambient urban particulate NO3- concentrations ranged from 0 to 11 microg/m3 and rural NO3- concentration was typically less than 1 microg/m3. Results of the intercomparisons of the semicontinuous measurements are presented, as are results of the comparisons between the semicontinuous and time-integrated filter-based measurements. The comparisons at both sites, in most cases, indicated similar performance characteristics. In addition, charge balance calculations, based on major soluble ionic components of atmospheric aerosol from the PILS-IC and the filter measurements, indicated slightly acidic aerosol at both locations.     

An operational assessment of the application of the relative reduction factors in the demonstration of attainment of the 8-hr ozone National Ambient Air Quality Standard

The U.S. Environmental Protection Agency in 1997 revised the 1-hr ozone (O3) National Ambient Air Quality Standard (NAAQS) to one based on an 8-hr average, resulting in potential nonattainment status for substantial portions of the eastern United States. The regulatory process provides for the development of a state implementation plan that includes a demonstration that the projected future O3 concentrations will be at or below the NAAQS based on photochemical modeling and analytical techniques. In this study, four photochemical modeling systems, based on two photochemical models, Community Model for Air Quality and the Comprehensive Air Quality Model with extensions, and two emissions processing models, Sparse Matrix Optimization Kernel for Emissions and Emissions Modeling System, were applied to the eastern United States, with emphasis on the northeastern Ozone Transport Region in terms of their response to oxides of nitrogen and volatile organic carbon-focused controls on the estimated design values. With the 8-hr O3 NAAQS set as a bright-line test, it was found that a given area could be termed as being in or out of attainment of the NAAQS depending upon the modeling system. This suggests the need to provide an estimate of model-to-model uncertainty in the relative reduction factor (RRF) for a better understanding of the uncertainty in projecting the status of an area's attainment. Results indicate that the model-to-model differences considered in this study introduce     

A combined model–observation approach to estimate historic gridded fields of PM2.5 mass and species concentrations

This paper introduces a methodology for estimating gridded fields of total and speciated fine particulate matter (PM_2.5) concentrations for time periods and regions not covered by observational data. The methodology is based on performing long-term regional scale meteorological and air quality simulations and then integrating these simulations with available observational data. To illustrate this methodology, we present an application in which year-round simulations with a meteorological model (the National Center for Atmospheric Research/Penn State Mesoscale Model, hereafter referred to as MM5) and a photochemical air quality model (the Community Multiscale Air Quality Model, hereafter referred to as CMAQ) have been performed over the northeastern United States for 1988–2005. Model evaluation results for total PM_2.5 mass and individual species for the time period from 2000 to 2005 show that model performance varies by species, season, and location. Therefore, an approach is developed to adjust CMAQ output with factors based on these three variables. The adjusted model values for total PM_2.5 mass for 2000–2005 are compared against independent measurements not utilized for the adjustment approach. This comparison reveals that the adjusted model values have a lower root mean square error (RMSE) and higher correlation coefficients than the original model values. Furthermore, the PM_2.5 estimates from these adjusted model values are compared against an alternate method for estimating historic PM_2.5 values that is based on PM_2.5/PM₁₀ ratios calculated at co-located monitors. Results reveal that both methods yield estimates of historic PM_2.5 mass that are broadly consistent; however, the adjusted CMAQ values provide greater spatial coverage and information for PM_2.5 species in addition to total PM_2.5 mass. Finally, strengths and limitations of the proposed approach are discussed in the context of potential uses of this method.     

Impact of biogenic emission uncertainties on the simulated response of ozone and fine particulate matter to anthropogenic emission reductions

The role of emissions of volatile organic compounds and nitric oxide from biogenic sources is becoming increasingly important in regulatory air quality modeling as levels of anthropogenic emissions continue to decrease and stricter health-based air quality standards are being adopted. However, considerable uncertainties still exist in the current estimation methodologies for biogenic emissions. The impact of these uncertainties on ozone and fine particulate matter (PM2.5) levels for the eastern United States was studied, focusing on biogenic emissions estimates from two commonly used biogenic emission models, the Model of Emissions of Gases and Aerosols from Nature (MEGAN) and the Biogenic Emissions Inventory System (BEIS). Photochemical grid modeling simulations were performed for two scenarios: one reflecting present day conditions and the other reflecting a hypothetical future year with reductions in emissions of anthropogenic oxides of nitrogen (NOx). For ozone, the use of MEGAN emissions resulted in a higher ozone response to hypothetical anthropogenic NOx emission reductions compared with BEIS. Applying the current U.S. Environmental Protection Agency guidance on regulatory air quality modeling in conjunction with typical maximum ozone concentrations, the differences in estimated future year ozone design values (DVF) stemming from differences in biogenic emissions estimates were on the order of 4 parts per billion (ppb), corresponding to approximately 5% of the daily maximum 8-hr ozone National Ambient Air Quality Standard (NAAQS) of 75 ppb. For PM2.5, the differences were 0.1-0.25 microg/m3 in the summer total organic mass component of DVFs, corresponding to approximately 1-2% of the value of the annual PM2.5 NAAQS of 15 microg/m3. Spatial variations in the ozone and PM2.5 differences also reveal that the impacts of different biogenic emission estimates on ozone and PM2.5 levels are dependent on ambient levels of anthropogenic emissions.