Engineering Papers⌕ Search

Engineering topics

Yang, Ping

Publications and source records attributed to Yang, Ping.

216 records · Page 12

Effect of the Inhomogeneity of Ice Crystals on Retrieving Ice Cloud Optical Thickness and Effective Particle Size

Spherical or spheroidal air bubbles are generally trapped in the formation of rapidly growing ice crystals. In this study the single-scattering properties of inhomogeneous ice crystals containing air bubbles are investigated. Specifically, a computational model based on an improved geometric-optics method (IGOM) has been developed to simulate the scattering of light by randomly oriented hexagonal ice crystals containing spherical or spheroidal air bubbles. A combination of the ray-tracing technique and the Monte Carlo method is used. The effect of the air bubbles within ice crystals is to smooth the phase functions, diminish the 22deg and 46deg halo peaks, and substantially reduce the backscatter relative to bubble-free particles. These features vary with the number, sizes, locations and shapes of the air bubbles within ice crystals. Moreover, the asymmetry factors of inhomogeneous ice crystals decrease as the volume of air bubbles increases. Cloud reflectance lookup tables were generated at wavelengths 0.65 m and 2.13 m with different air-bubble conditions to examine the impact of the bubbles on retrieving ice cloud optical thickness and effective particle size. The reflectances simulated for inhomogeneous ice crystals are slightly larger than those computed for homogenous ice crystals at a wavelength of 0.65 microns. Thus, the retrieved cloud optical thicknesses are reduced by employing inhomogeneous ice cloud models. At a wavelength of 2.13 microns, including air bubbles in ice cloud models may also increase the reflectance. This effect implies that the retrieved effective particle sizes for inhomogeneous ice crystals are larger than those retrieved for homogeneous ice crystals, particularly, in the case of large air bubbles.

Xie, Yu↗

Optical Properties of Ice Particles in Young Contrails

The single-scattering properties of four types of ice crystals (pure ice crystals, ice crystals with an internal mixture of ice and black carbon, ice crystals coated with black carbon, and soot coated with ice) in young contrails are investigated at wavelengths 0.65 and 2.13 micrometers using Mie codes from coated spheres. The four types of ice crystals have distinct differences in their single-scattering properties because of the embedded black carbon. The bulk scattering properties of young contrails consisting of the four types of ice crystals are further investigated by averaging their single-scattering properties over a typical ice particle size distribution found in young contrails. The effect of the radiative properties of the four types of ice particles on the Stokes parameters I, Q, U, and V is also investigated for different viewing zenith angles and relative azimuth angles with a solar zenith angle of 30 degrees using a vector radiative transfer model based on the adding-doubling technique. The Stokes parameters at a wavelength of 0.65 micrometers show pronounced differences for the four types of ice crystals. Those at a wavelength of 2.13 micrometers show similar variations with the viewing zenith angle and relative azimuth angle, but their values are noticeably different.

Hong, Gang↗

Atmospheric, Cloud, and Surface Parameters Retrieved from Satellite Ultra-spectral Infrared Sounder Measurements

An advanced retrieval algorithm with a fast radiative transfer model, including cloud effects, is used for atmospheric profile and cloud parameter retrieval. This physical inversion scheme has been developed, dealing with cloudy as well as cloud-free radiance observed with ultraspectral infrared sounders, to simultaneously retrieve surface, atmospheric thermodynamic, and cloud microphysical parameters. A fast radiative transfer model, which applies to the clouded atmosphere, is used for atmospheric profile and cloud parameter retrieval. A one-dimensional (1-d) variational multivariable inversion solution is used to improve an iterative background state defined by an eigenvector-regression-retrieval. The solution is iterated in order to account for non-linearity in the 1-d variational solution. This retrieval algorithm is applied to the MetOp satellite Infrared Atmospheric Sounding Interferometer (IASI) launched on October 19, 2006. IASI possesses an ultra-spectral resolution of 0.25 cm(exp -1) and a spectral coverage from 645 to 2760 cm(exp -1). Preliminary retrievals of atmospheric soundings, surface properties, and cloud optical/microphysical properties with the IASI measurements are obtained and presented.

Zhou, Daniel K.↗

Cloudy Sounding and Cloud-Top Height Retrieval From AIRS Alone Single Field-of-View Radiance Measurements

