Engineering PapersSearch

Engineering topics

Baum, Bryan A.

Publications and source records attributed to Baum, Bryan A..

At least 19 records

Toward Global Harmonization of Derived Cloud Products

Formerly known as the Cloud Retrieval Evaluation Workshop (CREW; see the list of acronyms used in this paper below) group (Roebeling et al. 2013, 2015), the International Cloud Working Group (ICWG) was created and endorsed during the 42nd Meeting of CGMS. The CGMS-ICWG provides a forum for space agencies to seek coherent progress in science and applications and also to act as a bridge between space agencies and the cloud remote sensing and applications community. The ICWG plans to serve as a forum to exchange and enhance knowledge on state-of-the-art cloud parameter retrievals algorithms, to stimulate support for training in the use of cloud parameters, and to encourage space agencies and the cloud remote sensing community to share knowledge. The ICWG plans to prepare recommendations to guide the direction of future research-for example, on observing severe weather events or on process studies-and to influence relevant programs of the WMO, WCRP, GCOS, and the space agencies.

cloud parameters

Degree of Ice Particle Surface Roughness Inferred from Polarimetric Observations

The degree of surface roughness of ice particles within thick, cold ice clouds is inferred from multidirectional, multi-spectral satellite polarimetric observations over oceans, assuming a column-aggregate particle habit. An improved roughness inference scheme is employed that provides a more noise-resilient roughness estimate than the conventional best-fit approach. The improvements include the introduction of a quantitative roughness parameter based on empirical orthogonal function analysis and proper treatment of polarization due to atmospheric scattering above clouds. A global 1-month data sample supports the use of a severely roughened ice habit to simulate the polarized reflectivity associated with ice clouds over ocean. The density distribution of the roughness parameter inferred from the global 1- month data sample and further analyses of a few case studies demonstrate the significant variability of ice cloud single-scattering properties. However, the present theoretical results do not agree with observations in the tropics. In the extra-tropics, the roughness parameter is inferred but 74% of the sample is out of the expected parameter range. Potential improvements are discussed to enhance the depiction of the natural variability on a global scale.

surface

Assessment of the Accuracy of the Conventional Ray-Tracing Technique: Implications in Remote Sensing and Radiative Transfer Involving Ice Clouds.

A fundamental problem in remote sensing and radiative transfer simulations involving ice clouds is the ability to compute accurate optical properties for individual ice particles. While relatively simple and intuitively appealing, the conventional geometric-optics method (CGOM) is used frequently for the solution of light scattering by ice crystals. Due to the approximations in the ray-tracing technique, the CGOM accuracy is not well quantified. The result is that the uncertainties are introduced that can impact many applications. Improvements in the Invariant Imbedding T-matrix method (II-TM) and the Improved Geometric-Optics Method (IGOM) provide a mechanism to assess the aforementioned uncertainties. The results computed by the II-TMþIGOM are considered as a benchmark because the IITM solves Maxwell's equations from first principles and is applicable to particle size parameters ranging into the domain at which the IGOM has reasonable accuracy. To assess the uncertainties with the CGOM in remote sensing and radiative transfer simulations, two independent optical property datasets of hexagonal columns are developed for sensitivity studies by using the CGOM and the II-TMþIGOM, respectively. Ice cloud bulk optical properties obtained from the two datasets are compared and subsequently applied to retrieve the optical thickness and effective diameter from Moderate Resolution Imaging Spectroradiometer (MODIS) measurements. Additionally, the bulk optical properties are tested in broadband radiative transfer (RT) simulations using the general circulation model (GCM) version of the Rapid Radiative Transfer Model (RRTMG) that is adopted in the National Center for Atmospheric Research (NCAR) Community Atmosphere Model (CAM, version 5.1). For MODIS retrievals, the mean bias of uncertainties of applying the CGOM in shortwave bands (0.86 and 2.13 micrometers) can be up to 5% in the optical thickness and as high as 20% in the effective diameter, depending on cloud optical thickness and effective diameter. In the MODIS infrared window bands centered at 8.5, 11, and 12 micrometers biases in the optical thickness and effective diameter are up to 12% and 10%, respectively. The CGOM-based simulation errors in ice cloud radiative forcing calculations are on the order of 10Wm(exp 2).

remote sensing

Influence of Ice Particle Surface Roughening on the Global Cloud Radiative Effect

