We report on the vertical distributions of Saharan dust aerosols over the N.E. Mediterranean region, which were obtained during a typical dust outbreak on August 2000, by two lidar systems located in Athens and Thessaloniki, Greece, in the frame of the European EARLINET project. MODIS and ground sun spectrophotometric data, as well as air-mass backward trajectories confirmed the existence of Saharan dust in the case examined, which was also successfully forecasted by the DREAM dust model. The lidar data analysis for the period 2000-2002 made possible, for the first time, an estimation of the vertical extent of free tropospheric dust layers [mean values of the aerosol backscatter and extinction coefficients and the extinction-to-backscatter ratio (lidar ratio, LR) at 355 nm], as well as a seasonal distribution of Saharan dust outbreaks over Greece, under cloud-free conditions. A mean value of the lidar ratio at 355 nm was obtained over Athens (53±1 sr) and over Thessaloniki (44±2 sr) during the Saharan dust outbreaks. The corresponding aerosol optical thickness (AOT) at 355 nm, in the altitude range 0-5 km, was 0.69±0.12 and 0.65±0.10 for Athens and Thessaloniki, respectively (within the dust layer the AOT was 0.23 and 0.21, respectively). Air-mass back-trajectory analysis performed in the period 2000-2002 for all Saharan dust outbreaks over the N.E. Mediterranean indicated the main pathways followed by the dust aerosols.

Air quality forecast systems need reliable and accurate representations of the planetary boundary layer (PBL) to perform well. An important question is how accurately numerical weather prediction models can reproduce conditions in diverse synoptic flow types. Here, observations from the summer 2014 HygrA-CD (Hygroscopic Aerosols to Cloud Droplets) experimental campaign are used to validate simulations from the Weather Research and Forecasting (WRF) model over the complex, urban terrain of the Greater Athens Area. Three typical atmospheric flow types were identified during the 39-day campaign based on 2-day backward trajectories: Continental, Etesians, and Saharan. It is shown that the numerical model simulations differ dramatically depending on the PBL scheme, atmospheric dynamics, and meteorological parameter (e.g., 2-m air temperature). Eight PBL schemes from WRF version 3.4 are tested with daily simulations on an inner domain at 1-km grid spacing. Near-surface observations of 2-m air temperature and relative humidity and 10-m wind speed are collected from multiple meteorological stations. Estimates of the PBL height come from measurements using a multiwavelength Raman lidar, with an adaptive extended Kalman filter technique. Vertical profiles of atmospheric variables are obtained from radiosonde launches, along with PBL heights calculated using bulk Richardson number. Daytime maximum PBL heights ranged from 2.57 km during Etesian flows, to as low as 0.37 km during Saharan flows. The largest differences between model and observations are found with simulated PBL height during Saharan synoptic flows. During the daytime, campaign-averaged near-surface variables show WRF tended to have a cool, moist bias with higher simulated wind speeds than the observations, especially near the coast. It is determined that non-local PBL schemes give the most agreeable solutions when compared with observations.Peer ReviewedPostprint (published version

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.

This paper presents a detailed description of LIRIC (LIdar-Radiometer Inversion Code) algorithm for simultaneous processing of coincident lidar and radiometric (sun photometric) observations for the retrieval of the aerosol concentration vertical profiles. As the lidar/radiometric input data we use measurements from European Aerosol Research Lidar Network (EARLINET) lidars and collocated sun-photometers of Aerosol Robotic Network (AERONET). The LIRIC data processing provides sequential inversion of the combined lidar and radiometric data. The algorithm starts with the estimations of column-integrated aerosol parameters from radiometric measurements followed by the retrieval of height dependent concentrations of fine and coarse aerosols from lidar signals using integrated column characteristics of aerosol layer as a priori constraints. The use of polarized lidar observations allows us to discriminate between spherical and non-spherical particles of the coarse aerosol mode.The LIRIC software package was implemented and tested at a number of EARLINET stations. Intercomparison of the LIRIC-based aerosol retrievals was performed for the observations by seven EARLINET lidars in Leipzig, Germany on 25 May 2009. We found close agreement between the aerosol parameters derived from different lidars that supports high robustness of the LIRIC algorithm. The sensitivity of the retrieval results to the possible reduction of the available observation data is also discussed.