Aerosol lidar intercomparison in the framework of the EARLINET project. 3. Raman lidar algorithm for aerosol extinction, backscatter, and lidar ratio

An intercomparison of the algorithms used to retrieve aerosol extinction and backscatter starting from Raman lidar signals has been performed by 11 groups of lidar scientists involved in the European Aerosol Research Lidar Network (EARLINET). This intercomparison is part of an extended quality assurance program performed on aerosol lidars in the EARLINET. Lidar instruments and aerosol backscatter algorithms were tested separately. The Raman lidar algorithms were tested by use of synthetic lidar data, simulated at 355, 532, 386, and 607 nm, with realistic experimental and atmospheric conditions taken into account. The intercomparison demonstrates that the data-handling procedures used by all the lidar groups provide satisfactory results. Extinction profiles show mean deviations from the correct solution within 10% in the planetary boundary layer (PBL), and backscatter profiles, retrieved by use of algorithms based on the combined Raman elastic-backscatter lidar technique, show mean deviations from solutions within 20% up to 2 km. The intercomparison was also carried out for the lidar ratio and produced profiles that show a mean deviation from the solution within 20% in the PBL. The mean value of this parameter was also calculated within a lofted aerosol layer at higher altitudes that is representative of typical layers related to special events such as Saharan dust outbreaks, forest fires, and volcanic eruptions. Here deviations were within 15%.     

Aerosol lidar intercomparison in the framework of the EARLINET project. 1. Instruments

In the framework of the European Aerosol Research Lidar Network to Establish an Aerosol Climatology (EARLINET), 19 aerosol lidar systems from 11 European countries were compared. Aerosol extinction or backscatter coefficient profiles were measured by at least two systems for each comparison. Aerosol extinction coefficients were derived from Raman lidar measurements in the UV (351 or 355 nm), and aerosol backscatter profiles were calculated from pure elastic backscatter measurements at 351 or 355, 532, or 1064 nm. The results were compared for height ranges with high and low aerosol content. Some systems were additionally compared with sunphotometers and starphotometers. Predefined maximum deviations were used for quality control of the results. Lidar systems with results outside those limits could not meet the quality assurance criterion. The algorithms for deriving aerosol backscatter profiles from elastic lidar measurements were tested separately, and the results are described in Part 2 of this series of papers [Appl. Opt. 43, 977-989 (2004)]. In the end, all systems were quality assured, although some had to be modified to improve their performance. Typical deviations between aerosol backscatter profiles were 10% in the planetary boundary layer and 0.1 x 10(-6) m(-1) sr(-1) in the free troposphere.     

Characteristics of aerosol optical properties in pollution and Asian dust episodes over Beijing, China

Aerosol optical properties were continuously measured with the National Institute for Environmental Studies (NIES) compact Raman lidar over Beijing, China, from 15 to 31 December 2007. The results indicated that in a moderate pollution episode, the averaged aerosol extinction below 1 km height was 0.39+/-0.15 km(-1) and the lidar ratio was 60.8+/-13.5 sr; in heavy pollution episode, they were 1.97+/-0.91 km(-1) and 43.7+/-8.3 sr; in an Asian dust episode, they were 0.33+/-0.11 km(-1) and 38.3+/-9.8 sr. The total depolarization ratio was mostly below 10% in the pollution episode, whereas it was larger than 20% in the Asian dust episode. The distinct characteristics of aerosol optical properties in moderate and heavy pollution episodes were attributed to the difference in air mass trajectory and the ambient atmospheric conditions such as relative humidity.     

Particle backscatter, extinction, and lidar ratio profiling with Raman lidar in south and north China

