The state-of-the-science Community Multiscale Air Quality (CMAQ) modeling system, which has recently been extended for hemispheric-scale modeling applications (referred to as H-CMAQ), is applied to ...study the trans-Pacific transport, a phenomenon recognized as a potential source of air pollution in the US, during April 2010. The results of this analysis are presented in two parts. In the previous paper (Part 1), model evaluation for tropospheric ozone (O
) was presented and an air mass characterization method was developed. Results from applying this newly established method pointed to the importance of emissions as the factor to enhance the surface O
mixing ratio over the US. In this subsequent paper (Part 2), emission impacts are examined based on mathematically rigorous sensitivity analysis using the higher-order decoupled direct method (HDDM) implemented in H-CMAQ. The HDDM sensitivity coefficients indicate the presence of a NO
-sensitive regime during April 2010 over most of the Northern Hemisphere. By defining emission source regions over the US and east Asia, impacts from these emission sources are examined. At the surface, during April 2010, the emission impacts of the US and east Asia are comparable over the western US with a magnitude of about 3ppbv impacts on monthly mean O
all-hour basis, whereas the impact of domestic emissions dominates over the eastern US with a magnitude of about 10ppbv impacts on monthly mean O
. The positive correlation (
= 0.63) between surface O
mixing ratios and domestic emission impacts is confirmed. In contrast, the relationship between surface O
mixing ratios and emission impacts from east Asia exhibits a flat slope when considering the entire US. However, this relationship has strong regional differences between the western and eastern US; the western region exhibits a positive correlation (
= 0
36-0.38), whereas the latter exhibits a flat slope (
0
1). Based on the comprehensive evaluation of H-CMAQ, we extend the sensitivity analysis for O
aloft. The results reveal the significant impacts of emissions from east Asia on the free troposphere (defined as 750 to 250hPa) over the US (impacts of more than 5ppbv) and the dominance of stratospheric air mass on upper model layer (defined as 250 to 50hPa) over the US (impacts greater than 10ppbv). Finally, we estimate changes of trans-Pacific transport by taking into account recent emission trends from 2010 to 2015 assuming the same meteorological condition. The analysis suggests that the impact of recent emission changes on changes in the contribution of trans-Pacific transport to US O
levels was insignificant at the surface level and was small (less than 1ppbv) over the free troposphere.
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
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
). For daily 24 h average PM
, 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
concentrations (0.1-0.25 μg m
) 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
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 design of emission abatement measures to effectively reduce high ground-level ozone (O3) concentrations in urban areas is very complex. In addition to the strongly non-linear chemistry of this ...secondary pollutant, precursors can be released by a variety of sources in different regions, and locally produced O3 is mixed with that transported from the regional or continental scales. All of these processes depend also on the specific meteorological conditions and topography of the study area. Consequently, high-resolution comprehensive modeling tools are needed to understand the drivers of photochemical pollution and to assess the potential of local strategies to reduce adverse impacts from high tropospheric O3 levels. In this study, we apply the Integrated Source Apportionment Method (ISAM) implemented in the Community Multiscale Air Quality (CMAQ v5.3.2) model to investigate the origin of summertime O3 in the Madrid region (Spain). Consistent with previous studies, our results confirm that O3 levels are dominated by non-local contributions, representing around 70 % of mean values across the region. Nonetheless, precursors emitted by local sources, mainly road traffic, play a more important role during O3 peaks, with contributions as high as 25 ppb. The potential impact of local measures is higher under unfavorable meteorological conditions associated with regional accumulation patterns. These findings suggest that this modeling system may be used in the future to simulate the potential outcomes of specific emission abatement measures to prevent high-O3 episodes in the Madrid metropolitan area.
The Community Multiscale Air Quality (CMAQ) model is a comprehensive multipollutant air quality modeling system developed and maintained by the US Environmental Protection Agency's (EPA) Office of ...Research and Development (ORD). Recently, version 5.1 of the CMAQ model (v5.1) was released to the public, incorporating a large number of science updates and extended capabilities over the previous release version of the model (v5.0.2). These updates include the following: improvements in the meteorological calculations in both CMAQ and the Weather Research and Forecast (WRF) model used to provide meteorological fields to CMAQ, updates to the gas and aerosol chemistry, revisions to the calculations of clouds and photolysis, and improvements to the dry and wet deposition in the model. Sensitivity simulations isolating several of the major updates to the modeling system show that changes to the meteorological calculations result in enhanced afternoon and early evening mixing in the model, periods when the model historically underestimates mixing. This enhanced mixing results in higher ozone (O
) mixing ratios on average due to reduced NO titration, and lower fine particulate matter (PM
) concentrations due to greater dilution of primary pollutants (e.g., elemental and organic carbon). Updates to the clouds and photolysis calculations greatly improve consistency between the WRF and CMAQ models and result in generally higher O
mixing ratios, primarily due to reduced cloudiness and attenuation of photolysis in the model. Updates to the aerosol chemistry result in higher secondary organic aerosol (SOA) concentrations in the summer, thereby reducing summertime PM
bias (PM
is typically underestimated by CMAQ in the summer), while updates to the gas chemistry result in slightly higher O
and PM
on average in January and July. Overall, the seasonal variation in simulated PM
generally improves in CMAQv5.1 (when considering all model updates), as simulated PM
concentrations decrease in the winter (when PM
is generally overestimated by CMAQ) and increase in the summer (when PM
is generally underestimated by CMAQ). Ozone mixing ratios are higher on average with v5.1 vs. v5.0.2, resulting in higher O
mean bias, as O
tends to be overestimated by CMAQ throughout most of the year (especially at locations where the observed O
is low); however, O
correlation is largely improved with v5.1. Sensitivity simulations for several hypothetical emission reduction scenarios show that v5.1 tends to be slightly more responsive to reductions in NO
(NO + NO
), VOC and SO
(SO
+ SO
) emissions than v5.0.2, representing an improvement as previous studies have shown CMAQ to underestimate the observed reduction in O
due to large, widespread reductions in observed emissions.
Numerous scientific upgrades to the representation of secondary organic aerosol (SOA) are incorporated into the Community Multiscale Air Quality (CMAQ) modeling system. Additions include several ...recently identified SOA precursors: benzene, isoprene, and sesquiterpenes; and pathways: in-cloud oxidation of glyoxal and methylglyoxal, particle-phase oligomerization, and acid enhancement of isoprene SOA. NO x -dependent aromatic SOA yields are also added along with new empirical measurements of the enthalpies of vaporization and organic mass-to-carbon ratios. For the first time, these SOA precursors, pathways and empirical parameters are included simultaneously in an air quality model for an annual simulation spanning the continental U.S. Comparisons of CMAQ-modeled secondary organic carbon (OCsec) with semiempirical estimates screened from 165 routine monitoring sites across the U.S. indicate the new SOA module substantially improves model performance. The most notable improvement occurs in the central and southeastern U.S. where the regionally averaged temporal correlations (r) between modeled and semiempirical OCsec increase from −0.5 to 0.8 and −0.3 to 0.8, respectively, when the new SOA module is employed. Wintertime OCsec results improve in all regions of the continental U.S. and the seasonal and regional patterns of biogenic SOA are better represented.
Satellite observations of tropospheric NO2 columns can provide top-down observational constraints on emissions estimates of nitrogen oxides (NOx). Mass-balance-based methods are often applied for ...this purpose but do not isolate near-surface emissions from those aloft, such as lightning emissions. Here, we introduce an inverse modeling framework that couples satellite chemical data assimilation to a chemical transport model. In the framework, satellite-constrained emissions totals are inferred using model simulations with and without data assimilation in the iterative finite-difference mass-balance method. The approach improves the finite-difference mass-balance inversion by isolating the near-surface emissions increment. We apply the framework to separately estimate lightning and anthropogenic NOx emissions over the Northern Hemisphere for 2019. Using overlapping observations from the Ozone Monitoring Instrument (OMI) and the Tropospheric Monitoring Instrument (TROPOMI), we compare separate NOx emissions inferences from these satellite instruments, as well as the impacts of emissions changes on modeled NO2 and O3. OMI inferences of anthropogenic emissions consistently lead to larger emissions than TROPOMI inferences, attributed to a low bias in TROPOMI NO2 retrievals. Updated lightning NOx emissions from either satellite improve the chemical transport model's low tropospheric O3 bias. The combined lighting and anthropogenic emissions updates improve the model's ability to reproduce measured ozone by adjusting natural, long-range, and local pollution contributions. Thus, the framework informs and supports the design of domestic and international control strategies.
Ammonia (NH3) emissions have large impacts on air quality and nitrogen
deposition, influencing human health and the well-being of sensitive
ecosystems. Large uncertainties exist in the “bottom-up” ...NH3 emission
inventories due to limited source information and a historical lack of
measurements, hindering the assessment of NH3-related environmental
impacts. The increasing capability of satellites to measure NH3
abundance and the development of modeling tools enable us to better
constrain NH3 emission estimates at high spatial resolution. In this
study, we constrain the NH3 emission estimates from the widely used
2011 National Emissions Inventory (2011 NEI) in the US using Infrared
Atmospheric Sounding Interferometer NH3 column density measurements
(IASI-NH3) gridded at a 36 km by 36 km horizontal resolution. With a
hybrid inverse modeling approach, we use the Community Multiscale Air Quality Modeling System (CMAQ) and its multiphase adjoint model to optimize NH3 emission estimates in April, July, and October.
Our optimized emission estimates suggest that the total NH3 emissions
are biased low by 26 % in 2011 NEI in April with overestimation in the Midwest
and underestimation in the Southern States. In July and October, the
estimates from NEI agree well with the optimized emission estimates, despite
a low bias in hotspot regions. Evaluation of the inversion performance using
independent observations shows reduced underestimation in simulated ambient
NH3 concentration in all 3 months and reduced underestimation in
NH4+ wet deposition in April. Implementing the optimized NH3
emission estimates improves the model performance in simulating PM2.5
concentration in the Midwest in April. The model results suggest that the
estimated contribution of ammonium nitrate would be biased high in a priori
NEI-based assessments. The higher emission estimates in this study also
imply a higher ecological impact of nitrogen deposition originating from
NH3 emissions.
Crop residue burning is a common land management practice that results in emissions of a variety of pollutants with negative health impacts. Modeling systems are used to estimate air quality impacts ...of crop residue burning to support retrospective regulatory assessments and also for forecasting purposes. Ground and airborne measurements from a recent field experiment in the Pacific Northwest focused on cropland residue burning was used to evaluate model performance in capturing surface and aloft impacts from the burning events. The Community Multiscale Air Quality (CMAQ) model was used to simulate multiple crop residue burns with 2 km grid spacing using field-specific information and also more general assumptions traditionally used to support National Emission Inventory based assessments. Field study specific information, which includes area burned, fuel consumption, and combustion completeness, resulted in increased biomass consumption by 123 tons (60% increase) on average compared to consumption estimated with default methods in the National Emission Inventory (NEI) process. Buoyancy heat flux, a key parameter for model predicted fire plume rise, estimated from fuel loading obtained from field measurements can be 30% to 200% more than when estimated using default field information. The increased buoyancy heat flux resulted in higher plume rise by 30% to 80%. This evaluation indicates that the regulatory air quality modeling system can replicate intensity and transport (horizontal and vertical) features for crop residue burning in this region when region-specific information is used to inform emissions and plume rise calculations. Further, previous vertical emissions allocation treatment of putting all cropland residue burning in the surface layer does not compare well with measured plume structure and these types of burns should be modeled more similarly to prescribed fires such that plume rise is based on an estimate of buoyancy.
Display omitted
•Detailed evaluation of fuel loading and emission factors of PM2.5 and CO for bluegrass and winter wheat fuel types.•30% to 200% underestimation in buoyancy heat flux with default field information underestimate plume height up to 80.•Improved plume structure modeling for crop residual burning with field measurements based buoyancy heat flux estimation.
The state-of-the-science Community Multiscale Air Quality (CMAQ) Modeling System has recently been extended for hemispheric-scale modeling applications (referred to as H-CMAQ). In this study, ...satellite-constrained estimation of the degassing SO2 emissions from 50 volcanoes over the Northern Hemisphere is incorporated into H-CMAQ, and their impact on tropospheric sulfate aerosol (SO42-) levels is assessed for 2010. The volcanic degassing improves predictions of observations from the Acid Deposition Monitoring Network in East Asia (EANET), the United States Clean Air Status and Trends Network (CASTNET), and the United States Integrated Monitoring of Protected Visual Environments (IMPROVE). Over Asia, the increased SO42- concentrations were seen to correspond to the locations of volcanoes, especially over Japan and Indonesia. Over the USA, the largest impacts that occurred over the central Pacific were caused by including the Hawaiian Kilauea volcano, while the impacts on the continental USA were limited to the western portion during summertime. The emissions of the Soufrière Hills volcano located on the island of Montserrat in the Caribbean Sea affected the southeastern USA during the winter season. The analysis at specific sites in Hawaii and Florida also confirmed improvements in regional performance for modeled SO42- by including volcanoes SO2 emissions. At the edge of the western USA, monthly averaged SO42- enhancements greater than 0.1 µgm-3 were noted within the boundary layer (defined as surface to 750 hPa) during June–September. Investigating the change on SO42- concentration throughout the free troposphere revealed that although the considered volcanic SO2 emissions occurred at or below the middle of free troposphere (500 hPa), compared to the simulation without the volcanic source, SO42- enhancements of more than 10 % were detected up to the top of the free troposphere (250 hPa). Our model simulations and comparisons with measurements across the Northern Hemisphere indicate that the degassing volcanic SO2 emissions are an important source and should be considered in air quality model simulations assessing background SO42- levels and their source attribution.
The state-of-the-science Community Multiscale Air Quality (CMAQ) Modeling System has recently been extended for hemispheric-scale modeling applications (referred to as H-CMAQ). In this study, ...satellite-constrained estimation of the degassing SO
emissions from 50 volcanoes over the Northern Hemisphere is incorporated into H-CMAQ, and their impact on tropospheric sulfate aerosol (
) levels is assessed for 2010. The volcanic degassing improves predictions of observations from the Acid Deposition Monitoring Network in East Asia (EANET), the United States Clean Air Status and Trends Network (CASTNET), and the United States Integrated Monitoring of Protected Visual Environments (IMPROVE). Over Asia, the increased
concentrations were seen to correspond to the locations of volcanoes, especially over Japan and Indonesia. Over the USA, the largest impacts that occurred over the central Pacific were caused by including the Hawaiian Kilauea volcano, while the impacts on the continental USA were limited to the western portion during summertime. The emissions of the Soufrière Hills volcano located on the island of Montserrat in the Caribbean Sea affected the southeastern USA during the winter season. The analysis at specific sites in Hawaii and Florida also confirmed improvements in regional performance for modeled
by including volcanoes SO
emissions. At the edge of the western USA, monthly averaged
enhancements greater than 0.1μgm
were noted within the boundary layer (defined as surface to 750hPa) during June- September. Investigating the change on
concentration throughout the free troposphere revealed that although the considered volcanic SO
emissions occurred at or below the middle of free troposphere (500hPa), compared to the simulation without the volcanic source,
enhancements of more than 10% were detected up to the top of the free troposphere (250hPa). Our model simulations and comparisons with measurements across the Northern Hemisphere indicate that the degassing volcanic SO
emissions are an important source and should be considered in air quality model simulations assessing background
levels and their source attribution.