Information content of data measured with a multiple-field-of-view lidar

Lidars with multiple fields of view (MFOVs) are promising tools for gaining information on cloud particle size. We perform a study of the information content of MFOV lidar data with the use of eigenvalue analysis. The approach we have developed permits an understanding of the main features of MFOV lidars and provides a way to relate the accuracy of particle size estimation with the measurement uncertainty and the scattering geometry such as the cloud-base height and the lidar sounding depth. Second-order scattering computations are performed for an extended range of particle sizes and for a wide range of lidar fields of view (FOVs). The results obtained allow us to specify the areas of possible applications of these lidars in cloud studies. Comparison of results obtained with polarized and cross-polarized scattered components demonstrate that the cross-polarized signal should provide a more stable retrieval and is preferable when double scattering is highly dominant. Our analysis allows for the estimation of the optimal number of FOVs in the system and their angular distribution, so this work can be a useful tool for practical MFOV lidar design.     

Standoff determination of the particle size and concentration of small optical depth clouds based on double scattering measurements: concept and experimental validation with bioaerosols

A multiple-field-of-view (MFOV) lidar is used to characterize size and optical depth of low concentration of bioaerosol clouds. The concept relies on the measurement of the forward scattered light by using the background aerosols at various distances at the back of a subvisible cloud. It also relies on the subtraction of the background aerosol forward scattering contribution and on the partial attenuation of the first-order backscattering. The validity of the concept developed to retrieve the effective diameter and the optical depth of low concentration bioaerosol clouds with good precision is demonstrated using simulation results and experimental MFOV lidar measurements. Calculations are also done to show that the method presented can be extended to small optical depth cloud retrieval.     

Retrieval of droplet-size density distribution from multiple-field-of-view cross-polarized lidar signals: theory and experimental validation

Multiple-field-of-view (MFOV) secondary-polarization lidar signals are used to calculate the particle-size density distribution (PSD) at the base of a cloud. At the cloud base, multiple scattering is weak and single backscattering is predominant by many orders of magnitude. Because secondary polarization is a direct measure of multiple scattering, it is therefore advantageous to use secondary polarization. A mathematical relation among the PSD, the lidar fields of view, the scattering angles, and the angular depolarization is derived to facilitate use of secondary polarization. The model is supported by experimental MFOV lidar measurements carried out in a controlled environment, and its limitations and restrictions are discussed.     

Error analysis of Raman differential absorption lidar ozone measurements in ice clouds

A formalism for the error treatment of lidar ozone measurements with the Raman differential absorption lidar technique is presented. In the presence of clouds wavelength-dependent multiple scattering and cloud-particle extinction are the main sources of systematic errors in ozone measurements and necessitate a correction of the measured ozone profiles. Model calculations are performed to describe the influence of cirrus and polar stratospheric clouds on the ozone. It is found that it is sufficient to account for cloud-particle scattering and Rayleigh scattering in and above the cloud; boundary-layer aerosols and the atmospheric column below the cloud can be neglected for the ozone correction. Furthermore, if the extinction coefficient of the cloud is ?0.1 km(-1), the effect in the cloud is proportional to the effective particle extinction and to a particle correction function determined in the limit of negligible molecular scattering. The particle correction function depends on the scattering behavior of the cloud particles, the cloud geometric structure, and the lidar system parameters. Because of the differential extinction of light that has undergone one or more small-angle scattering processes within the cloud, the cloud effect on ozone extends to altitudes above the cloud. The various influencing parameters imply that the particle-related ozone correction has to be calculated for each individual measurement. Examples of ozone measurements in cirrus clouds are discussed.     

Multiple-scattering lidar equation

A multiple-scattering lidar equation is derived from a phenomenological representation of the scattering processes. The contributions are separated into the unscattered, singly scattered, and multiply scattered illumination of the scattering volume, a single backscattering reflection from that volume, and the unscattered and multiply scattered propagation back to the receiver. The equation is obtained in the form of analytic expressions that explicitly show the signal dependence on the extinction coefficient, the effective particle size, the range, and the receiver field of view. Consistent agreement is found with Monte Carlo calculations and published laboratory measurements. Numerical simulations demonstrate that measurements made at three or more fields of view can be inverted to solve for the extinction coefficient and the effective particle radius. The multiple scatterings taken into account in the proposed equation are the small-angle diffraction scatterings; the wide-angle scatterings caused by refraction and reflection are considered lost, except for one backscattering at an angle close to 180°. Consequently, the equation is applicable to cases in which the projection of the lidar receiver field of view on the cloud is of the order of the angular width of the diffraction peak of the phase function times the penetration depth into the cloud.     

