Effects of Air Pollutants from Wildfires on Downwind Ecosystems: Observations, Knowledge Gaps, and Questions for Assessing Risk

Wildfires have increased in frequency and area burned, trends expected to continue with climate change. Among other effects, fires release pollutants into the atmosphere, representing a risk to human health and downwind terrestrial and aquatic ecosystems. While human health risks are well studied, the ecological impacts to downwind ecosystems are not, and this gap may present a constraint on developing an adequate assessment of the ecological risks associated with downwind wildfire exposure. Here, we first screened the scientific literature to assess general knowledge about pathways and end points of a conceptual model linking wildfire generated pollutants and other materials to downwind ecosystems. We found a substantial body of literature on the composition of wildfire derived pollution and materials in the atmosphere and subsequent transport, yet little observational or experimental work on their effects on downwind ecological end points. This dearth of information raises many questions related to adequately assessing the ecological risk of downwind exposure, especially given increasing wildfire trends. To guide future research, we pose eight questions within the well-established US EPA ecological risk assessment paradigm that if answered would greatly improve ecological risk assessment and, ultimately, management strategies needed to reduce potential wildfire impacts.

Modeling past and future spatiotemporal distributions of airborne allergenic pollen across the contiguous United States

Exposures to airborne allergenic pollen have been increasing under the influence of changing climate. A modeling system incorporating pollen emissions and atmospheric transport and fate processes has been developed and applied to simulate spatiotemporal distributions of two major aeroallergens, oak and ragweed pollens, across the contiguous United States (CONUS) for both historical (year 2004) and future (year 2047) conditions. The transport and fate of pollen presented here is simulated using our adapted version of the Community Multiscale Air Quality (CMAQ) model. Model performance was evaluated using observed pollen counts at monitor stations across the CONUS for 2004. Our analysis shows that there is encouraging consistency between observed seasonal mean concentrations and corresponding simulated seasonal mean concentrations (oak: Pearson = 0.35, ragweed: Pearson = 0.40), and that the model was able to capture the statistical patterns of observed pollen concentration distributions in 2004 for most of the pollen monitoring stations. Simulation of pollen levels for a future year (2047) considered conditions corresponding to the RCP8.5 scenario. Modeling results show substantial regional variability both in the magnitude and directionality of changes in pollen metrics. Ragweed pollen season is estimated to start earlier and last longer for all nine climate regions of the CONUS, with increasing average pollen concentrations in most regions. The timing and magnitude of oak pollen season vary across the nine climate regions, with the largest increases in pollen concentrations expected in the Northeast region.

Flexible Bayesian Ensemble Machine Learning Framework for Predicting Local Ozone Concentrations

3D-grid-based chemical transport models, such as the Community Multiscale Air Quality (CMAQ) modeling system, have been widely used for predicting concentrations of ambient air pollutants. However, typical horizontal resolutions of nationwide CMAQ simulations (12 × 12 km2) cannot capture local-scale gradients for accurately assessing human exposures and environmental justice disparities. In this study, a Bayesian ensemble machine learning (BEML) framework, which integrates 13 learning algorithms, was developed for downscaling CMAQ estimates of ozone daily maximum 8 h averages to the census tract level, across the contiguous US, and was demonstrated for 2011. Three-stage hyperparameter tuning and targeted validations were designed to ensure the ensemble model's ability to interpolate, extrapolate, and capture concentration peaks. The Shapley value metric from coalitional game theory was applied to interpret the drivers of subgrid gradients. The flexibility (transferability) of the 2011-trained BEML model was further tested by evaluating its ability to estimate fine-scale concentrations for other years (2012-2017) without retraining. To demonstrate the feasibility of using the BEML approach to strictly "data-limited" situations, the model was applied to downscale CMAQ outputs for a future-year scenario-based simulation that considers effects of variations in meteorology associated with climate change.

The Role of Temperature in Modifying the Risk of Ozone-Attributable Mortality under Future Changes in Climate: A Proof-of-Concept Analysis

Air pollution risk assessments typically estimate ozone-attributable mortality counts using concentration-response (C-R) parameters from epidemiologic studies that treat temperature as a potential confounder. However, some recent epidemiologic studies have indicated that temperature can modify the relationship between short-term ozone exposure and mortality, which has potentially important implications when considering the impacts of climate change on public health. This proof-of-concept analysis quantifies counts of temperature-modified ozone-attributable mortality using temperature-stratified C-R parameters from a multicity study in which the pooled ozone-mortality effect coefficients change in concert with daily temperature. Meteorology downscaled from two global climate models is used with a photochemical transport model to simulate ozone concentrations over the 21st century using two emission inventories: one holding air pollutant emissions constant at 2011 levels and another accounting for reduced emissions through the year 2040. The late century climate models project increased summer season temperatures, which in turn yields larger total counts of ozone-attributable deaths in analyses using temperature-stratified C-R parameters compared to the traditional temperature confounder approach. This analysis reveals substantial heterogeneity in the magnitude and distribution of the temperature-stratified ozone-attributable mortality results, which is a function of regional variability in both the C-R relationship and the model-predicted temperature and ozone.

Modeling future asthma attributable to fine particulate matter (PM2.5) in a changing climate: a health impact assessment

