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.

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.

Improved identification of the solution space of aerosol microphysical properties derived from the inversion of profiles of lidar optical data, part 2: simulations with synthetic optical data.

We developed a mathematical scheme that allows us to improve retrieval products obtained from the inversion of multiwavelength Raman/HSRL lidar data, commonly dubbed "3 backscatter+2 extinction" (3ß+2α) lidar. This scheme works independently of the automated inversion method that is currently being developed in the framework of the Aerosol-Cloud-Ecosystem (ACE) mission and which is successfully applied since 2012 [Atmos. Meas. Tech.7, 3487 (2014)10.5194/amt-7-3487-2014; "Comparison of aerosol optical and microphysical retrievals from HSRL-2 and in-situ measurements during DISCOVER-AQ 2013 (California and Texas)," in International Laser Radar Conference, July 2015, paper PS-C1-14] to data collected with the first airborne multiwavelength 3ß+2α high spectral resolution lidar (HSRL) developed at NASA Langley Research Center. The mathematical scheme uses gradient correlation relationships we presented in part 1 of our study [Appl. Opt.55, 9839 (2016)APOPAI0003-693510.1364/AO.55.009839] in which we investigated lidar data products and particle microphysical parameters from one and the same set of optical lidar profiles. For an accurate assessment of regression coefficients that are used in the correlation relationships we specially designed the proximate analysis method that allows us to search for a first-estimate solution space of particle microphysical parameters on the basis of a look-up table. The scheme works for any shape of particle size distribution. Simulation studies demonstrate a significant stabilization of the various solution spaces of the investigated aerosol microphysical data products if we apply this gradient correlation method in our traditional regularization technique. Surface-area concentration can be estimated with an uncertainty that is not worse than the measurement error of the underlying extinction coefficients. The retrieval uncertainty of the effective radius is as large as ±0.07 µm for fine mode particles and approximately 100% for particle size distributions composed of fine (submicron) and coarse (supermicron) mode particles. The volume concentration uncertainty is defined by the sum of the uncertainty of surface-area concentration and the uncertainty of the effective radius. The uncertainty of number concentration is better than 100% for any radius domain between 0.03 and 10 µm. For monomodal PSDs, the uncertainties of the real and imaginary parts of the CRI can be restricted to ±0.1 and ±0.01 on the domains [1.3; 1.8] and [0; 0.1], respectively.