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.     

Ozone differential absorption lidar algorithm intercomparison

An intercomparison of ozone differential absorption lidar algorithms was performed in 1996 within the framework of the Network for the Detection of Stratospheric Changes (NDSC) lidar working group. The objective of this research was mainly to test the differentiating techniques used by the various lidar teams involved in the NDSC for the calculation of the ozone number density from the lidar signals. The exercise consisted of processing synthetic lidar signals computed from simple Rayleigh scattering and three initial ozone profiles. Two of these profiles contained perturbations in the low and the high stratosphere to test the vertical resolution of the various algorithms. For the unperturbed profiles the results of the simulations show the correct behavior of the lidar processing methods in the low and the middle stratosphere with biases of less than 1% with respect to the initial profile to as high as 30 km in most cases. In the upper stratosphere, significant biases reaching 10% at 45 km for most of the algorithms are obtained. This bias is due to the decrease in the signal-to-noise ratio with altitude, which makes it necessary to increase the number of points of the derivative low-pass filter used for data processing. As a consequence the response of the various retrieval algorithms to perturbations in the ozone profile is much better in the lower stratosphere than in the higher range. These results show the necessity of limiting the vertical smoothing in the ozone lidar retrieval algorithm and questions the ability of current lidar systems to detect long-term ozone trends above 40 km. Otherwise the simulations show in general a correct estimation of the ozone profile random error and, as shown by the tests involving the perturbed ozone profiles, some inconsistency in the estimation of the vertical resolution among the lidar teams involved in this experiment.     

Vertical profiles of microphysical particle properties derived from inversion with two-dimensional regularization of multiwavelength Raman lidar data: experiment

Inversion with two-dimensional (2-D) regularization is a new methodology that can be used for the retrieval of profiles of microphysical properties, e.g., effective radius and complex refractive index of atmospheric particles from complete (or sections) of profiles of optical particle properties. The optical profiles are acquired with multiwavelength Raman lidar. Previous simulations with synthetic data have shown advantages in terms of retrieval accuracy compared to our so-called classical one-dimensional (1-D) regularization, which is a method mostly used in the European Aerosol Research Lidar Network (EARLINET). The 1-D regularization suffers from flaws such as retrieval accuracy, speed, and ability for error analysis. In this contribution, we test for the first time the performance of the new 2-D regularization algorithm on the basis of experimental data. We measured with lidar an aged biomass-burning plume over West/Central Europe. For comparison, we use particle in situ data taken in the smoke plume during research aircraft flights upwind of the lidar. We find good agreement for effective radius and volume, surface-area, and number concentrations. The retrieved complex refractive index on average is lower than what we find from the in situ observations. Accordingly, the single-scattering albedo that we obtain from the inversion is higher than what we obtain from the aircraft data. In view of the difficult measurement situation, i.e., the large spatial and temporal distances between aircraft and lidar measurements, this test of our new inversion methodology is satisfactory.     

Raman lidar monitoring of extinction and backscattering of African dust layers and dust characterization

Results on the monitoring of strong African dust outbreaks at Lecce in the southeastern corner of Italy (40 degrees 20' N, 18 degrees 6' E) during May 2001 are presented. This activity has been performed in the framework of the European Aerosol Research Lidar Network (EARLINET). The lidar station of Lecce is located on a flat rural area that is approximately 800 km from the northern Africa coast. So it is closer to Africa than most of all other EARLINET stations and allow monitoring African dust transport early in its life cycle, at all levels in the plume. An elastic-backscatter Raman lidar based on a XeF excimer laser (351 nm) has been used to monitor the time evolution and vertical structure of the dust layers and get independent measurements of the aerosol extinction and backscatter coefficients. The findings are presented in terms of vertical profiles of the extinction and backscatter coefficients and of the lidar ratio. A quite deep dust layer extending between 2 and 6 km and characterized by a backscatter coefficient of approximately 0.0016 (km sr)(-1), a lidar ratio of approximately 50 sr, and an aerosol optical depth of 0.26 was observed on 17 May 2001 between 18:55 and 20:07 UT. The layer persisted for approximately five days. Dust layers of lower optical thickness and shorter persistence time have generally been monitored at the lidar site during African dust outbreaks. Results on the chemical and morphological characterization of the dust collected at the lidar station are also given to further support the origin of the monitored aerosol layers.     

Aerosol observations by lidar in the nocturnal boundary layer