Exposure to fine particulate matter (PM_2.5) is associated with asthma development as well as asthma exacerbation in children. PM_2.5 can be directly emitted or can form in the atmosphere from pollutant precursors. PM_2.5 emitted and formed in the atmosphere is influenced by meteorology; future changes in climate may alter the concentration and distribution of PM_2.5. Our aim is to estimate the future burden of climate change and PM_2.5 on new and exacerbated cases of childhood asthma. Projected concentrations of PM_2.5 are based on the Geophysical Fluid Dynamics Laboratory Coupled Model version 3 climate model, the Representative Concentration Pathway 8.5 greenhouse gas scenario, and two air pollution emissions datasets: a 2011 emissions dataset and a 2040 emissions dataset that reflects substantial reductions in emissions of PM_2.5 as compared to the 2011 inventory. We estimate additional PM_2.5-attributable asthma as well as PM_2.5-attributable albuterol inhaler use for four future years (2030, 2050, 2075, and 2095) relative to the year 2000. Exacerbations, regardless of the trigger, are counted as attributable to PM_2.5 if the incident disease is attributable to PM_2.5. We project 38 thousand (95% CI 36, 39 thousand) additional PM_2.5-attributable incident childhood asthma cases and 29 million (95% CI 27, 31 million) additional PM_2.5-attributable albuterol inhaler uses per year in 2030, increasing to 200 thousand (95% CI 190, 210 thousand) additional incident cases and 160 million (95% CI 150, 160 million) inhaler uses per year by 2095 relative to 2000 under the 2011 emissions dataset. These additional PM_2.5-attributable incident asthma cases and albuterol inhaler use would cost billions of additional U.S. dollars per year by the late century. These outcomes could be mitigated by reducing air pollution emissions.

Evaluating long-term emission impacts of large-scale electric vehicle deployment in the US using a human-Earth systems model

While large-scale adoption of electric vehicles (EVs) globally would reduce carbon dioxide (CO2) and traditional air pollutant emissions from the transportation sector, emissions from the electric sector, refineries, and potentially other sources would change in response. Here, a multi-sector human-Earth systems model is used to evaluate the net long-term emission implications of large-scale EV adoption in the US over widely differing pathways of the evolution of the electric sector. Our results indicate that high EV adoption would decrease net CO2 emissions through 2050, even for a scenario where all electric sector capacity additions through 2050 are fossil fuel technologies. Greater net CO2 reductions would be realized for scenarios that emphasize renewables or decarbonization of electricity production. Net air pollutant emission changes in 2050 are relatively small compared to expected overall decreases from recent levels to 2050. States participating in the Regional Greenhouse Gas Initiative experience greater CO2 and air pollutant reductions on a percentage basis. These results suggest that coordinated, multi-sector planning can greatly enhance the climate and environmental benefits of EVs. Additional factors are identified that influence the net emission impacts of EVs, including the retirement of coal capacity, refinery operations under reduced gasoline demands, and price-induced fuel switching in residential heating and in the industrial sector.

Regional temperature-ozone relationships across the U.S. under multiple climate and emissions scenarios

The potential effects of 21st century climate change on ozone (O₃) concentrations in the United States are investigated using global climate simulations to drive higher-resolution regional meteorological and chemical transport models. Community Earth System Model (CESM) and Coupled Model version 3 (CM3) simulations of the Representative Concentration Pathway 8.5 scenario are dynamically downscaled using the Weather Research and Forecasting model, and the resulting meteorological fields are used to drive the Community Multiscale Air Quality model. Air quality is modeled for five 11-year periods using both a 2011 air pollutant emission inventory and a future projection accounting for full implementation of promulgated regulatory controls. Across the U.S., CESM projects daily maximum temperatures during summer to increase 1-4°C by 2050 and 2-7°C by 2095, while CM3 projects warming of 2-7°C by 2050 and 4-11°C by 2095. The meteorological changes have geographically varying impacts on O₃ concentrations. Using the 2011 emissions dataset, O₃ increases 1-5 ppb in the central Great Plains and Midwest by 2050 and more than 10 ppb by 2095, but it remains unchanged or even decreases in the Gulf Coast, Maine, and parts of the Southwest. Using the projected emissions, modeled increases are attenuated while decreases are amplified, indicating that planned air pollution control measures ameliorate the ozone climate penalty. The relationships between changes in maximum temperature and changes in O₃ concentrations are examined spatially and quantified to explore the potential for developing an efficient approach for estimating air quality impacts of other future climate scenarios.Implications: The effects of climate change on ozone air quality in the United States are investigated using two global climate model simulations of a high warming scenario for five decadal periods in the 21st century. Warming summer temperatures simulated under both models lead to higher ozone concentrations in some regions, with the magnitude of the change increasing with temperature over the century. The magnitude and spatial extent of the increases are attenuated under a future emissions projection that accounts for regulatory controls. Regional linear regression relationships are developed as a first step toward development of a reduced form model for efficient estimation of the health impacts attributable to changes in air quality resulting from a climate change scenario.

The Detailed Emissions Scaling, Isolation, and Diagnostic (DESID) module in the Community Multiscale Air Quality (CMAQ) modeling system version 5.3.2

Air quality modeling for research and regulatory applications often involves executing many emissions sensitivity cases to quantify impacts of hypothetical scenarios, estimate source contributions, or quantify uncertainties. Despite the prevalence of this task, conventional approaches for perturbing emissions in chemical transport models like the Community Multiscale Air Quality (CMAQ) model require extensive offline creation and finalization of alternative emissions input files. This workflow is often time-consuming, error-prone, inconsistent among model users, difficult to document, and dependent on increased hard disk resources. The Detailed Emissions Scaling, Isolation, and Diagnostic (DESID) module, a component of CMAQv5.3 and beyond, addresses these limitations by performing these modifications online during the air quality simulation. Further, the model contains an Emission Control Interface which allows users to prescribe both simple and highly complex emissions scaling operations with control over individual or multiple chemical species, emissions sources, and spatial areas of interest. DESID further enhances the transparency of its operations with extensive error-checking and optional gridded output of processed emission fields. These new features are of high value to many air quality applications including routine perturbation studies, atmospheric chemistry research, and coupling with external models (e.g., energy system models, reduced-form models).

