Variational methods of constructing monotone approximations for atmospheric chemistry models

A new method of constructing efficient monotone numerical schemes for solving direct, adjoint, and inverse atmospheric chemistry problems is presented. It is a synthesis of variational principles combined with splitting and decomposition methods and a constructive implementation of Euler integrating multipliers (EIM) bymeans of a local adjoint problem technique. To increase the efficiency of calculations, a method of decomposing the multicomponent substance transformation operators in terms of the mechanisms of reactions is also proposed. With analytical EIMs, the systems of stiff ODEs are decomposed and reduced to equivalent systems of integral equations solved by noniterative multistage algorithms of a given order of accuracy. An unconventional variational method of constructing mutually consistent algorithms for direct and adjoint problems and sensitivity studies for complex models with constraints is described.