Sergey L. Napelenok

Model Representation of Secondary Organic Aerosol in CMAQv4.7

Numerous scientific upgrades to the representation of secondary organic aerosol (SOA) are incorporated into the Community Multiscale Air Quality (CMAQ) modeling system. Additions include several recently identified SOA precursors: benzene, isoprene, and sesquiterpenes; and pathways: in-cloud oxidation of glyoxal and methylglyoxal, particle-phase oligomerization, and acid enhancement of isoprene SOA. NO(x)-dependent aromatic SOA yields are also added along with new empirical measurements of the enthalpies of vaporization and organic mass-to-carbon ratios. For the first time, these SOA precursors, pathways and empirical parameters are included simultaneously in an air quality model for an annual simulation spanning the continental U.S. Comparisons of CMAQ-modeled secondary organic carbon (OC(sec)) with semiempirical estimates screened from 165 routine monitoring sites across the U.S. indicate the new SOA module substantially improves model performance. The most notable improvement occurs in the central and southeastern U.S. where the regionally averaged temporal correlations (r) between modeled and semiempirical OC(sec) increase from 0.5 to 0.8 and 0.3 to 0.8, respectively, when the new SOA module is employed. Wintertime OC(sec) results improve in all regions of the continental U.S. and the seasonal and regional patterns of biogenic SOA are better represented.

Description and evaluation of the Community Multiscale Air Quality (CMAQ) modeling system version 5.1

The Community Multiscale Air Quality (CMAQ) model is a comprehensive multipollutant air quality modeling system developed and maintained by the US Environmental Protection Agency’s (EPA) Office of Research and Development (ORD). Recently, version 5.1 of the CMAQ model (v5.1) was released to the public, incorporating a large number of science updates and extended capabilities over the previous release version of the model (v5.0.2). These updates include the following: improvements in the meteorological calculations in both CMAQ and the Weather Research and Forecast (WRF) model used to provide meteorological fields to CMAQ, updates to the gas and aerosol chemistry, revisions to the calculations of clouds and photolysis, and improvements to the dry and wet deposition in the model. Sensitivity simulations isolating several of the major updates to the modeling system show that changes to the meteorological calculations result in enhanced afternoon and early evening mixing in the model, periods when the model historically underestimates mixing. This enhanced mixing results in higher ozone (O3) mixing ratios on average due to reduced NO titration, and lower fine particulate matter (PM2.5) concentrations due to greater dilution of primary pollutants (e.g., elemental and organic carbon). Updates to the clouds and photolysis calculations greatly improve consistency between the WRF and CMAQ models and result in generally higher O3 mixing ratios, primarily due to reduced cloudiness and attenuation of photolysis in the model. Updates to the aerosol chemistry result in higher secondary organic aerosol (SOA) concentrations in the summer, thereby reducing summertime PM2.5 bias (PM2.5 is typically underestimated by CMAQ in the summer), while updates to the gas chemistry result in slightly higher O3 and PM2.5 on average in January and July. Overall, the seasonal variation in simulated PM2.5 generally improves in CMAQv5.1 (when considering all model updates), as simulated PM2.5 concentrations decrease in the winter (when PM2.5 is generally overestimated by CMAQ) and increase in the summer (when PM2.5 is generally underestimated by CMAQ). Ozone mixing ratios are higher on average with v5.1 vs. v5.0.2, resulting in higher O3 mean bias, as O3 tends to be overestimated by CMAQ throughout most of the year (especially at locations where the observed O3 is low); however, O3 correlation is largely improved with v5.1. Sensitivity simulations for several hypothetical emission reduction scenarios show that v5.1 tends to be slightly more responsive to reductions in NOx (NO + NO2), VOC and SOx (SO2 + SO4) emissions than v5.0.2, representing an improvement as previous studies have shown CMAQ to underestimate the observed reduction in O3 due to large, widespread reductions in observed emissions.

A Direct Sensitivity Approach to Predict Hourly Ozone Resulting from Compliance with the National Ambient Air Quality Standard

In setting primary ambient air quality standards, the EPAs responsibility under the law is to establish standards that protect public health. As part of the current review of the ozone National Ambient Air Quality Standard (NAAQS), the US EPA evaluated the health exposure and risks associated with ambient ozone pollution using a statistical approach to adjust recent air quality to simulate just meeting the current standard level, without specifying emission control strategies. One drawback of this purely statistical concentration rollback approach is that it does not take into account spatial and temporal heterogeneity of ozone response to emissions changes. The application of the higher-order decoupled direct method (HDDM) in the community multiscale air quality (CMAQ) model is discussed here to provide an example of a methodology that could incorporate this variability into the risk assessment analyses. Because this approach includes a full representation of the chemical production and physical transport of ozone in the atmosphere, it does not require assumed background concentrations, which have been applied to constrain estimates from past statistical techniques. The CMAQ-HDDM adjustment approach is extended to measured ozone concentrations by determining typical sensitivities at each monitor location and hour of the day based on a linear relationship between first-order sensitivities and hourly ozone values. This approach is demonstrated by modeling ozone responses for monitor locations in Detroit and Charlotte to domain-wide reductions in anthropogenic NOx and VOCs emissions. As seen in previous studies, ozone response calculated using HDDM compared well to brute-force emissions changes up to approximately a 50% reduction in emissions. A new stepwise approach is developed here to apply this method to emissions reductions beyond 50% allowing for the simulation of more stringent reductions in ozone concentrations. Compared to previous rollback methods, this application of modeled sensitivities to ambient ozone concentrations provides a more realistic spatial response of ozone concentrations at monitors inside and outside the urban core and at hours of both high and low ozone concentrations.

