Wildfires Increase Concentrations of Hazardous Air Pollutants in Downwind Communities.

Due in part to climate change, wildfire activity is increasing, with the potential for greater public health impact from smoke in downwind communities. Studies examining the health effects of wildfire smoke have focused primarily on fine particulate matter (PM2.5), but there is a need to better characterize other constituents, such as hazardous air pollutants (HAPs). HAPs are chemicals known or suspected to cause cancer or other serious health effects that are regulated by the United States (US) Environmental Protection Agency. Here, we analyzed concentrations of 21 HAPs in wildfire smoke from 2006 to 2020 at 309 monitors across the western US. Additionally, we examined HAP concentrations measured in a major population center (San Jose, CA) affected by multiple fires from 2017 to 2020. We found that concentrations of select HAPs, namely acetaldehyde, acrolein, chloroform, formaldehyde, manganese, and tetrachloroethylene, were all significantly elevated on smoke-impacted versus nonsmoke days (P < 0.05). The largest median increase on smoke-impacted days was observed for formaldehyde, 1.3 µg/m3 (43%) higher than that on nonsmoke days. Acetaldehyde increased 0.73 µg/m3 (36%), and acrolein increased 0.14 µg/m3 (34%). By better characterizing these chemicals in wildfire smoke, we anticipate that this research will aid efforts to reduce exposures in downwind communities.

An Analysis of CMAQ Gas Phase Dry Deposition over North America Through Grid-Scale and Land-Use Specific Diagnostics in the Context of AQMEII4.

The fourth phase of the Air Quality Model Evaluation International Initiative (AQMEII4) is conducting a diagnostic intercomparison and evaluation of deposition simulated by regional-scale air quality models over North America and Europe. In this study, we analyze annual AQMEII4 simulations performed with the Community Multiscale Air Quality Model (CMAQ) version 5.3.1 over North America. These simulations were configured with both the M3Dry and Surface Tiled Aerosol and Gas Exchange (STAGE) dry deposition schemes available in CMAQ. A comparison of observed and modeled concentrations and wet deposition fluxes shows that the AQMEII4 CMAQ simulations perform similarly to other contemporary regional-scale modeling studies. During summer, M3Dry has higher ozone (O3) deposition velocities (Vd) and lower mixing ratios than STAGE for much of the eastern U.S. while the reverse is the case over eastern Canada and along the West Coast. In contrast, during winter STAGE has higher O3 Vd and lower mixing ratios than M3Dry over most of the southern half of the modeling domain while the reverse is the case for much of the northern U.S. and southern Canada. Analysis of the diagnostic variables defined for the AQMEII4 project, i.e. grid-scale and land-use (LU) specific effective conductances and deposition fluxes for the major dry deposition pathways, reveals generally higher summertime stomatal and wintertime cuticular grid-scale effective conductances for M3Dry and generally higher soil grid-scale effective conductances (for both vegetated and bare soil) for STAGE in both summer and winter. On a domain-wide basis, the stomatal grid-scale effective conductances account for about half of the total O3 Vd during daytime hours in summer for both schemes. Employing LU-specific diagnostics, results show that daytime Vd varies by a factor of 2 between LU categories. Furthermore, M3Dry vs. STAGE differences are most pronounced for the stomatal and vegetated soil pathway for the forest LU categories, with M3Dry estimating larger effective conductances for the stomatal pathway and STAGE estimating larger effective conductances for the vegetated soil pathway for these LU categories. Annual domain total O3 deposition fluxes differ only slightly between M3Dry (74.4 Tg/year) and STAGE (76.2 Tg/yr), but pathway-specific fluxes to individual LU types can vary more substantially on both annual and seasonal scales which would affect estimates of O3 damages to sensitive vegetation. A comparison of two simulations differing only in their LU classification scheme shows that the differences in LU cause seasonal mean O3 mixing ratio differences on the order of 1 ppb across large portions of the domain, with the differences generally largest during summer and in areas characterized by the largest differences in the fractional coverages of the forest, planted/cultivated, and grassland LU categories. These differences are generally smaller than the M3Dry vs. STAGE differences outside the summer season but have a similar magnitude during summer. Results indicate that the deposition impacts of LU differences are caused both by differences in the fractional coverages and spatial distributions of different LU categories as well as the characterization of these categories through variables like surface roughness and vegetation fraction in look-up tables used in the land-surface model and deposition schemes. Overall, the analyses and results presented in this study illustrate how the diagnostic grid-scale and LU-specific dry deposition variables adopted for AQMEII4 can provide insights into similarities and differences between the CMAQ M3Dry and STAGE dry deposition schemes that affect simulated pollutant budgets and ecosystem impacts from atmospheric pollution.