Ice clouds influence the climate system by changing the radiation budget and large-scale circulation. Therefore, climate models need to have an accurate representation of ice clouds and their radiative effects. In this paper, new broadband parameterizations for ice cloud bulk scattering properties are developed for severely roughened ice particles. The parameterizations are based on a general habit mixture that includes nine habits (droxtals, hollow/solid columns, plates, solid/hollow bullet rosettes, aggregate of solid columns, and small/large aggregates of plates). The scattering properties for these individual habits incorporate recent advances in light-scattering computations. The influence of ice particle surface roughness on the ice cloud radiative effect is determined through simulations with the Fu-Liou and the GCM version of the Rapid Radiative Transfer Model (RRTMG) codes and the National Center for Atmospheric Research Community Atmosphere Model (CAM, version 5.1). The differences in shortwave (SW) and longwave (LW) radiative effect at both the top of the atmosphere and the surface are determined for smooth and severely roughened ice particles. While the influence of particle roughening on the single-scattering properties is negligible in the LW, the results indicate that ice crystal roughness can change the SW forcing locally by more than 10 W m(exp −2) over a range of effective diameters. The global-averaged SW cloud radiative effect due to ice particle surface roughness is estimated to be roughly 1-2 W m(exp −2). The CAM results indicate that ice particle roughening can result in a large regional SW radiative effect and a small but nonnegligible increase in the global LW cloud radiative effect.

clouds

Spectrally Consistent Scattering, Absorption, and Polarization Properties of Atmospheric Ice Crystals at Wavelengths from 0.2 to 100 um

A data library is developed containing the scattering, absorption, and polarization properties of ice particles in the spectral range from 0.2 to 100 microns. The properties are computed based on a combination of the Amsterdam discrete dipole approximation (ADDA), the T-matrix method, and the improved geometric optics method (IGOM). The electromagnetic edge effect is incorporated into the extinction and absorption efficiencies computed from the IGOM. A full set of single-scattering properties is provided by considering three-dimensional random orientations for 11 ice crystal habits: droxtals, prolate spheroids, oblate spheroids, solid and hollow columns, compact aggregates composed of eight solid columns, hexagonal plates, small spatial aggregates composed of 5 plates, large spatial aggregates composed of 10 plates, and solid and hollow bullet rosettes. The maximum dimension of each habit ranges from 2 to 10,000 microns in 189 discrete sizes. For each ice crystal habit, three surface roughness conditions (i.e., smooth, moderately roughened, and severely roughened) are considered to account for the surface texture of large particles in the IGOM applicable domain. The data library contains the extinction efficiency, single-scattering albedo, asymmetry parameter, six independent nonzero elements of the phase matrix (P11, P12, P22, P33, P43, and P44), particle projected area, and particle volume to provide the basic single-scattering properties for remote sensing applications and radiative transfer simulations involving ice clouds. Furthermore, a comparison of satellite observations and theoretical simulations for the polarization characteristics of ice clouds demonstrates that ice cloud optical models assuming severely roughened ice crystals significantly outperform their counterparts assuming smooth ice crystals.

Optical properties

Comparison of PARASOL Observations with Polarized Reflectances Simulated Using Different Ice Habit Mixtures

Insufficient knowledge of the habit distribution and the degree of surface roughness of ice crystals within ice clouds is a source of uncertainty in the forward light scattering and radiative transfer simulations required in downstream applications involving these clouds. The widely used MODerate Resolution Imaging Spectroradiometer (MODIS) Collection 5 ice microphysical model assumes a mixture of various ice crystal shapes with smooth-facets except aggregates of columns for which a moderately rough condition is assumed. When compared with PARASOL (Polarization and Anisotropy of Reflectances for Atmospheric Sciences coupled with Observations from a Lidar) polarized reflection data, simulations of polarized reflectance using smooth particles show a poor fit to the measurements, whereas very rough-faceted particles provide an improved fit to the polarized reflectance. In this study a new microphysical model based on a mixture of 9 different ice crystal habits with severely roughened facets is developed. Simulated polarized reflectance using the new ice habit distribution is calculated using a vector adding-doubling radiative transfer model, and the simulations closely agree with the polarized reflectance observed by PARASOL. The new general habit mixture is also tested using a spherical albedo differences analysis, and surface roughening is found to improve the consistency of multi-angular observations. It is suggested that an ice model incorporating an ensemble of different habits with severely roughened surfaces would potentially be an adequate choice for global ice cloud retrievals.

Cole, Benjamin H.

Multilayer Cloud Detection with the MODIS Near-Infrared Water Vapor Absorption Band