Aerosol observations by lidar in the nocturnal boundary layer (NBL) were performed in Potenza, Southern Italy, from 20 January to 20 February 1997. Measurements during nine winter nights were considered, covering a variety of boundary-layer conditions. The vertical profiles of the aerosol backscattering coefficient at 355 and 723.37 nm were determined through a Klett-modified iterative procedure, assuming the extinction-to-backscattering ratio within the NBL has a constant value. Aerosol average size characteristics were retrieved from almost simultaneous profiles of the aerosol backscattering coefficient at 355 and 723.37 nm, the measurements being consistent with an accumulation mode radius not exceeding 0.4 mum. Similar results in terms of aerosol sizes were obtained from measurements of the extinction-to-backscattering ratio profile at 355 nm performed on six nights during the measurement campaign. Backscattering profiles at 723.37 nm were also converted into profiles of aerosol liquid water content.     

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.     

Determination of slope in lidar data using a duplicate of the inverted function

An iterative method for determining slope in noisy lidar data is considered based on the use of a corrected ('shaped') inverted function and an assumed behavior of the unknown function of interest (an 'image function'). The method is utilized for extracting extinction- coefficient profiles from data of multiangle measurements. The sequence and specifics of the retrieval procedure, results of simulations, and essentials of the practical retrieval of particulate extinction-coefficient profiles from signals of the elastic scanning lidar are considered. The methodology may be applicable when extracting the extinction-coefficient profiles from an elastic lidar operating in a multiangle scanning mode, a combined Raman elastic-backscatter lidar, or a high spectral resolution lidar operating in a fixed angular position.     

Simulations of the observation of clouds and aerosols with the Experimental Lidar in Space Equipment system

We carried out a simulation study for the observation of clouds and aerosols with the Japanese Experimental Lidar in Space Equipment (ELISE), which is a two-wavelength backscatter lidar with three detection channels. The National Space Development Agency of Japan plans to launch the ELISE on the Mission Demonstrate Satellite 2 (MDS-2). In the simulations, the lidar return signals for the ELISE are calculated for an artificial, two-dimensional atmospheric model including different types of clouds and aerosols. The signal detection processes are simulated realistically by inclusion of various sources of noise. The lidar signals that are generated are then used as input for simulations of data analysis with inversion algorithms to investigate retrieval of the optical properties of clouds and aerosols. The results demonstrate that the ELISE can provide global data on the structures and optical properties of clouds and aerosols. We also conducted an analysis of the effects of cloud inhomogeneity on retrievals from averaged lidar profiles. We show that the effects are significant for space lidar observations of optically thick broken clouds.     

Algorithm improvement and validation of National Institute for Environmental Studies ozone differential absorption lidar at the Tsukuba Network for Detection of Stratospheric Change complementary station

Recently, a data processing and retrieval algorithm (version 2) for ozone, aerosol, and temperature lidar measurements was developed for an ozone lidar system at the National Institute for Environmental Studies (NIES) in Tsukuba (36 degrees N,140 degrees E), Japan. A method for obtaining the aerosol boundary altitude and the aerosol extinction-to-backscatter ratio in the version 2 algorithm enables a more accurate determination of the vertical profiles of aerosols and a more accurate correction of the systematic errors caused by aerosols in the vertical profile of ozone. Improvements in signal processing are incorporated for the correction of systematic errors such as the signal-induced noise and the dead-time effect. The mean vertical ozone profiles of the NIES ozone lidar were compared with those of the Stratospheric Aerosol and Gas Experiment II (SAGE II); they agreed well within a 5% relative difference in the 20-40 km altitude range and within 10% up to 45 km. The long-term variations in the NIES ozone lidar also showed good coincidence with the ozonesonde and SAGE II at 20, 25, 30, and 35 km. The temperatures retrieved from the NIES ozone lidar and those given by the National Center for Environmental Prediction agreed within 7 K in the 35-50 km range.     

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.     

Laser Sensing of a Subsurface Oceanic Layer. II. Polarization Characteristics of Signals

In Part I of this paper we calculated depth profiles and polarization characteristics of airborne lidar return signals by the Monte Carlo method. Here we calculate the polarization characteristics of lidar return signals for different types of water. We demonstrate the feasibility of polarization lidar application to the detection of underwater inhomogeneities of different origins. It is shown that simultaneous analysis of depth profiles of the lidar return signal power and signal depolarization ratio substantially increases the information content of airborne lidar sensing of seawater. We compare calculated results with the data of airborne lidar measurements for lambda = 0.53 mum.     

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.     

Retrieval of aerosol extinction coefficient profiles from Raman lidar data by inversion method

We regard the problem of differentiation occurring in the retrieval of aerosol extinction coefficient profiles from inelastic Raman lidar signals by searching for a stable solution of the resulting Volterra integral equation. An algorithm based on a projection method and iterative regularization together with the L-curve method has been performed on synthetic and measured lidar signals. A strategy to choose a suitable range for the integration within the framework of the retrieval of optical properties is proposed here for the first time to our knowledge. The Monte Carlo procedure has been adapted to treat the uncertainty in the retrieval of extinction coefficients.     

Stratospheric Aerosol and Gas Experiment III cloud data product