Uncertainty Analysis of Ozone Formation and Response to Emission Controls Using Higher-Order Sensitivities

Abstract Understanding ozone response to its precursor emissions is crucial for effective air quality management practices. This nonlinear response is usually simulated using chemical transport models, and the modeling results are affected by uncertainties in emissions inputs. In this study, a high ozone episode in the southeastern United States is simulated using the Community Multiscale Air Quality (CMAQ) model. Uncertainties in ozone formation and response to emissions controls due to uncertainties in emission rates are quantified using the Monte Carlo method. Instead of propagating emissions uncertainties through the original CMAQ, a reduced form of CMAQ is formulated using directly calculated first- and second-order sensitivities that capture the nonlinear ozone concentration-emission responses. This modification greatly reduces the associated computational cost. Quantified uncertainties in modeled ozone concentrations and responses to various emissions controls are much less than the uncertainties in emissions inputs. Average uncertainties in modeled ozone concentrations for the Atlanta area are less than 10% (as measured by the inferred coefficient of variance [ICOV]) even when emissions uncertainties are assumed to vary between a factor of 1.5 and 2. Uncertainties in the ozone responses generally decrease with increased emission controls. Average uncertainties (ICOV) in emission-normalized ozone responses range from 4 to 22%, with the smaller being associated with controlling of the relatively certain point nitrogen oxide (NOx) emissions and the larger resulting from controlling of the less certain mobile NOx emissions. These small uncertainties provide confidence in the model applications, such as in performance evaluation, attainment demonstration, and control strategy development.

Extreme value analysis for evaluating ozone control strategies

Tropospheric ozone is one of six criteria pollutants regulated by the US EPA, and has been linked to respiratory and cardiovascular endpoints and adverse effects on vegetation and ecosystems. Regional photochemical models have been developed to study the impacts of emission reductions on ozone levels. The standard approach is to run the deterministic model under new emission levels and attribute the change in ozone concentration to the emission control strategy. However, running the deterministic model requires substantial computing time, and this approach does not provide a measure of uncertainty for the change in ozone levels. Recently, a reduced form model (RFM) has been proposed to approximate the complex model as a simple function of a few relevant inputs. In this paper, we develop a new statistical approach to make full use of the RFM to study the effects of various control strategies on the probability and magnitude of extreme ozone events. We fuse the model output with monitoring data to calibrate the RFM by modeling the conditional distribution of monitoring data given the RFM using a combination of flexible semiparametric quantile regression for the center of the distribution where data are abundant and a parametric extreme value distribution for the tail where data are sparse. Selected parameters in the conditional distribution are allowed to vary by the RFM value and the spatial location. Also, due to the simplicity of the RFM, we are able to embed the RFM in our Bayesian hierarchical framework to obtain a full posterior for the model input parameters, and propagate this uncertainty to the estimation of the effects of the control strategies. We use the new framework to evaluate three potential control strategies, and find that reducing mobile-source emissions has a larger impact than reducing point-source emissions or a combination of several emission sources.

Differences Between Magnitudes and Health Impacts of BC Emissions Across the United States Using 12 km Scale Seasonal Source Apportionment

Recent assessments have analyzed the health impacts of PM2.5 from emissions from different locations and sectors using simplified or reduced-form air quality models. Here we present an alternative approach using the adjoint of the Community Multiscale Air Quality (CMAQ) model, which provides source-receptor relationships at highly resolved sectoral, spatial, and temporal scales. While damage resulting from anthropogenic emissions of BC is strongly correlated with population and premature death, we found little correlation between damage and emission magnitude, suggesting that controls on the largest emissions may not be the most efficient means of reducing damage resulting from anthropogenic BC emissions. Rather, the best proxy for locations with damaging BC emissions is locations where premature deaths occur. Onroad diesel and nonroad vehicle emissions are the largest contributors to premature deaths attributed to exposure to BC, while onroad gasoline emissions cause the highest deaths per amount emitted. Emissions in fall and winter contribute to more premature deaths (and more per amount emitted) than emissions in spring and summer. Overall, these results show the value of the high-resolution source attribution for determining the locations, seasons, and sectors for which BC emission controls have the most effective health benefits.

Development and evaluation of a physics‐based windblown dust emission scheme implemented in the CMAQ modeling system