Wildfires in the western United States are mobilizing PM2.5-associated nutrients and may be contributing to downwind cyanobacteria blooms.

Wildfire activity is increasing in the continental U.S. and can be linked to climate change effects, including rising temperatures and more frequent drought conditions. Wildfire emissions and large fire frequency have increased in the western U.S., impacting human health and ecosystems. We linked 15 years (2006-2020) of particulate matter (PM2.5) chemical speciation data with smoke plume analysis to identify PM2.5-associated nutrients elevated in air samples on smoke-impacted days. Most macro- and micro-nutrients analyzed (phosphorus, calcium, potassium, sodium, silicon, aluminum, iron, manganese, and magnesium) were significantly elevated on smoke days across all years analyzed. The largest percent increase was observed for phosphorus. With the exception of ammonium, all other nutrients (nitrate, copper, and zinc), although not statistically significant, had higher median values across all years on smoke vs. non-smoke days. Not surprisingly, there was high variation between smoke impacted days, with some nutrients episodically elevated >10 000% during select fire events. Beyond nutrients, we also explored instances where algal blooms occurred in multiple lakes downwind from high-nutrient fires. In these cases, remotely sensed cyanobacteria indices in downwind lakes increased two to seven days following the occurrence of wildfire smoke above the lake. This suggests that elevated nutrients in wildfire smoke may contribute to downwind algal blooms. Since cyanobacteria blooms can be associated with the production of cyanotoxins and wildfire activity is increasing due to climate change, this finding has implications for drinking water reservoirs in the western United States, and for lake ecology, particularly alpine lakes with otherwise limited nutrient inputs.

2002-2017 anthropogenic emissions data for air quality modeling over the United States.

The United States Environmental Protection Agency (US EPA) has developed a set of annual North American emissions data for multiple air pollutants across 18 broad source categories for 2002 through 2017. The sixteen new annual emissions inventories were developed using consistent input data and methods across all years. When a consistent method or tool was not available for a source category, emissions were estimated by scaling data from the EPA's 2017 National Emissions Inventory with scaling factors based on activity data and/or emissions control information. The emissions datasets are designed to support regional air quality modeling for a wide variety of human health and ecological applications. The data were developed to support simulations of the EPA's Community Multiscale Air Quality model but can also be used by other regional scale air quality models. The emissions data are one component of EPA's Air Quality Time Series Project which also includes air quality modeling inputs (meteorology, initial conditions, boundary conditions) and outputs (e.g., ozone, PM2.5 and constituent species, wet and dry deposition) for the Conterminous US at a 12 km horizontal grid spacing.

Beyond Particulate Matter Mass: Heightened Levels of Lead and Other Pollutants Associated with Destructive Fire Events in California.

As the climate warms, wildfire activity is increasing, posing a risk to human health. Studies have reported on particulate matter (PM) in wildfire smoke, yet the chemicals associated with PM have received considerably less attention. Here, we analyzed 13 years (2006-2018) of PM2.5 chemical composition data from monitors in California on smoke-impacted days. Select chemicals (e.g., aluminum and sulfate) were statistically elevated on smoke-impacted days in over half of the years studied. Other chemicals, mostly trace metals harmful to human health (e.g., copper and lead), were elevated during particular fires only. For instance, in 2018, lead was more than 40 times higher on smoke days on average at the Point Reyes monitoring station, due mostly to the Camp Fire, burning approximately 200 km away. There was an association between these metals and the combustion of anthropogenic material (e.g., the burning of houses and vehicles). Although still currently rare, these infrastructure fires are likely becoming more common and can mobilize trace metals in smoke far downwind, at levels generally unseen except in the most polluted areas of the country. We hope a better understanding of the chemicals in wildfire smoke will assist in the communication and reduction of public health risks.

The Community Multiscale Air Quality (CMAQ) model versions 5.3 and 5.3.1: system updates and evaluation.