High-spectral resolution measurements from the Atmospheric Infrared Sounder (AIRS) onboard the EOS (Earth Observing System) Aqua satellite provide unique information about atmospheric state, surface and cloud properties. This paper presents an AIRS alone single field-of-view (SFOV) retrieval algorithm to simultaneously retrieve temperature, humidity and ozone profiles under all weather conditions, as well as cloud top pressure (CTP) and cloud optical thickness (COT) under cloudy skies. For optically thick cloud conditions the above-cloud soundings are derived, whereas for clear skies and optically thin cloud conditions the profiles are retrieved from 0.005 hPa down to the earth's surface. Initial validation has been conducted by using the operational MODIS (Moderate Resolution Imaging Spectroradiometer) product, ECMWF (European Center of Medium range Weather Forecasts) analysis fields and radiosonde observations (RAOBs). These inter-comparisons clearly demonstrate the potential of this algorithm to process data from 38 high-spectral infrared (IR) sounder instruments.

Weisz, Elisabeth↗

A Fast Infrared Radiative Transfer Model for Overlapping Clouds

A fast infrared radiative transfer model (FIRTM2) appropriate for application to both single-layered and overlapping cloud situations is developed for simulating the outgoing infrared spectral radiance at the top of the atmosphere (TOA). In FIRTM2 a pre-computed library of cloud reflectance and transmittance values is employed to account for one or two cloud layers, whereas the background atmospheric optical thickness due to gaseous absorption can be computed from a clear-sky radiative transfer model. FIRTM2 is applicable to three atmospheric conditions: 1) clear-sky, 2) single-layered ice or water cloud, and 3) two simultaneous cloud layers in a column (e.g., ice cloud overlying water cloud). Moreover, FIRTM2 outputs the derivatives (i.e., Jacobians) of the TOA brightness temperature with respect to cloud optical thickness and effective particle size. Sensitivity analyses have been carried out to assess the performance of FIRTM2 for two spectral regions, namely the longwave (LW) band (587.3 - 1179.5/cm) and the short-to-medium wave (SMW) band (1180.1 - 2228.9/cm). The assessment is carried out in terms of brightness temperature differences (BTD) between FIRTM2 and the well-known discrete ordinates radiative transfer model (DISORT), henceforth referred to as BTD (F-D). The BTD (F-D) values for single-layered clouds are generally less than 0.8 K. For the case of two cloud layers (specifically ice cloud over water cloud), the BTD(F-D) values are also generally less than 0.8 K except for the SMW band for the case of a very high altitude (>15 km) cloud comprised of small ice particles. Note that for clear-sky atmospheres, FIRTM2 reduces to the clear-sky radiative transfer model that is incorporated into FIRTM2, and the errors in this case are essentially those of the clear-sky radiative transfer model.

Niu, Jianguo↗

Potential Nighttime Contamination of CERES Clear-sky Field of View by Optically Thin Cirrus during the CRYSTAL-FACE Campaign

We investigate the outgoing broadband longwave (LW, 5 to approx. 200 microns) and window (WIN, 8 to approx. 12 microns) channel radiances at the top of atmosphere (TOA) under clear-sky conditions, using data acquired by the Cloud and the Earth s Radiant Energy System (CERES) and Moderate-Resolution Imaging Spectroradiometer (MODIS) instruments onboard the NASA Terra satellite platform. In this study, detailed analyses are performed on the CERES Single Scanner Footprint TOA/Surface Fluxes and Clouds product to understand the radiative effect of thin cirrus. The data are acquired over the Florida area during the Cirrus Regional Study of Tropical Anvils and Cirrus Layers Florida Area Cirrus Experiment (CRYSTAL-FACE) field program. Of particular interest is the anisotropy associated with the radiation field. Measured CERES broadband radiances are compared to those obtained from rigorous radiative transfer simulations. Analysis of results from this comparison indicates that the simulated radiances tend to be larger than their measured counterparts, with differences ranging from 2.1% to 8.3% for the LW band and from 1.7% to 10.6% for the WIN band. The averaged difference in radiance is approximately 4% for both the LW and WIN channels. A potential cause for the differences could be the presence of thin cirrus (i.e., optically thin ice clouds with visible optical thicknesses smaller than approximately 0.3). The detection and quantitative analysis of these thin cirrus clouds are challenging even with sophisticated multispectral instruments. While large differences in radiance between the CERES observations and the theoretical calculations are found, the corresponding difference in the anisotropic factors is very small (0.2%). Furthermore, sensitivity studies show that the influence due to a 1 K bias of the surface temperature on the errors of the LW and WIN channel radiances is of the same order as that associated with a 2% bias of the surface emissivity. The LW and WIN errors associated with a 5% bias of water vapor amount in the lower atmosphere in conjunction with a 50% bias of water vapor amount in the upper atmosphere is similar to that of a 1 K bias of the vertical temperature profile. Even with the uncertainties considered for these various factors, the simulated LW and WIN radiances are still larger than the observed radiances if thin cirrus clouds are excluded.

