Cloud tomography: Role of constraints and a new algorithm

[1] Retrieving spatial distributions of cloud liquid water content from limited-angle emission data (passive microwave cloud tomography) is ill-posed, and a small inaccuracy in the data and/or numerical treatments may result in a large error in the retrieval. Proper handling of the ill-posedness is an ongoing challenge to the atmospheric remote sensing community. In this paper we first analyze the major regularization methods that each apply a single but different constraint to their retrievals and extend these methods to allow for multiple constraints. We then develop a new iterative algorithm that can also incorporate complex physical constraints with great flexibility. To understand the influences of different constraints on the retrievals, we use the new iterative algorithm with various combinations of constraints to retrieve a stratocumulus cloud simulated with a large-eddy-simulation model. For this relatively homogeneous cloud case, the standard least squares method with no constraints, as expected, performs very poorly, and yields a very large retrieval error, making this method nearly useless. Adding a nonnegativity constraint reduces the mean retrieval error by a factor of 6 but the internal structure of the cloud is still not reproduced in the retrieval. Adding a smoothness constraint dramatically improves the retrieved spatial structure of the cloud, and brings the mean error down further, although the retrieved cloud top edges are still considerably blurred. Further adding the so-called double-side constraint (based on scaled adiabatic profiles) produces the best result; the retrieval faithfully reproduces the cloud water structure with a mean retrieval error of only one third of that of the nonnegativity and smoothness constrained method.

[1]  Dong Huang,et al.  Determination of cloud liquid water distribution using 3D cloud tomography , 2008 .

[2]  J. Drake,et al.  Field Tests of an Airborne Remote Sensing Technique for Measuring the Distribution of Liquid Water in Convective Cloud , 1988 .

[3]  J. Snider,et al.  Liquid Water Distribution Obtained from Coplanar Scanning Radiometers , 1986 .

[4]  P. Hobbs,et al.  A Model for Particle Microphysics,Turbulent Mixing, and Radiative Transfer in the Stratocumulus-Topped Marine Boundary Layer and Comparisons with Measurements , 1995 .

[5]  Dimitri P. Bertsekas,et al.  Nonlinear Programming , 1997 .

[6]  S. Twomey Iterafive Nonlinear Inversion Methods for Tomographic Problems , 1987 .

[7]  J. Verlinde,et al.  The Removal of Turbulent Broadening in Radar Doppler Spectra Using Linear Inversion with Double-Sided Constraints , 2000 .

[8]  J. Drake,et al.  A Theoretical Study of the Accuracy of Tomographic Retrieval of Cloud Liquid with an Airborne Radiometer , 1988 .

[9]  A Tikhonov,et al.  Solution of Incorrectly Formulated Problems and the Regularization Method , 1963 .

[10]  C. Lawson,et al.  Solving least squares problems , 1976, Classics in applied mathematics.

[11]  Bert W. Rust,et al.  Constrained Least Squares Interval Estimation , 1985 .

[12]  P. Krehbiel,et al.  Determination of Cloud Liquid Water Distribution by Inversion of Radiometric Data , 1985 .

[13]  W. Arnott,et al.  Particle size distribution retrieval from multispectral optical depth: Influences of particle nonsphericity and refractive index , 1999 .

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

[15]  S. Twomey Introduction to the Mathematics of Inversion in Remote Sensing and Indirect Measurements , 1997 .

[16]  Per Christian Hansen,et al.  Analysis of Discrete Ill-Posed Problems by Means of the L-Curve , 1992, SIAM Rev..

[17]  David L. Phillips,et al.  A Technique for the Numerical Solution of Certain Integral Equations of the First Kind , 1962, JACM.

[18]  Ben G. Fitzpatrick,et al.  Bayesian analysis in inverse problems , 1991 .