Chris A. Hostetler

The CALIPSO Automated Aerosol Classification and Lidar Ratio Selection Algorithm

Abstract Descriptions are provided of the aerosol classification algorithms and the extinction-to-backscatter ratio (lidar ratio) selection schemes for the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) aerosol products. One year of CALIPSO level 2 version 2 data are analyzed to assess the veracity of the CALIPSO aerosol-type identification algorithm and generate vertically resolved distributions of aerosol types and their respective optical characteristics. To assess the robustness of the algorithm, the interannual variability is analyzed by using a fixed season (June–August) and aerosol type (polluted dust) over two consecutive years (2006 and 2007). The CALIPSO models define six aerosol types: clean continental, clean marine, dust, polluted continental, polluted dust, and smoke, with 532-nm (1064 nm) extinction-to-backscatter ratios Sa of 35 (30), 20 (45), 40 (55), 70 (30), 65 (30), and 70 (40) sr, respectively. This paper presents the global distributions of the CALIPSO a...

Fully Automated Detection of Cloud and Aerosol Layers in the CALIPSO Lidar Measurements

Abstract Accurate knowledge of the vertical and horizontal extent of clouds and aerosols in the earth’s atmosphere is critical in assessing the planet’s radiation budget and for advancing human understanding of climate change issues. To retrieve this fundamental information from the elastic backscatter lidar data acquired during the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) mission, a selective, iterated boundary location (SIBYL) algorithm has been developed and deployed. SIBYL accomplishes its goals by integrating an adaptive context-sensitive profile scanner into an iterated multiresolution spatial averaging scheme. This paper provides an in-depth overview of the architecture and performance of the SIBYL algorithm. It begins with a brief review of the theory of target detection in noise-contaminated signals, and an enumeration of the practical constraints levied on the retrieval scheme by the design of the lidar hardware, the geometry of a space-based remote sensing pl...

The CALIPSO Lidar Cloud and Aerosol Discrimination: Version 2 Algorithm and Initial Assessment of Performance

Abstract The Cloud–Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) satellite was launched in April 2006 to provide global vertically resolved measurements of clouds and aerosols. Correct discrimination between clouds and aerosols observed by the lidar aboard the CALIPSO satellite is critical for accurate retrievals of cloud and aerosol optical properties and the correct interpretation of measurements. This paper reviews the theoretical basis of the CALIPSO lidar cloud and aerosol discrimination (CAD) algorithm, and describes the enhancements made to the version 2 algorithm that is used in the current data release (release 2). The paper also presents a preliminary assessment of the CAD performance based on one full day (12 August 2006) of expert manual classification and on one full month (July 2006) of the CALIOP 5-km cloud and aerosol layer products. Overall, the CAD algorithm works well in most cases. The 1-day manual verification suggests that the success rate is in the neighborh...

Airborne High Spectral Resolution Lidar for profiling aerosol optical properties

A compact, highly robust airborne High Spectral Resolution Lidar (HSRL) that provides measurements of aerosol backscatter and extinction coefficients and aerosol depolarization at two wavelengths has been developed, tested, and deployed on nine field experiments (over 650 flight hours). A unique and advantageous design element of the HSRL system is the ability to radiometrically calibrate the instrument internally, eliminating any reliance on vicarious calibration from atmospheric targets for which aerosol loading must be estimated. This paper discusses the design of the airborne HSRL, the internal calibration and accuracy of the instrument, data products produced, and observations and calibration data from the first two field missions: the Joint Intercontinental Chemical Transport Experiment--Phase B (INTEX-B)/Megacity Aerosol Experiment--Mexico City (MAX-Mex)/Megacities Impacts on Regional and Global Environment (MILAGRO) field mission (hereafter MILAGRO) and the Gulf of Mexico Atmospheric Composition and Climate Study/Texas Air Quality Study II (hereafter GoMACCS/TexAQS II).

Fully automated analysis of space-based lidar data: an overview of the CALIPSO retrieval algorithms and data products

The Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) satellite will be launched in April of 2005, and will make continuous measurements of the Earths atmosphere for the following three years. Retrieving the spatial and optical properties of clouds and aerosols from the CALIPSO lidar backscatter data will be confronted by a number of difficulties that are not faced in the analysis of ground-based data. Among these are the very large distance from the target, the high speed at which the satellite traverses the ground track, and the ensuing low signal-to-noise ratios that result from the mass and power restrictions imposed on space-based platforms. In this work we describe an integrated analysis scheme that employs a nested, multi-grid averaging technique designed to optimize tradeoffs between spatial resolution and signal-to-noise ratio. We present an overview of the three fundamental retrieval algorithms (boundary location, feature classification, and optical properties analysis), and illustrate their interconnections using data product examples that include feature top and base altitudes, feature type (i.e., cloud or aerosol), and layer optical depths.

CALIPSO Lidar Calibration Algorithms. Part I: Nighttime 532-nm Parallel Channel and 532-nm Perpendicular Channel