The Detailed Emissions Scaling, Isolation, and Diagnostic (DESID) module in the Community Multiscale Air Quality (CMAQ) modeling system version 5.3.2

Air quality modeling for research and regulatory applications often involves executing many emissions sensitivity cases to quantify impacts of hypothetical scenarios, estimate source contributions, or quantify uncertainties. Despite the prevalence of this task, conventional approaches for perturbing emissions in chemical transport models like the Community Multiscale Air Quality (CMAQ) model require extensive offline creation and finalization of alternative emissions input files. This workflow is often time-consuming, error-prone, inconsistent among model users, difficult to document, and dependent on increased hard disk resources. The Detailed Emissions Scaling, Isolation, and Diagnostic (DESID) module, a component of CMAQv5.3 and beyond, addresses these limitations by performing these modifications online during the air quality simulation. Further, the model contains an Emission Control Interface which allows users to prescribe both simple and highly complex emissions scaling operations with control over individual or multiple chemical species, emissions sources, and spatial areas of interest. DESID further enhances the transparency of its operations with extensive error-checking and optional gridded output of processed emission fields. These new features are of high value to many air quality applications including routine perturbation studies, atmospheric chemistry research, and coupling with external models (e.g., energy system models, reduced-form models).

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) O₃ and daily PM_2.5 values from several CMAQ531 simulations to a similarly configured CMAQ521 simulation encompassing 2016. For MDA8 O₃, CMAQ531 has higher O₃ in the winter versus CMAQ521, due primarily to reduced dry deposition to snow, which strongly reduces wintertime O₃ bias (2-4 ppbv monthly average). MDA8 O₃ is lower with CMAQ531 throughout the rest of the year, particularly in spring, due in part to reduced O₃ from the lateral boundary conditions (BCs), which generally increases MDA8 O₃ bias in spring and fall ( 0.5 μg m^-3). For daily 24 h average PM_2.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 O₃ mixing ratios in the regional CMAQ simulation in higher latitudes during winter (due to the decreased O₃ dry deposition to snow in CMAQ53) and lower O₃ mixing ratios in middle and lower latitudes year-round (due to reduced O₃ 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 O₃ mixing ratios and higher PM_2.5 concentrations (0.1-0.25 μg m^-3) throughout the year. Finally, comparisons of the M3Dry and STAGE deposition models showed that MDA8 O₃ is generally higher with M3Dry outside of summer, while PM_2.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 NH₃, 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.

The contribution of wildland fire emissions to deposition in the U S: implications for tree growth and survival in the Northwest

Ecosystems require access to key nutrients like nitrogen (N) and sulfur (S) to sustain growth and healthy function. However, excessive deposition can also damage ecosystems through nutrient imbalances, leading to changes in productivity and shifts in ecosystem structure. While wildland fires are a known source of atmospheric N and S, little has been done to examine the implications of wildland fire deposition for vulnerable ecosystems. We combine wildland fire emission estimates, atmospheric chemistry modeling, and forest inventory data to (a) quantify the contribution of wildland fire emissions to N and S deposition across the U S, and (b) assess the subsequent impacts on tree growth and survival rates in areas where impacts are likely meaningful based on the relative contribution of fire to total deposition. We estimate that wildland fires contributed 0.2 kg N ha-1 yr-1 and 0.04 kg S ha-1 yr-1 on average across the U S during 2008-2012, with maxima up to 1.4 kg N ha-1 yr-1 and 0.6 kg S ha-1 yr-1 in the Northwest representing over ~30% of total deposition in some areas. Based on these fluxes, exceedances of S critical loads as a result of wildland fires are minimal, but exceedances for N may affect the survival and growth rates of 16 tree species across 4.2 million hectares, with the most concentrated impacts occurring in Oregon, northern California, and Idaho. Understanding the broader environmental impacts of wildland fires in the U S will inform future decision making related to both fire management and ecosystem services conservation.

Associations Between Simulated Future Changes in Climate, Air Quality, and Human Health

IMPORTANCE

Future changes in climate are likely to adversely affect human health by affecting concentrations of particulate matter sized less than 2.5 μm (PM2.5) and ozone (O3) in many areas. However, the degree to which these outcomes may be mitigated by reducing air pollutant emissions is not well understood.

OBJECTIVE

To model the associations between future changes in climate, air quality, and human health for 2 climate models and under 2 air pollutant emission scenarios.

DESIGN, SETTING, AND PARTICIPANTS

This modeling study simulated meteorological conditions over the coterminous continental US during a 1995 to 2005 baseline and over the 21st century (2025-2100) by dynamically downscaling representations of a high warming scenario from the Community Earth System Model (CESM) and the Coupled Model version 3 (CM3) global climate models. Using a chemical transport model, PM2.5 and O3 concentrations were simulated under a 2011 air pollutant emission data set and a 2040 projection. The changes in PM2.5 and O3-attributable deaths associated with climate change among the US census-projected population were estimated for 2030, 2050, 2075, and 2095 for each of 2 emission inventories and climate models. Data were analyzed from June 2018 to June 2020.

MAIN OUTCOMES AND MEASURES

The main outcomes were simulated change in summer season means of the maximum daily 8-hour mean O3, annual mean PM2.5, population-weighted exposure, and the number of avoided or incurred deaths associated with these pollutants. Results are reported for 2030, 2050, 2075, and 2095, compared with 2000, for 2 climate models and 2 air pollutant emissions data sets.

RESULTS

The projected increased maximum daily temperatures through 2095 were up to 7.6 °C for the CESM model and 11.8 °C for the CM3 model. Under each climate model scenario by 2095, compared with 2000, an estimated additional 21 000 (95% CI, 14 000-28 000) PM2.5-attributable deaths and 4100 (95% CI, 2200-6000) O3-attributable deaths were projected to occur. These projections decreased to an estimated 15 000 (95% CI, 10 000-20 000) PM2.5-attributable deaths and 640 (95% CI, 340-940) O3-attributable deaths when simulated using a future emission inventory that accounted for reduced anthropogenic emissions.

CONCLUSIONS AND RELEVANCE

These findings suggest that reducing future air pollutant emissions could also reduce the climate-driven increase in deaths associated with air pollution by hundreds to thousands.

Ozone-related asthma emergency department visits in the US in a warming climate

Ozone exposure is associated with higher risk of asthma-related emergency department visits. The meteorological conditions that govern ozone concentration are projected to be more favorable to ozone formation over much of the United States due to continued climate change, even as emissions of anthropogenic ozone precursors are expected to decrease by 2050. Our goal is to quantify the health benefits of a climate change mitigation scenario versus a "business-as-usual" scenario, defined by the United Nations Intergovernmental Panel on Climate Change Representative Concentration Pathways (RCPs) 4.5 and 8.5, respectively, using the health impact analytical program Benefits Mapping and Analysis Program - Community Edition (BenMAP - CE) to project the number of asthma ED visits in 2045-2055. We project an annual average of 3100 averted ozone-related asthma ED visits during the 2045-2055 period under RCP4.5 versus RCP8.5, with all other factors held constant, which translates to USD $1.7 million in averted costs annually. We identify counties with tens to hundreds of avoided ozone-related asthma ED visits under RCP4.5 versus RCP8.5. Overall, we project a heterogeneous distribution of ozone-related asthma ED visits at different spatial resolutions, specifically national, regional, and county levels, and a substantial net health and economic benefit of climate change mitigation.

Air pollution control strategies directly limiting national health damages in the US

Exposure to fine particulate matter (PM2.5) from fuel combustion significantly contributes to global and US mortality. Traditional control strategies typically reduce emissions for specific air pollutants and sectors to maintain pollutant concentrations below standards. Here we directly set national PM2.5 mortality cost reduction targets within a global human-earth system model with US state-level energy systems, in scenarios to 2050, to identify endogenously the control actions, sectors, and locations that most cost-effectively reduce PM2.5 mortality. We show that substantial health benefits can be cost-effectively achieved by electrifying sources with high primary PM2.5 emission intensities, including industrial coal, building biomass, and industrial liquids. More stringent PM2.5 reduction targets expedite the phaseout of high emission intensity sources, leading to larger declines in major pollutant emissions, but very limited co-benefits in reducing CO2 emissions. Control strategies limiting health damages achieve the greatest emission reductions in the East North Central and Middle Atlantic states.

Projections of Atmospheric Nitrogen Deposition to the Chesapeake Bay Watershed

Atmospheric deposition is among the largest pathways of nitrogen loading to the Chesapeake Bay Watershed (CBW). The interplay between future climate and emission changes in and around the CBW will likely shift the future nutrient deposition abundance and chemical regime (e.g., oxidized vs. reduced nitrogen). In this work, a Representative Concentration Pathway (RCP) from the Community Earth System Model is dynamically downscaled using the Weather Research and Forecasting (WRF) and Community Multiscale Air Quality (CMAQ) model coupled to the agro-economic Environmental Policy Integrated Climate (EPIC) model. The relative impacts of emission and climate changes on atmospheric nutrient deposition are explored for a recent historical period and a period centered on 2050. The projected regional emissions in CMAQ reflect current federal and state regulations, which use baseline and projected emission years 2011 and 2040, respectively. The historical simulations of 2-m temperature and precipitation have cool and dry biases, and temperature and precipitation are projected to both increase. Ammonium wet deposition agrees well with observations, but nitrate wet deposition is underpredicted. Climate and deposition changes increase simulated future ammonium fertilizer application. In the CBW at 2050, these changes (along with widespread decreases in anthropogenic nitrogen oxide and sulfur oxide emissions, and relatively constant NH3 emissions) decrease total nitrogen deposition by 21%, decrease annual average oxidized nitrogen deposition by 44%, and increase reduced nitrogen deposition by 10%. These results emphasize the importance of decreased anthropogenic emissions on the control of future nitrogen loading to the Chesapeake Bay in a changing climate.

State-level drivers of future fine particulate matter mortality in the United States

Future fine particulate matter (PM₂.₅) concentrations and resulting health impacts will be largely determined by factors such as energy use, fuel choices, emission controls, state and national policies, and demographcs. In this study, a human-earth system model is used to estimate PM_2.5 mortality costs (PMMC) due to air pollutant emissions from each US state over the period 2015 to 2050, considering current major air quality and energy regulations. Contributions of various socioeconomic and energy factors to PMMC are quantified using the Logarithmic Mean Divisia Index. National PMMC are estimated to decrease 25% from 2015 to 2050, driven by decreases in energy intensity and PMMC per unit consumption of electric sector coal and transportation liquids. These factors together contribute 68% of the decrease, primarily from technology improvements and air quality regulations. States with greater population and economic growth, but with fewer clean energy resources, are more likely to face significant challenges in reducing future PMMC from their emissions. In contrast, states with larger projected decreases in PMMC have smaller increases in population and per capita GDP, and greater decreases in electric sector coal share and PMMC per unit fuel consumption.

Development of a semi-mechanistic allergenic pollen emission model

Modeling pollen emission processes is crucial for studying the spatiotemporal distributions of airborne allergenic pollen. A semi-mechanistic emission model was developed based on mass balance of pollen grain fluxes in the surroundings of allergenic plants. The emission model considers direct emission and resuspension and accounts for influences of temperature, wind velocity, and relative humidity. Modules of this emission model have been developed and parameterized with multiple years of pollen count observations to provide pollen season onset and duration, hourly flowering likelihood, and vegetation coverage for oak and ragweed, as two examples. The simulated spatiotemporal pattern of pollen emissions generally follows the corresponding pattern of area coverage of allergenic plants and diurnal pattern of hourly flowering likelihood. It is found that oak pollen emissions start from the Southern part of the Contiguous United States (CONUS) in March and then shift gradually toward the Northern CONUS, with a maximum emission flux of 5.8 × 10⁶ pollen/(m² h). On the other hand, ragweed pollen emissions start from the Northern CONUS in August and then shift gradually toward the Southern CONUS. The mean ragweed emission flux during August to September can increase up to 2.4 × 10⁶ pollen/(m² h). This emission model is robust with respect to the input parameters for oak and ragweed. Qualitative evaluations of the model performance indicated that the simulated pollen emission is strongly correlated with the plant coverages and observed pollen counts. This model could also be applied to other pollen species given the relevant parameters.

Projected Changes in Maternal Heat Exposure During Early Pregnancy and the Associated Congenital Heart Defect Burden in the United States

Background More intense and longer-lasting heat events are expected in the United States as a consequence of climate change. This study aimed to project the potential changes in maternal heat exposure during early pregnancy (3-8 weeks post conception) and the associated burden of congenital heart defects ( CHD s) in the future. Methods and Results This study expanded on a prior nationwide case-control study that evaluated the association between CHD s and maternal heat exposure during early pregnancy in summer and spring. We defined multiple indicators of heat exposure, and applied published odds ratios obtained for the matching season of the baseline (1995-2005) into the projection period (2025-2035) to estimate potential changes in CHD burden throughout the United States. Increases in maternal heat exposure were projected across the United States and to be larger in the summer. The Midwest will potentially have the highest increase in summer maternal exposure to excessively hot days (3.42; 95% CI, 2.99-3.88 per pregnancy), heat event frequency (0.52; 95% CI, 0.44-0.60) and heat event duration (1.73; 95% CI, 1.49-1.97). We also found large increases in specific CHD subtypes during spring, including a 34.0% (95% CI, 4.9%-70.8%) increase in conotruncal CHD in the South and a 38.6% (95% CI , 9.9%-75.1%) increase in atrial septal defect in the Northeast. Conclusions Projected increases in maternal heat exposure could result in an increased CHD burden in certain seasons and regions of the United States.

Estimating environmental co-benefits of U.S. low-carbon pathways using an integrated assessment model with state-level resolution

There are many technological pathways that can lead to reduced carbon dioxide emissions. However, these pathways can have substantially different impacts on other environmental endpoints, such as air quality and energy-related water demand. This study uses an integrated assessment model with state-level resolution of the energy system to compare environmental impacts of alternative low-carbon pathways for the United States. One set of pathways emphasizes nuclear energy and carbon capture and storage, while another set emphasizes renewable energy, including wind, solar, geothermal power, and bioenergy. These are compared with pathways in which all technologies are available. Air pollutant emissions, mortality costs attributable to particulate matter smaller than 2.5 μm in diameter, and energy-related water demands are evaluated for 50% and 80% carbon dioxide reduction targets in 2050. The renewable low-carbon pathways require less water withdrawal and consumption than the nuclear and carbon capture pathways. However, the renewable low-carbon pathways modeled in this study produce higher particulate matter-related mortality costs due to greater use of biomass in residential heating. Environmental co-benefits differ among states because of factors such as existing technology stock, resource availability, and environmental and energy policies.

Southeast Atmosphere Studies: learning from model-observation syntheses

Concentrations of atmospheric trace species in the United States have changed dramatically over the past several decades in response to pollution control strategies, shifts in domestic energy policy and economics, and economic development (and resulting emission changes) elsewhere in the world. Reliable projections of the future atmosphere require models to not only accurately describe current atmospheric concentrations, but to do so by representing chemical, physical and biological processes with conceptual and quantitative fidelity. Only through incorporation of the processes controlling emissions and chemical mechanisms that represent the key transformations among reactive molecules can models reliably project the impacts of future policy, energy and climate scenarios. Efforts to properly identify and implement the fundamental and controlling mechanisms in atmospheric models benefit from intensive observation periods, during which collocated measurements of diverse, speciated chemicals in both the gas and condensed phases are obtained. The Southeast Atmosphere Studies (SAS, including SENEX, SOAS, NOMADSS and SEAC4RS) conducted during the summer of 2013 provided an unprecedented opportunity for the atmospheric modeling community to come together to evaluate, diagnose and improve the representation of fundamental climate and air quality processes in models of varying temporal and spatial scales. This paper is aimed at discussing progress in evaluating, diagnosing and improving air quality and climate modeling using comparisons to SAS observations as a guide to thinking about improvements to mechanisms and parameterizations in models. The effort focused primarily on model representation of fundamental atmospheric processes that are essential to the formation of ozone, secondary organic aerosol (SOA) and other trace species in the troposphere, with the ultimate goal of understanding the radiative impacts of these species in the southeast and elsewhere. Here we address questions surrounding four key themes: gas-phase chemistry, aerosol chemistry, regional climate and chemistry interactions, and natural and anthropogenic emissions. We expect this review to serve as a guidance for future modeling efforts.

Influence of uncertainties in burned area estimates on modeled wildland fire PM2.5 and ozone pollution in the contiguous U.S

Wildland fires are a major source of fine particulate matter (PM_2.5), one of the most harmful ambient pollutants for human health globally. To represent the influence of wildland fire emissions on atmospheric composition, regional and global chemical transport models rely on emission inventories developed from estimates of burned area (i.e. fire size and location). While different methods of estimating annual burned area agree reasonably well in the western U.S. (within 20-30% for most years during 2002-2014), estimates for the southern U.S. can vary by more than a factor of 5. These differences in burned area lead to significant variability in the spatial and temporal allocation of emissions across fire emission inventory platforms. In this work, we implement wildland fire emission estimates for 2011 from three different products - the USEPA National Emission Inventory (NEI), the Fire Inventory of NCAR (FINN), and the Global Fire Emission Database (GFED4s) - into the Community Multiscale Air Quality (CMAQ) model to quantify and characterize differences in simulated PM and ozone concentrations across the contiguous U.S. (CONUS) due to the fire emission inventory used. The NEI is developed specifically for the U.S., while both FINN and GFED4s are available globally. We find that NEI emissions lead to the largest increases in modeled annual average PM_2.5 (0.85 μg m^-3) and April-September maximum daily 8-h ozone (0.28 ppb) nationally compared to a "no fire" baseline, followed by FINN (0.33 μg m^-3 and 0.22 ppb) and GFED4s (0.12 μg m^-3 and 0.17 ppb). Annual mean enhancements in wildland fire pollution are highest in the southern U.S. across all three inventories (over 4 μg m^-3 and 2 ppb in some areas), but show considerable spatial variability within these regions. We also examine the representation of five individual fire events during 2011 and find that of the two global inventories, FINN reproduces more of the acute changes in pollutant concentrations modeled with NEI and shown in surface observations during each of the episodes investigated compared to GFED4s. Understanding the sensitivity of modeling fire-related PM_2.5 and ozone in the U.S. to burned area estimation approaches will inform future efforts to assess the implications of present and future fire activity for air quality and human health at national and global scales.

The potential effects of climate change on air quality across the conterminous U.S. at 2030 under three Representative Concentration Pathways

The potential impacts of climate change on regional ozone (O3) and fine particulate (PM2.5) air quality in the United States are investigated by linking global climate simulations with regional scale meteorological and chemical transport models. Regional climate at 2000 and at 2030 under three Representative Concentration Pathways (RCPs) is simulated by using the Weather Research and Forecasting (WRF) model to downscale 11-year time slices from the Community Earth System Model (CESM). The downscaled meteorology is then used with the Community Multiscale Air Quality (CMAQ) model to simulate air quality during each of these 11-year periods. The analysis isolates the future air quality differences arising from climate-driven changes in meteorological parameters and specific natural emissions sources that are strongly influenced by meteorology. Other factors that will affect future air quality, such as anthropogenic air pollutant emissions and chemical boundary conditions, are unchanged across the simulations. The regional climate fields represent historical daily maximum and daily minimum temperatures well, with mean biases less than 2 K for most regions of the U.S. and most seasons of the year and good representation of variability. Precipitation in the central and eastern U.S. is well simulated for the historical period, with seasonal and annual biases generally less than 25%, with positive biases exceeding 25% in the western U.S. throughout the year and in part of the eastern U.S. during summer. Maximum daily 8-h ozone (MDA8 O3) is projected to increase during summer and autumn in the central and eastern U.S. The increase in summer mean MDA8 O3 is largest under RCP8.5, exceeding 4 ppb in some locations, with smaller seasonal mean increases of up to 2 ppb simulated during autumn and changes during spring generally less than 1 ppb. Increases are magnified at the upper end of the O3 distribution, particularly where projected increases in temperature are greater. Annual average PM2.5 concentration changes range from -1.0 to 1.0 μg m-3. Organic PM2.5 concentrations increase during summer and autumn due to increased biogenic emissions. Aerosol nitrate decreases during winter, accompanied by lesser decreases in ammonium and sulfate, due to warmer temperatures causing increased partitioning to the gas phase. Among meteorological factors examined to account for modeled changes in pollution, temperature and isoprene emissions are found to have the largest changes and the greatest impact on O3 concentrations.

A Maieutic Exploration of Nudging Strategies for Regional Climate Applications Using the WRF Model

The use of nudging in the Weather Research and Forecasting (WRF) Model to constrain regional climate downscaling simulations is gaining in popularity because it can reduce error and improve consistency with the driving data. While some attention has been paid to whether nudging is beneficial for downscaling, very little research has been performed to determine best practices. In fact, many published papers use the default nudging configuration (which was designed for numerical weather prediction), follow practices used by colleagues, or adapt methods developed for other regional climate models. Here, a suite of 45 three-year simulations is conducted with WRF over the continental United States to systematically and comprehensively examine a variety of nudging strategies. The simulations here use a longer test period than did previously published works to better evaluate the robustness of each strategy through all four seasons, through multiple years, and across nine regions of the United States. The analysis focuses on the evaluation of 2-m temperature and precipitation, which are two of the most commonly required downscaled output fields for air quality, health, and ecosystems applications. Several specific recommendations are provided to effectively use nudging in WRF for regional climate applications. In particular, spectral nudging is preferred over analysis nudging. Spectral nudging performs best in WRF when it is used toward wind above the planetary boundary layer (through the stratosphere) and temperature and moisture only within the free troposphere. Furthermore, the nudging toward moisture is very sensitive to the nudging coefficient, and the default nudging coefficient in WRF is too high to be used effectively for moisture.

Projecting state-level air pollutant emissions using an integrated assessment model: GCAM-USA

Integrated Assessment Models (IAMs) characterize the interactions among human and earth systems. IAMs typically have been applied to investigate future energy, land use, and emission pathways at global to continental scales. Recent directions in IAM development include enhanced technological detail, greater spatial and temporal resolution, and the inclusion of air pollutant emissions. These developments expand the potential applications of IAMs to include support for air quality management and for coordinated environmental, climate, and energy planning. Furthermore, these IAMs could help decision makers more fully understand tradeoffs and synergies among policy goals, identify important cross-sector interactions, and, via scenarios, consider uncertainties in factors such as population and economic growth, technology development, human behavior, and climate change. A version of the Global Change Assessment Model with U.S. state-level resolution (GCAM-USA) is presented that incorporates U.S.-specific emission factors, pollutant controls, and air quality and energy regulations. Resulting air pollutant emission outputs are compared to U.S. Environmental Protection Agency 2011 and projected inventories. A Quality Metric is used to quantify GCAM-USA performance for several pollutants at the sectoral and state levels. This information provides insights into the types of applications for which GCAM-USA is currently well suited and highlights where additional refinement may be warranted. While this analysis is specific to the U.S., the results indicate more generally the importance of enhanced spatial resolution and of considering national and sub-national regulatory constraints within IAMs.

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 (O₃) 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 O₃ 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 O₃ 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 O₃ 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 O₃ are much larger than the impacts of changing demographics. These results indicate the potential for future changes in O₃ exposure as a result of changes in climate that could impact human health.

Climate change impacts on projections of excess mortality at 2030 using spatially varying ozone-temperature risk surfaces

We project the change in ozone-related mortality burden attributable to changes in climate between a historical (1995-2005) and near-future (2025-2035) time period while incorporating a non-linear and synergistic effect of ozone and temperature on mortality. We simulate air quality from climate projections varying only biogenic emissions and holding anthropogenic emissions constant, thus attributing changes in ozone only to changes in climate and independent of changes in air pollutant emissions. We estimate non-linear, spatially varying, ozone-temperature risk surfaces for 94 US urban areas using observed data. Using the risk surfaces and climate projections we estimate daily mortality attributable to ozone exceeding 40 p.p.b. (moderate level) and 75 p.p.b. (US ozone NAAQS) for each time period. The average increases in city-specific median April-October ozone and temperature between time periods are 1.02 p.p.b. and 1.94 °F; however, the results varied by region. Increases in ozone because of climate change result in an increase in ozone mortality burden. Mortality attributed to ozone exceeding 40 p.p.b. increases by 7.7% (1.6-14.2%). Mortality attributed to ozone exceeding 75 p.p.b. increases by 14.2% (1.6 28.9%). The absolute increase in excess ozone mortality is larger for changes in moderate ozone levels, reflecting the larger number of days with moderate ozone levels.

The geographic distribution and economic value of climate change-related ozone health impacts in the United States in 2030

In this United States-focused analysis we use outputs from two general circulation models (GCMs) driven by different greenhouse gas forcing scenarios as inputs to regional climate and chemical transport models to investigate potential changes in near-term U.S. air quality due to climate change. We conduct multiyear simulations to account for interannual variability and characterize the near-term influence of a changing climate on tropospheric ozone-related health impacts near the year 2030, which is a policy-relevant time frame that is subject to fewer uncertainties than other approaches employed in the literature. We adopt a 2030 emissions inventory that accounts for fully implementing anthropogenic emissions controls required by federal, state, and/or local policies, which is projected to strongly influence future ozone levels. We quantify a comprehensive suite of ozone-related mortality and morbidity impacts including emergency department visits, hospital admissions, acute respiratory symptoms, and lost school days, and estimate the economic value of these impacts. Both GCMs project average daily maximum temperature to increase by 1-4°C and 1-5 ppb increases in daily 8-hr maximum ozone at 2030, though each climate scenario produces ozone levels that vary greatly over space and time. We estimate tens to thousands of additional ozone-related premature deaths and illnesses per year for these two scenarios and calculate an economic burden of these health outcomes of hundreds of millions to tens of billions of U.S. dollars (2010$).

IMPLICATIONS

Near-term changes to the climate have the potential to greatly affect ground-level ozone. Using a 2030 emission inventory with regional climate fields downscaled from two general circulation models, we project mean temperature increases of 1 to 4°C and climate-driven mean daily 8-hr maximum ozone increases of 1-5 ppb, though each climate scenario produces ozone levels that vary significantly over space and time. These increased ozone levels are estimated to result in tens to thousands of ozone-related premature deaths and illnesses per year and an economic burden of hundreds of millions to tens of billions of U.S. dollars (2010$).

Variation in estimated ozone-related health impacts of climate change due to modeling choices and assumptions

BACKGROUND

Future climate change may cause air quality degradation via climate-induced changes in meteorology, atmospheric chemistry, and emissions into the air. Few studies have explicitly modeled the potential relationships between climate change, air quality, and human health, and fewer still have investigated the sensitivity of estimates to the underlying modeling choices.

OBJECTIVES

Our goal was to assess the sensitivity of estimated ozone-related human health impacts of climate change to key modeling choices.

METHODS

Our analysis included seven modeling systems in which a climate change model is linked to an air quality model, five population projections, and multiple concentration-response functions. Using the U.S. Environmental Protection Agency's (EPA's) Environmental Benefits Mapping and Analysis Program (BenMAP), we estimated future ozone (O(3))-related health effects in the United States attributable to simulated climate change between the years 2000 and approximately 2050, given each combination of modeling choices. Health effects and concentration-response functions were chosen to match those used in the U.S. EPA's 2008 Regulatory Impact Analysis of the National Ambient Air Quality Standards for O(3).

