Development and evaluation of an advanced National Air Quality Forecasting Capability using the NOAA Global Forecast System version 16.

A new dynamical core, known as the Finite-Volume Cubed-Sphere (FV3) and developed at both NASA and NOAA, is used in NOAA's Global Forecast System (GFS) and in limited-area models for regional weather and air quality applications. NOAA has also upgraded the operational FV3GFS to version 16 (GFSv16), which includes a number of significant developmental advances to the model configuration, data assimilation, and underlying model physics, particularly for atmospheric composition to weather feedback. Concurrent with the GFSv16 upgrade, we couple the GFSv16 with the Community Multiscale Air Quality (CMAQ) model to form an advanced version of the National Air Quality Forecasting Capability (NAQFC) that will continue to protect human and ecosystem health in the US. Here we describe the development of the FV3GFSv16 coupling with a "state-of-the-science" CMAQ model version 5.3.1. The GFS-CMAQ coupling is made possible by the seminal version of the NOAA-EPA Atmosphere-Chemistry Coupler (NACC), which became a major piece of the next operational NAQFC system (i.e., NACC-CMAQ) on 20 July 2021. NACC-CMAQ has a number of scientific advancements that include satellite-based data acquisition technology to improve land cover and soil characteristics and inline wildfire smoke and dust predictions that are vital to predictions of fine particulate matter (PM2.5) concentrations during hazardous events affecting society, ecosystems, and human health. The GFS-driven NACC-CMAQ model has significantly different meteorological and chemical predictions compared to the previous operational NAQFC, where evaluation of NACC-CMAQ shows generally improved near-surface ozone and PM2.5 predictions and diurnal patterns, both of which are extended to a 72 h (3 d) forecast with this system.

Improving surface PM2.5 forecasts in the United States using an ensemble of chemical transport model outputs: 2. bias correction with satellite data for rural areas.

This work serves as the second of a two-part study to improve surface PM2.5 forecasts in the continental U.S. through the integrated use of multi-satellite aerosol optical depth (AOD) products (MODIS Terra/Aqua and VIIRS DT/DB), multi chemical transport model (CTM) (GEOS-Chem, WRF-Chem and CMAQ) outputs and ground observations. In part I of the study, a multi-model ensemble Kalman filter (KF) technique using three CTM outputs and ground observations was developed to correct forecast bias and generate a single best forecast of PM2.5 for next day over non-rural areas that have surface PM2.5 measurements in the proximity of 125 km. Here, with AOD data, we extended the bias correction into rural areas where the closest air quality monitoring station is at least 125 - 300 km away. First, we ensembled all of satellite AOD products to yield the single best AOD. Second, we corrected daily PM2.5 in rural areas from multiple models through the AOD spatial pattern between these areas and non-rural areas, referred to as "extended ground truth" or EGT, for today. Lastly, we applied the KF technique to update the bias in the forecast for next day using the EGT. Our results find that the ensemble of bias-corrected daily PM2.5 from three models for both today and next day show the best performance. Together, the two-part study develops a multi-model and multi-AOD bias correction technique that has the potential to improve PM2.5 forecasts in both rural and non-rural areas in near real time, and be readily implemented at state levels.

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.

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

The potential effects of 21st century climate change on ozone (O3) 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 O3 concentrations. Using the 2011 emissions dataset, O3 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 O3 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.

High-resolution dynamically downscaled rainfall and temperature projections for ecological life zones within Puerto Rico and for the US Virgin Islands.

The Weather Research and Forecasting (WRF) model and a combination of the Regional Spectral Model (RSM) and the Japanese Meteorological Agency Non-Hydrostatic Model (NHM) were used to dynamically downscale selected CMIP5 global climate models to provide 2-km projections with hourly model output for Puerto Rico and the U.S. Virgin Islands. Two 20-year time slices were downscaled for historical (1986-2005) and future (2041-2060) periods following RCP8.5. Projected changes to mean and extreme temperature and precipitation were quantified for Holdridge life zones within Puerto Rico and for the U.S. Virgin Islands. The evaluation reveals a persistent cold bias for all islands in the U.S. Caribbean, a dry bias across Puerto Rico, and a wet bias on the windward side of mountains within the U.S. Virgin Islands. Despite these biases, model simulations show a robust drying pattern for all islands that is generally larger for Puerto Rico (25% annual rainfall reduction for some life zones) than the U.S. Virgin Islands (12% island average). The largest precipitation reductions are found during the more convectively active afternoon and evening hours. Within Puerto Rico, the model uncertainty increases for the wetter life zones, especially for precipitation. Across the life zones, both models project unprecedented maximum and minimum temperatures that may exceed 200 days annually above the historical baseline with only small changes to the frequency of extreme rainfall. By contrast, in the U.S. Virgin Islands, there is no consensus on the location of the largest drying relative to the windward and leeward side of the islands. However, the models project the largest increases in maximum temperature on the southern side of St. Croix and in higher elevations of St. Thomas and St. John.

