WACCM-D Whole Atmosphere Community Climate Model with D-region ion chemistry

Energetic particle precipitation (EPP) and ion chemistry affect the neutral composition of the polar middle atmosphere. For example, production of odd nitrogen and odd hydrogen during strong events can decrease ozone by tens of percent. However, the standard ion chemistry parameterization used in atmospheric models neglects the effects on some important species, such as nitric acid. We present WACCM-D, a variant of the Whole Atmosphere Community Climate Model, which includes a set of lower ionosphere (D-region) chemistry: 307 reactions of 20 positive ions and 21 negative ions. We consider realistic ionization scenarios and compare the WACCM-D results to those from the Sodankylä Ion and Neutral Chemistry (SIC), a state-of-the-art 1-D model of the D-region chemistry. We show that WACCM-D produces well the main characteristics of the D-region ionosphere, as well as the overall proportion of important ion groups, in agreement with SIC. Comparison of ion concentrations shows that the WACCM-D bias is typically within ±10% or less below 70 km. At 70-90 km, when strong altitude gradients in ionization rates and/or ion concentrations exist, the bias can be larger for some groups but is still within tens of percent. Based on the good agreement overall and the fact that part of the differences are caused by different model setups, WACCM-D provides a state-of-the-art global representation of D-region ion chemistry and is therefore expected to improve EPP modeling considerably. These improvements are demonstrated in a companion paper by Andersson et al.