Retrieval of cloud optical parameters from space-based backscatter lidar data

We present an approach to estimating the multiple-scattering (MS) contribution to lidar return signals from clouds recorded from space that enables us to describe in more detail the return formation at the depth where first orders of scattering dominate. Estimates made have enabled us to propose a method for correcting solutions of single-scattering lidar equations for the MS contribution. We also describe an algorithm for reconstructing the profiles of the cloud scattering coefficient and the optical thickness tau under conditions of a priori uncertainties. The approach proposed is illustrated with results for optical parameters of cirrus and stratiform clouds determined from return signals calculated by the Monte Carlo method as well as from return signals acquired with the American spaceborne lidar during the Lidar In-Space Technology Experiment (LITE).     

New airborne scanning lidar system: applications for atmospheric remote sensing

A new scanning airborne-aerosol lidar system that has the potential to be a valuable atmospheric remote-sensing tool has been developed. The system has the ability to scan both parallel and perpendicular to an aircraft's flight path, and this ability permits both the three-dimensional rendering of the aerosol structure below the aircraft and the measurement of aerosol extinction and optical depth. The system has been integrated into a NASA P-3 aircraft and during a recent flight was used to acquire excellent data with both scanning modes. The system design, the application of the across-track scanning data to the study of the atmospheric boundary layer, and the computation of optical depth derived from along-track scan data are reported.     

Laser remote sensing of atmospheric ammonia using a CO2 lidar system

A CO2 differential-absorption lidar system has been used for the remote sensing of ammonia in the atmosphere. For CO2 lidar returns backscattered from topographic targets at ranges up to 2.7 km, the path-averaged sensitivity of the DIAL system was 5 ppb of NH3 . Concentrations of atmospheric ammonia were found to vary during the day from undetectable levels (<5 ppb) to as high as 20 ppb, depending on temperature and humidity conditions.     

Determination by spaceborne backscatter lidar of the structural parameters of atmospheric scattering layers

Spaceborne active lidar systems are under development to give new insight into the vertical distribution of clouds and aerosols in the atmosphere and to provide new information on variables required for improvement of forecast models and for understanding the radiative and dynamic processes that are linked to the dynamics of climate change. However, when they are operated from space, lidar systems are limited by atmospheric backscattered signals that have low signal-to-noise ratios (SNRs) on optically thin targets. Therefore specific methods of analysis have to be developed to ensure accurate determination of the geometric and optical properties of scattering layers in the atmosphere. A first approach to retrieving the geometric properties of semitransparent cloud and aerosol layers is presented as a function of false-alarm and no-detection probabilities for a given SNR. Simulations show that the geometric properties of thin cirrus clouds and the altitude of the top of the unstable atmospheric boundary layer can be retrieved with standard deviations smaller than 150 m for a vertical resolution of the lidar system in the 50-100-m range and a SNR of 3. The altitudes of the top of dense clouds are retrieved with a precision in altitude of better than 50 m, as this retrieval corresponds to a higher SNR value. Such methods have an important potential application to future spaceborne lidar missions.     

Fiber-based lidar for atmospheric water-vapor measurements

The design and evaluation of a prototype fiber-based lidar system for autonomous measurement of atmospheric water vapor are presented. The system components are described, along with current limitations and options for improvement. Atmospheric measurements show good agreement with modeled signal returns from 400 to 1000 m but are limited below 400 m as a result of errors in signal processing caused by violation of the assumptions used in the derivation of the differential absorption lidar equation.     

Microphysical aerosol parameters from multiwavelength lidar

The hybrid regularization technique developed at the Institute of Mathematics of Potsdam University (IMP) is used to derive microphysical properties such as effective radius, surface-area concentration, and volume concentration, as well as the single-scattering albedo and a mean complex refractive index, from multiwavelength lidar measurements. We present the continuation of investigations of the IMP method. Theoretical studies of the degree of ill-posedness of the underlying model, simulation results with respect to the analysis of the retrieval error of microphysical particle properties from multiwavelength lidar data, and a comparison of results for different numbers of backscatter and extinction coefficients are presented. Our analysis shows that the backscatter operator has a smaller degree of ill-posedness than the operator for extinction. This fact underlines the importance of backscatter data. Moreover, the degree of ill-posedness increases with increasing particle absorption, i.e., depends on the imaginary part of the refractive index and does not depend significantly on the real part. Furthermore, an extensive simulation study was carried out for logarithmic-normal size distributions with different median radii, mode widths, and real and imaginary parts of refractive indices. The errors of the retrieved particle properties obtained from the inversion of three backscatter (355, 532, and 1064 nm) and two extinction (355 and 532 nm) coefficients were compared with the uncertainties for the case of six backscatter (400, 710, 800 nm, additionally) and the same two extinction coefficients. For known complex refractive index and up to 20% normally distributed noise, we found that the retrieval errors for effective radius, surface-area concentration, and volume concentration stay below approximately 15% in both cases. Simulations were also made with unknown complex refractive index. In that case the integrated parameters stay below approximately 30%, and the imaginary part of the refractive index stays below 35% for input noise up to 10% in both cases. In general, the quality of the retrieved aerosol parameters depends strongly on the imaginary part owing to the degree of ill-posedness. It is shown that under certain constraints a minimum data set of three backscatter coefficients and two extinction coefficients is sufficient for a successful inversion. The IMP algorithm was finally tested for a measurement case.     