Lee, Yong-Keun↗

Impact of Cirrus Crystal Shape on Solar Spectral Irradiance: A Case Study for Subtropical Cirrus

Profiles of in situ measurements of ice crystal size distribution of subtropical cirrus were used to calculate solar spectral irradiances above and below the clouds. Spheres and nonspherical ice crystal habits (columns, hollows, plates, bullets, and aggregates) were assumed in the calculations. The simulation results were compared to irradiance measurements from the NASA Solar Spectral Flux Radiometer. The microphysical and radiation data were collected by three aircraft during CRYSTAL-FACE. Two cirrus cases (optical thickness of about 1 and 7) from two mission dates (26 and 23 July 2002) were investigated in detail. The measured downwelling and upwelling irradiance spectra above the cirrus could mostly be reproduced by the radiation model to within +/- 5-10% for most ice crystal habits. Below the cirrus the simulations disagreed with the measured irradiances due to surface albedo variability along the flight track, and nonoptimal colocation between the microphysical and irradiance measurements. The impact of shape characteristics of the crystals was important for the reflected irradiances above the optically thin cirrus, especially for small solar zenith angles, because in this case single-scattering dominated the solar radiation field. For the cirrus of moderate optical thickness the enhanced multiple scattering tended to diminish particular shape features caused by nonspherical single-scattering. Within the ice absorption bands the shape-related differences in the absorption characteristics of the individual nonspherical ice crystals were amplified if multiple scattering prevailed. Furthermore, it was found that below the cloud the shape sensitivity of the downwelling irradiance spectra is larger compared to the nonsphericity effects on reflected irradiances above the cirrus. Finally, it was shown that the calculated cirrus solar radiative forcing could vary by as much as 26% depending on the ice crystal habit.

Wendisch, Manfred↗

Inference of Ice Cloud Properties from High-spectral Resolution Infrared Observations

The theoretical basis is explored for inferring the microphysical properties of ice crystal from high-spectral resolution infrared observations. A radiative transfer model is employed to simulate spectral radiances to address relevant issues. The extinction and absorption efficiencies of individual ice crystals, assumed as hexagonal columns for large particles and droxtals for small particles, are computed from a combination of the finite- difference time-domain (FDTD) technique and a composite method. The corresponding phase functions are computed from a combination of FDTD and an improved geometric optics method (IGOM). Bulk scattering properties are derived by averaging the single- scattering properties of individual particles for 30 particle size distributions developed from in situ measurements and for additional four analytical Gamma size distributions for small particles. The non-sphericity of ice crystals is shown to have a significant impact on the radiative signatures in the infrared (IR) spectrum; the spherical particle approximation for inferring ice cloud properties may result in an overest&ation of the optical thickness and an inaccurate retrieval of effective particle size. Furthermore, we show that the error associated with the use of the Henyey-Greenstein phase function can be as larger as 1 K in terms of brightness temperature for larger particle effective size at some strong scattering wavenumbers. For small particles, the difference between the two phase functions is much less, with brightness temperatures generally differing by less than 0.4 K. The simulations undertaken in this study show that the slope of the IR brightness temperature spectrum between 790-960/cm is sensitive to the effective particle size. Furthermore, a strong sensitivity of IR brightness temperature to cloud optical thickness is noted within the l050-1250/cm region. Based on this spectral feature, a technique is presented for the simultaneous retrieval of the visible optical thickness and effective particle size from high spectral resolution infrared data under ice cloudy con&tion. The error analysis shows that the uncertainty of the retrieved optical thickness and effective particle size has a small range of variation. The error for retrieving particle size in conjunction with an uncertainty of 5 K in cloud'temperature, or a surface temperature uncertainty of 2.5 K, is less than 15%. The corresponding e m r in the uncertainty of optical thickness is within 5-2096, depending on the value of cloud optical thickness. The applicability of the technique is demonstrated using the aircraft-based High- resolution Interferometer Sounder (HIS) data from the Subsonic Aircraft: Contrail and Cloud Effects Special Study (SUCCESS) in 1996 and the First ISCCP Regional Experiment - Arctic Clouds Experiment (FIRE-ACE) in 1998.

