An Adaptive Auto‐Reduction Solver for Speeding Up Integration of Chemical Kinetics in Atmospheric Chemistry Models: Implementation and Evaluation in the Kinetic Pre‐Processor (KPP) Version 3.0.0

Abstract

AbstractKinetic integration of large and stiff chemical mechanisms is a computational bottleneck in models of atmospheric chemistry. It requires implicit solution of the coupled system of kinetic differential equations with time‐consuming construction and inversion of the Jacobian matrix. We present here a new version of the Kinetic Pre‐Processor (KPP 3.0.0) for fast integration of chemical kinetics featuring a range of improvements over previous versions in performance, diagnostics, versatility, and community openness. KPP 3.0.0 includes a new adaptive auto‐reduction solver to decrease the size of any mechanism locally and on the fly under conditions where full complexity is not needed, by partitioning species as “fast” or “slow” based on their local production and loss rates. Previous implementations of this adaptive solver suffered from excessive overhead in the repeated construction of the local Jacobian matrix or were hard‐wired to specific mechanisms. Here we retain the general applicability of the method to any mechanism and avoid overhead by using pre‐computed Jacobian matrix terms for the full mechanism and cropping the matrix locally to remove the slow species with no change in memory allocation. We apply this adaptive solver within KPP 3.0.0 to the GEOS‐Chem global 3‐D model of atmospheric chemistry and demonstrate a 32% reduction in solver time while maintaining a mean error lower than 1% for key species in the troposphere.