The Community Multiscale Air Quality (CMAQ) model version 5.3 (CMAQ53), released to the public in August 2019 and followed by version 5.3.1 (CMAQ531) in December 2019, contains numerous science updates, enhanced functionality, and improved computation efficiency relative to the previous version of the model, 5.2.1 (CMAQ521). Major science advances in the new model include a new aerosol module (AERO7) with significant updates to secondary organic aerosol (SOA) chemistry, updated chlorine chemistry, updated detailed bromine and iodine chemistry, updated simple halogen chemistry, the addition of dimethyl sulfide (DMS) chemistry in the CB6r3 chemical mechanism, updated M3Dry bidirectional deposition model, and the new Surface Tiled Aerosol and Gaseous Exchange (STAGE) bidirectional deposition model. In addition, support for the Weather Research and Forecasting (WRF) model's hybrid vertical coordinate (HVC) was added to CMAQ53 and the Meteorology-Chemistry Interface Processor (MCIP) version 5.0 (MCIP50). Enhanced functionality in CMAQ53 includes the new Detailed Emissions Scaling, Isolation and Diagnostic (DESID) system for scaling incoming emissions to CMAQ and reading multiple gridded input emission files. Evaluation of CMAQ531 was performed by comparing monthly and seasonal mean daily 8 h average (MDA8) O3 and daily PM2.5 values from several CMAQ531 simulations to a similarly configured CMAQ521 simulation encompassing 2016. For MDA8 O3, CMAQ531 has higher O3 in the winter versus CMAQ521, due primarily to reduced dry deposition to snow, which strongly reduces wintertime O3 bias (2-4 ppbv monthly average). MDA8 O3 is lower with CMAQ531 throughout the rest of the year, particularly in spring, due in part to reduced O3 from the lateral boundary conditions (BCs), which generally increases MDA8 O3 bias in spring and fall ( 0.5 µg m-3). For daily 24 h average PM2.5, CMAQ531 has lower concentrations on average in spring and fall, higher concentrations in summer, and similar concentrations in winter to CMAQ521, which slightly increases bias in spring and fall and reduces bias in summer. Comparisons were also performed to isolate updates to several specific aspects of the modeling system, namely the lateral BCs, meteorology model version, and the deposition model used. Transitioning from a hemispheric CMAQ (HCMAQ) version 5.2.1 simulation to a HCMAQ version 5.3 simulation to provide lateral BCs contributes to higher O3 mixing ratios in the regional CMAQ simulation in higher latitudes during winter (due to the decreased O3 dry deposition to snow in CMAQ53) and lower O3 mixing ratios in middle and lower latitudes year-round (due to reduced O3 over the ocean with CMAQ53). Transitioning from WRF version 3.8 to WRF version 4.1.1 with the HVC resulted in consistently higher (1.0-1.5 ppbv) MDA8 O3 mixing ratios and higher PM2.5 concentrations (0.1-0.25 µg m-3) throughout the year. Finally, comparisons of the M3Dry and STAGE deposition models showed that MDA8 O3 is generally higher with M3Dry outside of summer, while PM2.5 is consistently higher with STAGE due to differences in the assumptions of particle deposition velocities to non-vegetated surfaces and land use with short vegetation (e.g., grasslands) between the two models. For ambient NH3, STAGE has slightly higher concentrations and smaller bias in the winter, spring, and fall, while M3Dry has higher concentrations and smaller bias but larger error and lower correlation in the summer.

Variability in Observation-based Onroad Emission Constraints from a Near-road Environment.

This study uses Las Vegas near-road measurements of carbon monoxide (CO) and nitrogen oxides (NOx) to test the consistency of onroad emission constraint methodologies. We derive commonly used CO to NOx ratios (ΔCO:ΔNOx) from cross-road gradients and from linear regression using ordinary least squares (OLS) regression and orthogonal regression. The CO to NOx ratios are used to infer NOx emission adjustments for a priori emissions estimates from EPA's MOtor Vehicle Emissions Simulator (MOVES) model assuming unbiased CO. The assumption of unbiased CO emissions may not be appropriate in many circumstances but was implemented in this analysis to illustrate the range of NOx scaling factors that can be inferred based on choice of methods and monitor distance alone. For the nearest road estimates (25m), the cross-road gradient and ordinary least squares (OLS) agree with each other and are not statistically different from the MOVES-based emission estimate while ΔCO:ΔNOx from orthogonal regression is significantly higher than the emitted ratio from MOVES. Using further downwind measurements (i.e., 115m and 300m) increases OLS and orthogonal regression estimates of ΔCO:ΔNOx but not cross-road gradient ΔCO:ΔNOx. The inferred NOx emissions depend on the observation-based method, as well as the distance of the measurements from the roadway and can suggest either that MOVES NOx emissions are unbiased or that they should be adjusted downward by between 10% and 47%. The sensitivity of observation-based ΔCO:ΔNOx estimates to the selected monitor location and to the calculation method characterize the inherent uncertainty of these methods that cannot be derived from traditional standard-error based uncertainty metrics.