Huang, Hung-Lung↗

The Influence of Thermodynamic Phase on the Retrieval of Mixed-Phase Cloud Microphysical and Optical Properties in the Visible and Near Infrared Region

Cloud microphysical and optical properties are inferred from the bidirectional reflectances simulated for a single-layered cloud consisting of an external mixture of ice particles and liquid droplets. The reflectances are calculated with a rigorous discrete ordinates radiative transfer model and are functions of the cloud effective particle size, the cloud optical thickness, and the values of the ice fraction in the cloud (i.e., the ratio of ice water content to total water content). In the present light scattering and radiative transfer simulations, the ice fraction is assumed to be vertically homogeneous; the habit (shape) percentage as a function of ice particle size is consistent with that used for the Moderate Resolution Imaging Spectroradiometer (MODIS) operational (Collection 4 and earlier) cloud products; and the surface is assumed to be Lambertian with an albedo of 0.03. Furthermore, error analyses pertaining to the inference of the effective particle sizes and optical thicknesses of mixed-phase clouds are performed. Errors are calculated with respect to the assumption of a cloud containing solely liquid or ice phase particles. The analyses suggest that the effective particle size inferred for a mixed-phase cloud can be underestimated (or overestimated) if pure liquid phase (or pure ice phase) is assumed for the cloud, whereas the corresponding cloud optical thickness can be overestimated (or underestimated).

Lee, Joonsuk↗

Remote Sensing of Liquid Water and Ice Cloud Optical Thickness and Effective Radius in the Arctic: Application of Airborne Multispectral MAS Data

A multispectral scanning spectrometer was used to obtain measurements of the reflection function and brightness temperature of clouds, sea ice, snow, and tundra surfaces at 50 discrete wavelengths between 0.47 and 14.0 microns. These observations were obtained from the NASA ER-2 aircraft as part of the FIRE Arctic Clouds Experiment, conducted over a 1600 x 500 km region of the north slope of Alaska and surrounding Beaufort and Chukchi Seas between 18 May and 6 June 1998. Multispectral images of the reflection function and brightness temperature in 11 distinct bands of the MODIS Airborne Simulator (MAS) were used to derive a confidence in clear sky (or alternatively the probability of cloud), shadow, and heavy aerosol over five different ecosystems. Based on the results of individual tests run as part of the cloud mask, an algorithm was developed to estimate the phase of the clouds (water, ice, or undetermined phase). Finally, the cloud optical thickness and effective radius were derived for both water and ice clouds that were detected during one flight line on 4 June. This analysis shows that the cloud mask developed for operational use on MODIS, and tested using MAS data in Alaska, is quite capable of distinguishing clouds from bright sea ice surfaces during daytime conditions in the high Arctic. Results of individual tests, however, make it difficult to distinguish ice clouds over snow and sea ice surfaces, so additional tests were added to enhance the confidence in the thermodynamic phase of clouds over the Beaufort Sea. The cloud optical thickness and effective radius retrievals used 3 distinct bands of the MAS, with the newly developed 1.62 and 2.13 micron bands being used quite successfully over snow and sea ice surfaces. These results are contrasted with a MODIS-based algorithm that relies on spectral reflectance at 0.87 and 2.13 micron.

King, Michael D.↗

The Development of Midlatitude Cirrus Models for MODIS Using FIRE-I, FIRE-II, and ARM In Situ Data