Retrieval of atmospheric optical parameters from airborne flux measurements: application to the atmospheric correction of imagery

Methodologies that employ auxilliary flux data collected by upward- and downward-looking optical sensors to improve atmospheric corrections of airborne multispectral images are presented and evaluated. Such flux data often are collected in current airborne sensors to produce bidirectional reflectance factor (BRF) images and estimates of hemispherical-hemispherical reflectance. The fact that these images must then be corrected for atmospheric interference raises the question as to whether the auxilliary flux information can be employed to estimate some of the input parameters required by atmospheric correction models. Radiative transfer simulations are employed to demonstrate that the utilization of the downwelling and upwelling fluxes as a means of inferring intrinsic atmospheric optical information can be used to better characterize the local atmosphere and accordingly to improve the atmospheric corrections applied to the apparent BRF images.     

High spectral resolution lidar to measure optical scattering properties of atmospheric aerosols. 2: calibration and data analysis

The high spectral resolution lidar (HSRL) measures optical properties of atmospheric aerosols by interferometrically separating the elastic aerosol backscatter from the Doppler broadened molecular contribution. Calibration and data analysis procedures developed for the HSRL are described. Data obtained during flight evaluation testing of the HSRL system are presented with estimates of uncertainties due to instrument calibration. HSRL measurements of the aerosol scattering cross section are compared with in situ integrating nephelometer measurements.     

Multiple field of view lidar returns from atmospheric aerosols

As a laser pulse propagates into the atmosphere, it becomes broader in the lateral direction as a result of scattering by aerosols. The laser pulse may be described as the superposition of a central, unscattered component of reduced intensity and a surrounding scattered component. A multiple field of view lidar has been developed that makes simultaneous measurements of the backscattered power from the central pulse and multiply scattered power arising from the scattered component. Measurements from various types of atmospheric aerosols and precipitation are presented and compared with simulated returns. The results show how the multiply scattered signals are influenced by the distribution of the aerosols along the lidar path, the characteristic size of the aerosols, and the optical depth. It is shown that the multiple field of view lidar can provide meaningful, additional information about the aerosols that is not available from a conventional single field of view lidar.     

Ring-toric lens for focus-error sensing in optical data storage

We discuss the design and performance of diffractive ring-toric lenses for focus-error sensing in optical data storage. A ring-toric lens images a point source of light to a ring-shaped image. Focus-error sensing is accomplished by means of monitoring the change in ring radius: The ring expands in response to a diverging wave front, and the ring contracts in response to a converging wave front. We describe the use of a segmented phi detector to generate a focus-error signal (FES). We found that the FES slope, a measure of sensitivity to disk defocus, is higher for the ring-toric lenses described in this paper than for other techniques such as the astigmatic and the obscuration methods. We measured an FES slope of 0.7 per micrometer of disk defocus (mum(-1)). The corresponding theoretical FES slope is 0.96 mum(-1).     

Sensitivity of a lidar inversion algorithm to parameters relating atmospheric backscatter and extinction

A solution of the single-scattering lidar equation requires a relationship between the coefficients of backscatter beta(r) and extinction sigma(r) to be of the form beta(r) = C2sigma(r)k, where C2 and k are parameters independent of range r. The sensitivity of a particular lidar inversion algorithm to uncertainties in C2 and k is investigated using a measured lidar return which indicated the atmosphere to be essentially horizontally homogeneous during a reduced visibility condition. Starting with the measured power returned as a function of range, extinction coefficients and average visibilities are calculated using the inversion algorithm for different values of C2 and k and compared with those inferred from the lidar return using the slope method. The calculated extinction coefficients (and visibilities) were found to be extremely sensitive to uncertainties in C2. This questions the usefulness of the lidar inversion algorithm for aerosol extinction applications when the air mass characteristics change along the measurement path.     

Remote sensing of atmospheric turbulence and transverse winds from wave-front slope measurements from crossed optical paths

In the theory of atmospheric turbulence, the strength of the spatial variations of the index of refraction n is proportional to a parameter known as the atmospheric-structure constant. The atmosphericstructure constant is denoted C(2)(n)(z) and is a function of position along the optical path z. The characteristics of the temporal variations of the index of refraction are related to both C(2)(n)(z) and to the transverse wind velocity V(z). Current optical techniques for remotely sensing C(2)(n)(z) and V(z) rely primarily on the spatial or temporal cross-correlation properties of the intensity of the optical field. In the remote-sensing technique proposed here, we exploit the correlation properties of the wave-front slope measured from two point sources to obtain profiles of C(2)(n)(z) and V(z). The two sources are arranged to give crossed optical paths. The geometry of the crossed paths and the characteristics of the wave-front slope sensor determine the achievable resolution. The signal-to-noise ratio calculationsindicate the need for multiple measurements to obtain useful estimates of the desired quantities.     