Simulating lightning NO production in CMAQv5.2: performance evaluations.

This study assesses the impact of the lightning nitric oxide (LNO) production schemes in the Community Multiscale Air Quality (CMAQ) model on ground-level air quality as well as aloft atmospheric chemistry through detailed evaluation of model predictions of nitrogen oxides (NO x ) and ozone (O3) with corresponding observations for the US. For ground-level evaluations, hourly O3 and NO x values from the U.S. EPA Air Quality System (AQS) monitoring network are used to assess the impact of different LNO schemes on model prediction of these species in time and space. Vertical evaluations are performed using ozonesonde and P-3B aircraft measurements during the Deriving Information on Surface Conditions from Column and Vertically Resolved Observations Relevant to Air Quality (DISCOVER-AQ) campaign conducted in the Baltimore- Washington region during July 2011. The impact on wet deposition of nitrate is assessed using measurements from the National Atmospheric Deposition Program's National Trends Network (NADP NTN). Compared with the Base model (without LNO), the impact of LNO on surface O3 varies from region to region depending on the Base model conditions. Overall statistics suggest that for regions where surface O3 mixing ratios are already overestimated, the incorporation of additional NO from lightning generally increased model overestimation of mean daily maximum 8 h (DM8HR) O3 by 1-2 ppb. In regions where surface O3 is underestimated by the Base model, LNO can significantly reduce the underestimation and bring model predictions close to observations. Analysis of vertical profiles reveals that LNO can significantly improve the vertical structure of modeled O3 distributions by reducing underestimation aloft and to a lesser degree decreasing overestimation near the surface. Since the Base model underestimates the wet deposition of nitrate in most regions across the modeling domain with the exception of the Pacific Coast, the inclusion of LNO leads to reduction in biases and errors and an increase in correlation coefficients at almost all the NADP NTN sites. Among the three LNO schemes described in Kang et al. (2019), the hNLDN scheme, which is implemented using hourly observed lightning flash data from National Lightning Detection Network (NLDN), performs best for comparisons with ground-level values, vertical profiles, and wet deposition of nitrate; the mNLDN scheme (the monthly NLDN-based scheme) performed slightly better. However, when observed lightning flash data are not available, the linear regression-based parameterization scheme, pNLDN, provides an improved estimate for nitrate wet deposition compared to the base simulation that does not include LNO.

A Measurement-Model Fusion Approach for Improved Wet Deposition Maps and Trends.

Air quality models provide spatial fields of wet deposition (WD) and dry deposition that explicitly account for the transport and transformation of emissions from thousands of sources. However, many sources of uncertainty in the air quality model including errors in emissions and meteorological inputs (particularly precipitation) and incomplete descriptions of the chemical and physical processes governing deposition can lead to bias and error in the simulation of WD. We present an approach to bias correct Community Multiscale Air Quality model output over the contiguous United States using observation-based gridded precipitation data generated by the Parameter-elevation Regressions on Independent Slopes Model and WD observations at the National Atmospheric Deposition Program National Trends Network sites. A cross-validation analysis shows that the adjusted annual accumulated WD for NO3 -, NH4 +, and SO4 2- from 2002 to 2012 has less bias and higher correlation with observed values than the base model output without adjustment. Temporal trends in observed WD are captured well by the adjusted model simulations across the entire contiguous United States. Consistent with previous trend analyses, WD NO3 - and SO4 2- are shown to decrease during this period in the eastern half of the United States, particularly in the Northeast, while remaining nearly constant in the West. Trends in WD of NH4 + are more spatially and temporally heterogeneous, with some positive trends in the Great Plains and Central Valley of CA and slightly negative trends in the south.

Simulating lightning NO production in CMAQv5.2: evolution of scientific updates.

