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

2021 ◽  
Vol 14 (5) ◽  
pp. 2867-2897
Author(s):  
K. Wyat Appel ◽  
Jesse O. Bash ◽  
Kathleen M. Fahey ◽  
Kristen M. Foley ◽  
Robert C. Gilliam ◽  
...  

Abstract. 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.

2020 ◽  
Author(s):  
K. Wyat Appel ◽  
Jesse O. Bash ◽  
Kathleen M. Fahey ◽  
Kristen M. Foley ◽  
Robert C. Gilliam ◽  
...  

Abstract. 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/iodine chemistry; updated simple halogen chemistry; addition of dimethyl sulfide (DMS) chemistry in the CB6r3 chemical mechanism; updated M3Dry bi-directional deposition model; and the new Surface Tiled Aerosol and Gaseous Exchange (STAGE) bi-directional deposition model. In addition, support for the Weather Research and Forecasting (WRF) model’s hybrid vertical coordinate (HVC) system 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 CMAQ53 was performed by comparing monthly and seasonal mean daily 8-hr average (MDA8) O3 and 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. 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. For daily 24-hr average PM2.5, CMAQ531 has lower concentrations on average in spring and fall, higher concentrations in summer, and similar concentrations in winter as CMAQ521, which slightly increases bias in spring and fall and reduces bias in summer. Comparisons isolating updates to several specific aspects of the modeling system, namely the lateral BCs, meteorology model version, and the deposition model used, were also performed. 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 mid 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 system 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 some surfaces 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.


2010 ◽  
Vol 3 (4) ◽  
pp. 2291-2314
Author(s):  
G. Sarwar ◽  
K. W. Appel ◽  
A. G. Carlton ◽  
R. Mathur ◽  
K. Schere ◽  
...  

Abstract. A new condensed toluene mechanism is incorporated into the Community Multiscale Air Quality Modeling system. Model simulations are performed using the CB05 chemical mechanism containing the existing (base) and the new toluene mechanism for the western and eastern US for a summer month. With current estimates of tropospheric emission burden, the new toluene mechanism increases monthly mean daily maximum 8-h ozone by 1.0–3.0 ppbv in Los Angeles, Portland, Seattle, Chicago, Cleveland, northeastern US, and Detroit compared to that with the base toluene chemistry. It reduces model mean bias for ozone at elevated observed ozone mixing ratios. While the new mechanism increases predicted ozone, it does not enhance ozone production efficiency. Sensitivity study suggests that it can further enhance ozone if elevated toluene emissions are present. While changes in total fine particulate mass are small, predictions of in-cloud SOA increase substantially.


2018 ◽  
Author(s):  
Hervé Petetin ◽  
Bastien Sauvage ◽  
Herman G. J. Smit ◽  
François Gheusi ◽  
Fabienne Lohou ◽  
...  

Abstract. This paper investigates in an innovative way the climatological vertical stratification of relative humidity (RH) and ozone (O3) and carbon monoxide (CO) mixing ratios within the planetary boundary layer (PBL) and at the interface with the free troposphere (FT). The climatology includes all vertical profiles available at northern mid-latitudes over the period 1994–2016 in both IAGOS (In-service Aircraft for a Global Observing System) and WOUDC (World Ozone and Ultraviolet Radiation Data Centre) databases, which represents more than 90,000 vertical profiles. For all individual profiles, apart from the specific case of surface-based temperature inversions (SBIs), the PBL height is estimated following the elevated temperature inversion (EI) method. Several features of both SBIs and EIs are analysed, including their diurnal and seasonal variations. Based on these PBL height estimates (denoted h), the original approach introduced in this paper consists in building a so-called PBL-referenced vertical distribution of O3, CO and RH by averaging all individual profiles beforehand expressed as a function of z/h rather than z (with z the altitude). Using this vertical coordinate system allows to highlight the features existing at the PBL-FT interface that would have been smoothed otherwise. Results demonstrate that the frequently assumed well-mixed PBL remains an exception for both chemical species. Within the PBL, CO profiles are characterized by a mean vertical stratification (here defined as the standard deviation of the CO profile between the surface and the PBL top, normalized by the mean) of 11 %, with moderate seasonal and diurnal variations. A higher vertical stratification is observed for O3 mixing ratios (18 %), with stronger seasonal and diurnal variability (from ~ 10 % in spring/summer midday/afternoon to ~ 25 % in winter/fall night). This vertical stratification is distributed heterogeneously in the PBL with stronger vertical gradients observed at both the surface (due to dry deposition and titration by NO for O3; and due to surface emissions for CO) and the PBL-FT interface. These gradients vary with the season from lowest values in summer to highest ones in winter. Contrary to CO, the O3 vertical stratification was found to vary with the surface potential temperature following an interesting bell shape with weakest stratification for both lowest (typically negative) and highest temperatures, which could be due to a much lower O3 dry deposition under the presence of snow. Therefore, results demonstrate that EIs act as a geophysical interface separating air masses of distinct chemical composition and/or chemical regime. This is further supported by the analysis of the correlation of O3 and CO mixing ratios between the different altitude levels in the PBL and FT (the so-called vertical autocorrelation). Results indeed highlight lower correlations apart from the PBL-FT interface and higher correlations within each of the two atmospheric compartments (PBL and FT).