Influence of atmospheric and system parameters on multiple scattering in spaceborne backscatter lidar measurements

We have simulated backscatter signals of spaceborne lidar systems with the help of a Monte Carlo model. Calculations were performed for various combinations of system parameters. As typical examples of atmospheric observation targets, two kinds of cirrus cloud and two kinds of aerosol were considered. Both total multiple scattering and the significance of individual higher scattering orders are discussed. For all cases, an approximate multiple scattering factor F was calculated that can be used to correct the single-scattering lidar equation to account also for multiple scattering.     

Microphysical particle parameters from extinction and backscatter lidar data by inversion with regularization: experiment

We present effective radius, volume, surface-area, and number concentrations as well as mean complex refractive index of tropospheric particle size distributions based on lidar measurements at six wavelengths. The parameters are derived by means of an inversion algorithm that has been specifically designed for the inversion of available optical data sets. The data were taken on 20 June and on 20 July 1997 during the Aerosol Characterization Experiment ACE 2 (North Atlantic/Portugal) and on 9 August 1998 during the Lindenberg Aerosol Characterization Experiment LACE 98 (Lindenberg/Germany). Measurements on 20 June 1997 were taken in a clean-marine boundary layer, and a large value of 0.64 mum for the effective radius, a low value of 1.45 for the real part, and a negligible imaginary part of the complex refractive index were found. The single-scatter albedo was 0.98 at 532 nm. It was derived from the particle parameters with Mie-scattering calculations. In contrast, the particles were less than 0.2 mum in effective radius in a continental-polluted aerosol layer on 20 July 1997. The real part of the complex refractive index was ~1.6; the imaginary part showed values near 0.03i. The single-scatter albedo was 0.84. On 9 August 1998 an elevated particle layer located from 3000 to 6000 m was observed, which had originated from an area of biomass burning in northwestern Canada. Here the effective radius was ~0.24 mum, the real part of the complex refractive index was above 1.6, the imaginary part was ~0.04i, and the single-scatter albedo was 0.81. Excellent agreement has been found with results based on sunphotometer and in situ measurements that were performed during the field campaigns.     

Microphysical particle parameters from extinction and backscatter lidar data by inversion with regularization: theory

A method is proposed that permits one to retrieve physical parameters of tropospheric particle size distributions, e.g., effective radius, volume, surface-area, and number concentrations, as well as the mean complex refractive index on a routine basis from backscatter and extinction coefficients at multiple wavelengths. The optical data in terms of vertical profiles are derived from multiple-wavelength lidar measurements at 355, 400, 532, 710, 800, and 1064 nm for backscatter data and 355 and 532 nm for extinction data. The algorithm is based on the concept of inversion with regularization. Regularization is performed by generalized cross-validation. This method does not require knowledge of the shape of the particle size distribution and can handle measurement errors of the order of 20%. It is shown that at least two extinction data are necessary to retrieve the particle parameters to an acceptable accuracy. Simulations with monomodal and bimodal logarithmic-normal size distributions show that it is possible to derive effective radius, volume, and surface-area concentrations to an accuracy of +/-50%, the real part of the complex refractive index to +/-0.05, and the imaginary part to +/-50%. Number concentrations may have errors larger than +/-50%.     

Microphysical particle parameters from extinction and backscatter lidar data by inversion with regularization: simulation

A sensitivity study with an inversion scheme that permits one to retrieve physical parameters of tropospheric particle size distributions, e.g., effective radius, volume, surface-area, and number concentrations, as well as the mean complex refractive index from backscatter and extinction coefficients at multiple wavelengths is presented. The optical data for the analysis are derived from Mie-scattering calculations for monomodal and bimodal logarithmic-normal distributions in the particle size range between 0.01 and 10 microm. The complex refractive index is taken between 1.33 and 1.8 in the real part and between 0 and 0.1 in the imaginary part. The choice of these parameters takes account of properties of optically active atmospheric particles. The wavelengths were chosen at 355, 400, 532, 710, 800, and 1064 nm for the backscatter and at 355 and 532 nm for the extinction data, which are the available wavelengths of the two lidar systems at the Institute for Tropospheric Research. Cases of erroneous optical data of the order of as much as 20%, an unknown refractive index, which may also be wavelength and size dependent, as well as the a priori unknown modality of the particle size distribution were considered. It is shown that both extinction channels are necessary for determining the above-mentioned parameters within reasonable limits, i.e., effective radius, surface-area, and volume concentrations to an accuracy of +/-50%, the real part of the complex refractive index to +/-0.1, and the imaginary part to +/-50%. The number concentration may have errors larger than 50%. The overall performance of the inversion scheme permits the evaluation of experimental data on a routine basis.     

Geometric method of reconstructing systems from experimental data

A method of reconstructing model dynamic evolution equations reduced to the central invariant manifold is described, which is based on an analysis of experimental data from controlled objects. The proposed method takes into account the group transformations of phase trajectories, which retain the topological equivalence of local regions.