Aerosol Raman lidar observations of profiles of the particle extinction and backscatter coefficients and the respective extinction-to-backscatter ratio (lidar ratio) were performed under highly polluted conditions in the Pearl River Delta (PRD) in southern China in October 2004 and at Beijing during a clear period with moderately polluted to background aerosol conditions in January 2005. The anthropogenic haze in the PRD is characterized by volume light-extinction coefficients of particles ranging from approximately 200 to 800 Mm(-1) and lidar ratios mostly between 40 and 55 sr (average of 47+/-6 sr). Almost clean air masses were observed throughout the measurements of the Beijing campaign. These air masses originated from arid desert-steppe-like regions (greater Gobi area). Extinction values usually varied between 100 and 300 Mm(-1), and the lidar ratios were considerably lower (compared with PRD values) with values mostly from 30 to 45 sr (average of 38+/-7 sr). Gobi dust partly influenced the observations. Unexpectedly low lidar ratios of approximately 25 sr were found for a case of background aerosol with a low optical depth of 0.05. The low lidar ratios are consistent with Mie-scattering calculations applied to ground-based observations of particle size distributions.     

Ice clouds and Asian dust studied with lidar measurements of particle extinction-to-backscatter ratio, particle depolarization, and water-vapor mixing ratio over Tsukuba

The tropospheric particle extinction-to-backscatter ratio, the depolarization ratio, and the water-vapor mixing ratio were measured by use of a Raman lidar and a polarization lidar during the Asian dust seasons in 2001 and 2002 in Tsukuba, Japan. The apparent (not corrected for multiple-scattering effects) extinction-to-backscatter ratios (Sp) showed a dependence on the relative humidity with respect to ice (RHice) obtained from the lidar-derived water-vapor mixing ratio and radiosonde-derived temperature; they were mostly higher than 30 sr in dry air (RHice < 50%), whereas they were mostly lower than 30 sr in ice-supersaturated air (RHice > or = 100%), where the apparent extinction coefficients were larger than 0.036 km(-1). Both regions showed mean particle depolarization ratios of 20%-22%. Comparisons with theoretical calculations and the previous experiments suggest that the observed dependence of Sp on RHice is attributed to the difference in the predominant particles: nonspherical aerosols (mainly the Asian dust) in dry air and cloud particles in ice-supersaturated air.     

Particle backscatter and extinction profiling with the spaceborne high-spectral-resolution Doppler lidar ALADIN: methodology and simulations

The European Space Agency will launch the Atmospheric Laser Doppler Instrument (ALADIN) for global wind profile observations in the near future. The potential of ALADIN to measure the optical properties of aerosol and cirrus, as well, is investigated based on simulations. A comprehensive data analysis scheme is developed that includes (a) the correction of Doppler-shifted particle backscatter interference in the molecular backscatter channels (cross-talk effect), (b) a procedure that allows us to check the quality of the cross-talk correction, and (c) the procedures for the independent retrieval of profiles of the volume extinction and backscatter coefficients of particles considering the height-dependent ALADIN signal resolution. The error analysis shows that the particle backscatter and extinction coefficients, and the corresponding extinction-to-backscatter ratio (lidar ratio), can be obtained with an overall (systematic+statistical) error of 10%-15%, 15%-30%, and 20%-35%, respectively, in tropospheric aerosol and dust layers with extinction values from 50 to 200 Mm(-1); 700-shot averaging (50 km horizontal resolution) is required. Vertical signal resolution is 500 m in the lower troposphere and 1000 m in the free troposphere. In cirrus characterized by extinction coefficients of 200 Mm(-1) and an optical depth of >0.2, backscatter coefficients, optical depth, and column lidar ratios can be obtained with 25%-35% relative uncertainty and a horizontal resolution of 10 km (140 shots). In the stratosphere, only the backscatter coefficient of aerosol layers and polar stratospheric clouds can be retrieved with an acceptable uncertainty of 15%-30%. Vertical resolution is 2000 m.     

Two-wavelength lidar inversion algorithm for a two-component atmosphere