RESULTS

Different combinations of methodological choices produced a range of estimates of national O(3)-related mortality from roughly 600 deaths avoided as a result of climate change to 2,500 deaths attributable to climate change (although the large majority produced increases in mortality). The choice of the climate change and the air quality model reflected the greatest source of uncertainty, with the other modeling choices having lesser but still substantial effects.

CONCLUSIONS

Our results highlight the need to use an ensemble approach, instead of relying on any one set of modeling choices, to assess the potential risks associated with O(3)-related human health effects resulting from climate change.

Evaluating the effects of climate change on summertime ozone using a relative response factor approach for policymakers

The impact of climate change on surface-level ozone is examined through a multiscale modeling effort that linked global and regional climate models to drive air quality model simulations. Results are quantified in terms of the relative response factor (RRF(E)), which estimates the relative change in peak ozone concentration for a given change in pollutant emissions (the subscript E is added to RRF to remind the reader that the RRF is due to emission changes only). A matrix of model simulations was conducted to examine the individual and combined effects offuture anthropogenic emissions, biogenic emissions, and climate on the RRF(E). For each member in the matrix of simulations the warmest and coolest summers were modeled for the present-day (1995-2004) and future (2045-2054) decades. A climate adjustment factor (CAF(C) or CAF(CB) when biogenic emissions are allowed to change with the future climate) was defined as the ratio of the average daily maximum 8-hr ozone simulated under a future climate to that simulated under the present-day climate, and a climate-adjusted RRF(EC) was calculated (RRF(EC) = RRF(E) x CAF(C)). In general, RRF(EC) > RRF(E), which suggests additional emission controls will be required to achieve the same reduction in ozone that would have been achieved in the absence of climate change. Changes in biogenic emissions generally have a smaller impact on the RRF(E) than does future climate change itself The direction of the biogenic effect appears closely linked to organic-nitrate chemistry and whether ozone formation is limited by volatile organic compounds (VOC) or oxides of nitrogen (NO(x) = NO + NO2). Regions that are generally NO(x) limited show a decrease in ozone and RRF(EC), while VOC-limited regions show an increase in ozone and RRF(EC). Comparing results to a previous study using different climate assumptions and models showed large variability in the CAF(CB).