Inversion with regularization for the retrieval of tropospheric aerosol parameters from multiwavelength lidar sounding

We present an inversion algorithm for the retrieval of particle size distribution parameters, i.e., mean (effective) radius, number, surface area, and volume concentration, and complex refractive index from multiwavelength lidar data. In contrast to the classical Tikhonov method, which accepts only that solution for which the discrepancy reaches its global minimum, in our algorithm we perform the averaging of solutions in the vicinity of this minimum. This averaging stabilizes the underlying ill-posed inverse problem, particularly with respect to the retrieval of number concentration. Results show that, for typical tropospheric particles and 10% error in the optical data, the mean radius could be retrieved to better than 20% from a lidar on the basis of a Nd:YAG laser, which provides a combination of backscatter coefficients at 355, 532, and 1064 nm and extinction coefficients at 355 and 532 nm. The accuracy is improved if the lidar is also equipped with a hydrogen Raman shifter. In this case two additional backscatter coefficients at 416 and 683 nm are available. The combination of two extinction coefficients and five backscatter coefficients then allows one to retrieve not only averaged aerosol parameters but also the size distribution function. There was acceptable agreement between physical particle properties obtained from the evaluation of multiwavelength lidar data taken during the Lindenberg Aerosol Characterization Experiment in 1998 (LACE 98) and in situ data, which were taken aboard aircraft.     

Multiple-scattering-based lidar retrieval: method and results of cloud probings

Recent developments in the search for a practical method of exploiting the multiple-scattering contributions to lidar returns are consolidated in a robust retrieval algorithm. The theoretical basis is the small-angle diffusion approximation. This implies that the algorithm is limited to media of sufficient optical thickness to generate measurable multiple scattering and to geometries for which the receiver's footprint diameter is less than the scattering mean free path. The primary retrieval products are the range-resolved extinction coefficient and the effective particle diameter from which secondary products such as the particle volume mixing ratio and the extinction at other wavelengths can be calculated. We recall briefly earlier validation tests and present new data and analysis that demonstrate and quantify the solutions' accuracy. The results show that systematic lidar probings with the proposed multiple-scattering technique can provide valuable physical information on cloud formation and evolution.     

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.     

Inversion of multiwavelength Raman lidar data for retrieval of bimodal aerosol size distribution

We report on the feasibility of deriving microphysical parameters of bimodal particle size distributions from Mie-Raman lidar based on a triple Nd:YAG laser. Such an instrument provides backscatter coefficients at 355, 532, and 1064 nm and extinction coefficients at 355 and 532 nm. The inversion method employed is Tikhonov's inversion with regularization. Special attention has been paid to extend the particle size range for which this inversion scheme works to approximately 10 microm, which makes this algorithm applicable to large particles, e.g., investigations concerning the hygroscopic growth of aerosols. Simulations showed that surface area, volume concentration, and effective radius are derived to an accuracy of approximately 50% for a variety of bimodal particle size distributions. For particle size distributions with an effective radius of < 1 microm the real part of the complex refractive index was retrieved to an accuracy of +/- 0.05, the imaginary part was retrieved to 50% uncertainty. Simulations dealing with a mode-dependent complex refractive index showed that an average complex refractive index is derived that lies between the values for the two individual modes. Thus it becomes possible to investigate external mixtures of particle size distributions, which, for example, might be present along continental rims along which anthropogenic pollution mixes with marine aerosols. Measurement cases obtained from the Institute for Tropospheric Research six-wavelength aerosol lidar observations during the Indian Ocean Experiment were used to test the capabilities of the algorithm for experimental data sets. A benchmark test was attempted for the case representing anthropogenic aerosols between a broken cloud deck. A strong contribution of particle volume in the coarse mode of the particle size distribution was found.     

Distortions of the extinction coefficient profile caused by systematic errors in lidar data

The influence of lidar data systematic errors on the retrieved particulate extinction coefficient profile in clear atmospheres is investigated. Particularly, two sources of the extinction coefficient profile distortions are analyzed: (1) a zero-line offset remaining after subtraction of an inaccurately determined signal background component and (2) a far-end incomplete overlap due to poor adjustment of the lidar system optics. Inversion results for simulated lidar signals, obtained with the near- and far-end solutions, are presented that show advantages of the near-end solution for clear atmospheres.     

Multiply scattered aerosol lidar returns: inversion method and comparison with in situ measurements