A method for the boundary-value determination of aerosol extinction profiles from backscatter lidar measurements is presented. Artificially generated lidar signals from two-component inhomogeneous model atmospheres are inverted with the information from two wavelengths (532 and 1064 nm) simultaneously. The solution for the vertical aerosol extinction profile is formulated with Klett's far-end solution. The boundary value is expressed in terms of aerosol transmission along the lidar line according to Fernald's solution of the lidar equation. The aerosol transmission is determined iteratively with a transcendental equation on the assumption that a linear relationship exists between the extinction coefficients at both wavelengths. Inversion calculations are applied to model atmospheres with range-dependent lidar ratios representing the growth of aerosol particles caused by increasing relative humidity in the planetary boundary layer. For the inversion constant lidar ratios are assumed that vary between 40 and 70 sr. The numerical procedure turns out to be stable enough to provide meaningful results even in cases of misestimated lidar ratios. The application of the method is of less use for misestimated background radiation and low aerosol concentrations.     

Aerosol lidar intercomparison in the framework of the EARLINET project. 2. Aerosol backscatter algorithms

An intercomparison of aerosol backscatter lidar algorithms was performed in 2001 within the framework of the European Aerosol Research Lidar Network to Establish an Aerosol Climatology (EARLINET). The objective of this research was to test the correctness of the algorithms and the influence of the lidar ratio used by the various lidar teams involved in the EARLINET for calculation of backscatter-coefficient profiles from the lidar signals. The exercise consisted of processing synthetic lidar signals of various degrees of difficulty. One of these profiles contained height-dependent lidar ratios to test the vertical influence of those profiles on the various retrieval algorithms. Furthermore, a realistic incomplete overlap of laser beam and receiver field of view was introduced to remind the teams to take great care in the nearest range to the lidar. The intercomparison was performed in three stages with increasing knowledge on the input parameters. First, only the lidar signals were distributed; this is the most realistic stage. Afterward the lidar ratio profiles and the reference values at calibration height were provided. The unknown height-dependent lidar ratio had the largest influence on the retrieval, whereas the unknown reference value was of minor importance. These results show the necessity of making additional independent measurements, which can provide us with a suitable approximation of the lidar ratio. The final stage proves in general, that the data evaluation schemes of the different groups of lidar systems work well.     

Potential of lidar backscatter data to estimate solar aerosol radiative forcing

The potential to estimate solar aerosol radiative forcing (SARF) in cloudless conditions from backscatter data measured by widespread standard lidar has been investigated. For this purpose 132 days of sophisticated ground-based Raman lidar observations (profiles of particle extinction and backscatter coefficients at 532 nm wavelength) collected during two campaigns [the European Aerosol Research Lidar Network (EARLINET) and the Indian Ocean Experiment (INDOEX)] were analyzed. Particle extinction profiles were used as input for radiative transfer simulations with which to calculate the SARF, which then was plotted as a function of the column (i.e., height-integrated) particle backscatter coefficient (beta(c)). A close correlation between the SARF and beta(c) was found. SARF-beta(c) parameterizations in the form of polynomial fits were derived that exhibit an estimated uncertainty of +/-(10-30)%. These parameterizations can be utilized to analyze data of upcoming lidar satellite missions and for other purposes. The EARLINET-based parameterizations can be applied to lidar measurements at mostly continental, highly industrialized sites with limited maritime influence (Europe, North America), whereas the INDOEX parameterizations rather can be employed in polluted maritime locations, e.g., coastal regions of south and east Asia.     

Iterative method for the inversion of multiwavelength lidar signals to determine aerosol size distribution

Two iterative methods of inverting lidar backscatter signals to determine altitude profiles of aerosol extinction and altitude-resolved aerosol size distribution (ASD) are presented. The first method is for inverting two-wavelength lidar signals in which the shape of the ASD is assumed to be of power-law type, and the second method is for inverting multiwavelength lidar signals without assuming any a priori analytical form of ASD. An arbitrary value of the aerosol extinction-to-backscatter ratio (S(1)) is assumed initially to invert the lidar signals, and the ASD determined by use of the spectral dependence of the retrieved aerosol extinction coefficients is used to improve the value of S(1) iteratively. The methods are tested for different forms of altitude-dependent ASD's by use of simulated lidar-backscatter-signal profiles. The effect of random noise on the lidar backscatter signals is also studied.     

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.     

Aerosol and cloud backscatter at 1.06, 1.54, and 0.53 mum by airborne hard-target-calibrated Nd:YAG /methane Raman lidar

