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

Method for retrieving range-resolved aerosol microphysical properties from polarization lidar measurements

Aerosol microphysical properties, such as volume concentration (VC) and effective radius (ER), are of great importance to evaluate their radiative forcing and impacts on climate change. However, range-resolved aerosol VC and ER still cannot be obtained by remote sensing currently except for the column-integrated one from sun-photometer observation. In this study, a retrieval method of range-resolved aerosol VC and ER is firstly proposed based on the partial least squares regression (PLSR) and deep neural networks (DNN), combining polarization lidar and collocated AERONET (AErosol RObotic NETwork) sun-photometer observations. The results show that the measurement of widely-used polarization lidar can be reasonably used to derive the aerosol VC and ER, with the determination coefficient (R²) of 0.89 (0.77) for VC (ER) by use of the DNN method. Moreover, it is proven that the lidar-based height-resolved VC and ER at near-surface are well consistent with independent observations of collocated Aerodynamic Particle Sizer (APS). Additionally, we found that there are significant diurnal and seasonal variations of aerosol VC and ER in the atmosphere at Semi-Arid Climate and Environment Observatory of Lanzhou University (SACOL). Compared with columnar ones from the sun-photometer observation, this study provides a reliable and practical way to obtain full-day range-resolved aerosol VC and ER from widely-used polarization lidar observation, even under cloud conditions. Moreover, this study also can be applied to long-term observations by current ground-based lidar networks and spaceborne CALIPSO lidar, aiming to further evaluate aerosol climatic effects more accurately.

Correct(ed) Klett-Fernald algorithm for elastic aerosol backscatter retrievals: a sensitivity analysis

Retrieval algorithms for aerosol elastic backscatter lidars are most commonly based on a signal inversion known as the Klett solution. While often used, a sign error in the original publication has barely been mentioned or recognized in the scientific community. In this study, we present a corrected Klett inversion and a sensitivity analysis of its implementation under different atmospheric conditions. We show that the error's dimension depends on multiple factors, thus preventing trivial a posteriori corrections on the results calculated with the original, uncorrected Klett solution. Comparing the uncorrected with the corrected Klett solution, long integration pathways in combination with low aerosol concentrations lead to substantial relative deviations of more than 100%, whereas short integration pathways with high aerosol concentrations considerably reduce the relative deviations to magnitudes around 15%. The higher the altitude of layers with increased aerosol concentrations, the higher the deviations due to the incorrect usage of Klett's inversion, however on a slight level (≈1%).

Atmospheric aerosol particle size distribution from Lidar data based on the lognormal distribution mode

Analysis of the atmospheric Aerosol Particle Size Distribution (APSD) retrieved from Light Detection and Ranging (Lidar) data is one of the popular fields in atmospheric remote sensing. An APSD retrieval method based on the lognormal particle size distribution, the Mie theory, and the Second Simulation of the Satellite Signal in the Solar Spectrum (6S) were studied. According to the 6S, this method divides the main body of the aerosols into four basic components and calculates the APSD from Lidar data and the optical and microphysical characteristics of these components. Numerical simulation and experimental observations reveal that this method can obtain the APSD for particle sizes of >0.15 μm in the different vertical layers with good reliability.

Development of an Automatic Polarization Raman LiDAR for Aerosol Monitoring over Complex Terrain

High temporal and spatial resolution profiling of aerosol properties is required to study air pollution sources, aerosol transport, and features of atmospheric structures over complex terrain. A polarization Raman LiDAR with remote operation capability was developed for this purpose and deployed in the Vipava Valley, Slovenia, a location in the Alpine region where high concentrations of aerosols originating from a number of different local and remote sources were found. The system employs two high-power Nd:YAG pulsed lasers at 355 nm and 1064 nm as transmitters and provides the capability to extract the extinction coefficient, backscatter coefficients, depolarization ratio, Ångström exponent, and LiDAR ratio profiles. Automatized remote operation in an indoor environment provides a high duty cycle in all weather conditions. In addition to the detailed description of the device, an assessment of its potential and the retrieval uncertainties of the measured quantities is discussed. System optimization and performance studies include calibration of the depolarization ratio, merging of near-range (analog) and far-range (photon counting) data, determination of overlap functions, and validation of the retrieved observables with radiosonde data. Two cases for assessing LiDAR performance under specific weather conditions (during rain and in the presence of mineral dust) are also presented.

Automated, unsupervised inversion of multiwavelength lidar data with TiARA: assessment of retrieval performance of microphysical parameters using simulated data

