3+2 + X: what is the most useful depolarization input for retrieving microphysical properties of non-spherical particles from lidar measurements using the spheroid model of Dubovik et al. (2006)?

The typical multiwavelength aerosol lidar data set for inversion of optical to microphysical parameters is composed of three backscatter coefficients (β) at 355, 532, and 1064 nm and two extinction coefficients (α) at 355 and 532 nm. This data combination is referred to as a 3β+2α or 3+2 data set. This set of data is sufficient for retrieving some important microphysical particle parameters if the particles have spherical shape. Here, we investigate the effect of including the particle linear depolarization ratio (δ) as a third input parameter for the inversion of lidar data. The inversion algorithm is generally not used if measurements show values of δ that exceed 0.10 at 532 nm, i.e. in the presence of non-spherical particles such as desert dust, volcanic ash, and, under special circumstances, biomass-burning smoke. We use experimental data collected with instruments that are capable of measuring δ at all three lidar wavelengths with an inversion routine that applies the spheroidal light-scattering model of Dubovik et al. (2006) with a fixed axis-ratio distribution to replicate scattering properties of non-spherical particles. The inversion gives the fraction of spheroids required to replicate the optical data as an additional output parameter. This is the first systematic test of the effect of using all theoretically possible combinations of δ taken at 355, 532, and 1064 nm as input in the lidar data inversion. We find that depolarization information of at least one wavelength already provides useful information for the inversion of optical data that have been collected in the presence of non-spherical mineral dust particles. However, any choice of δλ will give lower values of the single-scattering albedo than the traditional 3+2 data set. We find that input data sets that include δ355 give a spheroid fraction that closely resembles the dust ratio we obtain from using β532 and δ532 in a methodology applied in aerosol-type separation. The use of δ355 in data sets of two or three δλ reduces the spheroid fraction that is retrieved when using δ532 and δ1064. Use of the latter two parameters without accounting for δ355 generally leads to high spheroid fractions that we consider not trustworthy. The use of three δλ instead of two δλ, including the constraint that one of these is measured at 355 nm does not provide any advantage over using 3+2+δ355 for the observations with varying contributions of mineral dust considered here. However, additional measurements at wavelengths different from 355 nm would be desirable for application to a wider range of aerosol scenarios that may include non-spherical smoke particles, which can have values of δ355 that are indistinguishable from those found for mineral dust. We therefore conclude that – depending on measurement capability – the future standard input for inversion of lidar data taken in the presence of mineral dust particles and using the spheroid model of Dubovik et al. (2006) might be 3+2+δ355 or 3+2+δ355+δ532.Peer reviewe

Polarization lidar: An extended three-signal calibration approach

We present a new formalism to calibrate a threesignal polarization lidar and to measure highly accurate height profiles of the volume linear depolarization ratios under realistic experimental conditions. The methodology considers elliptically polarized laser light, angular misalignment of the receiver unit with respect to the main polarization plane of the laser pulses, and cross talk among the receiver channels. A case study of a liquid-water cloud observation demonstrates the potential of the new technique. Long-term observations of the calibration parameters corroborate the robustness of the method and the long-term stability of the three-signal polarization lidar. A comparison with a second polarization lidar shows excellent agreement regarding the derived volume linear polarization ratios in different scenarios: A biomass burning smoke event throughout the troposphere and the lower stratosphere up to 16 km in height, a dust case, and also a cirrus cloud case. © Author(s) 2019

1064 nm rotational Raman lidar for particle extinction and lidar-ratio profiling: cirrus case study

For the first time, vertical profiles of the 1064 nm particle extinction coefficient obtained from Raman lidar observations at 1058 nm (nitrogen and oxygen rotational Raman backscatter) are presented. We applied the new technique in the framework of test measurements and performed several cirrus observations of particle backscatter and extinction coefficients, and corresponding extinction-to-backscatter ratios at the wavelengths of 355, 532, and 1064 nm. The cirrus backscatter coefficients were found to be equal for all three wavelengths keeping the retrieval uncertainties in mind. The multiple-scattering-corrected cirrus extinction coefficients at 355 nm were on average about 20–30 % lower than the ones for 532 and 1064 nm. The cirrus-mean extinction-to-backscatter ratio (lidar ratio) was 31 ± 5 sr (355 nm), 36 ± 5 sr (532 nm), and 38 ± 5 sr (1064 nm) in this single study. We further discussed the requirements needed to obtain aerosol extinction profiles in the lower troposphere at 1064 nm with good accuracy (20 % relative uncertainty) and appropriate temporal and vertical resolution

3+2 + X : what is the most useful depolarization input for retrieving microphysical properties of non-spherical particles from lidar measurements using the spheroid model of Dubovik et al. (2006)?