A novel aerosol lidar inversion method based on the use of multiple-scattering contributions measured by a multiple-field-of-view receiver is proposed. The method requires assumptions that restrict applications to aerosol particles large enough to give rise to measurable multiple scattering and depends on parameters that must be specified empirically but that have an uncertainty range of much less than the boundary value and the backscatter-to-extinction ratio of the conventional single-scattering inversion methods. The proposed method is applied to cloud measurements. The solutions obtained are the profiles of the scattering coefficient and the effective diameter of the cloud droplets. With mild assumptions on the form of the function, the full-size distribution is estimated at each range position from which the extinction coefficient at any visible and infrared wavelength and the liquid water content can be determined. Typical results on slant-path-integrated optical depth, vertical extinction profiles, and fluctuation statistics are compared with in situ data obtained in two field experiments. The inversion works well in all cases reported here, i.e., for water clouds at optical depths between ~0.1 and ~4.     

Analysis of lidar backscatter profiles in optically thin clouds

The solution of the lidar equation for profiles of backscatter and extinction in optically thin clouds is constrained by values of the cloud transmittance determined from the elastically scattered lidar signals below and above the cloud. The method is extended to those cases in which an aerosol layer lies below or above the cloud layer. Examples are given in both cases. An analytical expression for the average lidar ratio in the cloud is derived for those cases in which molecular scattering is significant.     

Raman lidar observations of cloud liquid water

We report the design and the performances of a Raman lidar for long-term monitoring of tropospheric aerosol backscattering and extinction coefficients, water vapor mixing ratio, and cloud liquid water. We focus on the system's capabilities of detecting Raman backscattering from cloud liquid water. After describing the system components, along with the current limitations and options for improvement, we report examples of observations in the case of low-level cumulus clouds. The measurements of the cloud liquid water content, as well as the estimations of the cloud droplet effective radii and number densities, obtained by combining the extinction coefficient and cloud water content within the clouds, are critically discussed.     

Determination of cloud microphysical properties by laser backscattering and extinction measurements

The extinction and backscattering of 514-nm laser radiation in polydisperse water droplet clouds have been studied in the laboratory. Three cloud size distributions with modal diameters of 0.02, 5, and 12 microm have been investigated. The relationships between the cloud optical parameters (attenuation coefficient sigma and volume backscattering coefficient beta(pi)) and the cloud water content C have been measured for each size distribution. It has been found that a linear relationship exists between sigma and C and between beta(pi) and C for cloud water content values up to 3 g/m3. The linear relationships obtained, however, have slopes which depend on the droplet size distribution. For a given water content both sigma and beta(pi) increase as the modal diameter decreases. The measured data are compared with existing theoretical analyses and discussed in terms of their application to lidar measurements of atmospheric clouds. It is concluded that the empirical information obtained can serve as a basis for quantitative lidar measurements.     

Significance of multiple scattering from tropospheric aerosols for ground-based backscatter lidar measurements

The influence of multiple scattering on the retrieval of extinction coefficients of tropospheric aerosols from ground-based backscatter lidar measurements is numerically modeled. In a first step, lidar returns are computed by means of a Monte Carlo code for model atmospheres with different aerosol types and different extinction coefficient profiles. In so doing, synthetic lidar signals with and without multiple scattering can be simulated. In a second step, both types of signal are inverted by the most frequently used analytical solution, which, however, is based on the single-scatter assumption. From a comparison of the results, the error of the retrieved aerosol-extinction profiles can be quantitatively determined. It was found that the contribution of multiply scattered photons to the lidar signals is typically below 10% and never exceeds 20%. The relative errors of the retrieved aerosol-extinction profile in the planetary boundary layer are still smaller; they were determined to be less than 3% for all aerosol types, even for extinction coefficients as large as 3.9 km(-1). Thus, for ground-based lidar measurements and typical meteorological conditions, errors caused by neglecting multiple scattering are by far less significant than other errors in lidar data evaluation.     

沙尘粒子的物理化学与光学特性

利用激光雷达、卫星、太阳光度计、地面能见度,结合分析化学等资料研究了长程运输到台湾地区的各种粒子的物理化学与光学特性。激光雷达与太阳光度计可以得到沙尘与云的高度分布、大气光学厚度、消光以及粒子大小分布等的信息。地面测站的化学分析等显示沙尘具有增强的水溶性离子:如Na+、K+、Mg2+、Ca2+及Cl-、NO3-、SO42-等。由于这些物质的亲水性,使它们成为云的凝结核。根据粒径分布随着湿度增加而变大,及在一些沙尘事件中观察到云的例案可以知道它们对气候的影响。     