Projecting changes in extreme rainfall from three tropical cyclones using the design-rainfall approach.

In the past quarter-century, Eastern North Carolina (ENC) experienced several devastating tropical cyclones that led to widespread flooding and damage. Historical climate records reflect an increasing trend in the frequency and intensity of extreme rainfall events across the eastern U.S., which is projected to continue to increase throughout the twenty-first century. Potential changes to extreme rainfall across ENC are explored and quantified for 2025-2100 for three tropical cyclones using an approach based on relative changes in future extreme rainfall frequencies (return periods) from dynamically downscaled projections. Maximum rainfall intensities at '2100' could increase locally by 168%, with widespread regional increases in total rainfall up to 44%. Although these magnitudes exceed the consensus in the literature, the values here are comparable to the most extreme rainfall events observed in the U.S. during the early twenty-first century, which suggests that the intensity of projected future events is already a present-day reality.

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.

Improving Surface PM2.5 Forecasts in the United States Using an Ensemble of Chemical Transport Model Outputs: 1. Bias Correction With Surface Observations in Nonrural Areas.

This work is the first of a two-part study that aims to develop a computationally efficient bias correction framework to improve surface PM2.5 forecasts in the United States. Here, an ensemble-based Kalman filter (KF) technique is developed primarily for nonrural areas with approximately 500 surface observation sites for PM2.5 and applied to three (GEOS-Chem, WRF-Chem, and WRF-CMAQ) chemical transport model (CTM) hindcast outputs for June 2012. While all CTMs underestimate daily surface PM2.5 mass concentration by 20-50%, KF correction is effective for improving each CTM forecast. Subsequently, two ensemble methods are formulated: (1) the arithmetic mean ensemble (AME) that equally weights each model and (2) the optimized ensemble (OPE) that calculates the individual model weights by minimizing the least-square errors. While the OPE shows superior performance than the AME, the combination of either the AME or the OPE with a KF performs better than the OPE alone, indicating the effectiveness of the KF technique. Overall, the combination of a KF with the OPE shows the best results. Lastly, the Successive Correction Method (SCM) was applied to spread the bias correction from model grids with surface PM2.5 observations to the grids lacking ground observations by using a radius of influence of 125 km derived from surface observations, which further improves the forecast of surface PM2.5 at the national scale. Our findings provide the foundation for the second part of this study that uses satellite-based aerosol optical depth (AOD) products to further improve the forecast of surface PM2.5 in rural areas by performing statistical analysis of model output.

Influence of bromine and iodine chemistry on annual, seasonal, diurnal, and background ozone: CMAQ simulations over the Northern Hemisphere.

Bromine and iodine chemistry has been updated in the Community Multiscale Air Quality (CMAQ) model to better capture the influence of natural emissions from the oceans on ozone concentrations. Annual simulations were performed using the hemispheric CMAQ model without and with bromine and iodine chemistry. Model results over the Northern Hemisphere show that including bromine and iodine chemistry in CMAQ not only reduces ozone concentrations within the marine boundary layer but also aloft and inland. Bromine and iodine chemistry reduces annual mean surface ozone over seawater by 25%, with lesser ozone reductions over land. The bromine and iodine chemistry decreases ozone concentration without changing the diurnal profile and is active throughout the year. However, it does not have a strong seasonal influence on ozone over the Northern Hemisphere. Model performance of CMAQ is improved by the bromine and iodine chemistry when compared to observations, especially at coastal sites and over seawater. Relative to bromine, iodine chemistry is approximately four times more effective in reducing ozone over seawater over the Northern Hemisphere (on an annual basis). Model results suggest that the chemistry modulates intercontinental transport and lowers the background ozone imported to the United States.

