Finite Difference of Adjoint or Adjoint of Finite Difference?Source: Monthly Weather Review:;1997:;volume( 125 ):;issue: 012::page 3373DOI: 10.1175/1520-0493(1997)125<3373:FDOAOA>2.0.CO;2Publisher: American Meteorological Society
Abstract: Adjoint models are used for atmospheric and oceanic sensitivity studies in order to efficiently evaluate the sensitivity of a cost function (e.g., the temperature or pressure at some target time tf, averaged over some region of interest) with respect to the three-dimensional model initial conditions. The time-dependent sensitivity, that is the sensitivity to initial conditions as function of the initial time ti, may be obtained directly and most efficiently from the adjoint model solution. There are two approaches to formulating an adjoint of a given model. In the first (?finite difference of adjoint?), one derives the continuous adjoint equations from the linearized continuous forward model equations and then formulates the finite-difference implementation of the continuous adjoint equations. In the second (?adjoint of finite difference?), one derives the finite-difference adjoint equations directly from the finite difference of the forward model. It is shown here that the time-dependent sensitivity obtained by using the second approach may result in a very strong nonphysical behavior such as a large-amplitude two-time-step leapfrog computational mode, which may prevent the solution from being used for time-dependent sensitivity studies. This is an especially relevant problem now, as this second approach is the one used by automatic adjoint compilers that are becoming widely used. The two approaches are analyzed in detail using both a simple model and the adjoint of a primitive equations ocean general circulation model. It is emphasized that both approaches are valid as long as they are used for obtaining the gradient or sensitivity at a single time, as needed in data assimilation, for example. Criteria are presented for the choice of the appropriate adjoint formulation for a given problem.
|
Collections
Show full item record
contributor author | Sirkes, Ziv | |
contributor author | Tziperman, Eli | |
date accessioned | 2017-06-09T16:11:42Z | |
date available | 2017-06-09T16:11:42Z | |
date copyright | 1997/12/01 | |
date issued | 1997 | |
identifier issn | 0027-0644 | |
identifier other | ams-63032.pdf | |
identifier uri | http://onlinelibrary.yabesh.ir/handle/yetl/4203991 | |
description abstract | Adjoint models are used for atmospheric and oceanic sensitivity studies in order to efficiently evaluate the sensitivity of a cost function (e.g., the temperature or pressure at some target time tf, averaged over some region of interest) with respect to the three-dimensional model initial conditions. The time-dependent sensitivity, that is the sensitivity to initial conditions as function of the initial time ti, may be obtained directly and most efficiently from the adjoint model solution. There are two approaches to formulating an adjoint of a given model. In the first (?finite difference of adjoint?), one derives the continuous adjoint equations from the linearized continuous forward model equations and then formulates the finite-difference implementation of the continuous adjoint equations. In the second (?adjoint of finite difference?), one derives the finite-difference adjoint equations directly from the finite difference of the forward model. It is shown here that the time-dependent sensitivity obtained by using the second approach may result in a very strong nonphysical behavior such as a large-amplitude two-time-step leapfrog computational mode, which may prevent the solution from being used for time-dependent sensitivity studies. This is an especially relevant problem now, as this second approach is the one used by automatic adjoint compilers that are becoming widely used. The two approaches are analyzed in detail using both a simple model and the adjoint of a primitive equations ocean general circulation model. It is emphasized that both approaches are valid as long as they are used for obtaining the gradient or sensitivity at a single time, as needed in data assimilation, for example. Criteria are presented for the choice of the appropriate adjoint formulation for a given problem. | |
publisher | American Meteorological Society | |
title | Finite Difference of Adjoint or Adjoint of Finite Difference? | |
type | Journal Paper | |
journal volume | 125 | |
journal issue | 12 | |
journal title | Monthly Weather Review | |
identifier doi | 10.1175/1520-0493(1997)125<3373:FDOAOA>2.0.CO;2 | |
journal fristpage | 3373 | |
journal lastpage | 3378 | |
tree | Monthly Weather Review:;1997:;volume( 125 ):;issue: 012 | |
contenttype | Fulltext |