Arrange and average algorithm for the retrieval of aerosol parameters from multiwavelength high-spectral-resolution lidar/Raman lidar data

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.

Improved identification of the solution space of aerosol microphysical properties derived from the inversion of profiles of lidar optical data, part 3: case studies

We conclude our series of publications on the development of the gradient correlation method (GCM), which can be used for an improved stabilization of the solution space of particle microphysical parameters derived from measurements with multiwavelength Raman and high-spectral-resolution lidar (3 backscatter +2 extinction coefficients). We show results of three cases studies. The data were taken with a ground-based multiwavelength Raman lidar during the Saharan Mineral Dust Experiment in the Cape Verde Islands (North Atlantic). These cases describe mixtures of dust with smoke. For our data analysis we separated the contribution of smoke to the total signal and only used these optical profiles for the test of GCM. The results show a significant stabilization of the solution space of the particle microphysical parameter retrieval on the particle radius domain from 0.03 to 10 μm, the real part of the complex refractive index domain from 1.3 to 1.8, and the imaginary part from 0 to 0.1. This new method will be included in the Tikhonov Advanced Regularization Algorithm, which is a fully automated, unsupervised algorithm that is used for the analysis of data collected with the worldwide first airborne 3 backscatter +2 extinction high-spectral-resolution lidar developed by NASA Langley Research Center.

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.

Improved identification of the solution space of aerosol microphysical properties derived from the inversion of profiles of lidar optical data, part 1: theory

Multiwavelength Raman/high spectral resolution lidars that measure backscatter coefficients at 355, 532, and 1064 nm and extinction coefficients at 355 and 532 nm can be used for the retrieval of particle microphysical parameters, such as effective and mean radius, number, surface-area and volume concentrations, and complex refractive index, from inversion algorithms. In this study, we carry out a correlation analysis in order to investigate the degree of dependence that may exist between the optical data taken with lidar and the underlying microphysical parameters. We also investigate if the correlation properties identified in our study can be used as a priori or a posteriori constraints for our inversion scheme so that the inversion results can be improved. We made the simplifying assumption of error-free optical data in order to find out what correlations exist in the best case situation. Clearly, for practical applications, erroneous data need to be considered too. On the basis of simulations with synthetic optical data, we find the following results, which hold true for arbitrary particle size distributions, i.e., regardless of the modality or the shape of the size distribution function: surface-area concentrations and extinction coefficients are linearly correlated with a correlation coefficient above 0.99. We also find a correlation coefficient above 0.99 for the extinction coefficient versus (1) the ratio of the volume concentration to effective radius and (2) the product of the number concentration times the sum of the squares of the mean radius and standard deviation of the investigated particle size distributions. Besides that, we find that for particles of any mode fraction of the particle size distribution, the complex refractive index is uniquely defined by extinction- and backscatter-related Ångström exponents, lidar ratios at two wavelengths, and an effective radius.

Retrieval of the particle size distribution function from the data of lidar sensing under the assumption of known refractive index

This paper presents a method to retrieve the particle size distribution function from the data of vertical lidar sensing. We have used 462 data models obtained at the Zvenigorod AERONET site obtained in 2011-2012. For each laser shot, we considered both fine (with particle sizes in a range from 0.05 to 0.6 μm) and coarse aerosol fractions (from 0.6 to 10 μm), with emphasize on the coarse fraction. Our suggested method is a modification of the Tikhonov method. The Tikhonov method is not optimal for coarse particles because its stabilizer does not and cannot account for the presence of the coarse mode, i.e., existence of more than one maximum of the size distribution function. The components of the matrix Wu-1 located in quadrants II and IV are sensitive to the change of these parameters. Neglecting this fact will lead again to arbitrary estimates of the contribution of the coarse particles even for exact values on the main diagonal and the two diagonals adjacent to it. Our method allows the coarse fraction up to 2.5 μm to be determined unambiguously. For larger particles (>2.5  μm) we recommend using the available sets of the coefficients, but with the level of values to be determined.

Retrieval of aerosol parameters from multiwavelength lidar: investigation of the underlying inverse mathematical problem

We present an investigation of some important mathematical and numerical features related to the retrieval of microphysical parameters [complex refractive index, single-scattering albedo, effective radius, total number, surface area, and volume concentrations] of ambient aerosol particles using multiwavelength Raman or high-spectral-resolution lidar. Using simple examples, we prove the non-uniqueness of an inverse solution to be the major source of the retrieval difficulties. Some theoretically possible ways of partially compensating for these difficulties are offered. For instance, an increase in the variety of input data via combination of lidar and certain passive remote sensing instruments will be helpful to reduce the error of estimation of the complex refractive index. We also demonstrate a significant interference between Aitken and accumulation aerosol modes in our inversion algorithm, and confirm that the solutions can be better constrained by limiting the particle radii. Applying a combination of an analytical approach and numerical simulations, we explain the statistical behavior of the microphysical size parameters. We reveal and clarify why the total surface area concentration is consistent even in the presence of non-unique solution sets and is on average the most stable parameter to be estimated, as long as at least one extinction optical coefficient is employed. We find that for selected particle size distributions, the total surface area and volume concentrations can be quickly retrieved with fair precision using only single extinction coefficients in a simple arithmetical relationship.