A lidar instrument was developed to make simultaneous measurements at three distinct wavelengths in the visible and near infrared at 0.532, 1.064, and 1.54 mum with high cross-sectional calibration accuracy. Aerosol and cloud backscatter cross sections were acquired during November and December 1989 and May and June 1990 by the NASA DC-8 aircraft as part of the Global Backscatter Experiment. The instrument, methodology, and measurement results are described. A Nd:YAG laser produced 1.064- and 0.532-mum energy. The 1.54-mum transmitted pulse was generated by Raman-shifted downconversion of the 1.064-mum pulse through a Raman cell pressured with methane gas. The lidar could be pointed in the nadir or zenith direction from the aircraft. A hard-target-based calibration procedure was used to obtain the ratio of the system calibration between the three wavelengths, and the absolute calibration was referenced to the 0.532-mum lidar molecular backscatter cross section for the clearest scattering regions. From the relative wavelength calibration, the aerosol backscatter cross sections at the longer wavelengths are resolved for values as small as 1% of the molecular cross section. Backscatter measurement accuracies are better than 10(-9) (m sr)(-1) at 1.064 and 1.54 mum. Results from the Pacific Ocean region of the multiwavelength backscatter dependence are presented. Results show extensive structure and variation for the aerosol cross sections. The range of observed aerosol cross section is over 4 orders of magnitude, from less than 10(-9) (m sr)(-1) to greater than 10(-5) (m sr)(-1).     

Vertical variability of aerosol backscatter from an airborne-focused continuous-wave CO2 lidar at 9.1-microm wavelength

Atmospheric aerosol backscatter measurements taken with a continuous-wave focused Doppler lidar at 9.1-mum wavelength were obtained over western North America and the Pacific Ocean from 13 to 26 September 1995 as part of a NASA airborne mission. Backscatter variability was measured for ~52 flight hours, covering an equivalent horizontal distance of ~30,000 km in the troposphere. Some quasi-vertical backscatter profiles were also obtained during various ascents and descents at altitudes that ranged from ~0.1 to 12 km. Similarities and differences for aerosol loading over land and ocean were observed. A midtropospheric aerosol backscatter background mode near 3 x 10(-11) to 1 x 10(-10) m(-1) sr(-1) was obtained, which is consistent with those of previous airborne and ground-based data sets.     

Synergetic technique combining elastic backscatter lidar data and sunphotometer AERONET inversion for retrieval by layer of aerosol optical and microphysical properties

We present a so-called lidar and almucantar (LidAlm) algorithm that combines information provided by standard elastic backscatter lidar (i.e., calibrated attenuated backscatter coefficient profile at one or two wavelengths) and sunphotometer AERONET inversion of almucantar like measurements (i.e., column-integrated aerosol size distribution and refractive index). The purpose of the LidAlm technique is to characterize the atmospheric column by its different aerosol layers. These layers may be distinct or partially mixed, and they may contain different aerosol species (e.g., urban, desert, or biomass burning aerosols). The LidAlm synergetic technique provides the extinction and backscatter coefficient profiles, particle size distributions, and backscatter-to-extinction ratios for each aerosol layer. We present the LidAlm procedure and sensitivity studies. The applications are illustrated with examples of actual atmospheric conditions encountered in the Paris area.     

Wavelength Dependence of Backscatter by use of Aerosol Microphysics and Lidar Data Sets: Application to 2.1- mum Wavelength for Space-Based and Airborne Lidars

An aerosol microphysics dataset was used to model backscatter in the 0.35-11-mum wavelength range, with the results validated by comparison with measured cw and pulsed lidar backscatter obtained during two NASA-sponsored airborne field experiments. Different atmospheric features were encountered, with aerosol backscatter ranging over 4 orders of magnitude. Modeled conversion functions were used to convert existing lidar backscatter datasets to 2.1 mum. Resulting statistical distribution shows the midtropospheric aerosol backscatter background mode of beta(2.1) to be between ~3.0 x 10(-10) and ~1.3 x 10(-9) m(-1) sr(-1), ~10-20 times higher than that for beta(9.1); and a beta(2.1) boundary layer mode of ~1.0 x 10(-7) to ~1.3 x 10(-6) m(-1) sr(-1), ~3-5 times higher than beta(9.1).     

