In this study we use a new dust product developed using CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation) observations and EARLINET (European Aerosol Research Lidar Network) measurements and methods to provide a 3-D multiyear analysis on the evolution of Saharan dust over North Africa and Europe. The product uses a CALIPSO L2 backscatter product corrected with a depolarization-based method to separate pure dust in external aerosol mixtures and a Saharan dust lidar ratio (LR) based on long-term EARLINET measurements to calculate the dust extinction profiles. The methodology is applied on a 9-year CALIPSO dataset (2007–2015) and the results are analyzed here to reveal for the first time the 3-D dust evolution and the seasonal patterns of dust over its transportation paths from the Sahara towards the Mediterranean and Continental Europe. During spring, the spatial distribution of dust shows a uniform pattern over the Sahara desert. The dust transport over the Mediterranean Sea results in mean dust optical depth (DOD) values up to 0.1. During summer, the dust activity is mostly shifted to the western part of the desert where mean DOD near the source is up to 0.6. Elevated dust plumes with mean extinction values between 10 and 75 Mm⁻¹ are observed throughout the year at various heights between 2 and 6 km, extending up to latitudes of 40° N. Dust advection is identified even at latitudes of about 60° N, but this is due to rare events of episodic nature. Dust plumes of high DOD are also observed above the Balkans during the winter period and above northwest Europe during autumn at heights between 2 and 4 km, reaching mean extinction values up to 50 Mm⁻¹. The dataset is considered unique with respect to its potential applications, including the evaluation of dust transport models and the estimation of cloud condensation nuclei (CCN) and ice nuclei (IN) concentration profiles. Finally, the product can be used to study dust dynamics during transportation, since it is capable of revealing even fine dynamical features such as the particle uplifting and deposition on European mountainous ridges such as the Alps and Carpathian Mountains.

We present LIVAS (LIdar climatology of Vertical Aerosol Structure for space-based lidar simulation studies), a 3-D multi-wavelength global aerosol and cloud optical database, optimized to be used for future space-based lidar end-to-end simulations of realistic atmospheric scenarios as well as retrieval algorithm testing activities. The LIVAS database provides averaged profiles of aerosol optical properties for the potential spaceborne laser operating wavelengths of 355, 532, 1064, 1570 and 2050 nm and of cloud optical properties at the wavelength of 532 nm. The global database is based on CALIPSO observations at 532 and 1064 nm and on aerosol-type-dependent backscatter- and extinction-related Ångström exponents, derived from EARLINET (European Aerosol Research Lidar Network) ground-based measurements for the UV and scattering calculations for the IR wavelengths, using a combination of input data from AERONET, suitable aerosol models and recent literature. The required spectral conversions are calculated for each of the CALIPSO aerosol types and are applied to CALIPSO backscatter and extinction data corresponding to the aerosol type retrieved by the CALIPSO aerosol classification scheme. A cloud optical database based on CALIPSO measurements at 532 nm is also provided, neglecting wavelength conversion due to approximately neutral scattering behavior of clouds along the spectral range of LIVAS. Averages of particle linear depolarization ratio profiles at 532 nm are provided as well. Finally, vertical distributions for a set of selected scenes of specific atmospheric phenomena (e.g., dust outbreaks, volcanic eruptions, wild fires, polar stratospheric clouds) are analyzed and spectrally converted so as to be used as case studies for spaceborne lidar performance assessments. The final global data set includes 4-year (1 January 2008–31 December 2011) time-averaged CALIPSO (Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations) data on a uniform grid of 1° × 1° with the original high vertical resolution of CALIPSO in order to ensure realistic simulations of the atmospheric variability in lidar end-to-end simulations.

Six months of stratospheric aerosol observations with the European Aerosol Research Lidar Network (EARLINET) from August 2017 to January 2018 are presented. The decay phase of an unprecedented, record-breaking stratospheric perturbation caused by wildfire smoke is reported and discussed in terms of geometrical, optical, and microphysical aerosol properties. Enormous amounts of smoke were injected into the upper troposphere and lower stratosphere over fire areas in western Canada on 12 August 2017 during strong thunderstorm–pyrocumulonimbus activity. The stratospheric fire plumes spread over the entire Northern Hemisphere in the following weeks and months. Twenty-eight European lidar stations from northern Norway to southern Portugal and the eastern Mediterranean monitored the strong stratospheric perturbation on a continental scale. The main smoke layer (over central, western, southern, and eastern Europe) was found at heights between 15 and 20 km since September 2017 (about 2 weeks after entering the stratosphere). Thin layers of smoke were detected at heights of up to 22–23 km. The stratospheric aerosol optical thickness at 532 nm decreased from values > 0.25 on 21–23 August 2017 to 0.005–0.03 until 5–10 September and was mainly 0.003–0.004 from October to December 2017 and thus was still significantly above the stratospheric background (0.001–0.002). Stratospheric particle extinction coefficients (532 nm) were as high as 50–200 Mm⁻¹ until the beginning of September and on the order of 1 Mm⁻¹ (0.5–5 Mm⁻¹) from October 2017 until the end of January 2018. The corresponding layer mean particle mass concentration was on the order of 0.05–0.5 µg m⁻³ over these months. Soot particles (light-absorbing carbonaceous particles) are efficient ice-nucleating particles (INPs) at upper tropospheric (cirrus) temperatures and available to influence cirrus formation when entering the tropopause from above. We estimated INP concentrations of 50–500 L⁻¹ until the first days in September and afterwards 5–50 L⁻¹ until the end of the year 2017 in the lower stratosphere for typical cirrus formation temperatures of −55 ^∘C and an ice supersaturation level of 1.15. The measured profiles of the particle linear depolarization ratio indicated a predominance of nonspherical smoke particles. The 532 nm depolarization ratio decreased slowly with time in the main smoke layer from values of 0.15–0.25 (August–September) to values of 0.05–0.10 (October–November) and < 0.05 (December–January). The decrease of the depolarization ratio is consistent with aging of the smoke particles, growing of a coating around the solid black carbon core (aggregates), and thus change of the shape towards a spherical form. We found ascending aerosol layer features over the most southern European stations, especially over the eastern Mediterranean at 32–35^∘ N, that ascended from heights of about 18–19 to 22–23 km from the beginning of October to the beginning of December 2017 (about 2 km per month). We discuss several transport and lifting mechanisms that may have had an impact on the found aerosol layering structures.

