Engineering Papers⌕ Search

Engineering topics

Mishchenko, Michael I.

Publications and source records attributed to Mishchenko, Michael I..

At least 73 records · Page 4

Direct Simulation of Multiple Scattering by Discrete Random Media Illuminated by Gaussian Beams

The conventional orientation-averaging procedure developed in the framework of the superposition T-matrix approach is generalized to include the case of illumination by a Gaussian beam (GB). The resulting computer code is parallelized and used to perform extensive numerically exact calculations of electromagnetic scattering by volumes of discrete random medium consisting of monodisperse spherical particles. The size parameters of the scattering volumes are 40, 50, and 60, while their packing density is fixed at 5%. We demonstrate that all scattering patterns observed in the far-field zone of a random multisphere target and their evolution with decreasing width of the incident GB can be interpreted in terms of idealized theoretical concepts such as forward-scattering interference, coherent backscattering (CB), and diffuse multiple scattering. It is shown that the increasing violation of electromagnetic reciprocity with decreasing GB width suppresses and eventually eradicates all observable manifestations of CB. This result supplements the previous demonstration of the effects of broken reciprocity in the case of magneto-optically active particles subjected to an external magnetic field.

Mackowski, Daniel W.↗

Comprehensive T-Matrix Reference Database: A 2007-2009 Update

The T-matrix method is among the most versatile, efficient, and widely used theoretical techniques for the numerically exact computation of electromagnetic scattering by homogeneous and composite particles, clusters of particles, discrete random media, and particles in the vicinity of an interface separating two half-spaces with different refractive indices. This paper presents an update to the comprehensive database of T-matrix publications compiled by us previously and includes the publications that appeared since 2007. It also lists several earlier publications not included in the original database.

Mishchenko, Michael I.↗

Toward Unified Satellite Climatology of Aerosol Properties. 3. MODIS Versus MISR Versus AERONET

We use the full duration of collocated pixel-level MODIS-Terra and MISR aerosol optical thickness (AOT) retrievals and level 2 cloud-screened quality-assured AERONET measurements to evaluate the likely individual MODIS and MISR retrieval accuracies globally over oceans and land. We show that the use of quality-assured MODIS AOTs as opposed to the use of all MODIS AOTs has little effect on the resulting accuracy. The MODIS and MISR relative standard deviations (RSTDs) with respect to AERONET are remarkably stable over the entire measurement record and reveal nearly identical overall AOT performances of MODIS and MISR over the entire suite of AERONET sites. This result is used to evaluate the likely pixel-level MODIS and MISR performances on the global basis with respect to the (unknown) actual AOTs. For this purpose, we use only fully compatible MISR and MODIS aerosol pixels. We conclude that the likely RSTDs for this subset of MODIS and MISR AOTs are 73% over land and 30% over oceans. The average RSTDs for the combined [AOT(MODIS)+AOT(MISR)]/2 pixel-level product are close to 66% and 27%, respectively, which allows us to recommend this simple blend as a better alternative to the original MODIS and MISR data. These accuracy estimates still do not represent the totality of MISR and quality-assured MODIS pixel-level AOTs since an unaccounted for and potentially significant source of errors is imperfect cloud screening. Furthermore, many collocated pixels for which one of the datasets reports a retrieval, whereas the other one does not may also be problematic.

Mishchenko, Michael I.↗

The Study of Biomass Emissions for Defining Radiative Forcing of Climate

Accurate quantification of the amounts of trace gases and particulate matter emitted from vegetation fires and other sources of biomass burning (agricultural waste and biofuels) on a regional and global basis is required by a number of users, including scientists studying a wide range of atmospheric processes, national governments who are required to report greenhouse gas emissions, and those interested in quantifying the sources of air pollution that affect human health at regional scales. Over the past decade, improvements in the ability to detect and map fires using a number of different satellite systems have been achieved, largely through efforts coordinated through working groups organized by the IGBP Data and Information System and Global Observation of Forest Cover (GOFC) projects. In addition, significant advances and improvement in or understanding of the emissions factors for biomass burning in different biomes has resulted through efforts by the Biomass Burning Experiment (BIBEX) organized through the International Global Atmospheric Chemistry project. A number of satellite-based fire data products have been generated, and a number of new products will shortly be available. These new data products will provide the basis for estimating emissions from biomass burning on a global basis. However, a number of issues remain concerning the availability of other data sets needed to generate these estimates. Recognizing these issues, the GOFC-Fire Satellite Validation Workshop (held in Lisbon, Portugal on 9-11 July 2001), recommended that a workshop focusing on Improving Global Estimates of Atmospheric Emissions from Biomass Burning be organized. This workshop was held from 17- 19 July 2002 on the campus of the University of Maryland, College Park, Maryland. This workshop served as the annual meeting of the GOFC/GOLD-Fire Program. The overall goals of the meeting were to review the information products generated from satellite imagery and other sources that are currently available for developing emission estimates from biomass burning, evaluate areas where improved or additional products would be beneficial, and recommend products for use by the atmospheric science community.