Detailed in situ data from cirrus clouds have been collected during dedicated field Campaigns, but the use of the size and habit distribution data has been lagging in the development of more realistic cirrus scattering models. In this study, the authors examine the use of in situ cirrus data collected during three field campaigns to develop more realistic midlatitude cirrus microphysical models. Data are used from the First International Satellite Cloud Climatology Project (ISCCP) Regional Experiment (FIRE)-I (1986) and FIRE-II (1991) campaigns and from a recent Atmospheric Radiation Measurement (ARM) Program campaign held in March-April of 2000. The microphysical models are based on measured vertical distributions of both particle size and particle habit and are used to develop new scattering models for a suite of moderate-resolution imaging spectroradiometer (MODIS) bands spanning visible. near-infrared, and infrared wavelengths. The sensitivity of the resulting scattering properties to the underlying assumptions of the assumed particle size and habit distributions are examined. It is found that the near-infrared bands are sensitive not only to the discretization of the size distribution but also to the assumed habit distribution. In addition. the results indicate that the effective diameter calculated from a given size distribution tends to be sensitive to the number of size bins that are used to discretize the data and also to the ice-crystal habit distribution.

Nasiri, Shaima L.↗

Parameterization of Cloud Optical Properties for a Mixture of Ice Particles for use in Atmospheric Models

Based on the single-scattering optical properties that are pre-computed using an improve geometric optics method, the bulk mass absorption coefficient, single-scattering albedo, and asymmetry factor of ice particles have been parameterized as a function of the mean effective particle size of a mixture of ice habits. The parameterization has been applied to compute fluxes for sample clouds with various particle size distributions and assumed mixtures of particle habits. Compared to the parameterization for a single habit of hexagonal column, the solar heating of clouds computed with the parameterization for a mixture of habits is smaller due to a smaller cosingle-scattering albedo. Whereas the net downward fluxes at the TOA and surface are larger due to a larger asymmetry factor. The maximum difference in the cloud heating rate is approx. 0.2 C per day, which occurs in clouds with an optical thickness greater than 3 and the solar zenith angle less than 45 degrees. Flux difference is less than 10 W per square meters for the optical thickness ranging from 0.6 to 10 and the entire range of the solar zenith angle. The maximum flux difference is approximately 3%, which occurs around an optical thickness of 1 and at high solar zenith angles.

Chou, Ming-Dah↗

Parameterization of Shortwave Cloud Optical Properties for a Mixture of Ice Particle Habits for use in Atmospheric Models

Based on the single-scattering optical properties pre-computed with an improved geometric optics method, the bulk absorption coefficient, single-scattering albedo, and asymmetry factor of ice particles have been parameterized as a function of the effective particle size of a mixture of ice habits, the ice water amount, and spectral band. The parameterization has been applied to computing fluxes for sample clouds with various particle size distributions and assumed mixtures of particle habits. It is found that flux calculations are not overly sensitive to the assumed particle habits if the definition of the effective particle size is consistent with the particle habits that the parameterization is based. Otherwise, the error in the flux calculations could reach a magnitude unacceptable for climate studies. Different from many previous studies, the parameterization requires only an effective particle size representing all ice habits in a cloud layer, but not the effective size of individual ice habits.

Chou, Ming-Dah↗

Seasonal and Global Variations of Water Vapor and High Clouds Observed with MODIS near-IR Channels

The NASA Moderate Resolution Imaging Spectrometer (MODIS) on the Terra Spacecraft has been collecting scientific data since February of 2000. MODIS is a major facility instrument for remote sensing of the atmosphere, land surfaces, and ocean color. On the MODIS instruments, there are five channels located within and around the .0.94 micron water vapor band absorption region for remote sensing of atmospheric water vapor. There is also a channel located at 1.375 micron for detecting thin cirrus clouds. We will describe the basic principles for using these near-IR channels for remote sensing of water vapor and high clouds. Based on our analysis of two years# measurements with these channels, we have found that reliable observations of water vapor and high clouds on regional and global scales can be made. We will present results on daily, seasonal and annual variations of water vapor and high clouds.

Gao, Bo-Cai↗

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↗

Retrieval of Polar Stratospheric Cloud Microphysical Properties From Lidar Measurements: Dependence on Particle Shape Assumptions

