On the accuracy of covariance matrix: Hessian versus Gauss-Newton methods in atmospheric remote sensing with infrared spectroscopy

Journal article

Authors/Editors

Research Areas

No matching items found.

Publication Details

Author list: Tsidu GM

Publisher: Elsevier

Place: OXFORD

Publication year: 2005

Journal: Journal of Quantitative Spectroscopy and Radiative Transfer (0022-4073)

Journal acronym: J QUANT SPECTROSC RA

Volume number: 96

Issue number: 1

Start page: 103

End page: 121

Number of pages: 19

ISSN: 0022-4073

Languages: English-Great Britain (EN-GB)

View in Web of Science | View on publisher site | View citing articles in Web of Science

Abstract

Most retrieval schemes use a linear approximation of the radiative transfer function within each iteration as well as for error analysis. Like most standard methods, the improved Hessian method relies on a quadratic form of the cost function and linear approximation in the error analysis. Often, there is no robust criterion in determining step size that can be used to calculate covariance matrix by discrete perturbation of the cost function in the Hessian approach. The Hessian method improved recently, however, overcomes this problem by employing adaptive algorithm which uses small step sizes in steep directions and large step sizes in flat directions of the cost function. The results of retrievals of atmospheric trace gases from simulated limb emission spectra show that Gauss-Newton algorithm and the improved Hessian generally give nearly identical volume mixing ratios and error covariance matrices in the original state vector space. Due to interlevel correlations, however, the agreement in the uncertainities in the original state vector coordinate system is partly lost in a space in which the elements of state vector are independent after orthogonal coordinate transformation. The significant discrepancies between the estimated uncertainities by the two methods are found to be related with elements of state vector that are dominantly controlled by flattest eigenvector directions of the inverse covariance matrix. The improved Hessian method determines the uncertainities in those shallowest directions with better accuracy than Gauss-Newton approach. The performance of the Hessian method is also found to be better in resolving structures related to the shallowest eigenvector directions as revealed by better vertical resolutions in the retrieved profiles of the trace species. (c) 2005 Elsevier Ltd. All rights reserved.

Keywords

Gauss-Newton method, Hessian method, inverse problem, IR spectroscopy

Documents

No matching items found.