2012 ◽  
Vol 12 (24) ◽  
pp. 11951-11974 ◽  
Author(s):  
V. Michoud ◽  
A. Kukui ◽  
M. Camredon ◽  
A. Colomb ◽  
A. Borbon ◽  
...  

Abstract. Chemical Ionisation Mass Spectrometer measurements of hydroxyl radical (OH) and the sum of hydroperoxy and organic peroxy (HO2+RO2) radicals were conducted during the MEGAPOLI summer field campaign at the SIRTA observatory near Paris, France, in July 2009. OH and (HO2+RO2) showed a typical diurnal variation with averaged daytime maxima values around 5×106 and 1.2×108 molecule cm−3, respectively. Simultaneously, a large number of ancillary measurements, such as NOx, O3, HONO, HCHO and other VOCs were also conducted. These data provide an opportunity to assess our understanding of the radical chemistry in a suburban environment by comparing the radical observations to calculations. First, OH mixing ratios were estimated by a simple Photo Stationary State (PSS) calculation. PSS calculations overestimate the OH mixing ratio by 50%, especially at NOx mixing ratios lower than 10 ppb, suggesting that some loss processes were missing in the calculation at low NOx. Then, a photochemical box model simulation based on the Master Chemical Mechanism (MCM) and constrained by ancillary measurements was run to calculate radical concentrations. Three different modelling procedures were tested, varying the way the unconstrained secondary species were estimated, to cope with the unavoidable lack of their measurements. They led to significant differences in simulated radical concentrations. OH and (HO2+RO2) concentrations estimated by two selected model version were compared with measurements. These versions of the model were chosen because they lead, respectively, to the higher and lower simulated radical concentrations and are thus the two extremes versions. The box model showed better results than PSS calculations, with a slight overestimation of 12% and 5%, for OH and (HO2+RO2) respectively, in average for the reference model, and an overestimation of approximately 20% for OH and an underestimation for (HO2+RO2) for the other selected model version. Thus, we can conclude from our study that OH and (HO2+RO2) radical levels agree on average with observations within the uncertainty range. Finally, an analysis of the radical budget, on a daily basis (06:00–18:00 UTC), indicates that HONO photolysis (~35%), O3 photolysis (~23%), and aldehydes and ketones photolysis (~16% for formaldehyde and 18% for others) are the main radical initiation pathways. According to the MCM modelling, the reactions of RO2 with NO2 (~19%), leading mainly to PAN formation, is a significant termination pathway in addition to the main net loss via reaction of OH with NO2 (~50%).


2012 ◽  
Vol 12 (10) ◽  
pp. 27173-27218 ◽  
Author(s):  
Y. Xie ◽  
F. Paulot ◽  
W. P. L. Carter ◽  
C. G. Nolte ◽  
D. J. Luecken ◽  
...  

