# Search Results

## You are looking at 11 - 17 of 17 items for

- Author or Editor: Robert F. Cahalan x

- Refine by Access: All Content x

## Abstract

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 ranges 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.

## Abstract

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 ranges 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.

## Abstract

Conventional wisdom is that lidar pulses do not significantly penetrate clouds having an optical thickness exceeding about *τ* = 2, and that no returns are detectible from more than a shallow skin depth. Yet optically thicker clouds of *τ* ≫ 2 reflect a larger fraction of visible photons and account for much of the earth’s global average albedo. As cloud-layer thickness grows, an increasing fraction of reflected photons are scattered multiple times within the cloud and return from a diffuse concentric halo that grows around the incident pulse, increasing in horizontal area with layer physical thickness. The reflected halo is largely undetected by narrow field-of-view (FOV) receivers commonly used in lidar applications. Cloud Thickness from Offbeam Returns (THOR) is an airborne wide-angle detection system with multiple FOVs, capable of observing the diffuse halo as a wide-angle signal, from which the physical thickness of optically thick clouds can be retrieved. This paper describes the THOR system, demonstrates that the halo signal is stronger for thicker clouds, and presents a validation of physical thickness retrievals for clouds having *τ* > 20, from NASA’s P-3B flights over the Department of Energy’s Atmospheric Radiation Measurement Southern Great Plains site, using the lidar, radar, and other ancillary ground-based data.

## Abstract

Conventional wisdom is that lidar pulses do not significantly penetrate clouds having an optical thickness exceeding about *τ* = 2, and that no returns are detectible from more than a shallow skin depth. Yet optically thicker clouds of *τ* ≫ 2 reflect a larger fraction of visible photons and account for much of the earth’s global average albedo. As cloud-layer thickness grows, an increasing fraction of reflected photons are scattered multiple times within the cloud and return from a diffuse concentric halo that grows around the incident pulse, increasing in horizontal area with layer physical thickness. The reflected halo is largely undetected by narrow field-of-view (FOV) receivers commonly used in lidar applications. Cloud Thickness from Offbeam Returns (THOR) is an airborne wide-angle detection system with multiple FOVs, capable of observing the diffuse halo as a wide-angle signal, from which the physical thickness of optically thick clouds can be retrieved. This paper describes the THOR system, demonstrates that the halo signal is stronger for thicker clouds, and presents a validation of physical thickness retrievals for clouds having *τ* > 20, from NASA’s P-3B flights over the Department of Energy’s Atmospheric Radiation Measurement Southern Great Plains site, using the lidar, radar, and other ancillary ground-based data.

## Abstract

Zonally averaged meteorological fields can have large variances in polar regions due to purely geometrical effects, because fewer statistically independent areas contribute to zonal means near the poles than near the equator. A model of a stochastic field with homogeneous statistics on the sphere is presented as an idealized example of the phenomenon. We suggest a quantitative method for isolating the geometrical effect and use it in examining the variance of the zonally averaged 500 mb geopotential height field.

## Abstract

Zonally averaged meteorological fields can have large variances in polar regions due to purely geometrical effects, because fewer statistically independent areas contribute to zonal means near the poles than near the equator. A model of a stochastic field with homogeneous statistics on the sphere is presented as an idealized example of the phenomenon. We suggest a quantitative method for isolating the geometrical effect and use it in examining the variance of the zonally averaged 500 mb geopotential height field.

## Abstract

