CMAQv5.4 Release Notes: Chemistry: Regional Atmospheric Chemistry Mechanism (RACM) - fisidi/CMAQ GitHub Wiki

Bug Fixes to Regional Atmospheric Chemistry Mechanism Version 2 (RACM2)

Golam Sarwar, U.S. Environmental Protection Agency

Type of update: Bug Fix

Release Version/Date: CMAQv5.4

Description:

Two different types of bugs are fixed in RACM2:

The first update:

CMAQv5.4 has two options for calculating in-line biogenic emissions: BEIS and MEGAN. For BEIS, CMAQ uses an emission profile for calculating in-line biogenic emissions. Xiaoyang Chen at Northeastern University notified that CMAQ (BEIS) does not generate any monoterpene emissions when in-line option is enabled with RACM2. Emission profile “B10RD” is currently used for RACM2; however, it does not contain correct mapping which results in no monoterpene emissions. An emission profile “B3V10” was generated when RACM2 was initially implemented in CMAQ which contains correct mapping of model species. The model is revised to remove emission profile “B10RD” and add the emission profile “B3V10”. Biogenic emissions calculation using MEGAN does not use this profile and works properly.

The second update:

RACM2 in CMAQv533 contains several bugs related to kinetics which are now corrected. (1) Rate constant for reaction #R40 is missing a negative sign (mech_racm2_ae6_aq.def) (2) Bill Stockwell identified that raw data files for HCHO photolysis contain errors and provided updated raw data files (HCHO_MOL_RACM2 and HCHO_RAD_RACM2) (3) CSQY_DATA.F contains an initialization problem (CSQY_DATA.F)

Significance and Impact:

The first update:

Model tests were completed with the revised emissions profile. It generates correct biogenic emissions for RACM2 using BEIS.

The second update:

Model sensitivity runs were completed using the existing and updated models (by correcting kinetics information) for 10 days in summer. Ratios of predicted 10-day mean HCHO photolysis rate coefficients with the existing and updated models are shown in Figure 1. For HCHO photolysis, the updated model increases photolysis rate coefficients for the molecular channel [Figure 1(a)] while reducing the photolysis rate coefficients for the radical channel [Figure 1(b)].

Figure 1: (a) Ratio of mean photolysis rate coefficient for the molecular channel of HCHO photolysis (values with updated model / values with existing model) (b) ratio of mean photolysis rate coefficients for the radical channel of HCHO photolysis (values with updated model / values with existing model)

Predicted 10-day mean ozone concentrations with the existing model are shown in Figure 2(a) and differences in predicted mean ozone concentrations with the updated and existing models are shown in Figure 2(b). It has mixed impacts on predicted ozone. The update increases ozone over some areas while reducing it over other areas.

Figure 2: (a) Predicted 10-day mean ozone concentrations with the existing model (b) differences in predicted mean ozone concentrations with the updated and existing models (predictions with updated model – predictions with existing model)

Impacts on other model species are also small. For example, differences in predicted mean sulfate concentrations with the updated and existing models are shown in Figure 3(a) and differences in predicted mean aerosol nitrate concentrations with updated and existing models are shown in Figure 3(b). Impacts on sulfate and aerosol nitrate concentrations are small.

Figure 3: (a) Differences in predicted mean sulfate concentrations with the updated and existing models (predictions with updated model – predictions with existing model) (b) differences in predicted mean aerosol nitrate concentrations with the updated and existing models (predictions with updated model – predictions with existing model)

Internal PRs PR#739, PR#790