ALBERT

Description: We review the basic multiple scattering theory of off-beam lidar returns from optically thick clouds using the diffusion approximation. The shape of the temporal signal - the stretched pulse - depends primarily on the physical thickness of the cloud whereas its spatial counterpart - the diffuse spot - conveys specific information on the cloud's optical thickness, as do the absolute returns. This makes observation of the weak off-beam lidar returns an attractive prospect in remote sensing of cloud properties. By estimating the signal-to-noise ratio, we show that night-time measurements can be performed with existing technology. By the same criterion, day-time operation is a challenge that can only be met with a combination of cutting-edge techniques in filtering and in laser sources.

Description: An analysis of nadir reflectivity Fourier spatial power spectra and autocorrelation functions at solar wavelengths and for cloudy conditions has been carried out. The data come from Landsat Thematic Mapper (TM) observations, while Monte Carlo (MC) simulations are used to aid the interpretation of the Landsat results. We show that radiative processes produce consistent signatures on power spectra and autocorrelation functions. The former take a variety of forms not shown or explained in previous observational studies. We demonstrate that the TM spectra can potentially be affected by both radiative "roughening" at intermediate scales (approx. 1 -5 km), being more prevalent at large solar zenith angles, and the already documented radiative "smoothing" at small scales (less than 1 km). These processes are wavelength dependent, as shown by systematic differences between conservative (for cloud droplets) TM band 4 (approx. 0.8 microns) and absorbing band 7 (approx. 2.2 microns): band 7 exhibits more roughening and less smoothing. This is confirmed quantitatively by comparing least-squared fitted power spectral slopes for the two bands. It is also corroborated by a slower decrease with distance of autocorrelation function values for band 4 compared to band 7. The appearance of roughening at large solar zenith angles is a result of side illumination and shadowing and adds an additional complexity to the power spectra. MC spectra are useful in illustrating that scale invariant optical depth fields can produce complex power spectra that take a variety of shapes under different conditions. We show that radiative roughening increases with the decrease of single scattering albedo and with the increase of solar zenith angle (as in the observations). For high Sun there is also a clear shift of the radiative smoothing scale to smaller values as droplet absorption increases. The shape of the power spectrum is sensitive to the magnitude and type of cloud top height variability, with the spectral signatures of decorrelation between reflectance and optical depth at large scales becoming stronger as the magnitude of cloud top variations increase. Finally, the usefulness of power spectral analysis in evaluating the skill of novel optical depth retrieval techniques in removing 3D radiative effects is demonstrated. New techniques using inverse Non-local Independent Pixel Approximation (NIPA) and Normalized Difference of Nadir Reflectivity (NDNR) yield optical depth fields which better match the scale-by-scale variability of the true optical depth field.

Description: Cloud radiative properties are sensitive to drop size and other parameters of cloud micro-structure, but also to cloud shape,spacing, and other parameters of cloud macro-structure, including internal fractal structure. New information on cloud structure is being derived from a variety of cloud radars. Ongoing field programs such as Department of Energy's Atmospheric Radiation Measurement (DoE/ARM) are improving the measurement and modelling of physical and radiative properties of clouds. A parallel effort is underway to improve cloud remote sensing, especially from the new suite of EOS-AM1 instruments which will provide higher spectral, spatial resolution, and/or angular resolution. Key parameters for improving pixel-scale retrievals are cloud thickness and photon mean-free-path, which together determine the scale of "radiative smoothing" of cloud fluxes and radiances. This scale has been observed as a change in the spatial spectrum of Landsat cloud radiances, and was also recently found with the Goddard micropulse lidar, by searching for returns from directions nonparallel to the incident beam. "Offbeam" Lidar returns are now being used to estimate the cloud "radiative Green's function", G,which depends on cloud thickness and may be used to retrieve that important quantity. G is also being applied to improving simple IPA estimates of cloud radiative properties. This and other measurements of 3D transfer in clouds, coupled with Monte Carlo and other 3D transfer methods, are beginning to provide a better understanding of the dependence of adiation on cloud inhomogeneity, and to suggest new retrieval and parameterization algorithms which take account of cloud inhomogeneity. An international "Intercomparison of 3D Radiation Codes" or I3RC, program is beginning to coordinate and evaluate the variety of 3D radiative transfer methods now available, and to make them more widely available. Information is on the Web at: http://climate.qsfc.nasa.crov/I3RC. Input consists of selected cloud fields derived from data sources such as radar, microwave and satellite, and from models involved in the GEWEX Cloud Systems Studies. Output is selected radiative quantities that characterize the large- scale properties of the fields of radiative fluxes and heating. Several example cloud fields will be used to illustrate.