Updates to the Noah Land Surface Model in WRF-CMAQ to Improve Simulated Meteorology, Air Quality, and Deposition.

Regional, state, and local environmental regulatory agencies often use Eulerian models to investigate the potential impacts on pollutant deposition and air quality from changes in land use, anthropogenic and natural emissions, and climate. The Noah land surface model (LSM) in the Weather Research and Forecasting (WRF) model is widely used with the Community Multiscale Air Quality (CMAQ) model for such investigations, but there are many inconsistencies that need to be changed so that they are consistent with dry deposition and emission processes. In this work, the Noah LSM in WRFv3.8.1 is improved in its linkage to CMAQv5.2 by adding important parameters to the WRF/Noah output, updating the WRF soil and vegetation reference tables that influence CMAQ wet and dry photochemical deposition processes, and decreasing WRF/Noah's top soil layer depth to be consistent with CMAQ processes (e.g., windblown dust and bidirectional ammonia exchange). The modified WRF/Noah-CMAQ system (both off-line and coupled) impacts meteorological predictions of 2-m temperature (T2; increases and decreases), 2-m mixing ratio (Q2; decreases), and 10-m wind speed (WSPD10; decreases) in the United States. These changes are mostly driven by leaf area index values and aerodynamic roughness lengths updated in the vegetation tables based on satellite data, with additional impacts from soil tables updated based on recent soil data. Improvements in the consistency in the treatment of land surface processes between CMAQ and WRF resulted in improvements in both estimated meteorological (e.g., T2, WSPD10, and latent heat fluxes) and chemical (e.g., ozone, sulfur dioxide, and windblown dust) model estimates.

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.

Developing PIDF Curves From Dynamically Downscaled WRF Model Fields to Examine Extreme Precipitation Events in Three Eastern U.S. Metropolitan Areas.

Extreme precipitation events influence watershed, agriculture, and urban management. The probability of extreme precipitation is estimated for storm water management using precipitation intensity-duration-frequency (PIDF) curves. This study explores developing PIDF curves from dynamically downscaled 36- and 12-km simulations using the Weather Research and Forecasting (WRF) model. Three modeled data sets are examined: 36-km WRF model forced with 2.5° (~275-km) NCEP-DOE AMIP-II Reanalysis (R2); 36-km WRF model forced with 0.75° (~80-km) ERA-Interim; and 12-km WRF model forced with ERA-Interim. The WRF outputs are verified against historical observations for three cities in the Eastern United States using a 23-year period (1988-2010). The 36-km WRF data set driven by R2 produced PIDF curves that were acceptable at the 12- to 24-hr durations, but those WRF data could not realistically simulate extremes represented by the high-intensity, short-duration precipitation events. Increasing the resolution of WRF's driving data from R2 to ERA-Interim did not improve WRF's representation of precipitation events. Using 12-km grid spacing enhances WRF's ability to reproduce PIDF curves developed from observations. Finer grid spacing dramatically improves the frequency and intensity of the 1- to 3-hr events and improves the 6- to 24-hr events. However, improvements with the 12-km WRF data did not apply equally to all study sites, suggesting further modifications to the WRF configuration and/or methodology are necessary. Although imperfect, the results here lend confidence to using modeled data to construct PIDF curves, which could be valuable for projecting changes to parameters used in urban and environmental planning.

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.

Effects of Mosaic Land Use on Dynamically Downscaled WRF Simulations of the Contiguous U.S.

The representation of land use (LU) in meteorological modeling strongly influences the simulation of fluxes of heat, moisture, and momentum; affecting the accuracy of 2-m temperature and precipitation. Here, the Weather Research and Forecasting (WRF) model is used with the Noah land surface model to compare a mosaic approach, which accounts for subgrid scale variability of LU types, to the default option which only considers the dominant category in each grid cell. Three-year historical dynamically downscaled WRF simulations are generated using a 12-km domain over the contiguous U.S. to assess the sensitivities to using mosaic LU and to changes to parameters associated with LU and soil categories. Compared to dominant LU, mosaic LU features decreased coverage of forest and agricultural types and increased low-density urban LU throughout much of the eastern and central U.S. However, highly urbanized areas show the opposite trend, as mosaic LU represents partial greenspace within areas that are exclusively urban within dominant LU. Mosaic LU results in widespread increases in sensible heat fluxes and 2-m temperatures, with reductions in latent heat flux, 2-m mixing ratio, and monthly precipitation across the central and eastern U.S. These changes exacerbate an existing warm bias found with dominant LU but reduce overestimations of precipitation. Highly urbanized areas in the eastern U.S. tend to have cooler, more realistic temperatures with mosaic LU relative to dominant LU. A pair of runs with updated surface parameters corroborates these results. Overall, differences between the simulations are largely attributable to their representations of urban LU.

