Implementation and evaluation of an array of chemical solvers in the Global Chemical Transport Model GEOS-Chem

Abstract. This paper discusses the implementation and performance of an array of gas-phase chemistry solvers for the state-of-the-science GEOS-Chem global chemical transport model. The implementation is based on the Kinetic PreProcessor (KPP). Two perl parsers automatically generate the needed interfaces between GEOS-Chem and KPP, and allow access to the chemical simulation code without any additional programming effort. This work illustrates the potential of KPP to positively impact global chemical transport modeling by providing additional functionality as follows. (1) The user can select a highly efficient numerical integration method from an array of solvers available in the KPP library. (2) KPP offers a wide variety of user options for studies that involve changing the chemical mechanism (e.g., a set of additional reactions is automatically translated into efficient code and incorporated into a modified global model). (3) This work provides access to tangent linear, continuous adjoint, and discrete adjoint chemical models, with applications to sensitivity analysis and data assimilation.

[1]  Florian A. Potra,et al.  The kinetic preprocessor KPP*/a software environment for solving chemical kinetics , 2002 .

[2]  John H. Seinfeld,et al.  Global modeling of secondary organic aerosol formation from aromatic hydrocarbons: high- vs. low-yield pathways , 2007 .

[3]  R. Turco,et al.  SMVGEAR: A sparse-matrix, vectorized gear code for atmospheric models , 1994 .

[4]  G. Carmichael,et al.  Direct and adjoint sensitivity analysis of chemical kinetic systems with KPP: Part I—theory and software tools , 2003, Atmospheric Environment.

[5]  G. Socher Conclusions References , 2000 .

[6]  D. Fonteyn,et al.  Four‐dimensional variational chemical assimilation of CRISTA stratospheric measurements , 2001 .

[7]  J. Seinfeld,et al.  Development of the adjoint of GEOS-Chem , 2006 .

[8]  J. Blom,et al.  Benchmarking Stiff ODE Solvers for Atmospheric Chemistry Problems II: Rosenbrock Methods , 1996 .

[9]  Adrian Sandu,et al.  Benchmarking stiff ode solvers for atmospheric chemistry problems II: Rosenbrock solvers , 1997 .

[10]  G. Carmichael,et al.  Direct and adjoint sensitivity analysis of chemical kinetic systems with KPP: II—numerical validation and applications , 2003 .

[11]  Adrian Sandu,et al.  Benchmarking Stiff ODE Solvers for Atmospheric Chemistry Problems I: Implicit versus Explicit , 1996 .

[12]  Adrian Sandu,et al.  The adjoint of CMAQ. , 2007, Environmental science & technology.

[13]  Adrian Sandu,et al.  Technical note: Simulating chemical systems in Fortran90 and Matlab with the Kinetic PreProcessor KPP-2.1 , 2005 .

[14]  Emil M. Constantinescu,et al.  Predicting air quality: Improvements through advanced methods to integrate models and measurements , 2008, J. Comput. Phys..

[15]  M. Jacobson Improvement of SMVGEAR II on vector and scalar machines through absolute error tolerance control , 1998 .

[16]  Simon Chabrillat,et al.  4D-Var assimilation of MIPAS chemical observations: ozone and nitrogen dioxide analyses , 2008 .

[17]  Alan C. Hindmarsh,et al.  Description and use of LSODE, the Livermore Solver for Ordinary Differential Equations , 1993 .