Abstract. The CMAQ model in combination with observations for INTEX-NA/ICARTT 2004 are used to evaluate recent advances in isoprene oxidation chemistry and provide constraints on isoprene nitrate yields, isoprene nitrate lifetimes, and NOx recycling rates. We incorporate recent advances in isoprene oxidation chemistry into the SAPRC-07 chemical mechanism within the US EPA Community Multiscale Air Quality (CMAQ) model. The results show improved model performance for a range of species compared against aircraft observations from the INTEX-NA/ICARTT 2004 field campaign. We further investigate the key processes in isoprene nitrate chemistry and evaluate the impact of uncertainties in the isoprene nitrate yield, NOx (NOx = NO + NO2) recycling efficiency, dry deposition velocity, and RO2 + HO2 reaction rates. We focus our examination in the Southeastern United States, which is impacted by both abundant isoprene emissions and high levels of anthropogenic pollutants. We find that NOx concentrations increase by 4–9% as a result of reduced removal by isoprene nitrate chemistry. O3 increases by 2 ppbv as a result of changes in NOx. OH concentrations increase by 30%, which can be primarily attributed to greater HOx production. We find that the model can capture observed total alkyl and multifunctional nitrates (∑ANs) and their relationship with O3, by assuming either an isoprene nitrate yield of 6% and daytime lifetime of 6 h or a yield of 12% and lifetime of 4 h. Uncertainties in the isoprene nitrates can impact ozone production by 10% and OH concentrations by 6%. The uncertainties in NOx recycling efficiency appear to have larger effects than uncertainties in isoprene nitrate yield and dry deposition velocity. Further progress depends on improved understanding of isoprene oxidation pathways, the rate of NOx recycling from isoprene nitrates, and the fate of the secondary, tertiary, and further oxidation products of isoprene.


2016 ◽  
Author(s):  
K. Wyat Appel ◽  
Sergey L. Napelenok ◽  
Kristen M. Foley ◽  
Havala O. T. Pye ◽  
Christian Hogrefe ◽  
...  

Abstract. The Community Multiscale Air Quality (CMAQ) model is a comprehensive multi-pollutant air quality modeling system developed and maintained by the U.S. 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 which incorporates a large number of science updates and extended capabilities over the previous release version of the model (v5.0.2). These updates include 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 generally result in greater afternoon and early evening mixing in the model, times when the model historically underestimates mixing. The result is higher ozone (O3) mixing ratios on average due to reduced NO titration and lower fine particulate matter (PM2.5) concentrations due to greater dilution of primary pollutants (e.g. elemental and organic carbon). Updates to the clouds and photolysis calculations greatly improve consistency between the WRF and CMAQ models and result in generally higher O3 mixing ratios, primarily due to reduced cloudiness and reduced attenuation of photolysis in the model. Updates to the aerosol chemistry results in higher secondary organic aerosol (SOA) concentrations in the summer, thereby reducing PM2.5 bias, while updates to the gas chemistry result in generally increased O3 in January and July (small) and slightly higher PM2.5 concentrations on average in both January and July. Overall, seasonal variation in simulated PM2.5 generally improves in the new model version, as concentrations decrease in the winter (when PM2.5 is overestimated by CMAQ v5.0.2) and increase in the summer (when PM2.5 is underestimated by CMAQ v5.0.2). Ozone mixing ratios are higher on average with v5.1 versus v5.0.2, resulting in higher O3 mean bias, as O3 tends to be overestimated by CMAQ throughout most of the year (especially at locations where the observed O3 is low), however both the error and correlation are largely improved with v5.1. Sensitivity simulations for several hypothetical emission reduction scenarios showed that v5.1 tends to be slightly more responsive to reductions in NOx (NO + NO2), VOC and SOx (SO2 + SO4) emissions than v5.0.2, representing an improvement as previous studies have shown CMAQ to underestimate the observed reduction in O3 due to large, widespread reductions in observed emissions. Finally, the computational efficiency of the model was significantly improved in v5.1, which keeps runtimes similar to v5.0.2 despite the added complexity to the model.