Atmospheric deposition and exceedances of critical loads from 1800-2025 for the conterminous United States.

Atmospheric deposition of nitrogen (N) and sulfur (S) has increased dramatically over pre-industrial levels, with many potential impacts on terrestrial and aquatic ecosystems. Quantitative thresholds, termed "critical loads" (CLs), have been developed to estimate the deposition rate above which damage is thought to occur. However, there remains no comprehensive comparison of when, where, and over what time periods individual CLs have been exceeded. We addressed this knowledge gap by combining several published data sources for historical and contemporary deposition, and overlaying these on six CL types from the National Critical Loads Database (NCLDv2.5; terrestrial acidification, aquatic acidification, lichen, nitrate leaching, plant community composition, and forest-tree health) to examine exceedances from 1800 to 2011. We expressed CLs as the minimum, 10th, and 50th percentiles within 12-km grid cells. Minimum CLs were relatively uniform across the country (200-400 eq·ha-1 ·yr-1 ), and have been exceeded for decades beginning in the early 20th century. The area exceeding minimum CLs peaked in the 1970s and 1980s, exposing 300,000 to 3 million km2 (depending on the CL type) to harmful levels of deposition, with a total area exceeded of 5.8 million km2 (~70% of the conterminous United States). Since then, deposition levels have dropped, especially for S, with modest reductions in exceedance by 2011 for all CL types, totaling 5.2 million km2 in exceedance. The 10th and 50th percentile CLs followed similar trends, but were not consistently available at the 12-km grid scale. We also examined near-term future deposition and exceedances in 2025 under current air quality regulations, and under various scenarios of climate change and additional nitrogen management controls. Current regulations were projected to reduce exceedances of any CL from 5.2 million km2 in 2011 to 4.8 million km2 in 2025. None of the additional N management or climate scenarios significantly affected areal exceedances, although exceedance severity declined. In total, it is clear that many CLs have been exceeded for decades, and are likely to remain so in the short term under current policies. Additionally, we suggest many areas for improvement to enhance our understanding of deposition and its effects to support informed decision making.

Examining WRF's Sensitivity to Contemporary Land-Use Datasets across the Contiguous United States Using Dynamical Downscaling.

Land-use (LU) representation plays a critical role in simulating air-surface interactions that affect meteorological conditions and regional climate. In the Noah LSM within the WRF Model, LU categories are used to set the radiative properties of the surface and to influence exchanges of heat, moisture, and momentum between the air and land surface. Previous literature examined the sensitivity of WRF simulations to LU using short-term meteorological modeling approaches. Here, the sensitivity to LU representation is studied using continental-scale dynamical downscaling, which typically uses longer temporal and larger spatial scales. Two LU datasets, the U.S. Geological Survey (USGS) dataset and the 2006 National Land Cover Dataset (NLCD), are utilized in 3-yr dynamically downscaled WRF simulations over a historical period. Precipitation and 2-m air temperature are evaluated against observation-based datasets for simulations covering the contiguous United States. The WRF-NLCD simulation tends to produce lower precipitation than the WRF-USGS run, with slightly warmer mean monthly temperatures. However, WRF-NLCD results in more notable increases in the frequency of hot days [i.e., days with temperature >90°F (32.2°C)]. These changes are attributable to reductions in forest and agricultural area in the NLCD relative to USGS. There is also subtle but important sensitivity to the method of interpolating LU data to the WRF grid in the model preprocessing. In all cases, the sensitivity resulting from changes in the LU is smaller than model error. Although this sensitivity is small, it persists across spatial and temporal scales.

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.

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.

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.

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.