Characterization of smoke and dust episode over West Africa: comparison of MERRA-2 modeling with multiwavelength Mie-Raman lidar observations.

Observations of multiwavelength Mie-Raman lidar taken during the SHADOW field campaign are used to analyze a smoke-dust episode over West Africa on 24-27 December 2015. For the case considered, the dust layer extended from the ground up to approximately 2000 m while the elevated smoke layer occurred in the 2500-4000 m range. The profiles of lidar measured backscattering, extinction coefficients, and depolarization ratios are compared with the vertical distribution of aerosol parameters provided by the Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2). The MERRA-2 model simulated the correct location of the near-surface dust and elevated smoke layers. The values of modeled and observed aerosol extinction coefficients at both 355 and 532 nm are also rather close. In particular, for the episode reported, the mean value of difference between the measured and modeled extinction coefficients at 355 nm is 0.01 km-1 with SD of 0.042 km-1. The model predicts significant concentration of dust particles inside the elevated smoke layer, which is supported by an increased depolarization ratio of 15 % observed in the center of this layer. The modeled at 355 nm the lidar ratio of 65 sr in the near-surface dust layer is close to the observed value (70 ± 10) sr. At 532 nm, however, the simulated lidar ratio (about 40 sr) is lower than measurements (55 ± 8 sr). The results presented demonstrate that the lidar and model data are complimentary and the synergy of observations and models is a key to improve the aerosols characterization.

Application of the lamp mapping technique for overlap function for Raman lidar systems.

Traditionally, the lidar water vapor mixing ratio (WVMR) is corrected for overlap using data from another instrument, such as a radiosonde. Here we introduce a new experimental method to determine the overlap function using the lamp mapping technique (LMT), which relies on the lidar optics and detection system. The LMT discussed here involves a standard halogen lamp being scanned over the aperture of a Raman lidar telescope in synchronization with the lidar detection system [Appl. Opt.50, 4622 (2011)APOPAI0003-693510.1364/AO.50.004622, Appl. Opt.53, 8538 (2014)APOPAI0003-693510.1364/AO.53.008535]. In this paper, we show results for a LMT-determined overlap function for individual channels, as well as a WVMR overlap function. We found that the LMT-determined WVMR overlap functions deviate within 5% of the traditional radiosonde-determined overlap.

Gluing for Raman lidar systems using the lamp mapping technique.

In the context of combined analog and photon counting (PC) data acquisition in a Lidar system, glue coefficients are defined as constants used for converting an analog signal into a virtual PC signal. The coefficients are typically calculated using Lidar profile data taken under clear, nighttime conditions since, in the presence of clouds or high solar background, it is difficult to obtain accurate glue coefficients from Lidar backscattered data. Here we introduce a new method in which we use the lamp mapping technique (LMT) to determine glue coefficients in a manner that does not require atmospheric profiles to be acquired and permits accurate glue coefficients to be calculated when adequate Lidar profile data are not available. The LMT involves scanning a halogen lamp over the aperture of a Lidar receiver telescope such that the optical efficiency of the entire detection system is characterized. The studies shown here involve two Raman lidar systems; the first from Howard University and the second from NASA/Goddard Space Flight Center. The glue coefficients determined using the LMT and the Lidar backscattered method agreed within 1.2% for the water vapor channel and within 2.5% for the nitrogen channel for both Lidar systems. We believe this to be the first instance of the use of laboratory techniques for determining the glue coefficients for Lidar data analysis.

Assessing the temperature dependence of narrow-band Raman water vapor lidar measurements: a practical approach.

Narrow-band detection of the Raman water vapor spectrum using the lidar technique introduces a concern over the temperature dependence of the Raman spectrum. Various groups have addressed this issue either by trying to minimize the temperature dependence to the point where it can be ignored or by correcting for whatever degree of temperature dependence exists. The traditional technique for performing either of these entails accurately measuring both the laser output wavelength and the water vapor spectral passband with combined uncertainty of approximately 0.01 nm. However, uncertainty in interference filter center wavelengths and laser output wavelengths can be this large or larger. These combined uncertainties translate into uncertainties in the magnitude of the temperature dependence of the Raman lidar water vapor measurement of 3% or more. We present here an alternate approach for accurately determining the temperature dependence of the Raman lidar water vapor measurement. This alternate approach entails acquiring sequential atmospheric profiles using the lidar while scanning the channel passband across portions of the Raman water vapor Q-branch. This scanning is accomplished either by tilt-tuning an interference filter or by scanning the output of a spectrometer. Through this process a peak in the transmitted intensity can be discerned in a manner that defines the spectral location of the channel passband with respect to the laser output wavelength to much higher accuracy than that achieved with standard laboratory techniques. Given the peak of the water vapor signal intensity curve, determined using the techniques described here, and an approximate knowledge of atmospheric temperature, the temperature dependence of a given Raman lidar profile can be determined with accuracy of 0.5% or better. A Mathematica notebook that demonstrates the calculations used here is available from the lead author.