2012 ◽  
Vol 12 (5) ◽  
pp. 2567-2585 ◽  
Author(s):  
Y. Kanaya ◽  
A. Hofzumahaus ◽  
H.-P. Dorn ◽  
T. Brauers ◽  
H. Fuchs ◽  
...  

Abstract. A photochemical box model constrained by ancillary observations was used to simulate OH and HO2 concentrations for three days of ambient observations during the HOxComp field campaign held in Jülich, Germany in July 2005. Daytime OH levels observed by four instruments were fairly well reproduced to within 33% by a base model run (Regional Atmospheric Chemistry Mechanism with updated isoprene chemistry adapted from Master Chemical Mechanism ver. 3.1) with high R2 values (0.72–0.97) over a range of isoprene (0.3–2 ppb) and NO (0.1–10 ppb) mixing ratios. Daytime HO2(*) levels, reconstructed from the base model results taking into account the sensitivity toward speciated RO2 (organic peroxy) radicals, as recently reported from one of the participating instruments in the HO2 measurement mode, were 93% higher than the observations made by the single instrument. This also indicates an overprediction of the HO2 to OH recycling. Together with the good model-measurement agreement for OH, it implies a missing OH source in the model. Modeled OH and HO2(*) could only be matched to the observations by addition of a strong unknown loss process for HO2(*) that recycles OH at a high yield. Adding to the base model, instead, the recently proposed isomerization mechanism of isoprene peroxy radicals (Peeters and Müller, 2010) increased OH and HO2(*) by 28% and 13% on average. Although these were still only 4% higher than the OH observations made by one of the instruments, larger overestimations (42–70%) occurred with respect to the OH observations made by the other three instruments. The overestimation in OH could be diminished only when reactive alkanes (HC8) were solely introduced to the model to explain the missing fraction of observed OH reactivity. Moreover, the overprediction of HO2(*) became even larger than in the base case. These analyses imply that the rates of the isomerization are not readily supported by the ensemble of radical observations. One of the measurement days was characterized by low isoprene concentrations (∼0.5 ppb) and OH reactivity that was well explained by the observed species, especially before noon. For this selected period, as opposed to the general behavior, the model tended to underestimate HO2(*). We found that this tendency is associated with high NOx concentrations, suggesting that some HO2 production or regeneration processes under high NOx conditions were being overlooked; this might require revision of ozone production regimes.


2009 ◽  
Vol 9 (11) ◽  
pp. 3731-3743 ◽  
Author(s):  
M. Mena-Carrasco ◽  
G. R. Carmichael ◽  
J. E. Campbell ◽  
D. Zimmerman ◽  
Y. Tang ◽  
...  

Abstract. The impact of Mexico City (MCMA) emissions is examined by studying its effects on air quality, photochemistry, and on ozone production regimes by combining model products and aircraft observations from the MILAGRO experiment during March 2006. The modeled influence of MCMA emissions to enhancements in surface level NOx, CO, and O3 concentrations (10–30% increase) are confined to distances <200 km, near surface. However, the extent of the influence is significantly larger at higher altitudes. Broader MCMA impacts (some 900 km Northeast of the city) are shown for specific outflow conditions in which enhanced ozone, NOy, and MTBE mixing ratios over the Gulf of Mexico are linked to MCMA by source tagged tracers and sensitivity runs. This study shows that the "footprint" of MCMA on average is fairly local, with exception to reactive nitrogen, which can be transported long range in the form of PAN, acting as a reservoir and source of NOx with important regional ozone formation implications. The simulated effect of MCMA emissions of anthropogenic aerosol on photochemistry showed a maximum regional decrease of 40% in J[NO2→NO+O], and resulting in the reduction of ozone production by 5–10%. Observed ozone production efficiencies are evaluated as a function of distance from MCMA, and by modeled influence from MCMA. These tend to be much lower closer to MCMA, or in those points where modeled contribution from MCMA is large. This research shows that MCMA emissions do effect on regional air quality and photochemistry, both contributing large amounts of ozone and its precursors, but with caveat that aerosol concentrations hinder formation of ozone to its potential due to its reduction in photolysis rates.


Sign in / Sign up

Export Citation Format

Share Document