Particle size distribution retrieval from multispectral optical depth: Influences of particle nonsphericity and refractive index

Abstract

Abstract. Retrieval of size distributions from multispectral optical depth measurements requires solution of an ill-posed inverse problem. The ill-posedness causes problems such as solution ambiguity. Size distribution retrieval becomes more complicated in the presence of nonspherical particles and/or refractive index errors. A new retrieval algorithm is first developed which allows for both smoothing and nonnegativity constraint along with the L-curve method for choosing the Lagrange multiplier that controls the degree of the imposed smoothing. This new algorithm is compared to an iterative algorithm and the method of truncated singular value decomposition, demonstrating that the new algorithm outperforms the other two. With the new algorithm to perform the size distribution retrieval and with the T-matrix method to calculate optical depth for a given cloud consisting of randomly oriented finite circular cylinders, the influence of particle nonsphericity on the size distribution retrieval is investigated by use of the Mie theory for spheres as well as the anomalous diffraction theory for cylinders. The results show that spurious particles and even spurious particle modes occur for both approximations. The effect of refractive index errors are also investigated, showing that even a small perturbation of refractive indices can cause serious distortions of retrieved size distributions. A further examination reveals that either applying an approximate lightscattering theory (Mie theory or anomalous diffraction theory) to nonspherical particles or using incorrect refractive indices results in systematic errors in the model which in turn conspire with the ill-posedness inherent in the retrieval to cause the distortions of retrieved size distributions. The retrieval process essentially transforms the model error into the error in retrieved size distribution, yet improves the agreement between true and retrieved optical depth. IntroductionVarious properties (e.g., dynamical, optical, electrical, and chemical) of atmospheric particles depend on particle size. Number distribution of particles with respect to size (size distribution) is a fundamental subject in aerosol and cloud physics [Pruppacher and Klett, 1978]. Its importance has been increasingly recognized with the concern over the effects of aerosols

Probing ice clouds by broadband mid-infrared extinction spectroscopy: case studies from ice nucleation experiments in the AIDA aerosol and cloud chamber

“…For example, Baran (2004) has recently summarised concepts for representing the various shapes of ice crystals that exist in cirrus clouds as well as computational methods for calculating their single-scattering properties. As clearly shown by Liu et al (1999), the infrared retrieval results are notably influenced by different assumptions on the shape of the ice crystals. This was demonstrated by calculating a synthetic infrared extinction spectrum for a given size distribution of finite circular cylinders of aspect ratio φ=1 using the T-matrix method and inverting the spectrum applying Mie theory.…”

“…In natural cirrus clouds, there is a considerable variability in the size of the ice crystals, ranging from <10 µm to >1 mm (Lynch et al, 2002). In the AIDA chamber studies, the size of the ice crystals is mainly determined by the number concentration of ice activated seed aerosol particles and the prevailing water vapour pressure, which, due to the ice-coated chamber walls, exponentially increases with increasing temperature.…”

“…The basic procedure for the quantitative analysis of the extinction spectra of the ice crystals has been adopted from previous laboratory studies on the characterisation of aerosols from infrared optical depth measurements (Arnott et al, 1997b;Liu et al, 1999;Zasetsky et al, 2004). It is based on the Lambert Beer equation, relating the measured optical depth τ ν j at a specific wavenumberν j according to…”

“…Due to the important impact of ice crystal clouds (cirrus) on the earth-atmosphere radiation balance, numerous field and laboratory studies have been conducted to explore and to quantify the two distinct freezing mechanisms that lead to the formation of the ice phase: (i) homogeneous freezing of aqueous solution droplets, and (ii) heterogeneous ice nucleation by solid materials, so-called ice nuclei (Lynch et al, 2002). For this purpose, field missions, like the CRYSTAL-FACE measurement campaign on tropical cirrus cloud physical properties and formation processes, are featuring extensive instrumentation for the characterisation of both the aerosol and the ice particles (http:// cloud1.arc.nasa.gov/crystalface/instruments.html).…”

Cloud tomography: Role of constraints and a new algorithm

“…[22] Liu et al [1999] developed an approach by combining the Tikhonov regularization with the NNLS in retrieval of the particle size distribution from measurements of multispectral optical depths, and dubbed this approach as smoothness-constrained NNLS (SCNNLS). Essentially, the SCNNLS requires solving the following minimization problem:…”