This work describes the lightning nitric oxide (LNO) production schemes in the Community Multiscale Air Quality (CMAQ) model. We first document the existing LNO production scheme and vertical distribution algorithm. We then describe updates that were made to the scheme originally based on monthly National Lightning Detection Network (mNLDN) observations. The updated scheme uses hourly NLDN (hNLDN) observations. These NLDN-based schemes are good for retrospective model applications when historical lightning data are available. For applications when observed data are not available (i.e., air quality forecasts and climate studies that assume similar climate conditions), we have developed a scheme that is based on linear and log-linear parameters derived from regression of multiyear historical NLDN (pNLDN) observations and meteorological model simulations. Preliminary assessment for total column LNO production reveals that the mNLDN scheme overestimates LNO by over 40% during summer months compared with the updated hNLDN scheme that reflects the observed lightning activity more faithfully in time and space. The pNLDN performance varies with year, but it generally produced LNO columns that are comparable to hNLDN and mNLDN, and in most cases it outperformed mNLDN. Thus, when no observed lightning data are available, pNLDN can provide reasonable estimates of LNO emissions over time and space for this important natural NO source that influences air quality regulations.

Characterizing CO and NO y Sources and Relative Ambient Ratios in the Baltimore Area Using Ambient Measurements and Source Attribution Modeling.

Modeled source attribution information from the Community Multiscale Air Quality model was coupled with ambient data from the 2011 Deriving Information on Surface conditions from Column and Vertically Resolved Observations Relevant to Air Quality Baltimore field study. We assess source contributions and evaluate the utility of using aircraft measured CO and NO y relationships to constrain emission inventories. We derive ambient and modeled ΔCO:ΔNO y ratios that have previously been interpreted to represent CO:NO y ratios in emissions from local sources. Modeled and measured ΔCO:ΔNO y are similar; however, measured ΔCO:ΔNO y has much more daily variability than modeled values. Sector-based tagging shows that regional transport, on-road gasoline vehicles, and nonroad equipment are the major contributors to modeled CO mixing ratios in the Baltimore area. In addition to those sources, on-road diesel vehicles, soil emissions, and power plants also contribute substantially to modeled NO y in the area. The sector mix is important because emitted CO:NO x ratios vary by several orders of magnitude among the emission sources. The model-predicted gasoline/diesel split remains constant across all measurement locations in this study. Comparison of ΔCO:ΔNO y to emitted CO:NO y is challenged by ambient and modeled evidence that free tropospheric entrainment, and atmospheric processing elevates ambient ΔCO:ΔNO y above emitted ratios. Specifically, modeled ΔCO:ΔNO y from tagged mobile source emissions is enhanced 5-50% above the emitted ratios at times and locations of aircraft measurements. We also find a correlation between ambient formaldehyde concentrations and measured ΔCO:ΔNO y suggesting that secondary CO formation plays a role in these elevated ratios. This analysis suggests that ambient urban daytime ΔCO:ΔNO y values are not reflective of emitted ratios from individual sources.

Characterizing the impact of projected changes in climate and air quality on human exposures to ozone.