Abstract The Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) mission was launched in April 2006 and has continuously acquired collocated multisensor observations of the spatial and optical properties of clouds and aerosols in the earth’s atmosphere. The primary payload aboard CALIPSO is the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP), which makes range-resolved measurements of elastic backscatter at 532 and 1064 nm and linear depolarization ratios at 532 nm. CALIOP measurements are important in reducing uncertainties that currently limit understanding of the global climate system, and it is essential that these measurements be accurately calibrated. This work describes the procedures used to calibrate the 532-nm measurements acquired during the nighttime portions of the CALIPSO orbits. Accurate nighttime calibration of the 532-nm parallel-channel data is fundamental to the success of the CALIOP measurement scheme, because the nighttime calibration is used to infer...

Status and performance of the CALIOP lidar

The Cloud-Aerosol LIdar with Orthogonal Polarization (CALIOP) is the primary instrument on the CALIPSO satellite, which is scheduled to launch in 2005. CALIOP will provide profiles of total backscatter at two wavelengths, from which aerosol and cloud profiles will be derived. The instrument also measures the linear depolarization of the backscattered return, allowing discrimination of cloud phase and the identification of the presence of non-spherical aerosols. CALIOP is complete and has been tested in a ground-based configuration. This paper provides information on basic characteristics and performance of CALIOP.

Gravity wave models for the horizontal wave number spectra of atmospheric velocity and density fluctuations

Gravity wave models for the horizontal wave number spectra of atmospheric velocity and density fluctuations are derived by assuming that both saturated and unsaturated waves obey the polarization and dispersion relations and that the joint vertical wave number (m) and temporal frequency (ω) spectrum is separable. In the wave number region h* = ƒm*/N ≤ h ≤ m*, where h = (k2 + l2)1/2, k is the zonal wave number, l is the meridional wave number, ƒ is the inertial frequency, N is the buoyancy frequency, and m* is vertical wave number of the largest scale saturated wave, the model predicts that the two-dimensional isotropic spectrum is proportional to h−(p+1) and the one-dimensional spectra are proportional to k−p and l−p when the temporal spectrum is proportional to ω−p. In the saturation region h, k, l, m > m*, the spectra are proportional to h−(q+1), k−q and l−q when the saturated m spectrum is proportional to m−q. The joint (k, l, m) and (k, l, ω) model spectra are not separable. The models are consistent with existing observations of horizontal wave number spectra in the lower stratosphere and upper mesosphere.

Estimating random errors due to shot noise in backscatter lidar observations

We discuss the estimation of random errors due to shot noise in backscatter lidar observations that use either photomultiplier tube (PMT) or avalanche photodiode (APD) detectors. The statistical characteristics of photodetection are reviewed, and photon count distributions of solar background signals and laser backscatter signals are examined using airborne lidar observations at 532 nm using a photon-counting mode APD. Both distributions appear to be Poisson, indicating that the arrival at the photodetector of photons for these signals is a Poisson stochastic process. For Poisson- distributed signals, a proportional, one-to-one relationship is known to exist between the mean of a distribution and its variance. Although the multiplied photocurrent no longer follows a strict Poisson distribution in analog-mode APD and PMT detectors, the proportionality still exists between the mean and the variance of the multiplied photocurrent. We make use of this relationship by introducing the noise scale factor (NSF), which quantifies the constant of proportionality that exists between the root mean square of the random noise in a measurement and the square root of the mean signal. Using the NSF to estimate random errors in lidar measurements due to shot noise provides a significant advantage over the conventional error estimation techniques, in that with the NSF, uncertainties can be reliably calculated from or for a single data sample. Methods for evaluating the NSF are presented. Algorithms to compute the NSF are developed for the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations lidar and tested using data from the Lidar In-space Technology Experiment.

Observations of horizontal and vertical wave number spectra of gravity wave motions in the stratosphere and mesosphere over the mid‐Pacific

We present vertical and horizontal wave number spectra of density perturbations in the upper stratosphere (25–40 km) and the upper mesosphere (∼80–105 km) inferred from ∼45 hours of airborne Na/Rayleigh lidar observations in the vicinity of Hawaii. The mean density variance is (0.92%)2 for the upper stratosphere and (6.1%)2 and for the upper mesosphere. The mean vertical shear variance and Richardsons number for the upper mesosphere are (6.1%/km)2 and 0.69, respectively. The observed vertical wave number m and horizontal wave number k spectra exhibit power law dependencies. The mean slopes of the vertical wave number spectra are −2.5 for the upper stratosphere and −3.1 for the upper mesosphere. The mean slope of the horizontal wave number spectra is −1.9 for the upper mesosphere. The mean characteristic vertical and horizontal wavelengths inferred from the spectra are 13.4 and 5810 km, respectively. In all cases the magnitudes of the upper stratosphere m spectra are more than an order of magnitude smaller than those of the upper mesosphere at all observable scales. In the stratosphere the m spectra exhibit significant energy at low wave numbers that are less than the values expected for the characteristic wave numbers. This suggests the presence of gravity wave sources capable of generating significant energy at long vertical scales. The source is believed to be geostrophic adjustment of the jet stream. A large-amplitude, highly coherent, quasi-monochromatic, short horizontal scale (λh ∼ 16 km) oscillation in the Na layer was observed on the March 22 flight. The oscillation resembled short horizontal scale structures often observed in OH airglow images and was found to be consistent with a gravity wave interpretation.