Description: Measurements of the distribution of reflected light from a laser beam incident on an aqueous suspension of particles or "cloud" with known thickness and particle size distribution are reported. The distribution is referred to as the "cloud radiative Green's function", G. In the diffusion domain, G is sensitive to cloud thickness, allowing that important quantity to be retrieved. The goal of the laboratory simulation is to provide preliminary estimates of sensitivity of G to cloud thickness,for use in the optimal design of an offbeam Lidar instrument for remote sensing of cloud thickness (THOR, Thickness from Offbeam Returns). These clouds of polystyrene microspheres suspended in water are analogous to real clouds of water droplets suspended in air. The microsphere size distribution is roughly lognormal, from 0.5 microns to 25 microns, similar to real clouds. Density of suspended spheres is adjusted so mean-free-path of visible photons is about 10 cm, approximately 1000 times smaller than in real clouds. The light source is a ND:YAG laser at 530 nm. Detectors are flux and photon-counting Photomultiplier Tube (PMTS), with a glass probe for precise positioning. A Labview 5 VI controls positioning, and data acquisition, via an NI Motion Control board connected to a stepper motor driving an Edmund linear slider, and a 16-channel 16-bit NI-DAQ board. The stepper motor is accurate to 10 microns, and step size is selectable from the VI software. Far from the incident beam, the rate of exponential increase as the direction of the incident beam is approached scales as expected from diffusion theory, linearly with the cloud thickness, and inversely as the square root of the reduced optical thickness, and is independent of particle size. Near the beam the signal begins to increase faster than exponential, due to single and low-order scattering near the backward direction, and here the distribution depends on particle size. Results are being used to verify 3D Monte Carlo radiative transfer simulations, used to estimate signal-to-noise ratios for remotely sensed off beam returns, for both homogeneous and inhomogeneous clouds. Signal-to-noise estimates show that unfiltered observations are straight forward at night, while narrow band pass filters are being studied for day.

Description: Cloud radiative properties are sensitive to drop size and other parameters of cloud micro-structure, but also to cloud shape,spacing, and other parameters of cloud macro-structure, including internal fractal structure. New information on cloud structure is being derived from a variety of cloud radars, and ongoing field programs such as DoE/ARM. These programs are improving the measurement and modelling of physical and radiative properties of clouds. A parallel effort is underway to improve cloud remote sensing, especially from the new suite of EOS-AM1 instruments which will provide higher spectral, spatial resolution, and/or angular resolution. Key parameters for improving pixel-scale retrievals are cloud thickness and photon mean-free-path, which together determine the scale of "radiative smoothing" of cloud fluxes and radiances. This scale has been observed as a change in the spatial spectrum of Landsat cloud radiances, and was also recently found with the Goddard micropulse lidar, by searching for returns from directions nonparallel to the incident beam. "Offbeam" Lidar returns are now being used to estimate the cloud "radiative Green's function", G,which depends on cloud thickness and may be used to retrieve that important quantity. G is also being applied to improving simple IPA estimates of cloud radiative properties. This and other measurements of 3D transfer in clouds, coupled with Monte Carlo and other 3D transfer methods, are beginning to provide a better understanding of the dependence of radiation on cloud inhomogeneity, and to suggest new retrieval and parameterization algorithms which take account of cloud inhomogeneity.

Description: To improve radiative transfer calculations for inhomogeneous clouds, a consistent means of modeling inhomogeneity is needed. One current method of modeling cloud inhomogeneity is through the use of fractal parameters. This method is based on the supposition that cloud inhomogeneity over a large range of scales is related. An analysis technique named wavelet analysis provides a means of studying the multiscale nature of cloud inhomogeneity. In this paper, the authors discuss the analysis and modeling of cloud inhomogeneity through the use of wavelet analysis. Wavelet analysis as well as other windowed analysis techniques are used to study liquid water path (LWP) measurements obtained during the marine stratocumulus phase of the First ISCCP (International Satellite Cloud Climatology Project) Regional Experiment. Statistics obtained using analysis windows, which are translated to span the LWP dataset, are used to study the local (small scale) properties of the cloud field as well as their time dependence. The LWP data are transformed onto an orthogonal wavelet basis that represents the data as a number of times series. Each of these time series lies within a frequency band and has a mean frequency that is half the frequency of the previous band. Wavelet analysis combined with translated analysis windows reveals that the local standard deviation of each frequency band is correlated with the local standard deviation of the other frequency bands. The ratio between the standard deviation of adjacent frequency bands is 0.9 and remains constant with respect to time. This ratio defined as the variance coupling parameter is applicable to all of the frequency bands studied and appears to be related to the slope of the data's power spectrum. Similar analyses are performed on two cloud inhomogeneity models, which use fractal-based concepts to introduce inhomogeneity into a uniform cloud field. The bounded cascade model does this by iteratively redistributing LWP at each scale using the value of the local mean. This model is reformulated into a wavelet multiresolution framework, thereby presenting a number of variants of the bounded cascade model. One variant introduced in this paper is the 'variance coupled model,' which redistributes LWP using the local standard deviation and the variance coupling parameter. While the bounded cascade model provides an elegant two- parameter model for generating cloud inhomogeneity, the multiresolution framework provides more flexibility at the expense of model complexity. Comparisons are made with the results from the LWP data analysis to demonstrate both the strengths and weaknesses of these models.