Sensitivity of the lidar solution to errors of the aerosol backscatter-to-extinction ratio: influence of a monotonic change in the aerosol extinction coefficient

Unlike other errors in the lidar equation solution for the two-component atmosphere, the error of the measured aerosol extinction coefficient caused by inaccuracies in the assumed aerosol backscatter-to-extinction ratios significantly depends on the aerosol spatial inhomogeneity. In a slightly nonhomogeneous atmosphere, an incorrect value in the assumed aerosol backscatter-to-extinction ratio does not significantly corrupt the measurement result, whereas in an atmosphere with a large monotonic change of the aerosol extinction [e.g., in the lower troposphere], the incorrect value yields a large distortion of the retrieved extinction-coefficient profile. In the latter case, even the far-end solution can produce a large error in the retrieved extinction coefficient. The analytical formulas for the determination of the range errors, obtained for the Klett and the optical-depth solutions, show that these errors significantly depend on the method of the boundary-condition determination. Distortions of the retrieved aerosol extinction profiles are, in general, larger if the assumed aerosol backscatter-to-extinction ratio is underestimated in relation to the real value.     

Spaceborne profiling of atmospheric temperature and particle extinction with pure rotational Raman lidar and of relative humidity in combination with differential absorption lidar: performance simulations

The performance of a spaceborne temperature lidar based on the pure rotational Raman (RR) technique in the UV has been simulated. Results show that such a system deployed onboard a low-Earth-orbit satellite would provide global-scale clear-sky temperature measurements in the troposphere and lower stratosphere with precisions that satisfy World Meteorological Organization (WMO) threshold observational requirements for numerical weather prediction and climate research applications. Furthermore, nighttime temperature measurements would still be within the WMO threshold observational requirements in the presence of several cloud structures. The performance of aerosol extinction measurements from space, which can be carried out simultaneously with temperature measurements by RR lidar, is also assessed. Furthermore, we discuss simulations of relative humidity measurements from space obtained from RR temperature measurements and water-vapor data measured with the differential absorption lidar (DIAL) technique.     

Lidar inelastic multiple-scattering parameters of cirrus particle ensembles determined with geometrical-optics crystal phase functions

Multiple-scattering correction factors for cirrus particle extinction coefficients measured with Raman and high spectral resolution lidars are calculated with a radiative-transfer model. Cirrus particle-ensemble phase functions are computed from single-crystal phase functions derived in a geometrical-optics approximation. Seven crystal types are considered. In cirrus clouds with height-independent particle extinction coefficients the general pattern of the multiple-scattering parameters has a steep onset at cloud base with values of 0.5-0.7 followed by a gradual and monotonic decrease to 0.1-0.2 at cloud top. The larger the scattering particles are, the more gradual is the rate of decrease. Multiple-scattering parameters of complex crystals and of imperfect hexagonal columns and plates can be well approximated by those of projected-area equivalent ice spheres, whereas perfect hexagonal crystals show values as much as 70% higher than those of spheres. The dependencies of the multiple-scattering parameters on cirrus particle spectrum, base height, and geometric depth and on the lidar parameters laser wavelength and receiver field of view, are discussed, and a set of multiple-scattering parameter profiles for the correction of extinction measurements in homogeneous cirrus is provided.     

Corona-producing ice clouds: a case study of a cold mid-latitude cirrus layer

A high (14.0-km), cold (-71.0 degrees C) cirrus cloud was studied by ground-based polarization lidar and millimeter radar and aircraft probes on the night of 19 April 1994 from the Cloud and Radiation Testbed site in northern Oklahoma. A rare cirrus cloud lunar corona was generated by this 1-2-km-deep cloud, thus providing an opportunity to measure the composition in situ, which had previously been assumed only on the basis of lidar depolarization data and simple diffraction theory for spheres. In this case, corona ring analysis indicated an effective particle diameter of ~22 mum. A variety of in situ data corroborates the approximate ice-particle size derived from the passive retrieval method, especially near the cloud top, where impacted cloud samples show simple solid crystals. The homogeneous freezing of sulfuric acid droplets of stratospheric origin is assumed to be the dominant ice-particle nucleation mode acting in corona-producing cirrus clouds. It is speculated that this process results in a previously unrecognized mode of acid-contaminated ice-particle growth and that such small-particle cold cirrus clouds are potentially a radiatively distinct type of cloud.