We present an automatic aerosol classification method based solely on the European Aerosol Research Lidar Network (EARLINET) intensive optical parameters with the aim of building a network-wide classification tool that could provide near-real-time aerosol typing information. The presented method depends on a supervised learning technique and makes use of the Mahalanobis distance function that relates each unclassified measurement to a predefined aerosol type. As a first step (training phase), a reference dataset is set up consisting of already classified EARLINET data. Using this dataset, we defined 8 aerosol classes: clean continental, polluted continental, dust, mixed dust, polluted dust, mixed marine, smoke, and volcanic ash. The effect of the number of aerosol classes has been explored, as well as the optimal set of intensive parameters to separate different aerosol types. Furthermore, the algorithm is trained with literature particle linear depolarization ratio values. As a second step (testing phase), we apply the method to an already classified EARLINET dataset and analyze the results of the comparison to this classified dataset. The predictive accuracy of the automatic classification varies between 59 % (minimum) and 90 % (maximum) from 8 to 4 aerosol classes, respectively, when evaluated against pre-classified EARLINET lidar. This indicates the potential use of the automatic classification to all network lidar data. Furthermore, the training of the algorithm with particle linear depolarization values found in the literature further improves the accuracy with values for all the aerosol classes around 80 %. Additionally, the algorithm has proven to be highly versatile as it adapts to changes in the size of the training dataset and the number of aerosol classes and classifying parameters. Finally, the low computational time and demand for resources make the algorithm extremely suitable for the implementation within the single calculus chain (SCC), the EARLINET centralized processing suite.

The CALIPSO Level 3 (CL3) product is the most recent data set produced by the observations of the Cloud–Aerosol Lidar with Orthogonal Polarization (CALIOP) instrument onboard the Cloud–Aerosol Lidar and Pathfinder Satellite Observations (CALIPSO) space platform. The European Aerosol Research Lidar Network (EARLINET), based mainly on multi-wavelength Raman lidar systems, is the most appropriate ground-based reference for CALIPSO calibration/validation studies on a continental scale. In this work, CALIPSO data are compared against EARLINET monthly averaged profiles obtained by measurements performed during CALIPSO overpasses. In order to mitigate uncertainties due to spatial and temporal differences, we reproduce a modified version of CL3 data starting from CALIPSO Level 2 (CL2) data. The spatial resolution is finer and nearly 2° × 2° (latitude × longitude) and only simultaneous measurements are used for ease of comparison. The CALIPSO monthly mean profiles following this approach are called CALIPSO Level 3*, CL3*. We find good agreement on the aerosol extinction coefficient, yet in most of the cases a small CALIPSO underestimation is observed with an average bias of 0.02 km⁻¹ up to 4 km and 0.003 km⁻¹ higher above. In contrast to CL3 standard product, the CL3* data set offers the possibility to assess the CALIPSO performance also in terms of the particle backscatter coefficient keeping the same quality assurance criteria applied to extinction profiles. The mean relative difference in the comparison improved from 25 % for extinction to 18 % for backscatter, showing better performances of CALIPSO backscatter retrievals. Additionally, the aerosol typing comparison yielded a robust identification of dust and polluted dust. Moreover, the CALIPSO aerosol-type-dependent lidar ratio selection is assessed by means of EARLINET observations, so as to investigate the performance of the extinction retrievals. The aerosol types of dust, polluted dust, and clean continental showed noticeable discrepancy. Finally, the potential improvements of the lidar ratio assignment have been examined by adjusting it according to EARLINET-derived values.

In this paper, we report the first systematic comparison of 12-year modeled dust extinction profiles vs. Raman lidar measurements. We use the BSC-DREAM8b model, one of the most widely used dust regional models in the Mediterranean, and Potenza EARLINET lidar profiles for Saharan dust cases, the largest one-site database of dust extinction profiles. A total of 310 dust cases were compared for the May 2000–July 2012 period. The model reconstructs the measured layers well: profiles are correlated within 5% of significance for 60% of the cases and the dust layer center of mass as measured by lidar and modeled by BSC-DREAM8b differ on average 0.3 ± 1.0 km. Events with a dust optical depth lower than 0.1 account for 70% of uncorrelated profiles. Although there is good agreement in terms of profile shape and the order of magnitude of extinction values, the model overestimates the occurrence of dust layer top above 10 km. Comparison with extinction profiles measured by the Raman lidar shows that BSC-DREAM8b typically underestimates the dust extinction coefficient, in particular below 3 km. Lowest model–observation differences (below 17%) correspond to a lidar ratio at 532 nm and Ångström exponent at 355/532 nm of 60 ± 13 and 0.1 ± 0.6 sr, respectively. These are in agreement with values typically observed and modeled for pure desert dust. However, the highest differences (higher than 85%) are typically related to greater Ångström values (0.5 ± 0.6), denoting smaller particles. All these aspects indicate that the level of agreement decreases with an increase in mixing/modification processes.