Empirical Orthogonal Functions (EOF's), eigenvectors of the spatial cross-covariance matrix of a meteorological field, are reviewed with special attention given to the necessary weighting factors for gridded data and the sampling errors incurred when too small a sample is available. The geographical shape of an EOF shows large intersample variability when its associated eigenvalue is “close” to a neighboring one. A rule of thumb indicating when an EOF is likely to be subject to large sampling fluctuations is presented. An explicit example, based on the statistics of the 500 mb geopotential height field, displays large intersample variability in the EOF's for sample sizes of a few hundred independent realizations, a size seldom exceeded by meteorological data sets.

## Abstract

Empirical Orthogonal Functions (EOF's), eigenvectors of the spatial cross-covariance matrix of a meteorological field, are reviewed with special attention given to the necessary weighting factors for gridded data and the sampling errors incurred when too small a sample is available. The geographical shape of an EOF shows large intersample variability when its associated eigenvalue is “close” to a neighboring one. A rule of thumb indicating when an EOF is likely to be subject to large sampling fluctuations is presented. An explicit example, based on the statistics of the 500 mb geopotential height field, displays large intersample variability in the EOF's for sample sizes of a few hundred independent realizations, a size seldom exceeded by meteorological data sets.

## Abstract

An increase in the planetary albedo of the earth-atmosphere system by only 10% can decrease the equilibrium surface temperature to that of the last ice age. Nevertheless, albedo biases of 10% or greater would be introduced into large regions of current climate models if clouds were given their observed liquid water amounts, because of the treatment of clouds as plane parallel. Past work has addressed the effect of cloud shape on albedo; here the focus is on the within-cloud variability of the vertically integrated liquid water. The main result is an estimate of the “plane-parallel albedo bias” using the “independent pixel approximation,” which ignores net horizontal photon transport, from a simple fractal model of marine stratocumulus clouds that ignores the cloud shape. The use of the independent pixel approximation in this context will be justified in a separate Monte Carlo study.

The focus on marine stratocumulus clouds is due to their important role in cloud radiative forcing and also that, of the wide variety of earth's cloud types, they are most nearly plane parallel, so that they have the least albedo bias. The fractal model employed here reproduces both the probability distribution and the wavenumber spectrum of the stratocumulus liquid water path, as observed during the First ISCCP Regional Experiment (FIRE). The model distributes the liquid water by a cascade process, related to the upscale cascade of energy transferred from the cloud thickness scale to the mesoscale by approximately 2D motions. For simplicity, the cloud microphysical parameters are assumed homogeneous, as is the geometrical cloud thickness; and the mesoscale-averaged vertical optical thickness is kept fixed at each step of the cascade. A single new fractal parameter, 0 ≤ *f* ≤ 1, is introduced and determined empirically by the variance of the logarithm of the vertically integrated liquid water. In the case of conservative scattering, the authors are able to estimate the albedo bias analytically as a function of the fractal parameter *f*, mean vertical optical thickness *T _{ν}
*, and sun angle

*θ*. Typical observed values are

*f*= 0.5,

*T*= 15, and

_{ν}*θ*= 60°, which give an absolute bias of 0.09, or a relative bias equal to 15% of the plane-parallel albedo of 0.60. The reduced reflectivity of fractal stratocumulus clouds is approximately given by the plane-parallel reflectivity evaluated at a reduced “effective optical thickness,” which when

*f*= 0.5 is

*T*

_{eff}≈ 10.

Study of the diurnal cycle of stratocumulus liquid water during FIRE leads to a key unexpected result: the plane-parallel albedo bias is largest when the cloud fraction reaches 100%, that is, when any bias associated with the cloud fraction vanishes. This is primarily due to the variability increase with cloud fraction. Thus, the within-cloud fractal structure of stratocumulus has a more significant impact on estimates of its mesoscale-average albedo than does the cloud fraction.

## Abstract

An increase in the planetary albedo of the earth-atmosphere system by only 10% can decrease the equilibrium surface temperature to that of the last ice age. Nevertheless, albedo biases of 10% or greater would be introduced into large regions of current climate models if clouds were given their observed liquid water amounts, because of the treatment of clouds as plane parallel. Past work has addressed the effect of cloud shape on albedo; here the focus is on the within-cloud variability of the vertically integrated liquid water. The main result is an estimate of the “plane-parallel albedo bias” using the “independent pixel approximation,” which ignores net horizontal photon transport, from a simple fractal model of marine stratocumulus clouds that ignores the cloud shape. The use of the independent pixel approximation in this context will be justified in a separate Monte Carlo study.

The focus on marine stratocumulus clouds is due to their important role in cloud radiative forcing and also that, of the wide variety of earth's cloud types, they are most nearly plane parallel, so that they have the least albedo bias. The fractal model employed here reproduces both the probability distribution and the wavenumber spectrum of the stratocumulus liquid water path, as observed during the First ISCCP Regional Experiment (FIRE). The model distributes the liquid water by a cascade process, related to the upscale cascade of energy transferred from the cloud thickness scale to the mesoscale by approximately 2D motions. For simplicity, the cloud microphysical parameters are assumed homogeneous, as is the geometrical cloud thickness; and the mesoscale-averaged vertical optical thickness is kept fixed at each step of the cascade. A single new fractal parameter, 0 ≤ *f* ≤ 1, is introduced and determined empirically by the variance of the logarithm of the vertically integrated liquid water. In the case of conservative scattering, the authors are able to estimate the albedo bias analytically as a function of the fractal parameter *f*, mean vertical optical thickness *T _{ν}
*, and sun angle

*θ*. Typical observed values are

*f*= 0.5,

*T*= 15, and

_{ν}*θ*= 60°, which give an absolute bias of 0.09, or a relative bias equal to 15% of the plane-parallel albedo of 0.60. The reduced reflectivity of fractal stratocumulus clouds is approximately given by the plane-parallel reflectivity evaluated at a reduced “effective optical thickness,” which when

*f*= 0.5 is

*T*

_{eff}≈ 10.

Study of the diurnal cycle of stratocumulus liquid water during FIRE leads to a key unexpected result: the plane-parallel albedo bias is largest when the cloud fraction reaches 100%, that is, when any bias associated with the cloud fraction vanishes. This is primarily due to the variability increase with cloud fraction. Thus, the within-cloud fractal structure of stratocumulus has a more significant impact on estimates of its mesoscale-average albedo than does the cloud fraction.

## Abstract

Liquid clouds play a profound role in the global radiation budget, but it is difficult to retrieve their vertical profile remotely. Ordinary narrow-field-of-view (FOV) lidars receive a strong return from such clouds, but the information is limited to the first few optical depths. Wide-angle multiple-FOV lidars can isolate radiation that is scattered multiple times before returning to the instrument, often penetrating much deeper into the cloud than does the single-scattered signal. These returns potentially contain information on the vertical profile of the extinction coefficient but are challenging to interpret because of the lack of a fast radiative transfer model for simulating them. This paper describes a variational algorithm that incorporates a fast forward model that is based on the time-dependent two-stream approximation, and its adjoint. Application of the algorithm to simulated data from a hypothetical airborne three-FOV lidar with a maximum footprint width of 600 m suggests that this approach should be able to retrieve the extinction structure down to an optical depth of around 6 and a total optical depth up to at least 35, depending on the maximum lidar FOV. The convergence behavior of Gauss–Newton and quasi-Newton optimization schemes are compared. Results are then presented from an application of the algorithm to observations of stratocumulus by the eight-FOV airborne Cloud Thickness from Off-Beam Lidar Returns (THOR) lidar. It is demonstrated how the averaging kernel can be used to diagnose the effective vertical resolution of the retrieved profile and, therefore, the depth to which information on the vertical structure can be recovered. This work enables more rigorous exploitation of returns from spaceborne lidar and radar that are subject to multiple scattering than was previously possible.

## Abstract

Liquid clouds play a profound role in the global radiation budget, but it is difficult to retrieve their vertical profile remotely. Ordinary narrow-field-of-view (FOV) lidars receive a strong return from such clouds, but the information is limited to the first few optical depths. Wide-angle multiple-FOV lidars can isolate radiation that is scattered multiple times before returning to the instrument, often penetrating much deeper into the cloud than does the single-scattered signal. These returns potentially contain information on the vertical profile of the extinction coefficient but are challenging to interpret because of the lack of a fast radiative transfer model for simulating them. This paper describes a variational algorithm that incorporates a fast forward model that is based on the time-dependent two-stream approximation, and its adjoint. Application of the algorithm to simulated data from a hypothetical airborne three-FOV lidar with a maximum footprint width of 600 m suggests that this approach should be able to retrieve the extinction structure down to an optical depth of around 6 and a total optical depth up to at least 35, depending on the maximum lidar FOV. The convergence behavior of Gauss–Newton and quasi-Newton optimization schemes are compared. Results are then presented from an application of the algorithm to observations of stratocumulus by the eight-FOV airborne Cloud Thickness from Off-Beam Lidar Returns (THOR) lidar. It is demonstrated how the averaging kernel can be used to diagnose the effective vertical resolution of the retrieved profile and, therefore, the depth to which information on the vertical structure can be recovered. This work enables more rigorous exploitation of returns from spaceborne lidar and radar that are subject to multiple scattering than was previously possible.

The interaction of clouds with solar and terrestrial radiation is one of the most important topics of climate research. In recent years it has been recognized that only a full three-dimensional (3D) treatment of this interaction can provide answers to many climate and remote sensing problems, leading to the worldwide development of numerous 3D radiative transfer (RT) codes. The international Intercomparison of 3D Radiation Codes (I3RC), described in this paper, sprung from the natural need to compare the performance of these 3D RT codes used in a variety of current scientific work in the atmospheric sciences. I3RC supports intercomparison and development of both exact and approximate 3D methods in its effort to 1) understand and document the errors/limits of 3D algorithms and their sources; 2) provide “baseline” cases for future code development for 3D radiation; 3) promote sharing and production of 3D radiative tools; 4) derive guidelines for 3D radiative tool selection; and 5) improve atmospheric science education in 3D RT. Results from the two completed phases of I3RC have been presented in two workshops and are expected to guide improvements in both remote sensing and radiative energy budget calculations in cloudy atmospheres.

The interaction of clouds with solar and terrestrial radiation is one of the most important topics of climate research. In recent years it has been recognized that only a full three-dimensional (3D) treatment of this interaction can provide answers to many climate and remote sensing problems, leading to the worldwide development of numerous 3D radiative transfer (RT) codes. The international Intercomparison of 3D Radiation Codes (I3RC), described in this paper, sprung from the natural need to compare the performance of these 3D RT codes used in a variety of current scientific work in the atmospheric sciences. I3RC supports intercomparison and development of both exact and approximate 3D methods in its effort to 1) understand and document the errors/limits of 3D algorithms and their sources; 2) provide “baseline” cases for future code development for 3D radiation; 3) promote sharing and production of 3D radiative tools; 4) derive guidelines for 3D radiative tool selection; and 5) improve atmospheric science education in 3D RT. Results from the two completed phases of I3RC have been presented in two workshops and are expected to guide improvements in both remote sensing and radiative energy budget calculations in cloudy atmospheres.