Lamp mapping technique for independent determination of the water vapor mixing ratio calibration factor for a Raman lidar system.

We have investigated a technique that allows for the independent determination of the water vapor mixing ratio calibration factor for a Raman lidar system. This technique utilizes a procedure whereby a light source of known spectral characteristics is scanned across the aperture of the lidar system's telescope and the overall optical efficiency of the system is determined. Direct analysis of the temperature-dependent differential scattering cross sections for vibration and vibration-rotation transitions (convolved with narrowband filters) along with the measured efficiency of the system, leads to a theoretical determination of the water vapor mixing ratio calibration factor. A calibration factor was also obtained experimentally from lidar measurements and radiosonde data. A comparison of the theoretical and experimentally determined values agrees within 5%. We report on the sensitivity of the water vapor mixing ratio calibration factor to uncertainties in parameters that characterize the narrowband transmission filters, the temperature-dependent differential scattering cross section, and the variability of the system efficiency ratios as the lamp is scanned across the aperture of the telescope used in the Howard University Raman Lidar system.

Comments on "Accuracy of Raman lidar water vapor calibration and its applicability to long-term measurements".

In a recent publication, Leblanc and McDermid [Appl. Opt., 47, 5592 (2008)]APOPAI0003-693510.1364/AO.47.005592 proposed a hybrid calibration technique for Raman water vapor lidar involving a tungsten lamp and radiosondes. Measurements made with the lidar telescope viewing the calibration lamp were used to stabilize the lidar calibration determined by comparison with radiosonde. The technique provided a significantly more stable calibration constant than radiosondes used alone. The technique involves the use of a calibration lamp in a fixed position in front of the lidar receiver aperture. We examine this configuration and find that such a configuration likely does not properly sample the full lidar system optical efficiency. While the technique is a useful addition to the use of radiosondes alone for lidar calibration, it is important to understand the scenarios under which it will not provide an accurate quantification of system optical efficiency changes. We offer examples of these scenarios. Scanning of the full telescope aperture with the calibration lamp can circumvent most of these limitations. Based on the work done to date, it seems likely that the use of multiple calibration lamps in different fixed positions in front of the telescope may provide sufficient redundancy for long-term calibration needs. Further full-aperture scanning experiments, performed over an extended period of time, are needed to determine a "best practice" for the use of multiple calibration lamps in the hybrid technique.

Validation of the Raman lidar algorithm for quantifying aerosol extinction.

To calculate aerosol extinction from Raman lidar data, it is necessary to evaluate the derivative of a molecular Raman signal with respect to range. The typical approach taken in the lidar community is to make an a priori assumption about the functional behavior of the data to calculate the derivative. It has previously been shown that the use of the chi-squared technique to determine the most likely functional behavior of the data prior to actually calculating the derivative eliminates the need for making a priori assumptions. Here that technique is validated through numerical simulation and by application to a significant body of Raman lidar measurements. In general, we show that the chi-squared approach for evaluating extinction yields lower extinction uncertainty than traditional techniques. We also use the technique to study the feasibility of developing a general characterization of the extinction uncertainty that could permit the uncertainty in Raman lidar aerosol extinction measurements to be estimated accurately without the need of the chi-squared technique.

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.

Information content of multiwavelength lidar data with respect to microphysical particle properties derived from eigenvalue analysis.

The multiwavelength Raman lidar technique in combination with sophisticated inversion algorithms has been recognized as a new tool for deriving information about the microphysical properties of atmospheric aerosols. The input optical parameter sets, provided by respective aerosol Raman lidars, are at the theoretical lower limit at which these inversion algorithms work properly. For that reason there is ongoing intense discussion of the accuracy of these inversion methods and the possibility of simultaneous retrieval of the particle size distribution and the complex refractive index. We present results of the eigenvalue analysis, used to study the information content of multiwavelength lidar data with respect to microphysical particle properties. Such an analysis provides, on a rather mathematical basis, more insight into the limitations of these inversion algorithms regarding the accuracy of the retrieved parameters. We show that the effective radius may be retrieved to 50% accuracy and the real and imaginary part of the complex refractive index to +/- 0.05 and +/- 0.005i, if the imaginary part is < 0.02i. These results are in accordance with the classic approach of simulation studies with synthetic particle size distributions. Major difficulties are found with a particle effective radius of < 0.15 microm. In that case the complex refractive index may not be derived with sufficient accuracy. The eigenvalue analysis also shows that the accuracy of the derived parameters degrades if the imaginary part is > 0.02i. Furthermore it shows the importance of the simultaneous use of backscatter and extinction coefficients for the retrieval of microphysical parameters.

Active Raman sounding of the earth's water vapor field.