Penner, Joyce↗

Vector Radiative Transfer Equation for Arbitrarily Shaped and Arbitrarily Oriented Particles: A Microphysical Derivation from Statistical Electromagnetics

We use the concepts of statistical electromagnetics to derive the general radiative transfer equation (RTE) describing multiple scattering of polarized light by sparse discrete random media consisting of arbitrarily shaped and arbitrarily oriented particles. The derivation starts from the volume integral and Lippmann-Schwinger equations for the electric field scattered by a fixed N-particle system and proceeds to the vector form of the Foldy-Lax equations and their approximate far-field version. We then assume that particle positions are completely random and derive the vector RTE by applying the Twersky approximation to the coherent electric field and the Twersky and ladder approximations to the coherency dyad of the diffuse field in the limit N -> infinity. The concluding section discusses the physical meaning of the quantities entering the general vector RTE and the assumptions made in its derivation.

Mishchenko, Michael I.↗

Asymptotic Solutions for Optical Properties of Large Particles with Strong Absorption

The transverse wave condition is not applicable to the refracted electromagnetic wave within the context of geometrical optics when absorption is involved. Either the transverse magnetic (TM) or the transverse electric (TE) wave condition can be assumed for the wave to locally satisfy the electromagnetic boundary condition in a ray-tracing calculation. The assumed wave mode affects both the reflection and the refraction coefficients. As a result, nonunique solutions for these coefficients are inevitable. In this study the appropriate solutions for the Fresnel reflection-refraction coefficients are identified in light-scattering calculations based on the ray-tracing technique. In particular, a 3x2 refraction or transmission matrix is derived to account for the inhomogeneity of the refracted wave in an absorbing medium. An asymptotic solution that completely includes the effect of medium absorption on Fresnel coefficients is obtained for the scattering properties of a general polyhedral particle. Numerical results are presented for hexagonal plates and columns with both preferred and random orientations.

Yang, Ping↗

Microwave Polarized Signatures Generated within Cloud Systems: SSM/I Observations Interpreted with Radiative Transfer Simulations

Special Sensor Microwave /Imager (SSM/I) observations in cloud systems are studied over the tropics. Over optically thick cloud systems, presence of polarized signatures at 37 and 85 GHz is evidenced and analyzed with the help of cloud top temperature and optical thickness extracted from visible and IR satellite observations. Scattering signatures at 85 GHz (TbV(85) less than or = 250 K) are associated with polarization differences greater than or = 6 K, approx. 50%, of the time over ocean and approx. 40% over land. In addition. over thick clouds the polarization difference at 37 GHz is rarely negligible. The polarization differences at 37 and 85 GHz do not stem from the surface but are generated in regions of relatively homogeneous clouds having high liquid water content. To interpret the observations, a radiative transfer model that includes the scattering by non-spherical particles is developed. based on the T-matrix approach and using the doubling and adding method. In addition to handling randomly and perfectly oriented particles, this model can also simulate the effect of partial orientation of the hydrometeors. Microwave brightness temperatures are simulated at SSM/I frequencies and are compared with the observations. Polarization differences of approx. 2 K can be simulated at 37 GHz over a rain layer, even using spherical drops. The polarization difference is larger for oriented non-spherical particles. The 85 GHz simulations are very sensitive to the ice phase of the cloud. Simulations with spherical particles or with randomly oriented non-spherical ice particles cannot replicate the observed polarization differences. However, with partially oriented non-spherical particles, the observed polarized signatures at 85 GHz are explained, and the sensitivity of the scattering characteristics to the particle size, asphericity, and orientation is analyzed. Implications on rain and ice retrievals are discussed.

Prigent, Catherine↗

Constraints on PSC Particle Microphysics Derived From Lidar Observations