The impact of climate change on human and environmental health is of critical concern. Population exposures to air pollutants both indoors and outdoors are influenced by a wide range of air quality, meteorological, behavioral, and housing-related factors, many of which are also impacted by climate change. An integrated methodology for modeling changes in human exposures to tropospheric ozone (O3) owing to potential future changes in climate and demographics was implemented by linking existing modeling tools for climate, weather, air quality, population distribution, and human exposure. Human exposure results from the Air Pollutants Exposure Model (APEX) for 12 US cities show differences in daily maximum 8-h (DM8H) exposure patterns and levels by sex, age, and city for all scenarios. When climate is held constant and population demographics are varied, minimal difference in O3 exposures is predicted even with the most extreme demographic change scenario. In contrast, when population is held constant, we see evidence of substantial changes in O3 exposure for the most extreme change in climate. Similarly, we see increases in the percentage of the population in each city with at least one O3 exposure exceedance above 60 p.p.b and 70 p.p.b thresholds for future changes in climate. For these climate and population scenarios, the impact of projected changes in climate and air quality on human exposure to O3 are much larger than the impacts of changing demographics. These results indicate the potential for future changes in O3 exposure as a result of changes in climate that could impact human health.

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 NO x (NO + NO2), VOC and SO x (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 spectral method for spatial downscaling.

Complex computer models play a crucial role in air quality research. These models are used to evaluate potential regulatory impacts of emission control strategies and to estimate air quality in areas without monitoring data. For both of these purposes, it is important to calibrate model output with monitoring data to adjust for model biases and improve spatial prediction. In this article, we propose a new spectral method to study and exploit complex relationships between model output and monitoring data. Spectral methods allow us to estimate the relationship between model output and monitoring data separately at different spatial scales, and to use model output for prediction only at the appropriate scales. The proposed method is computationally efficient and can be implemented using standard software. We apply the method to compare Community Multiscale Air Quality (CMAQ) model output with ozone measurements in the United States in July 2005. We find that CMAQ captures large-scale spatial trends, but has low correlation with the monitoring data at small spatial scales.

Sensitivity of continental United States atmospheric budgets of oxidized and reduced nitrogen to dry deposition parametrizations.

Reactive nitrogen (Nr) is removed by surface fluxes (air-surface exchange) and wet deposition. The chemistry and physics of the atmosphere result in a complicated system in which competing chemical sources and sinks exist and impact that removal. Therefore, uncertainties are best examined with complete regional chemical transport models that simulate these feedbacks. We analysed several uncertainties in regional air quality model resistance analogue representations of air-surface exchange for unidirectional and bi-directional fluxes and their effect on the continental Nr budget. Model sensitivity tests of key parameters in dry deposition formulations showed that uncertainty estimates of continental total nitrogen deposition are surprisingly small, 5 per cent or less, owing to feedbacks in the chemistry and rebalancing among removal pathways. The largest uncertainties (5%) occur with the change from a unidirectional to a bi-directional NH3 formulation followed by uncertainties in bi-directional compensation points (1-4%) and unidirectional aerodynamic resistance (2%). Uncertainties have a greater effect at the local scale. Between unidirectional and bi-directional formulations, single grid cell changes can be up to 50 per cent, whereas 84 per cent of the cells have changes less than 30 per cent. For uncertainties within either formulation, single grid cell change can be up to 20 per cent, but for 90 per cent of the cells changes are less than 10 per cent.

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.

A comparison of statistical techniques for combining modeled and observed concentrations to create high-resolution ozone air quality surfaces.

Air quality surfaces representing pollutant concentrations across space and time are needed for many applications, including tracking trends and relating air quality to human and ecosystem health. The spatial and temporal characteristics of these surfaces may reveal new information about the associations between emissions, pollution levels, and human exposure and health outcomes that may not have been discernable before. This paper presents four techniques, ranging from simple to complex, to statistically combine observed and modeled daily maximum 8-hr ozone concentrations for a domain covering the greater New York State area for the summer of 2001. Cross-validation results indicate that, for the domain and time period studied, the simpler techniques (additive and multiplicative bias adjustment) perform as well as or better than the more complex techniques. However, the spatial analyses of the resulting ozone concentration surfaces revealed some problems with these simpler techniques in limited areas where the model exhibits difficulty in simulating the complex features such as those observed in the New York City area.

Combat-related guilt mediates the relations between exposure to combat-related abusive violence and psychiatric diagnoses.

BACKGROUND: This study examined the degree to which combat-related guilt mediated the relations between exposure to combat-related abusive violence and both Posttraumatic Stress Disorder (PTSD) and Major Depressive Disorder (MDD) in Vietnam Veterans. METHODS: Secondary analyses were conducted on data collected from 1,323 male Vietnam Veterans as part of a larger, multisite study. RESULTS: Results revealed that combat-related guilt partially mediated the association between exposure to combat-related abusive violence and PTSD, but completely mediated the association with MDD, with overall combat exposure held constant in the model. Follow-up analyses showed that, when comparing those participants who actually participated in combat-related abusive violence with those who only observed it, combat-related guilt completely mediated the association between participation in abusive violence and both PTSD and MDD. Moreover, when comparing those participants who observed combat-related abusive violence with those who had no exposure at all to it, combat-related guilt completely mediated the association between observation of combat-related abusive violence and MDD, but only partially mediated the association with PTSD. CONCLUSIONS: These findings suggest that guilt may be a mechanism through which abusive violence is related to PTSD and MDD among combat-deployed Veterans. These findings also suggest the importance of assessing abusive-violence related guilt among combat-deployed Veterans and implementing relevant interventions for such guilt whenever indicated.

Efficient probabilistic estimates of surface ozone concentration using an ensemble of model configurations and direct sensitivity calculations.

Because all models are a simplification of the phenomenon they aim to represent, it is often more useful to estimate the probability of an event rather than a single "best" model result. Previous air quality ensemble approaches have used computationally expensive simulations of separately developed modeling systems. We present an efficient method to generate ensembles with hundreds of members based on several structural configurations of a single air quality modeling system. We use the Decoupled Direct Method in three dimensions to directly calculate how ozone concentrations change as a result of changes in input parameters. The modeled probability estimate is compared to observations and is shown to have a high level of skill and improved resolution and sharpness. This approach can help resolve the practical limits of incorporating uncertainty estimation into deterministic air quality management modeling applications.