The typically weak cross-sections characteristic of Raman processes has historically limited their use in atmospheric remote sensing to nighttime application. However, with advances in instrumentation and techniques, it is now possible to apply Raman lidar to the monitoring of atmospheric water vapor, aerosols and clouds throughout the diurnal cycle. Upper tropospheric and lower stratospheric measurements of water vapor using Raman lidar are also possible but are limited to nighttime and require long integration times. However, boundary layer studies of water vapor variability can now be performed with high temporal and spatial resolution. This paper will review the current state-of-the-art of Raman lidar for high-resolution measurements of the atmospheric water vapor, aerosol and cloud fields. In particular, we describe the use of Raman lidar for mapping the vertical distribution and variability of atmospheric water vapor, aerosols and clouds throughout the evolution of dynamic meteorological events. The ability of Raman lidar to detect and characterize water in the region of the tropopause and the importance of high-altitude water vapor for climate-related studies and meteorological satellite performance are discussed.

Numerical simulation of light backscattering by spheres with off-center inclusion. Application to the lidar case.

A Mie backscattering model for spherical particles with off-center inclusion has been developed and tested. The program is capable of dealing with size parameter values up to approximately 1000, thus allowing one to simulate the optical behavior of a large variety of atmospheric aerosols, as well as cloud and precipitation particles. On the basis of this model, we simulated the optical properties of polydisperse composite atmospheric particles as observed by ground-based and airborne lidar systems. We have characterized optical properties in terms of host and inclusion radii, considering water particles with different composition inclusions. The performed modeling provides some insight into the so-called lidar bright- and dark-band phenomenon.

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.

Examination of the traditional Raman lidar technique. I. Evaluating the temperature-dependent lidar equations.

The essential information required for the analysis of Raman lidar water vapor and aerosol data acquired by use of a single laser wavelength is compiled here and in a companion paper [Appl. Opt. 42, 2593 (2003)]. Various details concerning the evaluation of the lidar equations when Raman scattering is measured are covered. These details include the influence of the temperature dependence of both pure rotational and vibrational-rotational Raman scattering on the lidar profile. The full temperature dependence of the Rayleigh-Mie and Raman lidar equations are evaluated by use of a new form of the lidar equation where all the temperature dependence is carried in a single term. The results indicate that, for the range of temperatures encountered in the troposphere, the magnitude of the temperature-dependent effect can reach 10% or more for narrowband Raman water-vapor measurements. Also, the calculation of atmospheric transmission, including the effects of depolarization, is examined carefully. Various formulations of Rayleigh cross-section determination commonly used in the lidar field are compared and reveal differences of as much as 5% among the formulations. The influence of multiple scattering on the measurement of aerosol extinction with the Raman lidar technique is considered, as are several photon pulse pileup-correction techniques.

Examination of the traditional raman lidar technique. II. Evaluating the ratios for water vapor and aerosols.

In a companion paper [Appl. Opt. 42, 2571 (2003)] the temperature dependence of Raman scattering and its influence on the Raman and Rayleigh-Mie lidar equations were examined. New forms of the lidar equation were developed to account for this temperature sensitivity.Here those results are used to derive the temperature-dependent forms of the equations for the water vapor mixing ratio, the aerosol scattering ratio, the aerosol backscatter coefficient, and the extinction-to-backscatter ratio. The error equations are developed, the influence of differential transmission is studied, and several laser sources are considered in the analysis. The results indicate that the temperature functions become significant when narrowband detection is used. Errors of 5% and more can be introduced into the water-vapor mixing ratio calculation at high altitudes, and errors larger than 10% are possible for calculations of aerosol scattering ratio and thus of aerosol backscatter coefficient and of extinction-to-backscatter ratio.

Angle- and size-dependent characteristics of incoherent Raman and fluorescent scattering by microspheres. 1. General expressions.

A dipole model is used to simulate incoherent Raman and fluorescent scattering by microspheres. The use of the addition theorem for spherical harmonics circumvents the need to evaluate double sums in the final formulas, thereby drastically reducing computational effort. Special attention is paid to consideration of backscattering geometry, which is important for lidar applications. The formulas derived for backscattering geometry decrease the computation time for size parameter x approdximately 100 by a factor of 200 compared with the time for calculations performed at other angles.

Angle- and size-dependent characteristics of incoherent Raman and fluorescent scattering by microspheres. 2. Numerical simulation.

The results of numerical simulation of inelastic scattering by microspheres with the use of a dipole model are presented. The formulas that are derived speed up the computation, thereby permitting larger-sized microspheres to be studied. The angular scattering cross section and depolarization are calculated for a wide range of size parameters as well as for different orientations of incident wave polarization. Calculations performed with small incremental changes in size permit the influence of morphology-dependent resonance (MDR) on the power and angular distribution of scattered radiation to be studied. TM and TE types of MDR produce enhanced scattering of the incident wave with vertical and horizontal polarization; the corresponding shape of the phase function becomes oscillatory. Special attention is paid to the simulation of backward scattering by water droplets, which is important for Raman lidar applications.

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.