Based on extensive T-matrix computations of light scattering by polydispersions of randomly oriented, rotationally symmetric nonspherical particles, we analyze existing lidar observations of polar stratospheric clouds (PSCs) and derive several constraints on PSC particle microphysical properties. We show that sharp-edged nonspherical particles (finite circular cylinders) exhibit less variability of lidar backscattering characteristics with particle size and aspect ratio than particles with smooth surfaces (spheroids). For PSC particles significantly smaller than the wavelength, the backscatter color index Alpha and the depolarization color index Beta are essentially shape-independent. Observations for type Ia PSCs can be reproduced by spheroids with aspect ratios larger than 1.2, oblate cylinders with diameter-to-length ratios greater than 1.6, and prolate cylinders with length-to-diameter ratios greater than 1.4. The effective equal-volume-sphere radius for type la PSCs is about 0.8 microns or larger. Type Ib PSCs are likely to be composed of spheres or nearly spherical particles with effective radii smaller than 0.8 microns. Observations for type II PSCs are consistent with large ice crystals (effective radius greater than 1 micron modeled as cylinders or prolate spheroids.

Liu, Li↗

Radiation Force Caused by Scattering, Absorption, and Emission of Light by Nonspherical Particles

General formulas for computing the radiation force exerted on arbitrarily oriented and arbitrarily shaped nonspherical particles due to scattering, absorption, and emission of electromagnetic radiation are derived. For randomly oriented particles with a plane of symmetry, the formula for the average radiation force caused by the particle response to external illumination reduces to the standard Debye formula derived from the Lorenz-Mie theory, whereas the average radiation force caused by emission vanishes.

Mishchenko, Michael I.↗

Scattering by Randomly Oriented Thin Ice Disks with Moderate Equivalent-Sphere Size Parameters

We use the T-matrix method to compute the scattering matrix for randomly oriented circular ice cylinders with diameter-to-length ratios 1 and 20 and surface-equivalent-sphere size parameters up to 12. We show that wavelength-sized, sharp-edged ice plates with extreme diameter-to-length ratios possess the same scattering properties as smooth plate-like spheroids: their phase functions are similar to those of surface-equivalent compact particles, whereas all other elements of the scattering matrix are typical of Rayleigh scattering.

Zakharova, Nadia T.↗

Scattering Properties of Needle-Like and plate-like Ice Spheroids with Moderate Size Parameters

We use the current advanced version of the T-matrix method to compute the optical cross sections, the asymmetry parameter of the phase function, and the scattering matrix elements of ice spheroids with aspect ratios up to 20 and surface-equivalent-sphere size parameters up to 12. We demonstrate that plate-like and needle-like particles with moderate size parameters possess unique scattering properties: their asymmetry parameters and phase functions are similar to those of surface-equivalent spheres, whereas all other elements of the scattering matrix are typical of particles much smaller than the wavelength (Rayleigh scatterers). This result may have important implications for optical particle sizing and remote sensing of the terrestrial and planetary atmospheres.

Zakharova, Nadia T.↗

Asymptotic Solutions for Optical Properties of Large Particles with Strong Absorption

For scattering calculations involving nonspherical particles such as ice crystals, we show that the transverse wave condition is not applicable to the refracted electromagnetic wave in the context of geometric optics when absorption is involved. Either the TM wave condition (i.e., where the magnetic field of the refracted wave is transverse with respect to the wave direction) or the TE wave condition (i.e., where the electric field is transverse with respect to the propagating direction of the wave) may be assumed for the refracted wave in an absorbing medium to locally satisfy the electromagnetic boundary condition in the ray tracing calculation. The wave mode assumed for the refracted wave affects both the reflection and refraction coefficients. As a result, a nonunique solution for these coefficients is derived from the electromagnetic boundary condition. In this study we have identified the appropriate solution for the Fresnel reflection/refraction coefficients in light scattering calculation based on the ray tracing technique. We present the 3 x 2 refraction or transmission matrix that completely accounts for the inhomogeneity of the refracted wave in an absorbing medium. Using the Fresnel coefficients for an absorbing medium, we derive an asymptotic solution in an analytical format for the scattering properties of a general polyhedral particle. Numerical results are presented for hexagonal plates and columns with both preferred and random orientations. The asymptotic theory can produce reasonable accuracy in the phase function calculations in the infrared window region (wavelengths near 10 micron) if the particle size (in diameter) is on the order of 40 micron or larger. However, since strong absorption is assumed in the computation of the single-scattering albedo in the asymptotic theory, the single scattering albedo does not change with variation of the particle size. As a result, the asymptotic theory can lead to substantial errors in the computation of single-scattering albedo for small and moderate particle sizes. However, from comparison of the asymptotic results with the FDTD solution, it is expected that a convergence between the FDTD results and the asymptotic theory results can be reached when the particle size approaches 200 micron. We show that the phase function at side-scattering and backscattering angles is insensitive to particle shape if the random orientation condition is assumed. However, if preferred orientations are assumed for particles, the phase function has a strong dependence on scattering azimuthal angle. The single-scattering albedo also shows very strong dependence on the inclination angle of incident radiation with respect to the rotating axis for the preferred particle orientations.

Yang, Ping↗

Toward Creating A Global Retrospective Climatology of Aerosol Properties

Tropospheric aerosols are thought to cause a significant direct and indirect climate forcing, but the magnitude of this forcing remains highly uncertain because of poor knowledge of global aerosol characteristics and their temporal changes. The standard long-term global product, the one-channel Advanced Very-High-Resolution Radiometer (AVHRR) aerosol optical thickness over the ocean, relies on a single predefined aerosol model and can be inaccurate in many cases. Furthermore, it provides no information on aerosol column number density, thus making it impossible to estimate the indirect aerosol effect on climate. Total Ozone Mapping Spectrometer (TOMS) data can be used to detect absorbing aerosols over land, but are insensitive to aerosols located below one kilometer. It is thus clear that innovative approaches must be employed in order to extract a more quantitative and accurate aerosol climatology from available satellite and other measurements, thus enabling more reliable estimates of the direct and indirect aerosol forcings. The Global Aerosol Climatology Project (GACP) was established in 1998 as part of the Global Energy and Water Cycle Experiment (GEWEX). Its main objective is to analyze satellite radiance measurements and field observations to infer the global distribution of aerosols, their properties, and their seasonal and interannual variations. The overall goal is to develop advanced global aerosol climatologies for the period of satellite data and to make the aerosol climatologies broadly available through the GACP web site.

Curran, Robert J.↗

Sensitivity Analysis of Aerosol Retrieval Algorithms Using Two-Channel Satellite Radiance Data

We discuss the methodology of interpreting channel 1 and 2 AVHRR radiance data to retrieve tropospheric aerosol properties over the ocean and describe a detailed analysis of the sensitivity of monthly average retrievals to the assumed aerosol models. We use real AVHRR data and accurate numerical techniques for computing single and multiple scattering and spectral absorption of light in the vertically inhomogeneous atmosphere-ocean system. Our analysis shows that two-channel algorithms can provide significantly more accurate retrievals of the aerosol optical thickness than one-channel algorithms and that imperfect cloud screening is the largest source of errors in the retrieved optical thickness. Both underestimating and overestimating aerosol absorption as well as strong variability of the aerosol refractive index may lead to regional and/or seasonal biases in optical thickness retrievals. The Angstrom exponent appears to be the most invariant aerosol size characteristic and should be retrieved along with optical thickness as the second aerosol parameter.

Geogdzhayev, Igor V.↗

Bidirectional Reflectance of Flat, Optically Thick Particulate Layers: An Efficient Radiative Transfer Solution and Applications to Snow and Soil Surfaces

We describe a simple and highly efficient and accurate radiative transfer technique for computing bidirectional reflectance of a macroscopically flat scattering layer composed of nonabsorbing or weakly absorbing, arbitrarily shaped, randomly oriented and randomly distributed particles. The layer is assumed to be homogeneous and optically semi-infinite, and the bidirectional reflection function (BRF) is found by a simple iterative solution of the Ambartsumian's nonlinear integral equation. As an exact Solution of the radiative transfer equation, the reflection function thus obtained fully obeys the fundamental physical laws of energy conservation and reciprocity. Since this technique bypasses the computation of the internal radiation field, it is by far the fastest numerical approach available and can be used as an ideal input for Monte Carlo procedures calculating BRFs of scattering layers with macroscopically rough surfaces. Although the effects of packing density and coherent backscattering are currently neglected, they can also be incorporated. The FORTRAN implementation of the technique is available on the World Wide Web at http://ww,,v.giss.nasa.gov/-crmim/brf.html and can be applied to a wide range of remote sensing, engineering, and biophysical problems. We also examine the potential effect of ice crystal shape on the bidirectional reflectance of flat snow surfaces and the applicability of the Henyey-Greenstein phase function and the 6-Eddington approximation in calculations for soil surfaces.

Mishchenko, Michael I.↗