The typical multiwavelength aerosol lidar data set for inversion of optical to microphysical parameters is composed of three backscatter coefficients (β) at 355, 532, and 1064 nm and two extinction coefficients (α) at 355 and 532 nm. This data combination is referred to as a 3β C 2α or 3 + 2 data set. This set of data is sufficient for retrieving some important microphysical particle parameters if the particles have spherical shape. Here, we investigate the effect of including the particle linear depolarization ratio (δ) as a third input parameter for the inversion of lidar data. The inversion algorithm is generally not used if measurements show values of d that exceed 0.10 at 532 nm, i.e. in the presence of nonspherical particles such as desert dust, volcanic ash, and, under special circumstances, biomass-burning smoke. We use experimental data collected with instruments that are capable of measuring d at all three lidar wavelengths with an inversion routine that applies the spheroidal light-scattering model of Dubovik et al. (2006) with a fixed axis-ratio distribution to replicate scattering properties of non-spherical particles. The inversion gives the fraction of spheroids required to replicate the optical data as an additional output parameter. This is the first systematic test of the effect of using all theoretically possible combinations of d taken at 355, 532, and 1064 nm as input in the lidar data inversion. We find that depolarization information of at least one wavelength already provides useful information for the inversion of optical data that have been collected in the presence of non-spherical mineral dust particles. However, any choice of d will give lower values of the single-scattering albedo than the traditional 3 + 2 data set. We find that input data sets that include d355 give a spheroid fraction that closely resembles the dust ratio we obtain from using β532 and d532 in a methodology applied in aerosol-type separation. The use of d355 in data sets of two or three d? reduces the spheroid fraction that is retrieved when using d532 and d1064. Use of the latter two parameters without accounting for d355 generally leads to high spheroid fractions that we consider not trustworthy. The use of three d instead of two δ, including the constraint that one of these is measured at 355 nm does not provide any advantage over using 3 + 2 + d355 for the observations with varying contributions of mineral dust considered here. However, additional measurements at wavelengths different from 355 nm would be desirable for application to a wider range of aerosol scenarios that may include non-spherical smoke particles, which can have values of d355 that are indistinguishable from those found for mineral dust. We therefore conclude that - depending on measurement capability - the future standard input for inversion of lidar data taken in the presence of mineral dust particles and using the spheroid model of Dubovik et al. (2006) might be 3+2Cδ355 or 3 + 2 + δ355 + δ532. © 2019 The Author(s)

Depolarization and lidar ratios at 355, 532, and 1064 nm and microphysical properties of aged tropospheric and stratospheric Canadian wildfire smoke

We present spectrally resolved optical and microphysical properties of western Canadian wildfire smoke observed in a tropospheric layer from 5-6.5 km height and in a stratospheric layer from 15-16 km height during a recordbreaking smoke event on 22 August 2017. Three polarization/ Raman lidars were run at the European Aerosol Research Lidar Network (EARLINET) station of Leipzig, Germany, after sunset on 22 August. For the first time, the linear depolarization ratio and extinction-to-backscatter ratio (lidar ratio) of aged smoke particles were measured at all three important lidar wavelengths of 355, 532, and 1064 nm. Very different particle depolarization ratios were found in the troposphere and in the stratosphere. The obviously compact and spherical tropospheric smoke particles caused almost no depolarization of backscattered laser radiation at all three wavelengths ( 500nm). The stratospheric smoke particles formed a pronounced accumulation mode (in terms of particle volume or mass) centered at a particle radius of 350-400 nm. The effective particle radius was 0.32 μm. The tropospheric smoke particles were much smaller (effective radius of 0.17 μm). Mass concentrations were of the order of 5.5 μgm-3 (tropospheric layer) and 40 μgm-3 (stratospheric layer) in the night of 22 August 2017. The single scattering albedo of the stratospheric particles was estimated to be 0.74, 0.8, and 0.83 at 355, 532, and 1064 nm, respectively

Investigations of boundary layer structure, cloud characteristics and vertical mixing of aerosols at Barbados with large eddy simulations

Large eddy simulations (LESs) are performed for the area of the Caribbean island Barbados to investigate island effects on boundary layer modification, cloud generation and vertical mixing of aerosols. Due to the presence of a topographically structured island surface in the domain center, the model setup has to be designed with open lateral boundaries. In order to generate inflow turbulence consistent with the upstream marine boundary layer forcing, we use the cell perturbation method based on finite amplitude potential temperature perturbations. In this work, this method is for the first time tested and validated for moist boundary layer simulations with open lateral boundary conditions. Observational data obtained from the SALTRACE field campaign is used for both model initialization and a comparison with Doppler wind and Raman lidar data. Several numerical sensitivity tests are carried out to demonstrate the problems related to “gray zone modeling” when using coarser spatial grid spacings beyond the inertial subrange of three-dimensional turbulence or when the turbulent marine boundary layer flow is replaced by laminar winds. Especially cloud properties in the downwind area west of Barbados are markedly affected in these kinds of simulations. Results of an additional simulation with a strong trade-wind inversion reveal its effect on cloud layer depth and location. Saharan dust layers that reach Barbados via long-range transport over the North Atlantic are included as passive tracers in the model. Effects of layer thinning, subsidence and turbulent downward transport near the layer bottom at z ≈ 1800 m become apparent. The exact position of these layers and strength of downward mixing is found to be mainly controlled atmospheric stability (especially inversion strength) and wind shear. Comparisons of LES model output with wind lidar data show similarities in the downwind vertical wind structure. Additionally, the model results accurately reproduce the development of the daytime convective boundary layer measured by the Raman lidar