A new windblown dust emission treatment was incorporated in the Community Multiscale Air Quality (CMAQ) modeling system. This new model treatment has been built upon previously developed physics-based parameterization schemes from the literature. A distinct and novel feature of this scheme, however, is the incorporation of a newly developed dynamic relation for the surface roughness length relevant to small-scale dust generation processes. Through this implementation, the effect of nonerodible elements on the local flow acceleration, drag partitioning, and surface coverage protection is modeled in a physically based and consistent manner. Careful attention is paid in integrating the new windblown dust treatment in the CMAQ model to ensure that the required input parameters are correctly configured. To test the performance of the new dust module in CMAQ, the entire year 2011 is simulated for the continental United States, with particular emphasis on the southwestern United States (SWUS) where windblown dust concentrations are relatively large. Overall, the model shows good performance with the daily mean bias of soil concentrations fluctuating in the range of ±1 μg m−3 for the entire year. Springtime soil concentrations are in quite good agreement (normalized mean bias of 8.3%) with observations, while moderate to high underestimation of soil concentration is seen in the summertime. The latter is attributed to the issue of representing the convective dust storms in summertime. Evaluations against observations for seven elevated dust events in the SWUS indicate that the new windblown dust treatment is capable of capturing spatial and temporal characteristics of dust outbreaks.

Premature deaths attributed to source-specific BC emissions in six urban US regions

Recent studies have shown that exposure to particulate black carbon (BC) has significant adverse health effects and may be more detrimental to human health than exposure to PM2.5 as a whole. Mobile source BC emission controls, mostly on diesel-burning vehicles, have successfully decreased mobile source BC emissions to less than half of what they were 30 years ago. Quantification of the benefits of previous emissions controls conveys the value of these regulatory actions and provides a method by which future control alternatives could be evaluated. In this study we use the adjoint of the Community Multiscale Air Quality (CMAQ) model to estimate highly-resolved spatial distributions of benefits related to emission reductions for six urban regions within the continental US. Emissions from outside each of the six chosen regions account for between 7% and 27% of the premature deaths attributed to exposure to BC within the region. While we estimate that nonroad mobile and onroad diesel emissions account for the largest number of premature deaths attributable to exposure to BC, onroad gasoline is shown to have more than double the benefit per unit emission relative to that of nonroad mobile and onroad diesel. Within the region encompassing New York City and Philadelphia, reductions in emissions from large industrial combustion sources that are not classified as EGUs (i.e., non-EGU) are estimated to have up to triple the benefits per unit emission relative to reductions to onroad diesel sectors, and provide similar benefits per unit emission to that of onroad gasoline emissions in the region. While onroad mobile emissions have been decreasing in the past 30 years and a majority of vehicle emission controls that regulate PM focus on diesel emissions, our analysis shows the most efficient target for stricter controls is actually onroad gasoline emissions.

Diagnostic Air Quality Model Evaluation of Source-Specific Primary and Secondary Fine Particulate Carbon

Ambient measurements of 78 source-specific tracers of primary and secondary carbonaceous fine particulate matter collected at four midwestern United States locations over a full year (March 2004-February 2005) provided an unprecedented opportunity to diagnostically evaluate the results of a numerical air quality model. Previous analyses of these measurements demonstrated excellent mass closure for the variety of contributing sources. In this study, a carbon-apportionment version of the Community Multiscale Air Quality (CMAQ) model was used to track primary organic and elemental carbon emissions from 15 independent sources such as mobile sources and biomass burning in addition to four precursor-specific classes of secondary organic aerosol (SOA) originating from isoprene, terpenes, aromatics, and sesquiterpenes. Conversion of the source-resolved model output into organic tracer concentrations yielded a total of 2416 data pairs for comparison with observations. While emission source contributions to the total model bias varied by season and measurement location, the largest absolute bias of -0.55 μgC/m(3) was attributed to insufficient isoprene SOA in the summertime CMAQ simulation. Biomass combustion was responsible for the second largest summertime model bias (-0.46 μgC/m(3) on average). Several instances of compensating errors were also evident; model underpredictions in some sectors were masked by overpredictions in others.

Bayesian analysis of a reduced-form air quality model.

Numerical air quality models are being used for assessing emission control strategies for improving ambient pollution levels across the globe. This paper applies probabilistic modeling to evaluate the effectiveness of emission reduction scenarios aimed at lowering ground-level ozone concentrations. A Bayesian hierarchical model is used to combine air quality model output and monitoring data in order to characterize the impact of emissions reductions while accounting for different degrees of uncertainty in the modeled emissions inputs. The probabilistic model predictions are weighted based on population density in order to better quantify the societal benefits/disbenefits of four hypothetical emission reduction scenarios in which domain-wide NO(x) emissions from various sectors are reduced individually and then simultaneously. Cross validation analysis shows the statistical model performs well compared to observed ozone levels. Accounting for the variability and uncertainty in the emissions and atmospheric systems being modeled is shown to impact how emission reduction scenarios would be ranked, compared to standard methodology.