Ice Cloud Optical Thickness, Effective Radius, And Ice Water Path Inferred From Fused MISR and MODIS Measurements Based on a Pixel‐Level Optimal Ice Particle Roughness Model
Explore the source record for details and available documents.
SEARCH · Search NASA
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Explore the source record for details and available documents.
The aircraft microphysics probe, PVM-100A, was tested in the Colorado State University dynamic cloud chamber to establish its ability to measure ice water content (IWC), PSA, and Re in ice clouds. Its response was compared to other means of measuring those ice-cloud parameters that included using FSSP-100 and 230-X 1-D optical probes for ice-crystal concentrations, a film-loop microscope for ice-crystal habits and dimensions, and an in-situ microscope for determining ice-crystal orientation. Intercomparisons were made in ice clouds containing ice crystals ranging in size from about 10 microns to 150 microns diameter, and ice crystals with plate, columnar, dendritic, and spherical shapes. It was not possible to determine conclusively that the PVM accurately measures IWC, PSA, and Re of ice crystals, because heat from the PVM evaporated in part the crystals in its vicinity in the chamber thus affecting its measurements. Similarities in the operating principle of the FSSP and PVM, and a comparison between Re measured by both instruments, suggest, however, that the PVM can make those measurements. The resolution limit of the PVM for IWC measurements was found to be on the order of 0.001 g/cubic m. Algorithms for correcting IWC measured by FSSP and PVM were developed.
The effective radius(r(sub e)) is a crucial variable in representing the radiative properties of cloud layers in general circulation models. This parameter is proportional to the condensed water content (CWC) divided by the extinction (sigma). For ice cloud layers, parameterizations for r(sub e), have been developed from aircraft in-situ measurements 1) indirectly, using data obtained from particle spectrometer probes and assumptions or observations about particle shape and mass to get the ice water content (IWC) and area to get sigma, and recently 2) from probes that measure IWC and sigma directly. This study compares [IWC/sigma] derived from the two methods using data sets acquired from comparable instruments on two aircraft, one sampling clouds at mid-levels and the other at upper-levels during the CRYSTAL-FACE field program in Florida in 2002. The sigma and IWC derived by each method are compared and evaluated in different ways for each aircraft data set. Direct measurements of sigma exceed those derived indirectly by a factor of two to two and a half. The IWC probes, relying on ice sublimation, appear to measure accurately except when the IWC is high or the particles too large to sublimate completely during the short transit time through the probe. The IWC estimated from the particle probes are accurate when direct measurements are available to provide constraints and useful information in high IWC/large particle situations. Because of the discrepancy in sigma estimates between the direct and indirect approaches, there is a factor of 2 to 3 difference in [IWC/sigma] between them. Although there are significant uncertainties involved in its use, comparisons with several independent data sources suggest that the indirect method is the more accurate of the two approaches. However, experiments are needed to resolve the source of the discrepancy in sigma.
Satellite observations show that ice cloud effective radius (r(sub e)) increases with ice water content (IWC) but decreases with aerosol optical thickness (AOT). Using least-squares fitting to the observed data, we obtain an analytical formula to describe the variations of r(sub e) with IWC and AOT for several regions with distinct characteristics of r(sub e) -IWC-AOT relationships. As IWC directly relates to convective strength and AOT represents aerosol loading, our empirical formula provides a means to quantify the relative roles of dynamics and aerosols in controlling r(sub e) in different geographical regions, and to establish a framework for parameterization of aerosol effects on r(sub e) in climate models.
A novel technique is used to estimate derivatives of ice effective radius with respect to height near convective cloud tops (dr(sub e)/dz) from airborne shortwave reflectance measurements and lidar. Values of dr(sub e)/dz are about −6 micrometer/km for cloud tops below the homogeneous freezing level, increasing to near 0 micrometer/km above the estimated level of neutral buoyancy. Retrieved dr(sub e)/dz compares well with previously documented remote sensing and in situ estimates. Effective radii decrease with increasing cloud top height, while cloud top extinction increases. This is consistent with weaker size sorting in high, dense cloud tops above the level of neutral buoyancy where fewer large particles are present and with stronger size sorting in lower cloud tops that are less dense. The results also confirm that cloud top trends of effective radius can generally be used as surrogates for trends with height within convective cloud tops. These results provide valuable observational targets for model evaluation.
Ice crystals in clouds are highly complex. Their sizes, macroscale shape (i.e., habit), mesoscale shape (i.e., aspect ratio of components) and microscale shape (i.e., surface roughness) determine optical properties and affect physical properties such as fall speeds, growth rates and aggregation efficiency. Our current understanding on the formation and evolution of ice crystals under various conditions can be considered poor. Commonly, ice crystal size and shape are related to ambient temperature and humidity, but global observational statistics on the variation of ice crystal size and particularly shape have not been available. Here we show results of a project aiming to infer ice crystal size, shape and scattering properties from a combination of MODIS measurements and POLDER-PARASOL multi-angle polarimetry. The shape retrieval procedure infers the mean aspect ratios of components of ice crystals and the mean microscale surface roughness levels, which are quantifiable parameters that mostly affect the scattering properties, in contrast to a habit. We present global statistics on the variation of ice effective radius, component aspect ratio, microscale surface roughness and scattering asymmetry parameter as a function of cloud top temperature, latitude, location, cloud type, season, etc. Generally, with increasing height, sizes decrease, roughness increases, asymmetry parameters decrease and aspect ratios increase towards unity. Some systematic differences are observed for clouds warmer and colder than the homogeneous freezing level. Uncertainties in the retrievals will be discussed. These statistics can be used as observational targets for modeling efforts and to better constrain other satellite remote sensing applications and their uncertainties.
The solar reflectance bi-spectral (SRBS) and infrared split-window (IRSpW) methods are two of the most popular techniques for passive ice cloud property retrievals from multispectral imagers. Ice clouds are usually assumed to be vertically homogeneous in global operational algorithms based on these methods, although significant vertical variations of ice particle size are typically observed in ice clouds. In this Study we investigate uncertainties in retrieved optical thickness, effective particle size, and ice water path introduced by a homogeneous cloud assumption in both the SRBS and IRSpW methods, and focus on whether the assumption can lead to significant discrepancies between the two methods. The study simulates the upwelling spectral radiance associated with vertically structured clouds and passes the results through representative SRBS and IRSpW retrieval algorithms. Cloud optical thickness is limited to values for which IRSpW retrievals are possible (optical thickness less than about 7). When the ice cloud is optically thin and yet has a significant ice particle size vertical variation, it is found that both methods tend to underestimate the effective radius and ice water path. The reason for the underestimation is the nonlinear dependence of ice particle scattering properties (extinction and single scattering albedo) on the effective radius. Because the nonlinearity effect is stronger in the IRSpW than the SRBS method, the IRSpW-based IWP tends to be smaller than the SRBS counterpart. When the ice cloud is moderately optically thick, the IRSpW method is relatively insensitive to cloud vertical structure and effective radius retrieval is weighted toward smaller ice particle size, while the weighting function makes the SRBS method more sensitive to the ice particle size in the upper portion of the cloud. As a result, when ice particle size increases monotonically toward cloud base, the two methods are in qualitative agreement; in the event that ice particle size decreases toward cloud base, the effective radius and ice water path retrievals based on the SRBS method are substantially larger than those from the IRSpW. The main findings of this Study Suggest that the homogenous cloud assumption can affect the SRBS and IRSpW methods to different extents and, consequently, can lead to significantly different retrievals. Therefore caution should be taken when comparing and combining the ice cloud property retrievals from these two methods.
MICROBASE is a baseline retrieval of cloud microphysical properties. It uses a combination of observations from the cloud radar, ceilometer, micropulse lidar, microwave radiometer, and balloon-borne radiosonde soundings in order to produce instantaneous vertical profiles of cloud liquid water content (LWC), cloud ice water content (IWC), liquid cloud particle effective radius (LIQRE), and ice cloud particle effective radius (ICERE). Uncertainites are also produced by the VAP. The inputs are by the following VAPs: Active Remote Sensing of Clouds (ARSCL) Merged Sounding (MERGESONDE) Microwave Radiometer Retrievals (MWRRET) The output are daily files.
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.
Cloud optical thickness and effective radius retrievals from solar reflectance measurements traditionally depend on a combination of spectral channels that are absorbing and non-absorbing for liquid water droplets. Reflectances in non-absorbing channels (e.g., 0.67, 0.86 micrometer bands) are largely dependent on cloud optical thickness, while longer wavelength absorbing channels (1.6, 2.1, and 3.7 micrometer window bands) provide cloud particle size information. Retrievals are complicated by the presence of an underlying ice/snow surface. At the shorter wavelengths, sea ice is both bright and highly variable, significantly increasing cloud retrieval uncertainty. However, reflectances at the longer wavelengths are relatively small and may be comparable to that of dark open water. Sea ice spectral albedos derived from Cloud Absorption Radiometer (CAR) measurements during April 1992 and June 1995 Arctic field deployments are used to illustrate these statements. A modification to the traditional retrieval technique is devised. The new algorithm uses a combination of absorbing spectral channels for which the snow/ice albedo is relatively small. Using this approach, preliminary retrievals have been made with the MODIS Airborne Simulator (MAS) imager flown aboard the NASA ER-2 during FIRE-ACE. Data from coordinated ER-2 and University of Washington CV-580 aircraft observations of liquid water stratus clouds on June 3 and June 6, 1998 have been examined. Size retrievals are compared with in situ cloud profile measurements of effective radius made with the CV-580 PMS FSSP probe, and optical thickness retrievals are compared with extinction profiles derived from the Gerber Scientific "g-meter" probe. MAS retrievals are shown to be in good agreement with the in situ measurements.
We have performed one dimensional radiative transfer calculations to evaluate the impact of cirrus clouds on the tropical radiation budget. We investigate the sensitivity of solar and infrared fluxes to cloud optical depth, particle size distributions, and cloud height. If the observed solar cloud forcing in excess of 100 W/sq m is to be attributed to cirrus anvils alone, then the optical depth of these anvils must be at least 5 (assuming 50% cloud cover and an ice crystal effective radius of 15 microns). The net radiative forcing of cirrus near the tropical tropopause is positive (heating) for cloud optical depths less than about 16 and negative (cooling) for larger optical depths. If cirrus clouds alone are responsible for the equal and opposite shortwave and longwave cloud forcing in excess of 100 W/sq m observed by Earth Radiation Budget Experiment (ERBE), then the cirrus must typically take the form of deep, optically thick clouds with relatively small particles (radii of 10-20 microns) and cloud-tops well below the tropopause. The maintenance of this balance on monthly time scales can be attributed to a variety of correlations: The cloud cover of optically thick cirrus or thin cirrus overlying low-level stratus clouds could vary; or cirrus anvil height cloud increase along with a decrease in the ice crystal effective radius and an increase in optical depth. It would be of great interest to determine observationally which of these correlations is responsible for the observed lack of variation in cloud forcing.
A data set containing radiation, meteorological, and cloud sensor observations is documented. It was prepared for use by the Department of Energy's Atmospheric Radiation Measurement (ARM) program and other interested scientists. These data are a precursor of the types of data that ARM Cloud And Radiation Testbed (CART) sites will provide. The data are from the Cloud Lidar And Radar Exploratory Test (CLARET) conducted by the Wave Propagation Laboratory during autumn 1989 in the Denver-Boulder area of Colorado primarily for the purpose of developing new cloud-sensing techniques on cirrus. After becoming aware of this experiment, ARM scientists requested archival of subsets or the data to assist in the developing ARM program. Five CLARET cases were selected: two with cirrus, one with stratus, one with mixed-phase clouds, and one with clear skies. The cases range from 2 to 9.5 h in length. A pyranometer, pyrgeometer, pyrheliometer, and an infrared radiometer constituted the ensemble of instruments that provided surface radiation data. A lidar, radar, and ceilometer observed the cloud geometrical structure, and visual reports and all-sky camera observations were assimilated to provide cloud cover data. Radiosondes, wind profiler, RASS (profiling virtual temperature), microwave radiometers (observing column integrated liquid water and water vapor), and standard surface measurements provided meteorological data. Satellite data from the stratus case and one cirrus case were analyzed for statistics on cloud cover and top height. The main body of the selected data are available on diskette from the Wave Propagation Laboratory or Los Alamos National Laboratory. In addition to documenting the data set, this report describes CLARET and gives a bibliography of publications associated with the project. Some preliminary results of CLARET' research are also summarized. Simultaneous CO 2 lidar and radar backscatter measurements were shown to provide estimates of the effective radius of ice particles. Simultaneous radar and infrared radiometer data appear useful for estimating column-integrated numbers and average sizes of ice cloud particles. Ice water content obtained with this method compared favorably with values from another empirical technique using radar data alone. Depolarization of the CO 2 lidar signal from ice clouds was surprisingly small, suggesting that calculation of backscatter from nonspherical particles for this lidar is a tractable problem. Examples are also cited of CO 2 lidar measurements of the effective radius of water cloud drop size distributions and of inference of the size of pristine ice crystals that assume a particular orientation in the air. These parameters are all important to radiative transfer through clouds.
In this study, cloud properties measured from the CALIPSO, CloudSat, and MODIS (CCM) are used for top-of-atmosphere (TOA) shortwave (SW) broadband (BB) irradiance computations. The CALIPSO and CloudSat active sensors provide detailed cloud vertical profiles, but these occasionally miss parts of the cloud columns due to the full attenuation of sensor signals, surface clutter, or insensitivity to a certain range of cloud particle sizes. As a result, the CCM-merged cloud extinction coefficient profiles can be underestimated. Therefore, we compare the column-integrated visible scaled cloud optical depth (VSCOD) of the CCM-merged cloud extinction coefficient profile with the MODIS-estimated VSCOD and apply a scaling factor to the CCM-merged cloud extinction profile. The VSCOD is defined as a visible cloud optical depth multiplied by (1¬–asymmetry parameter). The SW irradiances are computed using the scaled CCM-merged cloud extinction coefficient and effective radius profiles. It is shown that the multi-sensor-combined cloud profiles significantly reduce positive TOA SW BB biases, compared to those with MODIS-derived cloud properties only. The improvement is more pronounced for optically thick clouds, where MODIS ice particle effective radius is largely underestimated.
These investigations focused essentially on the macroscale organization of cloud and precipitation which occurred during the 4th Intensive Observing Period (IOP-4) of the Experiment for Rapidly Intensifying Cyclones over the Atlantic (ERICA). This experiment, held off the East Coast of the United States and Canada during the winter of 1989, documented several episodes of rapid cyclonic storm development. Also playing a major role as validation and ground truth in these studies are Special Sensor Microwave Imager (SSM/I) retrievals of precipitable water, total liquid water and ice, generated by other Marshall Space Flight Center (MSFC) supported investigations. Model simulations produced to date suggest that, while the large-scale atmospheric dynamics was an essential driving mechanism, the role of condensation was crucial in facilitating the exceptionally rapid spinup of the cyclone and the low surface pressure. A model simulation of the precipitation rate at the time of most rapid storm intensification is shown in the accompanying figure. Heavier precipitation rates in the crescent shaped region are associated with deep convection along the leading edge of a dry intrusion behind the surface low. The majority of precipitation in the stratiform region to the northeast involved the production of ice with deposition from vapor to ice being the dominant process of growth. Some small amount of mixed phase cloudiness was present. Model condensate distributions matched well with SSM/I observations. The accompanying SSM/I imagery which delineates areas of large (greater than several hundred micron effective radius) precipitating ice over the ocean suggests that the model has done well in capturing the essential mechanisms responsible for the horizontal distribution of precipitation.
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.
Since multilayer cloud scenes are common in the atmosphere and can be an important source of uncertainty in passive satellite sensor cloud retrievals, the MODIS MOD06 and MYD06 standard cloud optical property products include a multilayer cloud detection algorithm to assist with data quality assessment. This paper presents an evaluation of the Aqua MODIS MYD06 Collection 6 multilayer cloud detection algorithm through comparisons with active Cloud Profiling Radar (CPR) and Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) products that have the ability to provide cloud vertical distributions and directly classify multilayer cloud scenes and layer properties. To compare active sensor products with an imager such as MODIS, it is first necessary to define multilayer clouds in the context of their radiative impact on cloud retrievals. Three main parameters have thus been considered in this evaluation: (1) the maximum separation distance between two cloud layers, (2) the thermodynamic phase of those layers and (3) the upper-layer cloud optical thickness. The impact of including the Pavolonis–Heidinger multilayer cloud detection algorithm, introduced in Collection 6, to assist with multilayer cloud detection has also been assessed. For the year 2008, the MYD06 C6 multilayer cloud detection algorithm identifies roughly 20% of all cloudy pixels as multilayer (decreasing to about 13% if the Pavolonis–Heidinger algorithm output is not used). Evaluation against the merged CPR and CALIOP2B-CLDCLASS-lidar product shows that the MODIS multilayer detection results are quite sensitive to how multilayer clouds are defined in the radar and lidar product and that the algorithm performs better when the optical thickness of the upper cloud layer is greater than about 1.2 with a minimum layer separation distance of 1 km. Finally, we find that filtering the MYD06 cloud optical properties retrievals using the multilayer cloud flag improves aggregated statistics, particularly for ice cloud effective radius.
The one-dimensional cirrus model described in part 1 of this issue has been used to study the sensitivity of simulated cirrus microphysical and radiative properties to poorly known model parameters, poorly understood physical processes, and environmental conditions. Model parameters and physical processes investigated include nucleation rate, mode of nucleation (e.g., homogeneous freezing of aerosols and liquid droplets or heterogeneous deposition), ice crystal shape, and coagulation. These studies suggest that the leading sources of uncertainty in the model are the phase change (liquid-solid) energy barrier and the ice-water surface energy which dominate the homogeneous freezing nucleation rate and the coagulation sticking efficiency at low temperatures which controls the production of large ice crystals (radii greater than 100 mcirons). Environmental conditions considered in sensitivity tests were CN size distribution, vertical wind speed, and cloud height. We found that (unlike stratus clouds) variations in the total number of condensation nuclei (NC) have little effect on cirrus microphysical and radiative properties, since nucleation occurs only on the largest CN at the tail of the size distribution. The total number of ice crystals which nucleate has little or no relationship to the number of CN present and depends primarily on the temperature and the cooling rate. Stronger updrafts (more rapid cooling) generate higher ice number densities, ice water content, cloud optical depth, and net radiative forcing. Increasing the height of the clouds in the model leads to an increase in ice number density, a decrease in effective radius, and a decrease in ice water content. The most prominent effect of increasing cloud height was a rapid increase in the net cloud radiative forcing which can be attributed to the change in cloud temperature as well as change in cloud ice size distributions. It has long been recognized that changes in cloud height or cloud area have the greatest potential for causing feedbacks on climate change. Our results suggest that variations in vertical velocity or cloud microphysical changes associatd with cloud height changes may also be important.
Single-layer mixed-phase stratiform (MPS) Arctic clouds, which formed under conditions of large surface heat flux combined with general subsidence during a subperiod of the Atmospheric Radiation Measurement (ARM) Program Mixed-Phase Arctic Cloud Experiment (M-PACE), are simulated with a cloud resolving model (CRM). The CRM is implemented with either an advanced two-moment (M05) or a commonly used one-moment (L83) bulk microphysics scheme and a state-of-the-art radiative transfer scheme. The CONTROL simulation, that uses the M05 scheme and observed aerosol size distribution and ice nulei (IN) number concentration, reproduces the magnitudes and vertical structures of cloud liquid water content (LWC), total ice water content (IWC), number concentration and effective radius of cloud droplets as suggested by the M-PACE observations. It underestimates ice crystal number concentrations by an order of magnitude and overestimates effective radius of ice crystals by a factor of 2-3. The OneM experiment, that uses the L83 scheme, produces values of liquid water path (LWP) and ice plus snow water path (ISWP) that were about 30% and 4 times, respectively, of those produced by the CONTROL. Its vertical profile of IWC exhibits a bimodal distribution in contrast to the constant distribution of IWC produced in the CONTROL and observations.