Calculation of Jacobians for inverse radiative transfer: An efficient hybrid method

Abstract We present an accurate and numerically efficient procedure of calculating Jacobians by finite difference that consists of two components: (1) a method employing the saving of atmospheric layers that accelerates the solution to the equation of radiative transfer for solvers that use the Discrete Space formulation and (2) a method of perturbing the eigenmatrix spectrum associated with a reduced attenuation matrix. The procedure eliminates the need to call the eigenmatrix package, here, LAPACK a second time and provides insights into the fundamental properties of the attenuation matrix, useful for characterizing the accuracy of the derivatives calculated by finite difference methods. The computational complexity of the perturbation method is 8 n 3 + 22 n 2 , where n is one half the number of streams in the radiance field as opposed to 16 n 3 using LAPACK. The method is not limited to the calculation of base state radiances I ( ω ) and those associated with an ‘infinitesimal’ perturbation I ( ω + δ ω ) (from which the numerical derivative of I ( ω + δ ω ) with respect to δ ω may be approximated), but is also useful in the calculation of radiances associated with a ‘finite’ perturbation I ( ω + Δ ω ) from which a sensitivity can be calculated.

[1]  Vladimir V. Rozanov,et al.  GOMETRAN: A radiative transfer model for the satellite project GOME, the plane-parallel version , 1997 .

[2]  Clive D Rodgers,et al.  Inverse Methods for Atmospheric Sounding: Theory and Practice , 2000 .

[3]  V. N. Bogaevski,et al.  Matrix Perturbation Theory , 1991 .

[4]  Thomas Trautmann,et al.  A linearized radiative transfer model for ozone profile retrieval using the analytical forward-adjoint perturbation theory approach , 2001 .

[5]  James Demmel,et al.  Applied Numerical Linear Algebra , 1997 .

[6]  G. Stephens,et al.  Adjoint perturbation and selection rule methods for solar broadband two-stream fluxes in multi-layer media , 2000 .

[7]  David W. Lewis,et al.  Matrix theory , 1991 .

[8]  Knut Stamnes,et al.  A New Look at the Discrete Ordinate Method for Radiative Transfer Calculations in Anisotropically Scattering Atmospheres , 1981 .

[9]  Graeme L. Stephens,et al.  On the fundamental solution of the radiative transfer equation , 1988 .

[10]  Jean-François Mahfouf,et al.  Evaluation of the Jacobians of Infrared Radiation Models for Variational Data Assimilation , 2001 .

[11]  G. E. Hunt,et al.  Discrete space theory of radiative transfer I. Fundamentals , 1969, Proceedings of the Royal Society of London. A. Mathematical and Physical Sciences.

[12]  A. Denning,et al.  Global observations of the carbon budget: 1. Expected satellite capabilities for emission spectroscopy in the EOS and NPOESS eras , 2001 .

[13]  Steven,et al.  Objective Assessment of the Information Content of Visible and Infrared Radiance Measurements for Cloud Microphysical Property Retrievals over the Global Oceans. Part I: Liquid Clouds , 2006 .

[14]  Thomas P. Kurosu,et al.  A linearized discrete ordinate radiative transfer model for atmospheric remote-sensing retrieval , 2001 .

[15]  Laurence S. Rothman,et al.  Journal of Quantitative Spectroscopy & Radiative Transfer 96 , 2005 .

[16]  Robert Spurr,et al.  Simultaneous derivation of intensities and weighting functions in a general pseudo-spherical discrete ordinate radiative transfer treatment , 2002 .

[17]  Angela Benedetti,et al.  Properties of reflected sunlight derived from a Green's function method , 2002 .

[18]  Tristan L'Ecuyer,et al.  Objective Assessment of the Information Content of Visible and Infrared Radiance Measurements for Cloud Microphysical Property Retrievals over the Global Oceans. Part II: Ice Clouds , 2006 .