The latest in a series of solar occultation satellite instruments, Stratospheric Aerosol and Gas Experiment (SAGE) III, was placed into orbit in December 2001, and data were obtained until March 2006. Measurements were made of the extinction attributable to aerosols and cloud at a number of wavelengths between 290 and 1550 nm. The analysis of data obtained by its predecessor, SAGE II, has shown that an intercomparison of such data at two or more wavelengths may be used to separate the effects of cloud and aerosol. This analysis has been done on a routine basis for many years using SAGE II data at 525 and 1020 nm and applied extensively to global studies of tropospheric cloud and aerosol. Here we describe the aerosol-cloud separation algorithm developed for use with the SAGE III data, which uses the extinction at 525, 1020, and 1550 nm. This algorithm is now being used to produce vertical profiles of cloud presence as a standard SAGE III data product. These profiles have a vertical resolution of 0.5 km and cover the altitude range from 6.0 to 30.0 km, and data are presently available from March 2002 onward. An outline is given of the development of this algorithm, the nature of the SAGE III data, and the algorithm performance. To maintain continuity with SAGE II cloud data, the relative performances of the SAGE II and SAGE III algorithms are also examined. An example of the application of the algorithm to SAGE III tropospheric data is shown and discussed.     

III-posed retrieval of aerosol extinction coefficient profiles from Raman lidar data by regularization

In the analysis of Raman lidar measurements of aerosol extinction, it is necessary to calculate the derivative of the logarithm of the ratio between the atmospheric number density and the range-corrected lidar-received power. The statistical fluctuations of the Raman signal can produce large fluctuations in the derivative and thus in the aerosol extinction profile. To overcome this difficult situation we discuss three methods: Tikhonov regularization, variational, and the sliding best-fit (SBF). Three methods are performed on the profiles taken from the European Aerosol Research Lidar Network lidar database simulated at the Raman shifted wavelengths of 387 and 607 nm associated with the emitted signals at 355 and 532 nm. Our results show that the SBF method does not deliver good results for low fluctuation in the profile. However, Tikhonov regularization and the variational method yield very good aerosol extinction coefficient profiles for our examples. With regard to, e.g., the 532 nm wavelength, the L2 errors of the aerosol extinction coefficient profile by using the SBF, Tikhonov, and variational methods with respect to synthetic noisy data are 0.0015(0.0024), 0.00049(0.00086), and 0.00048(0.00082), respectively. Moreover, the L2 errors by using the Tikhonov and variational methods with respect to a more realistic noisy profile are 0.0014(0.0016) and 0.0012(0.0016), respectively. In both cases the L2 error given in parentheses concerns the second example.     

Upper tropospheric temperature measurements with the use of a Raman lidar

Upper tropospheric temperature profiles were measured with the NASA Goddard Space Flight Center scanning Raman lidar five months after the eruption of Mt. Pinatubo. To derive temperatures in regions of high aerosol content, the aerosol transmission is calculated for the Raman N(2) return signals under cloud-free conditions. The lidar-derived aerosol backscattering ratio and an estimate of the aerosol extinction-to-backscatter ratio were used to compute the aerosol transmission. With a model reference temperature at 25 km, temperature profiles with a root-mean-square difference between the lidar and radiosonde temperatures of <2 K were obtained over an altitude range of 5-10 km for a 10-min integrated measurement with 300-m resolution.     

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.     

High-pulse-repetition-freqyency lidar system using a single telescope for transmission and reception

The design, construction, and operation of a stratospheric Rayleigh lidar system is outlined. The lidar system was designed to operate as a Doppler lidar; however, for the first stage of the project it was set up to operate in a manner similar to a more conventional stratospheric Rayleigh lidar. This system includes a number of unique design features, including a high-pulse-repetition-frequency laser and the use of a single 1-m-diameter telescope for transmission of the laser pulse and reception of the backscattered light. An associated high-speed rotating shutter system switches the optical system from the transmission to the reception mode. The system was operated at Adelaide, Australia (35° S, 138° E). Scattering ratio and temperature profiles are calculated for data collected during the period from 10 March 1992 to 11 May 1993. The scattering ratio profiles clearly show the reduction in the scattering from the stratospheric aerosol layer. This is due to the removal of the aerosol injected by the eruption of Mt. Pinatubo. The measured relative density profiles show very good agreement with the Cospar international reference atmosphere model densities, as do the temperature profiles calculated from these.     

Laser sensing of a subsurface oceanic layer. I. Effect of the atmosphere and wind-driven sea waves

Depth profiles and polarization characteristics of airborne lidar return signals have been calculated by the Monte Carlo method. We analyze some peculiarities of depth profiles of lidar return signals for a rough air-water interface. The distorting effect of the atmosphere on the lidar return signal structure is evaluated as a function of the geometry of the observations. Calculated results are compared with the data of airborne lidar measurements for lambda = 0.53 mum.