We evaluate the retrieval performance of the automated, unsupervised inversion algorithm, Tikhonov Advanced Regularization Algorithm (TiARA), which is used for the autonomous retrieval of microphysical parameters of anthropogenic and natural pollution particles. TiARA (version 1.0) has been developed in the past 10 years and builds on the legacy of a data-operator-controlled inversion algorithm used since 1998 for the analysis of data from multiwavelength Raman lidar. The development of TiARA has been driven by the need to analyze in (near) real time large volumes of data collected with NASA Langley Research Center's high-spectral-resolution lidar (HSRL-2). HSRL-2 was envisioned as part of the NASA Aerosols-Clouds-Ecosystems mission in response to the National Academy of Sciences (NAS) Decadal Study mission recommendations 2007. TiARA could thus also serve as an inversion algorithm in the context of a future space-borne lidar. We summarize key properties of TiARA on the basis of simulations with monomodal logarithmic-normal particle size distributions that cover particle radii from approximately 0.05 μm to 10 μm. The real and imaginary parts of the complex refractive index cover the range from non-absorbing to highly light-absorbing pollutants. Our simulations include up to 25% measurement uncertainty. The goal of our study is to provide guidance with respect to technical features of future space-borne lidars, if such lidars will be used for retrievals of microphysical data products, absorption coefficients, and single-scattering albedo. We investigate the impact of two different measurement-error models on the quality of the data products. We also obtain for the first time, to the best of our knowledge, a statistical view on systematic and statistical uncertainties, if a large volume of data is processed. Effective radius is retrieved to 50% accuracy for 58% of cases with an imaginary part up to 0.01i and up to 100% of cases with an imaginary part of 0.05i. Similarly, volume concentration, surface-area concentration, and number concentrations are retrieved to 50% accuracy in 56%-100% of cases, 99%-100% of cases, and 54%-87% of cases, respectively, depending on the imaginary part. The numbers represent measurement uncertainties of up to 15%. If we target 20% retrieval accuracy, the numbers of cases that fall within that threshold are 36%-76% for effective radius, 36%-73% for volume concentration, 98%-100% for surface-area concentration, and 37%-61% for number concentration. That range of numbers again represents a spread in results for different values of the imaginary part. At present, we obtain an accuracy of (on average) 0.1 for the real part. A case study from the ORCALES field campaign is used to illustrate data products obtained with TiARA.

Current Research in Lidar Technology Used for the Remote Sensing of Atmospheric Aerosols

Lidars are active optical remote sensing instruments with unique capabilities for atmospheric sounding. A manifold of atmospheric variables can be profiled using different types of lidar: concentration of species, wind speed, temperature, etc. Among them, measurement of the properties of aerosol particles, whose influence in many atmospheric processes is important but is still poorly stated, stands as one of the main fields of application of current lidar systems. This paper presents a review on fundamentals, technology, methodologies and state-of-the art of the lidar systems used to obtain aerosol information. Retrieval of structural (aerosol layers profiling), optical (backscatter and extinction coefficients) and microphysical (size, shape and type) properties requires however different levels of instrumental complexity; this general outlook is structured following a classification that attends these criteria. Thus, elastic systems (detection only of emitted frequencies), Raman systems (detection also of Raman frequency-shifted spectral lines), high spectral resolution lidars, systems with depolarization measurement capabilities and multi-wavelength instruments are described, and the fundamentals in which the retrieval of aerosol parameters is based is in each case detailed.

Vertical profiles of pure dust and mixed smoke-dust plumes inferred from inversion of multiwavelength Raman/polarization lidar data and comparison to AERONET retrievals and in situ observations

We present for the first time vertical profiles of microphysical properties of pure mineral dust (largely unaffected by any other aerosol types) on the basis of the inversion of optical data collected with multiwavelength polarization Raman lidar. The data were taken during the Saharan Mineral Dust Experiment (SAMUM) in Morocco in 2006. We also investigated two cases of mixed dust-smoke plumes on the basis of data collected during the second SAMUM field campaign that took place in the Republic of Cape Verde in 2008. Following the experience of the Aerosol Robotic Network (AERONET), the dust is modeled as a mixture of spherical particles and randomly oriented spheroids. The retrieval is performed from the full set of lidar input data (three backscatter coefficients, two extinction coefficients, and one depolarization ratio) and from a reduced set of data in which we exclude the depolarization ratio. We find differences of the microphysical properties depending on what kind of optical data combination we use. For the case of pure mineral dust, the results from these two sets of optical data are consistent and confirm the validity of the spheroid particle model for data inversion. Our results indicate that in the case of pure mineral dust we do not need depolarization information in the inversion. For the mixture of dust and biomass burning, there seem to be more limitations in the retrieval accuracy of the various data products. The evaluation of the quality of our data products is done by comparing our lidar-derived data products (vertically resolved) to results from AERONET Sun photometer observations (column-averaged) carried out at the lidar field site. Our results for dust effective radius show agreement with the AERONET observations within the retrieval uncertainties. Regarding the complex refractive index a comparison is difficult, as AERONET provides this parameter as wavelength-dependent quantity. In contrast, our inversion algorithm provides this parameter as a wavelength-independent quantity. We also show some comparison to results from airborne in situobservation. A detailed comparison to in situ results will be left for a future contribution.

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.

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.

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.