Data Collection 5 processing for the Moderate Resolution Imaging Spectroradiometer (MODIS) onboard the NASA Earth Observing System EOS Terra and Aqua spacecraft includes an algorithm for detecting multilayered clouds in daytime. The main objective of this algorithm is to detect multilayered cloud scenes, specifically optically thin ice cloud overlying a lower-level water cloud, that presents difficulties for retrieving cloud effective radius using single layer plane-parallel cloud models. The algorithm uses the MODIS 0.94 micron water vapor band along with CO2 bands to obtain two above-cloud precipitable water retrievals, the difference of which, in conjunction with additional tests, provides a map of where multilayered clouds might potentially exist. The presence of a multilayered cloud results in a large difference in retrievals of above-cloud properties between the CO2 and the 0.94 micron methods. In this paper the MODIS multilayered cloud algorithm is described, results of using the algorithm over example scenes are shown, and global statistics for multilayered clouds as observed by MODIS are discussed. A theoretical study of the algorithm behavior for simulated multilayered clouds is also given. Results are compared to two other comparable passive imager methods. A set of standard cloudy atmospheric profiles developed during the course of this investigation is also presented. The results lead to the conclusion that the MODIS multilayer cloud detection algorithm has some skill in identifying multilayered clouds with different thermodynamic phases

Wind, Galina

Introduction to MODIS Cloud Products

The Earth's radiative energy balance and hydrological cycle are fundamentally coupled with the distribution and properties of clouds. Therefore, the ability to remotely infer cloud properties and their variation in space and time is crucial for establishing climatologies as a reference for validation of present-day climate models and in assessing future climate change. Remote cloud observations also provide data sets useful for testing and improving cloud model physics, and for assimilation into numerical weather prediction models. The MODerate Resolution Imaging Spectroradiometer (MODIS) imagers on the Terra and Aqua Earth Observing System (EOS) platforms provide the capability for globally retrieving these properties using passive solar reflectance and infrared techniques. In addition to providing measurements similar to those offered on a suite of historical operational weather platforms such as the Advanced Very High Resolution Radiometer (AVHRR), the High-resolution Infrared Radiation Sounder (HIRS), and the Geostationary Operational Environmental Satellite (GOES), MODIS provides additional spectral and/or spatial resolution in key atmospheric bands, along with on-board calibration, to expand the capability for global cloud property retrievals. The core MODIS operational cloud products include cloud top pressure, thermodynamic phase, optical thickness, particle size, and water path, and are derived globally at spatial resolutions of either 1- or 5-km (referred to as Level-2 or pixel-level products). In addition, the MODIS atmosphere team (collectively providing cloud, aerosol, and clear sky products) produces a combined gridded product (referred to as Level-3) aggregated to a 1 equal-angle grid, available for daily, eight-day, and monthly time periods. The wealth of information available from these products provides critical information for climate studies as well as the continuation and improved understanding of existing satellite-based cloud climatologies derived from heritage instruments. This chapter provides an overview of the MODIS Level-2 and -3 operational cloud products.

Baum, Bryan A.

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

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

Combining Satellite Data, Trajectory Modeling and Surface Insolation Measurements to Deduce the Direct Radiative Effect of Smoke Aerosol

In this paper, we have introduced a method of inferring the radiative effect of smoke aerosols using a technique that combines satellite remote sensing with trajectory modeling. The results shown here clearly show large flux biases between theoretical and measured radiative fluxes correlate with the arrival of smoke aerosol to the area. Further analysis is required to convincingly demonstrate that the reason for these differences is the radiative effect of the smoke aerosols. To do this, the estimated fluxes taken from the ERA-15 will be recomputed every 3 hours using International Satellite Cloud Climatology Project (ISCCP) data set entitled DX gridded to a 1o equal angle resolution (see paper 7B.2 for details). Surface radiometric and ancillary data for several more Canadian surface sites are being obtained at minute temporal resolution. The ultimate purpose of this research is to derive aerosol smoke maps for fire events such as this to be included in an aerosol climatology and be incorporated in the computation of the earth's surface radiation budget to better understand the radiative effect of aerosols.

Stackhouse, Paul W., Jr.

Global Multispectral Cloud Retrievals from MODIS

The Moderate Resolution Imaging Spectroradiometer (MODIS) was developed by NASA and launched onboard the Terra spacecraft on December 18,1999 and Aqua spacecraft on May 4,2002. It achieved its final orbit and began Earth observations on February 24, 2000 for Terra and June 24, 2002 for Aqua. A comprehensive set of remote sensing algorithms for cloud masking and the retrieval of cloud physical and optical properties has been developed by members of the MODIS atmosphere science team. The archived products from these algorithms have applications in climate change studies, climate modeling, numerical weather prediction, as well as fundamental atmospheric research. In addition to an extensive cloud mask, products include cloud-top properties (temperature, pressure, effective emissivity), cloud thermodynamic phase, cloud optical and microphysical parameters (optical thickness, effective particle radius, water path), as well as derived statistics. We will describe the various cloud properties being analyzed on a global basis from both Terra and Aqua, and will show characteristics of cloud optical and microphysical properties as a function of latitude for land and ocean separately, and contrast the statistical properties of similar cloud types in various parts of the world.

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.

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

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