Discrete second order adjoints in atmospheric chemical transport modeling

  • Authors:

  • Adrian Sandu;Lin Zhang

  • Affiliations:

  • Department of Computer Science, Virginia Polytechnic Institute and State University, 2244 Knowledgeworks II, 2202 Kraft Drive, Blacksburg, VA 24061, United States;Department of Computer Science, Virginia Polytechnic Institute and State University, 2244 Knowledgeworks II, 2202 Kraft Drive, Blacksburg, VA 24061, United States

  • Venue:
  • Journal of Computational Physics
  • Year:
  • 2008

Quantified Score

Hi-index 31.45

Visualization

Abstract

Atmospheric chemical transport models (CTMs) are essential tools for the study of air pollution, for environmental policy decisions, for the interpretation of observational data, and for producing air quality forecasts. Many air quality studies require sensitivity analyses, i.e., the computation of derivatives of the model output with respect to model parameters. The derivatives of a cost functional (defined on the model output) with respect to a large number of model parameters can be calculated efficiently through adjoint sensitivity analysis. While the traditional (first order) adjoint models give the gradient of the cost functional with respect to parameters, second order adjoint models give second derivative information in the form of products between the Hessian of the cost functional and a vector (representing a perturbation in sensitivity analysis, a search direction in optimization, an eigenvector, etc.). In this paper we discuss the mathematical foundations of the discrete second order adjoint sensitivity method and present a complete set of computational tools for performing second order sensitivity studies in three-dimensional atmospheric CTMs. The tools include discrete second order adjoints of Runge-Kutta and of Rosenbrock time stepping methods for stiff equations together with efficient implementation strategies. Numerical examples illustrate the use of these computational tools in important applications like sensitivity analysis, optimization, uncertainty quantification and the calculation of directions of maximal error growth in three-dimensional atmospheric CTMs.