ACPD Atmospheric Chemistry and Physics Discussions ACPD 1680-7375 Copernicus GmbH Göttingen, Germany 10.5194/acpd-5-6215-2005 Technical note: application of α-QSS to the numerical integration of kinetic equations in tropospheric chemistry Liu F. 1 Schaller E. 1 Mott D. R. 2 Department of Environmental Meteorology, BTU Cottbus, Germany Laboratory for Computational Physics and Fluid Dynamics, Naval Research Laboratory, Washington, DC 20375-5320, USA 18 08 2005 5 4 6215 6262 This is an open-access article ditributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. This article is available from http://www.atmos-chem-phys-discuss.net/5/6215/2005/acpd-5-6215-2005.html The full text article is available as a PDF file from http://www.atmos-chem-phys-discuss.net/5/6215/2005/acpd-5-6215-2005.pdf

A major task in many applications of atmospheric chemistry transport problems is the numerical integration of stiff systems of Ordinary Differential Equations (ODEs) describing the chemical transformations. A faster solver that is easier to couple to the other physics in the problem is still needed. The integration method, α-QSS, corresponding to the solver CHEMEQ2 aims at meeting the demands of a process-split, reacting-flow simulation (Mott 2000; Mott and Oran, 2001). However, this integrator has yet to be applied to the numerical integration of kinetic equations in tropospheric chemistry. A zero-dimensional (box) model is developed to test how well CHEMEQ2 works on the tropospheric chemistry equations. This paper presents the testing results. The reference chemical mechanisms herein used are Regional Atmospheric Chemistry Mechanism (RACM) (Stockwell et al., 1997) and its secondary lumped successor Regional Lumped Atmospheric Chemical Scheme (ReLACS) (Crassier et al., 2000). The box model is forced and initialized by the DRY scenarios of Protocol Ver. 2 developed by EUROTRAC (Poppe et al., 2001). The accuracy of CHEMEQ2 is evaluated by comparing the results to solutions obtained with VODE. This comparison is made with parameters of the error tolerance, relative difference with respect to VODE scheme, trade off between accuracy and efficiency, global time step for integration etc. The study based on the comparison concludes that the single-point α-QSS approach is fast and moderately accurate as well as easy to couple to reacting flow simulation models, which makes CHEMEQ2 one of the best candidates for three-dimensional atmospheric Chemistry Transport Modelling (CTM) studies. In addition the RACM mechanism may be replaced by ReLACS mechanism for tropospheric chemistry transport modelling. The testing results also imply that the accuracy for chemistry numerical simulations is highly different from species to species. Therefore ozone is not the good choice for testing numerical ODE solvers or for evaluation of mechanisms because current tropospheric chemistry mechanisms are mainly designed for troposphere ozone prediction.