Cloud top heights and aerosol columnar properties from combined EarthCARE lidar and imager observations: the AM-CTH and AM-ACD products

The Earth Cloud, Aerosol and Radiation Explorer (EarthCARE) is a combination of multiple active and passive instruments on a single platform. The Atmospheric Lidar (ATLID) provides vertical information of clouds and aerosol particles along the satellite track. In addition, the Multi-Spectral Imager (MSI) collects multi-spectral information from the visible to the infrared wavelengths over a swath width of 150 km across the track. The ATLID–MSI Column Products processor (AM-COL) described in this paper combines the high vertical resolution of the lidar along track and the horizontal resolution of the imager across track to better characterize a three-dimensional scene. ATLID Level 2a (L2a) data from the ATLID Layer Products processor (A-LAY), MSI L2a data from the MSI Cloud Products processor (M-CLD) and the MSI Aerosol Optical Thickness processor (M-AOT), and MSI Level 1c (L1c) data are used as input to produce the synergistic columnar products: the ATLID–MSI Cloud Top Height (AM-CTH) and the ATLID–MSI Aerosol Column Descriptor (AM-ACD). The coupling of ATLID (measuring at 355 nm) and MSI (at ≥670 nm) provides multi-spectral observations of the aerosol properties. In particular, the Ångström exponent from the spectral aerosol optical thickness (AOT 355/670 nm) adds valuable information for aerosol typing. The AOT across track, the Ångström exponent and the dominant aerosol type are stored in the AM-ACD product. The accurate detection of the cloud top height (CTH) with lidar is limited to the ATLID track. The difference in the CTH detected by ATLID and retrieved by MSI is calculated along track. The similarity of MSI pixels across track with those along track is used to transfer the calculated CTH difference to the entire MSI swath. In this way, the accuracy of the CTH is increased to achieve the EarthCARE mission's goal of deriving the radiative flux at the top of the atmosphere with an accuracy of 10 W m−2 for a 100 km2 snapshot view of the atmosphere. The synergistic CTH difference is stored in the AM-CTH product. The quality status is provided with the products. It depends, e.g., on day/night conditions and the presence of multiple cloud layers. The algorithm was successfully tested using the common EarthCARE test scenes. Two definitions of the CTH from the model truth cloud extinction fields are compared: an extinction-based threshold of 20 Mm−1 provides the geometric CTH, and a cloud optical thickness threshold of 0.25 describes the radiative CTH. The first CTH definition was detected with ATLID and the second one with MSI. The geometric CTH is always higher than or equal to the radiative CTH

Measurement of the linear depolarization ratio of aged dust of three wavelengths (355, 532 and 1064 nm) simultaneously over Barbados

A ground-based polarization Raman lidar is presented, that is able to measure the depolarization ratio at three wavelengths (355, 532 and 1064 nm) simultaneously. This new feature is implemented for the first time in a Raman lidar. It provides a full dataset of 3 backscatter coefficients, two extinction coefficients and 3 depolarization ratios (3+2+3 lidar system). To ensure the data quality, it has been compared to the well characterized two wavelength polarization lidar POLIS. Measurements of long-range transported dust have been performed in the framework of the Saharan Aerosol Long-Range Transport and Aerosol-Cloud-Interaction Experiment (SALTRACE) in the Caribbean

Estimation of dust related ice nucleating particles in the atmosphere: Comparison of profiling and in-situ measurements

Vertical profiles of number concentrations of dust particles relevant for ice nucleation in clouds are derived from lidar measurements. The results are compared to coincidental airborne in-situ measurements of particle number and surface area concentrations in the dust layer. The observations were performed in long-range transported Saharan dust at Barbados and Asian dust at Cyprus. The Asian dust data analysis is ongoing. A comparison of Asian and Saharan dust will be given at the conference. Concentrations of ice nucleating particles in the order of 10 to 1000 per cm-3 in the dust layer are derived for a temperature of-25°C at Barbados. The method can be used to continuously monitor the concentration of ice nucleating dust particles vertically resolved from lidar measurements. © 2019 The Authors, published by EDP Sciences