Extinction-to-backscatter ratio of Asian dust observed with high-spectral-resolution lidar and Raman lidar

Extinction-to-backscatter ratio or lidar ratio is a key parameter in the issue of backscatter-lidar inversions. The lidar ratio of Asian dust was observed with a high-spectral-resolution lidar and a combined Raman elastic-backscatter lidar during the springs of 1998 and 1999. The measured values range from 42 to55 sr in most cases, with a mean of 51 sr. These values are significantly larger than those predicted by the Mie computations that incorporate measured Asian dust size distributions and a range of refractive index with a typical value of 1.55-0.005i. The enhancement of lidar ratio is mostly due to the nonsphericity of dust particles, as indicated by the T-matrix calculations for spheroid particles and a number of other theoretical studies. In addition, possible contamination of urban aerosols may also contribute somewhat in optically thin cases. Mie theory, although it can well describe spherical particle scattering, will not be sufficient to represent the scattering characteristics of irregular particles such as Asian dust, especially in directions larger than approximately 90 degrees when the size parameter is large.     

Imaginary refractive-index effects on desert-aerosol extinction versus backscatter relationships at 351 nm: numerical computations and comparison with Raman lidar measurements

A numerical model is used to investigate the dependence at 351 nm of desert-aerosol extinction and backscatter coefficients on particle imaginary refractive index (mi). Three ranges (-0.005 < or = mi < or = -0.001, -0.01 < or = mi < or = -0.001, and -0.02 < or = mi < or = -0.001) are considered, showing that backscatter coefficients are reduced as /mi/ increases, whereas extinction coefficients are weakly dependent on mi. Numerical results are compared with extinction and backscatter coefficients retrieved by elastic Raman lidar measurements performed during Saharan dust storms over the Mediterranean Sea. The comparison indicates that a range of -0.01 to -0.001 can be representative of Saharan dust aerosols and that the nonsphericity of mineral particles must be considered.     

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.     

Principal component analysis applied to multiwavelength lidar aerosol backscatter and extinction measurements

The use of powerful Raman backscatter lidars enables one to measure the stratospheric aerosol extinction profile independently of the backscatter, thereby obtaining additional information to aid in retrieving the physical characteristics of the sampled aerosol. We used principal component analysis to construct a self-consistent method for the retrieval of aerosol bulk physical and optical properties from multiwavelength elastic and/or inelastic Raman backscatter lidar signals. The procedure is applied to synthetic and actual lidar signals. We found that aerosol surface area and volume can be usefully estimated and that the use of Raman-derived extinction data leads to a notable improvement in the accuracy of the estimations.     

Tropospheric aerosol extinction coefficient profiles derived from scanning lidar measurements over Tsukuba, Japan, from 1990 to 1993

Mie scattering lidar was used to observe aerosol extinction coefficient profiles in the troposphere over Tsukuba (140 E, 36 N), Japan, for three years from March 1990 to February 1993, and data obtained in fair weather were analyzed. The lidar measurements were made by a vertical scanning mode to generate profiles of extinction coefficients from the lidar level to a 12-km altitude. The extinction coefficients were derived from the lidar signals using a two-component (air molecule and aerosol) lidar equation, in which the ratio of aerosol extinction to backscattering was assumed to be constant. Seasonal average profiles were derived from individual profiles. Three-year average profiles were also calculated and modeled using mathematical expressions. The model profile assumed (1) a constant extinction ratio in the atmospheric boundary layer (ABL), (2) an exponentially decreasing extinction ratio above the ABL, and (3) a constant extinction ratio in the upper troposphere where the extinction ratio can be defined as the ratio of the aerosol extinction coefficient to the air molecule extinction coefficient. The extinction ratios both in the ABL and in the upper troposphere and the scale height that was used to express the exponential decrease were used as three unknown parameters. Seasonal variation of optical thickness that was obtained by integrating extinction coefficients with height was also investigated.