Knowledge of particle sizes and number densities of polar stratospheric clouds (PSCs) is highly important, because they are critical parameters for the modeling of the ozone chemistry of the stratosphere. In situ measurements of PSC particles are rare. the main instrument for the accumulation of PSC data are lidar systems. Therefore the derivation of some microphysical properties of PSCS from the optical parameters measured by lidars would be highly beneficial for ozone research. Inversion of lidar data obtained in the presence of PSCs formed from crystalline particles type 11 and the various nitric acid tri Ydrrate (NAT) types cannot be easily accomplished, because a suitable scattering theory for small faceted crystals has not been readily available tip to now. As a consequence, the T-matrix method is commonly used for the interpretation of these PSC lidar data. Here the assumption is made that the optical properties of an ensemble of spheroids resemble those of crystalline PSCs, and microphysical properties of the PSC are inferred from the optical signatures of the PSC at two or more wavelengths. The problem with the T-matrix approach is that the assumption of spheroidal instead of faceted particles can lead to dramatically wrong results: Usually cloud particle properties are deduced from analysis of lidar profiles of backscatter ratio and depolarization ratio. The particle contribution to the backscatter ratio is given by the product of the particle number density and the backscattering cross section. The latter is proportional to the value of the particle's scattering phase function at 180 degrees scattering angle. At 180 degrees however, the phase functions of rough, faceted crystals and of spheroids with same maximum dimension differ by a factor of 6. From this it follows that for a PSC consisting of faceted crystals, the particle number density is underestimated by roughly the same factor if spheroidal particles are unrealistically assumed. We are currently developing a retrieval technique for determining the microphysical parameters of crystalline PSCs that takes into account the faceted shape of the PSC particles. This approach utilizes finite-difference time-domain (FDTD) calculations of particle optical properties. The accuracy and the free choice of the shape of the scattering particle make the FDTD technique a promising tool for the inversion of PSC lidar data. A first comparison of FDTD and T-matrix calculations will be presented.

Reichardt, Susanne↗

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↗

Radiative Properties of Cirrus Clouds in the Infrared (8-13 microns) Spectral Region

Atmospheric radiation in the infrared (IR) 8-13 microns spectral region contains a wealth of information that is very useful for the retrieval of ice cloud properties from aircraft or space-borne measurements. To provide the scattering and absorption properties of nonspherical ice crystals that are fundamental to the IR retrieval implementation, we use the finite-difference time domain (FDTD) method to solve for the extinction efficiency, single-scattering albedo, and the asymmetry parameter of the phase function for ice crystals smaller than 40 microns. For particles larger than this size, the improved geometric optics method (IGOM) can be employed to calculate the asymmetry parameter with an acceptable accuracy, provided that we properly account for the inhomogeneity of the refracted wave due to strong absorption inside the ice particle. A combination of the results computed from the two methods provides the asymmetry parameter for the entire practical range of particle sizes between 1 micron and 10000 microns over wavelengths ranging from 8 microns to 13 microns. For the extinction and absorption efficiency calculations, several methods including the IGOM, Mie solution for equivalent spheres (MSFES), and the anomalous diffraction theory (ADT) can lead to a substantial discontinuity in comparison with the FDTD solutions for particle sizes on the order of 40 microns. To overcome this difficulty, we have developed a novel approach called the stretched scattering potential method (SSPM). For the IR 8-13 microns spectral region, we show that SSPM is a more accurate approximation than ADT, MSFES, and IGOM. The SSPM solution can be further refined numerically. Through a combination of the FDTD and SSPM, we have computed the extinction and absorption efficiency for hexagonal ice crystals with sizes ranging from 1 to 10000 microns at 12 wavelengths between 8 and 13 microns Calculations of the cirrus bulk scattering and absorption properties are performed for 30 size distributions obtained from various field campaigns for midlatitude and tropical cirrus cloud systems. Parameterization of these bulk scattering properties is carried out by using second-order polynomial functions for the extinction efficiency and the single-scattering albedo and the power law expression for the asymmetry parameter. We note that the volume-normalized extinction coefficient can be separated into two parts: one is inversely proportional to effective size and is independent of wavelength, and the other is the wavelength-dependent effective extinction efficiency. Unlike conventional parameterization efforts, the present parameterization scheme is more accurate because only the latter part of the volume-normalized extinction coefficient is approximated in terms of an analytical expression. After averaging over size distribution, the single-scattering albedo is shown to decrease with an increase in effective size for wavelengths shorter than 10.0 microns whereas the opposite behavior is observed for longer wavelengths. The variation of the asymmetry parameter as a function of effective size is substantial when the effective size is smaller than 50 microns. For effective sizes larger than 100 microns, the asymmetry parameter approaches its asymptotic value. The results derived in this study can be useful to remote sensing applications involving IR window bands under cirrus cloud conditions.

Yang, Ping↗