IMPLICATIONS

We present a methodology for adjusting the RRF to account for the influence of climate change on ozone. The findings of this work suggest that in some geographic regions, climate change has the potential to negate decreases in surface ozone concentrations that would otherwise be achieved through ozone mitigation strategies. In regions of high biogenic VOC emissions relative to anthropogenic NO(x) emissions, the impact of climate change is somewhat reduced, while the opposite is true in regions of high anthropogenic NO(x) emissions relative to biogenic VOC emissions. Further, different future climate realizations are shown to impact ozone in different ways.

Trimethylsilyl derivatives of organic compounds in source samples and in atmospheric fine particulate matter

Source sample extracts of vegetative detritus, motor vehicle exhaust, tire dust paved road dust, and cigarette smoke have been silylated and analyzed by GC-MS to identify polar organic compounds that may serve as tracers for those specific emission sources of atmospheric fine particulate matter. Candidate molecular tracers were also identified in atmospheric fine particle samples collected in the San Joaquin Valley of California. A series of normal primary alkanols, dominated by even carbon-numbered homologues from C26 to C32, the secondary alcohol 10-nonacosanol, and some phytosterols are prominent polar compounds in the vegetative detritus source sample. No new polar organic compounds are found in the motor vehicle exhaust samples. Several hydrogenated resin acids are present in the tire dust sample, which might serve as useful tracers for those sources in areas that are heavily impacted by motor vehicle traffic. Finally, the alcohol and sterol emission profiles developed for all the source samples examined in this project are scaled according to the ambient fine particle mass concentrations attributed to those sources by a chemical mass balance receptor model that was previously applied to the San Joaquin Valley to compute the predicted atmospheric concentrations of individual alcohols and sterols. The resulting underprediction of alkanol concentrations at the urban sites suggests that alkanols may be more sensitive tracers for natural background from vegetative emissions (i.e., waxes) than the high molecular weight alkanes, which have been the best previously available tracers for that source.

Highly polar organic compounds present in wood smoke and in the ambient atmosphere

Fine particulate matter emitted during wood combustion is known to contribute a significant fraction of the total fine aerosol concentration in the atmosphere of both urban and rural areas. In the present study, additional organic compounds that may act as wood smoke tracers in the atmosphere are sought. Polar organic compounds in wood smoke fine particulate matter are converted to their trimethylsilyl derivatives and analyzed by gas chromatography/mass spectrometry. Silylation enables the detection of n-alkanols, plant sterols, and a number of compounds derived from wood lignin that have not previously been reported in wood smoke samples, as well as levoglucosan and related sugar anhydrides formed during the combustion of cellulose. The concentrations of these compounds measured in source emissions are compared to the concentrations in atmospheric fine particle samples collected at a rural background site and at two urban sites in California's San Joaquin Valley. On the basis of this analysis, the sugar anhydrides galactosan and mannosan can be listed along with levoglucosan as being among the most abundant organic compounds detected in all samples.