Search

Aerosol–cloud interactions in mixed-phase clouds (MPCs) are one of the most uncertain drivers of the hydrological cycle and climate change. A synergy of in situ, remote-sensing and modelling experiments were used to determine the source of ice-nucleating particles (INPs) for MPCs at Mount Helmos in the eastern Mediterranean. The influences of boundary layer turbulence, vertical aerosol distributions and meteorological conditions were also examined. When the observation site is in the free troposphere (FT), approximately 1 in ×106 aerosol particles serve as INPs around −25 °C. The INP abundance spans 3 orders of magnitude and increases in the following order: marine aerosols; continental aerosols; and, finally, dust plumes. Biological particles are important INPs observed in continental and marine aerosols, whereas they play a secondary, although important, role during Saharan dust events. Air masses in the planetary boundary layer (PBL) show both enriched INP concentrations and a higher proportion of INPs to total aerosol particles, compared with cases in the FT. The presence of precipitation/clouds enriches INPs in the FT but decreases INPs in the PBL. Additionally, new INP parameterizations are developed that incorporate the ratio of fluorescent-to-nonfluorescent or coarse-to-fine particles and predict >90 % of the observed INPs within an uncertainty range of a factor of 10; these new parameterizations exhibit better performance than current widely used parameterizations and allow ice formation in models to respond to variations in dust and biological particles. The improved parameterizations can help MPC formation simulations in regions with various INP sources or different regions with prevailing INP sources.

The purpose of this study is to understand the drivers of cloud droplet formation in orographic clouds. We used a combination of modeling, in situ, and remote sensing measurements at the high-altitude Helmos Hellenic Atmospheric Aerosol and Climate Change ((HAC)2) station, which is located at the top of Mt. Helmos (1314 m above sea level), Greece, during the Cloud–AerosoL InteractionS in the Helmos Background TropOsphere (CALISHTO) campaign in fall 2021 (https://calishto.panacea-ri.gr/, last access: 1 August 2024) to examine the origins of the aerosols (i.e., local aerosol from the planetary boundary layer (PBL) or long-range-transported aerosol from the free-tropospheric layer (FTL) contributing to the cloud condensation nuclei (CCN)), their characteristics (hygroscopicity, size distribution, and mixing state), and the vertical velocity distributions and resulting supersaturations. We found that the characteristics of the PBL aerosol were considerably different from FTL aerosol and use the aerosol particle number and equivalent mass concentration of the black carbon (eBC) in order to determine when (HAC)2 was within the FTL or PBL based on time series of the height of the PBL. During the (HAC)2 cloud events we sample a mixture of interstitial aerosol and droplet residues, which we characterize using a new approach that utilizes the in situ droplet measurements to determine time periods when the aerosol sample is purely interstitial. From the dataset we determine the properties (size distribution and hygroscopicity) of the pre-cloud, activated, and interstitial aerosol. The hygroscopicity of activated aerosol is found to be higher than that of the interstitial or pre-cloud aerosol. A series of closure studies with the droplet parameterization shows that cloud droplet concentration (Nd) and supersaturation can be predicted to within 25 % of observations when the aerosol size distributions correspond to pre-cloud conditions. The analysis of the characteristic supersaturation of each aerosol population indicates that droplet formation in clouds is aerosol-limited when formed in FTL air masses – hence droplet formation is driven by aerosol variations, while clouds formed in the PBL tend to be velocity-limited and droplet variations are driven by fluctuations in vertical velocity. Given that the cloud dynamics do not vary significantly between air masses, the variation in aerosol concentration and type is mostly responsible for these shifts in cloud microphysical state and sensitivity to aerosol. With these insights, the remote sensing of cloud droplets in such clouds can be used to infer either CCN spectra (when in the FTL) or vertical velocity (when in the PBL). In conclusion, we show that a coordinated measurement of aerosol and cloud properties, together with the novel analysis approaches presented here, allows for the determination of the drivers of droplet formation in orographic clouds and their sensitivity to aerosol and vertical velocity variations.

The Cloud-AerosoL InteractionS in the Helmos background TropOsphere (CALISHTO) campaign took place in autumn 2021 at the NCSR Demokritos background high-altitude Helmos Hellenic Atmospheric Aerosol and Climate Change station (HAC)2 to study the interactions between aerosols and clouds. The current study presents the chemical characterization of the non-refractory (NR) PM1 aerosol fraction using a time-of-flight aerosol chemical speciation monitor (ToF-ACSM). A comparative offline aerosol filter analysis by a high-resolution time-of-flight aerosol mass spectrometer (HR-ToF-AMS) showed consistent results regarding the species determined. Source apportionment applied on both datasets (ACSM-ToF and offline AMS analysis on filter extracts) yielded the same factors for the organic aerosol (one primary and two secondary factors). Additionally, the positive matrix factorization (PMF) model was applied on the total PM1 fraction by the ToF-ACSM (including both organic and inorganic ions). Five different types were identified, including a primary organic factor; ammonium nitrate; ammonium sulfate; and two secondary organic aerosols, one more oxidized and one less oxidized. The prevailing atmospheric conditions at the station, i.e., cloud presence, influence of emissions from the planetary boundary layer (PBL), and air mass origin, were also incorporated in the study. The segregation between PBL and free-troposphere (FT) conditions was made by combining data from remote sensing and in situ measurement techniques. The types of air masses arriving at the site were grouped as continental, marine, dust, and marine–dust based on back-trajectory data. Significant temporal variability in the aerosol characteristics was observed throughout the campaign; in September, air masses from within the PBL were sampled most of the time, resulting in much higher mass concentrations compared to October and November when concentrations were reduced by a factor of 5. Both in-cloud and FT measurement periods resulted in much lower concentration levels, while a similar composition was observed in PBL and FT conditions. We take advantage of using a recently developed “virtual-filtering” technique to separate interstitial from activated aerosol sampled from a PM10 inlet during cloudy periods. This allows the determination of the chemical composition of the interstitial aerosol during in-cloud periods. Ammonium sulfate, the dominant PMF factor in all conditions, contributed more when air masses were arriving at (HAC)2 during dust events, while a higher secondary organic aerosol contribution was observed when air masses arrived from continental Europe.

In August 2021, a historic heatwave was recorded in Greece which resulted in extreme wildfire events that strongly affected the air quality over the city of Athens. Saharan dust was also transferred over Greece on certain days of the same period due to the prevailing southern winds. The impact of these events on air quality and surface solar radiation is investigated in this study. Event characterization based on active and passive remote sensing instrumentation has been performed. The study shows that significantly increased levels of air pollution were recorded from the end of July to the first week of August. The smoke led to unusually high aerosol optical depth (AOD) values (up to 3.6 at 500 nm), high Ångström exponent (AE) (up to 2.4 at 440–870 nm), and a strong and negative dependence of single-scattering albedo (SSA) on wavelength that was observed to decrease from 0.93 at 440 nm to 0.86 at 1020 nm, while the dust event led to high AOD (up to 0.7 at 500 nm), low AE (up to 0.9 at 440–870 nm), and a positive dependence of SSA on wavelength that was observed to increase from 0.89 at 440 nm to 0.95 at 1020. Furthermore, the smoke plume was also detected over the PANhellenic GEophysical observatory of Antikythera on 7 August, which is about 240 km away from Athens. Increased AOD values (up to ∼ 0.90 at 500 nm) associated with a high fine-mode AOD (up to ∼ 0.85 at 500 nm) and decrease in SSA with wavelength suggested the dominance of fine biomass burning aerosols. The impact of dust and smoke on solar irradiance revealed significant differences in the spectral dependence of the attenuation caused by the two different aerosol types. The attenuation of solar irradiance in the ultraviolet (UV-B) spectrum was found to be much lower in the case of dust compared to smoke for similar AOD500 values. Differences were less pronounced in the near-infrared and visible spectral regions. The large AODs during the wildfires resulted in a decrease in the noon UV index by up to 53 %, as well as in the daily effective doses for the production of vitamin D (up to 50 %), in the daily photosynthetically active radiation (up to 21 %) and in the daily global horizontal irradiance (up to 17 %), with serious implications for health, agriculture, and energy. This study highlights the wider impacts of wildfires that are part of the wider problem for Mediterranean countries, whose frequency is predicted to increase in view of the projected increasing occurrence of summer heatwaves.

Since 2018, the Aeolus satellite of the European Space Agency (ESA) has acquired wind HLOS (horizontal line-of-sight) profiles throughout the troposphere and up to the lower stratosphere, filling a critical gap in the Global Observing System (GOS). Aeolus, carrying ALADIN (Atmospheric LAser Doppler INstrument), the first UV HSRL (High Spectral Resolution Lidar) Doppler lidar ever placed in space, provides also vertically resolved optical properties of particulates (aerosols and clouds). The present study focuses on the assessment of Aeolus L2A particulate backscatter coefficient (baseline 2A11), retrieved by the Standard Correct Algorithm (SCA), in the eastern Mediterranean, a region hosting a variety of aerosol species. Ground-based retrievals acquired by lidar instruments operating in Athens (central Greece), Thessaloniki (northern Greece) and Antikythera (southwestern Greece) serve as reference. All lidar stations provide routine measurements to the PANACEA (PANhellenic infrastructure for Atmospheric Composition and climatE chAnge) network. A set of ancillary data, including sun-photometric observations (AERONET), reanalysis products (CAMS and MERRA-2), satellite observations (MSG-SEVIRI and MODIS Aqua) and backward-trajectories modelling (FLEXPART), is utilized towards an optimum characterization of the probed atmospheric conditions under the absence of a classification scheme in Aeolus SCA profiles. First, emphasis is given on the assessment of Aeolus SCA backscatter coefficient under specific aerosol scenarios over Antikythera island. Due to the misdetection of the cross-polar component of the backscattered lidar signal, Aeolus underestimates the aerosol backscatter coefficient by up to 33 % when non-spherical mineral particles are recorded (10 July 2019). A good performance is revealed on 3 July 2019, when horizontally homogeneous loads of fine spherical particles are confined below 4 km. For other two cases (8 July 2020 and 5 August 2020), due to noise issues, the SCA performance degrades in terms of depicting the stratification of aerosol layers composed of particles of different origin. According to the statistical assessment analysis of 43 identified cases, a poor-to-moderate performance is revealed for the unfiltered (aerosols plus clouds) SCA profiles, which improves substantially when cloud-contaminated profiles are excluded from the collocated sample. This improvement is evident at both Aeolus vertical scales (regular scales have 24 bins and mid-bin scales have 23 bins), and it is justified by the drastic reduction in the bias (from 0.45 to 0.27 Mm−1 sr−1 for SCA and from 0.69 to 0.37 Mm−1 sr−1 for SCA mid-bin) and root mean square error (from 2.00 to 1.65 Mm−1 sr−1 for SCA and from 1.88 to 1.00 Mm−1 sr−1 for SCA mid-bin) scores. In the vertical, the SCA performance degrades at the lowermost bins due to either the contamination from surface signals or the increased noise levels for the aerosol retrievals. Among the three PANACEA stations, the best agreement is found at the remote site of Antikythera with respect to the urban sites of Athens and Thessaloniki. Finally, all key Cal/Val (calibration and validation) aspects necessary for future relevant studies, the recommendations for a possible Aeolus follow-on mission and an overview of the ongoing related activities are thoroughly discussed.

The purpose of this study is to investigate the ability of the Sentinel-5P TROPOspheric Monitoring Instrument (TROPOMI) to derive accurate geometrical features of lofted aerosol layers, selecting the Mediterranean Basin as the study area. Comparisons with ground-based correlative measurements constitute a key component in the validation of passive and active satellite aerosol products. For this purpose, we use ground-based observations from quality-controlled lidar stations reporting to the European Aerosol Research Lidar Network (EARLINET). An optimal methodology for validation purposes has been developed and applied using the EARLINET optical profiles and TROPOMI aerosol products, aiming at the in-depth evaluation of the TROPOMI aerosol layer height (ALH) product for the period 2018 to 2022 over the Mediterranean Basin. Seven EARLINET stations were chosen, taking into consideration their proximity to the sea, which provided 63 coincident aerosol cases for the satellite retrievals. In the following, we present the first validation results for the TROPOMI/S5P ALH using the optimized EARLINET lidar products employing the automated validation chain designed for this purpose. The quantitative validation at pixels over the selected EARLINET stations illustrates that the TROPOMI ALH product is consistent with the EARLINET lidar products, with a high correlation coefficient R=0.82 (R=0.51) and a mean bias of -0.51±0.77 km and -2.27±1.17 km over ocean and land, respectively. Overall, it appears that aerosol layer altitudes retrieved from TROPOMI are systematically lower than altitudes from the lidar retrievals. High-albedo scenes, as well as low-aerosol-load scenes, are the most challenging for the TROPOMI retrieval algorithm, and these results testify to the need to further investigate the underlying cause. This work provides a clear indication that the TROPOMI ALH product can under certain conditions achieve the required threshold accuracy and precision requirements of 1 km, especially when only ocean pixels are included in the comparison analysis. Furthermore, we describe and analyse three case studies in detail, one dust and two smoke episodes, in order to illustrate the strengths and limitations of the TROPOMI ALH product and demonstrate the presented validation methodology. The present analysis provides important additions to the existing validation studies that have been performed so far for the TROPOMI S5P ALH product, which were based only on satellite-to-satellite comparisons.

We introduce an automated aerosol type classification method, called Source Classification Analysis (SCAN). SCAN is based on predefined and characterized aerosol source regions, the time that the air parcel spends above each geographical region, and a number of additional criteria. The output of SCAN is compared with two independent aerosol classification methods, which use the intensive optical parameters from lidar data: (1) the Mahalanobis distance automatic aerosol type classification (MD) and (2) a neural network aerosol typing algorithm (NATALI). In this paper, data from the European Aerosol Research Lidar Network (EARLINET) have been used. A total of 97 free tropospheric aerosol layers from four typical EARLINET stations (i.e., Bucharest, Kuopio, Leipzig, and Potenza) in the period 2014–2018 were classified based on a 3β+2α+1δ lidar configuration. We found that SCAN, as a method independent of optical properties, is not affected by overlapping optical values of different aerosol types. Furthermore, SCAN has no limitations concerning its ability to classify different aerosol mixtures. Additionally, it is a valuable tool to classify aerosol layers based on even single (elastic) lidar signals in the case of lidar stations that cannot provide a full data set (3β+2α+1δ) of aerosol optical properties; therefore, it can work independently of the capabilities of a lidar system. Finally, our results show that NATALI has a lower percentage of unclassified layers (4 %), while MD has a higher percentage of unclassified layers (50 %) and a lower percentage of cases classified as aerosol mixtures (5 %).

In this study, we report on the retrieval of aerosol extinction profiles from ground-based scattered sunlight multi-axis differential optical absorption spectroscopy (MAX-DOAS) measurements, carried out at Athens, Greece. It is the first time that aerosol profiles are retrieved from MAX-DOAS measurements in Athens. The reported aerosol vertical distributions at 477 nm are derived from the oxygen dimer (O4) differential-slant-column-density observations at different elevation angles by applying the BRemen Optimal estimation REtrieval for Aerosol and trace gaseS (BOREAS) retrieval algorithm. Four case studies have been selected for validation purposes; the retrieved aerosol profiles and the corresponding aerosol optical depths (AODs) from the MAX-DOAS are compared with lidar extinction profiles and with sun-photometric measurements (Aerosol Robotic Network, AERONET, observations), respectively. Despite the different approach of each method regarding the retrieval of the aerosol information, the comparison with the lidar measurements at 532 nm reveals a very good agreement in terms of vertical distribution, with r>0.90 in all cases. The AODs from the MAX-DOAS and the sun photometer (the latter at 500 nm) show a satisfactory correlation (with 0.45 < r < 0.7 in three out of the four cases). The comparison indicates that the MAX-DOAS systematically underestimates the AOD in the cases of large particles (small Ångström exponent) and for measurements at small relative azimuthal angles between the viewing direction and the sun. Better agreement is achieved in the morning, at large relative azimuthal angles. Overall, the aerosol profiles retrieved from MAX-DOAS measurements are of good quality; thus, new perspectives are opened up for assessing urban aerosol pollution on a long-term basis in Athens from continuous and uninterrupted MAX-DOAS measurements.

Remote sensing measurements of aerosols using depolarization Raman lidar systems from four EARLINET (European Aerosol Research Lidar Network) stations are used for a comprehensive analysis of Saharan dust events over the Mediterranean basin in the period 2014–2017. In this period, 51 dust events regarding the geometrical, optical and microphysical properties of dust were selected, classified and assessed according to their radiative forcing effect on the atmosphere. From west to east, the stations of Granada, Potenza, Athens and Limassol were selected as representative Mediterranean cities regularly affected by Saharan dust intrusions. Emphasis was given on lidar measurements in the visible (532 nm) and specifically on the consistency of the particle linear depolarization ratio (δp532), the extinction-to-backscatter lidar ratio (LR532) and the aerosol optical thickness (AOT532) within the observed dust layers. We found mean δp532 values of 0.24±0.05, 0.26±0.06, 0.28±0.05 and 0.28±0.04, mean LR532 values of 52±8, 51±9, 52±9 and 49±6 sr and mean AOT532 values of 0.40±0.31, 0.11±0.07, 0.12±0.10 and 0.32±0.17, for Granada, Potenza, Athens and Limassol, respectively. The mean layer thickness values were found to range from ∼ 1700 to ∼ 3400 m a.s.l. Additionally, based also on a previous aerosol type classification scheme provided by airborne High Spectral Resolution Lidar (HSRL) observations and on air mass backward trajectory analysis, a clustering analysis was performed in order to identify the mixing state of the dusty layers over the studied area. Furthermore, a synergy of lidar measurements and modeling was used to analyze the solar and thermal radiative forcing of airborne dust in detail. In total, a cooling behavior in the solar range and a significantly lower heating behavior in the thermal range was estimated. Depending on the dust optical and geometrical properties, the load intensity and the solar zenith angle (SZA), the estimated solar radiative forcing values range from −59 to −22 W m−2 at the surface and from −24 to −1 W m−2 at the top of the atmosphere (TOA). Similarly, in the thermal spectral range these values range from +2 to +4 W m−2 for the surface and from +1 to +3 W m−2 for the TOA. Finally, the radiative forcing seems to be inversely proportional to the dust mixing ratio, since higher absolute values are estimated for less mixed dust layers.

The methodology of analysing the biomass burning events recorded in the database of the European Aerosol Research Lidar Network in the framework of the Aerosol, Clouds and Trace Gases Research Infrastructure is presented. The period of 2008–2017 was chosen to analyse all of the events stored in the database under the Forest Fire category for a total of 14 stations available. The data provided ranged from complete datasets (particle backscatter, extinction and linear depolarization ratio profiles) to single profiles (particle backscatter coefficient profile). Smoke layers geometry was evaluated and the mean optical properties within each layer were computed. The back-trajectory technique was used to double-check the source of all pollution layers. The biomass burning layers were identified by taking into account the presence of the fires along the back trajectory. The biomass burning events are analysed by the means of the intensive parameters. The analysis was structured in three directions: (I) common biomass burning source (fire) recorded by at least two stations, (II) long-range transport from North America, and (III) analysis over four geographical regions (south-eastern Europe, north-eastern Europe, central Europe, and south-western Europe). Based on back-trajectory calculations and fire locations, the lidar measurements can be labelled either as measurements of a “single fire” or “mixed fires” (case I), measurements of North American fires, or measurements of mixed North American and local fires (case II). The histogram of the fire locations reveals the smoke sources for each region. For each region, statistics on intensive parameters are performed. The source origin of the intensive parameters is categorized based on the continental origin of the air mass (European, African, Asian, North American, or a combination of them). The methodology presented here is meant to provide a perspective to explore a large number of lidar data and deliver novel approaches to analyse the intensive parameters based on the assigned biomass burning sources. A thorough consideration of all potential fire sources reveals that most of the time the lidar measurements characterize the smoke from a mixture of fires. A comprehensive discussion of all the results (based on the intensive parameters and the source locations) will be given in a companion paper submitted to the ACP EARLINET special issue.

A stand-alone lidar-based method for detecting airborne hazards for aviation in near real time (NRT) is presented. A polarization lidar allows for the identification of irregular-shaped particles such as volcanic dust and desert dust. The Single Calculus Chain (SCC) of the European Aerosol Research Lidar Network (EARLINET) delivers high-resolution preprocessed data: the calibrated total attenuated backscatter and the calibrated volume linear depolarization ratio time series. From these calibrated lidar signals, the particle backscatter coefficient and the particle depolarization ratio can be derived in temporally high resolution and thus provide the basis of the NRT early warning system (EWS). In particular, an iterative method for the retrieval of the particle backscatter is implemented. This improved capability was designed as a pilot that will produce alerts for imminent threats for aviation. The method is applied to data during two diverse aerosol scenarios: first, a record breaking desert dust intrusion in March 2018 over Finokalia, Greece, and, second, an intrusion of volcanic particles originating from Mount Etna, Italy, in June 2019 over Antikythera, Greece. Additionally, a devoted observational period including several EARLINET lidar systems demonstrates the network's preparedness to offer insight into natural hazards that affect the aviation sector.

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−1 until the beginning of September and on the order of 1 Mm−1 (0.5–5 Mm−1) 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−3 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−1 until the first days in September and afterwards 5–50 L−1 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 the evaluation activity of the European Aerosol Research Lidar Network (EARLINET) for the quantitative assessment of the Level 2 aerosol backscatter coefficient product derived by the Cloud-Aerosol Transport System (CATS) aboard the International Space Station (ISS; Rodier et al., 2015). The study employs correlative CATS and EARLINET backscatter measurements within a 50 km distance between the ground station and the ISS overpass and as close in time as possible, typically with the starting time or stopping time of the EARLINET performed measurement time window within 90 min of the ISS overpass, for the period from February 2015 to September 2016. The results demonstrate the good agreement of the CATS Level 2 backscatter coefficient and EARLINET. Three ISS overpasses close to the EARLINET stations of Leipzig, Germany; Évora, Portugal; and Dushanbe, Tajikistan, are analyzed here to demonstrate the performance of the CATS lidar system under different conditions. The results show that under cloud-free, relative homogeneous aerosol conditions, CATS is in good agreement with EARLINET, independent of daytime and nighttime conditions. CATS low negative biases are observed, partially attributed to the deficiency of lidar systems to detect tenuous aerosol layers of backscatter signal below the minimum detection thresholds; these are biases which may lead to systematic deviations and slight underestimations of the total aerosol optical depth (AOD) in climate studies. In addition, CATS misclassification of aerosol layers as clouds, and vice versa, in cases of coexistent and/or adjacent aerosol and cloud features, occasionally leads to non-representative, unrealistic, and cloud-contaminated aerosol profiles. Regarding solar illumination conditions, low negative biases in CATS backscatter coefficient profiles, of the order of 6.1 %, indicate the good nighttime performance of CATS. During daytime, a reduced signal-to-noise ratio by solar background illumination prevents retrievals of weakly scattering atmospheric layers that would otherwise be detectable during nighttime, leading to higher negative biases, of the order of 22.3 %.

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.

In this study we investigate the climatological behavior of the aerosol optical properties over Thessaloniki during the years 2003–2017. For this purpose, measurements of two independent instruments, a lidar and a sunphotometer, were used. These two instruments represent two individual networks, the European Lidar Aerosol Network (EARLINET) and the Aerosol Robotic Network (AERONET). They include different measurement schedules. Fourteen years of lidar and sunphotometer measurements were analyzed, independently of each other, in order to obtain the annual cycles and trends of various optical and geometrical aerosol properties in the boundary layer, in the free troposphere, and for the whole atmospheric column. The analysis resulted in consistent statistically significant and decreasing trends of aerosol optical depth (AOD) at 355 nm of −23.2 and −22.3 % per decade in the study period over Thessaloniki for the EARLINET and the AERONET datasets, respectively. Therefore, the analysis indicates that the EARLINET sampling schedule can be quite effective in producing data that can be applied to long-term climatological studies. It is also shown that the observed decreasing trend is mainly attributed to changes in the aerosol load inside the boundary layer. Seasonal profiles of the most dominant aerosol mixture types observed over Thessaloniki have been generated from the lidar data. The higher values of the vertically resolved extinction coefficient at 355 nm appear in summer, while the lower ones appear in winter. The dust component is more dominant in the free troposphere than in the boundary layer during summer. The biomass burning layers tend to arrive in the free troposphere during spring and summer. This kind of information can be quite useful for applications that require a priori aerosol profiles. For instance, they can be utilized in models that require aerosol climatological data as input, in the development of algorithms for satellite products, and also in passive remote-sensing techniques that require knowledge of the aerosol vertical distribution.

Accurate continuous measurements of relative humidity (RH) vertical profiles in the lower troposphere have become a significant scientific challenge. In recent years a synergy of various ground-based remote sensing instruments have been successfully used for RH vertical profiling, which has resulted in the improvement of spatial resolution and, in some cases, of the accuracy of the measurement. Some studies have also suggested the use of high-resolution model simulations as input datasets into RH vertical profiling techniques. In this paper we apply two synergetic methods for RH profiling, including the synergy of lidar with a microwave radiometer and high-resolution atmospheric modeling. The two methods are employed for RH retrieval between 100 and 6000 m with increased spatial resolution, based on datasets from the HygrA-CD (Hygroscopic Aerosols to Cloud Droplets) campaign conducted in Athens, Greece from May to June 2014. RH profiles from synergetic methods are then compared with those retrieved using single instruments or as simulated by high-resolution models. Our proposed technique for RH profiling provides improved statistical agreement with reference to radiosoundings by 27 % when the lidar–radiometer (in comparison with radiometer measurements) approach is used and by 15 % when a lidar model is used (in comparison with WRF-model simulations). Mean uncertainty of RH due to temperature bias in RH profiling was ∼ 4.34 % for the lidar–radiometer and ∼ 1.22 % for the lidar–model methods. However, maximum uncertainty in RH retrievals due to temperature bias showed that lidar-model method is more reliable at heights greater than 2000 m. Overall, our results have demonstrated the capability of both combined methods for daytime measurements in heights between 100 and 6000 m when lidar–radiometer or lidar–WRF combined datasets are available.

Particle depolarization ratio retrieved from lidar measurements are commonly used for aerosol-typing studies, microphysical inversion, or mass concentration retrievals. The particle depolarization ratio is one of the primary parameters that can differentiate several major aerosol components but only if the measurements are accurate enough. The accuracy related to the retrieval of particle depolarization ratios is the driving factor for assessing and improving the uncertainties of the depolarization products. This paper presents different depolarization calibration procedures used to improve the quality of the depolarization data. The results illustrate a significant improvement of the depolarization lidar products for all the selected lidar stations that have implemented depolarization calibration procedures. The calibrated volume and particle depolarization profiles at 532 nm show values that fall within a range that is generally accepted in the literature.

Lidar depolarization measurements distinguish between spherical and non-spherical aerosol particles based on the change of the polarization state between the emitted and received signal. The particle shape information in combination with other aerosol optical properties allows the characterization of different aerosol types and the retrieval of aerosol particle microphysical properties. Regarding the microphysical inversions, the lidar depolarization technique is becoming a key method since particle shape information can be used by algorithms based on spheres and spheroids, optimizing the retrieval procedure. Thus, the identification of the depolarization error sources and the quantification of their effects are crucial. This work presents a new tool to assess the systematic error of the volume linear depolarization ratio (δ), combining the Stokes–Müller formalism and the complete sampling of the error space using the lidar model presented in Freudenthaler (2016a). This tool is applied to a synthetic lidar system and to several EARLINET lidars with depolarization capabilities at 355 or 532 nm. The lidar systems show relative errors of δ larger than 100 % for δ values around molecular linear depolarization ratios (∼ 0.004 and up to ∼ 10 % for δ = 0.45). However, one system shows only relative errors of 25 and 0.22 % for δ = 0.004 and δ = 0.45, respectively, and gives an example of how a proper identification and reduction of the main error sources can drastically reduce the systematic errors of δ. In this regard, we provide some indications of how to reduce the systematic errors.

The simultaneous analysis of aerosol microphysical properties profiles at different European stations is made in the framework of the ChArMEx/EMEP 2012 field campaign (9–11 July 2012). During and in support of this campaign, five lidar ground-based stations (Athens, Barcelona, Bucharest, Évora, and Granada) performed 72 h of continuous lidar measurements and collocated and coincident sun-photometer measurements. Therefore it was possible to retrieve volume concentration profiles with the Lidar Radiometer Inversion Code (LIRIC). Results indicated the presence of a mineral dust plume affecting the western Mediterranean region (mainly the Granada station), whereas a different aerosol plume was observed over the Balkans area. LIRIC profiles showed a predominance of coarse spheroid particles above Granada, as expected for mineral dust, and an aerosol plume composed mainly of fine and coarse spherical particles above Athens and Bucharest. Due to the exceptional characteristics of the ChArMEx database, the analysis of the microphysical properties profiles' temporal evolution was also possible. An in-depth analysis was performed mainly at the Granada station because of the availability of continuous lidar measurements and frequent AERONET inversion retrievals. The analysis at Granada was of special interest since the station was affected by mineral dust during the complete analyzed period. LIRIC was found to be a very useful tool for performing continuous monitoring of mineral dust, allowing for the analysis of the dynamics of the dust event in the vertical and temporal coordinates. Results obtained here illustrate the importance of having collocated and simultaneous advanced lidar and sun-photometer measurements in order to characterize the aerosol microphysical properties in both the vertical and temporal coordinates at a regional scale. In addition, this study revealed that the use of the depolarization information as input in LIRIC in the stations of Bucharest, Évora, and Granada was crucial for the characterization of the aerosol types and their distribution in the vertical column, whereas in stations lacking depolarization lidar channels, ancillary information was needed. Results obtained were also used for the validation of different mineral dust models. In general, the models better forecast the vertical distribution of the mineral dust than the column-integrated mass concentration, which was underestimated in most of the cases.

This study investigates the concentration, cloud condensation nuclei (CCN) activity and hygroscopic properties of particles influenced by biomass burning in the eastern Mediterranean and their impacts on cloud droplet formation. Air masses sampled were subject to a range of atmospheric processing (several hours up to 3 days). Values of the hygroscopicity parameter, κ, were derived from CCN measurements and a Hygroscopic Tandem Differential Mobility Analyzer (HTDMA). An Aerosol Chemical Speciation Monitor (ACSM) was also used to determine the chemical composition and mass concentration of non-refractory components of the submicron aerosol fraction. During fire events, the increased organic content (and lower inorganic fraction) of the aerosol decreases the values of κ, for all particle sizes. Particle sizes smaller than 80 nm exhibited considerable chemical dispersion (where hygroscopicity varied up to 100 % for particles of same size); larger particles, however, exhibited considerably less dispersion owing to the effects of condensational growth and cloud processing. ACSM measurements indicate that the bulk composition reflects the hygroscopicity and chemical nature of the largest particles (having a diameter of ∼ 100 nm at dry conditions) sampled. Based on positive matrix factorization (PMF) analysis of the organic ACSM spectra, CCN concentrations follow a similar trend as the biomass-burning organic aerosol (BBOA) component, with the former being enhanced between 65 and 150 % (for supersaturations ranging between 0.2 and 0.7 %) with the arrival of the smoke plumes. Using multilinear regression of the PMF factors (BBOA, OOA-BB and OOA) and the observed hygroscopicity parameter, the inferred hygroscopicity of the oxygenated organic aerosol components is determined. We find that the transformation of freshly emitted biomass burning (BBOA) to more oxidized organic aerosol (OOA-BB) can result in a 2-fold increase of the inferred organic hygroscopicity; about 10 % of the total aerosol hygroscopicity is related to the two biomass-burning components (BBOA and OOA-BB), which in turn contribute almost 35 % to the fine-particle organic water of the aerosol. Observation-derived calculations of the cloud droplet concentrations that develop for typical boundary layer cloud conditions suggest that biomass burning increases droplet number, on average by 8.5 %. The strongly sublinear response of clouds to biomass-burning (BB) influences is a result of strong competition of CCN for water vapor, which results in very low maximum supersaturation (0.08 % on average). Attributing droplet number variations to the total aerosol number and the chemical composition variations shows that the importance of chemical composition increases with distance, contributing up to 25 % of the total droplet variability. Therefore, although BB may strongly elevate CCN numbers, the impact on droplet number is limited by water vapor availability and depends on the aerosol particle concentration levels associated with the background.

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.

In the framework of ACTRIS (Aerosols, Clouds, and Trace Gases Research Infrastructure Network) summer 2012 measurement campaign (8 June–17 July 2012), EARLINET organized and performed a controlled exercise of feasibility to demonstrate its potential to perform operational, coordinated measurements and deliver products in near-real time. Eleven lidar stations participated in the exercise which started on 9 July 2012 at 06:00 UT and ended 72 h later on 12 July at 06:00 UT. For the first time, the single calculus chain (SCC) – the common calculus chain developed within EARLINET for the automatic evaluation of lidar data from raw signals up to the final products – was used. All stations sent in real-time measurements of a 1 h duration to the SCC server in a predefined netcdf file format. The pre-processing of the data was performed in real time by the SCC, while the optical processing was performed in near-real time after the exercise ended. 98 and 79 % of the files sent to SCC were successfully pre-processed and processed, respectively. Those percentages are quite large taking into account that no cloud screening was performed on the lidar data. The paper draws present and future SCC users' attention to the most critical parameters of the SCC product configuration and their possible optimal value but also to the limitations inherent to the raw data. The continuous use of SCC direct and derived products in heterogeneous conditions is used to demonstrate two potential applications of EARLINET infrastructure: the monitoring of a Saharan dust intrusion event and the evaluation of two dust transport models. The efforts made to define the measurements protocol and to configure properly the SCC pave the way for applying this protocol for specific applications such as the monitoring of special events, atmospheric modeling, climate research and calibration/validation activities of spaceborne observations.

Systematic measurements of dust concentration profiles at a continental scale were recently made possible by the development of synergistic retrieval algorithms using combined lidar and sun photometer data and the establishment of robust remote-sensing networks in the framework of Aerosols, Clouds, and Trace gases Research InfraStructure Network (ACTRIS)/European Aerosol Research Lidar Network (EARLINET). We present a methodology for using these capabilities as a tool for examining the performance of dust transport models. The methodology includes considerations for the selection of a suitable data set and appropriate metrics for the exploration of the results. The approach is demonstrated for four regional dust transport models (BSC-DREAM8b v2, NMMB/BSC-DUST, DREAMABOL, DREAM8-NMME-MACC) using dust observations performed at 10 ACTRIS/EARLINET stations. The observations, which include coincident multi-wavelength lidar and sun photometer measurements, were processed with the Lidar-Radiometer Inversion Code (LIRIC) to retrieve aerosol concentration profiles. The methodology proposed here shows advantages when compared to traditional evaluation techniques that utilize separately the available measurements such as separating the contribution of dust from other aerosol types on the lidar profiles and avoiding model assumptions related to the conversion of concentration fields to aerosol extinction values. When compared to LIRIC retrievals, the simulated dust vertical structures were found to be in good agreement for all models with correlation values between 0.5 and 0.7 in the 1–6 km range, where most dust is typically observed. The absolute dust concentration was typically underestimated with mean bias values of -40 to -20 μg m−3 at 2 km, the altitude of maximum mean concentration. The reported differences among the models found in this comparison indicate the benefit of the systematic use of the proposed approach in future dust model evaluation studies.

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.

This paper presents a new application of assimilating lidar signals to aerosol forecasting. It aims at investigating the impact of a ground-based lidar network on the analysis and short-term forecasts of aerosols through a case study in the Mediterranean basin. To do so, we employ a data assimilation (DA) algorithm based on the optimal interpolation method developed in the Polair3D chemistry transport model (CTM) of the Polyphemus air quality modelling platform. We assimilate hourly averaged normalised range-corrected lidar signals (PR2) retrieved from a 72 h period of intensive and continuous measurements performed in July 2012 by ground-based lidar systems of the European Aerosol Research Lidar Network (EARLINET) integrated into the Aerosols, Clouds, and Trace gases Research InfraStructure (ACTRIS) network and an additional system in Corsica deployed in the framework of the pre-ChArMEx (Chemistry-Aerosol Mediterranean Experiment)/TRAQA (TRAnsport à longue distance et Qualité de l'Air) campaign. This lidar campaign was dedicated to demonstrating the potential operationality of a research network like EARLINET and the potential usefulness of assimilation of lidar signals to aerosol forecasts. Particles with an aerodynamic diameter lower than 2.5 μm (PM2.5) and those with an aerodynamic diameter higher than 2.5 μm but lower than 10 μm (PM10–2.5) are analysed separately using the lidar observations at each DA step. First, we study the spatial and temporal influences of the assimilation of lidar signals on aerosol forecasting. We conduct sensitivity studies on algorithmic parameters, e.g. the horizontal correlation length (Lh) used in the background error covariance matrix (50 km, 100 km or 200 km), the altitudes at which DA is performed (0.75–3.5 km, 1.0–3.5 km or 1.5–3.5 km a.g.l.) and the assimilation period length (12 h or 24 h). We find that DA with Lh = 100 km and assimilation from 1.0 to 3.5 km a.g.l. during a 12 h assimilation period length leads to the best scores for PM10 and PM2.5 during the forecast period with reference to available measurements from surface networks. Secondly, the aerosol simulation results without and with lidar DA using the optimal parameters (Lh = 100 km, an assimilation altitude range from 1.0 to 3.5 km a.g.l. and a 12 h DA period) are evaluated using the level 2.0 (cloud-screened and quality-assured) aerosol optical depth (AOD) data from AERONET, and mass concentration measurements (PM10 or PM2.5) from the French air quality (BDQA) network and the EMEP-Spain/Portugal network. The results show that the simulation with DA leads to better scores than the one without DA for PM2.5, PM10and AOD. Additionally, the comparison of model results to evaluation data indicates that the temporal impact of assimilating lidar signals is longer than 36 h after the assimilation period.

We examine sunsets painted by famous artists as proxy information for the aerosol optical depth after major volcanic eruptions. Images derived from precision colour protocols applied to the paintings were compared to online images, and found that the latter, previously analysed, provide accurate information. Aerosol optical depths (AODs) at 550 nm, corresponding to Northern Hemisphere middle latitudes, calculated by introducing red-to-green (R / G) ratios from a large number of paintings to a radiative transfer model, were significantly correlated with independent proxies from stratospheric AOD and optical extinction data, the dust veil index, and ice core volcanic indices. AODs calculated from paintings were grouped into 50-year intervals from 1500 to 2000. The year of each eruption and the 3 following years were defined as "volcanic". The remaining "non-volcanic" years were used to provide additional evidence of a multidecadal increase in the atmospheric optical depths during the industrial "revolution". The increase of AOD at 550 nm calculated from the paintings grows from 0.15 in the middle 19th century to about 0.20 by the end of the 20th century. To corroborate our findings, an experiment was designed in which a master painter/colourist painted successive sunsets during and after the passage of Saharan aerosols over the island of Hydra in Greece. Independent solar radiometric measurements confirmed that the master colourist's R / G ratios which were used to model his AODs, matched the AOD values measured in situ by co-located sun photometers during the declining phase of the Saharan aerosol. An independent experiment was performed to understand the difference between R / G ratios calculated from a typical volcanic aerosol and those measured from the mineral aerosol during the Hydra experiment. It was found that the differences in terms of R / G ratios were small, ranging between −2.6% and +1.6%. Also, when analysing different parts of cloudless skies of paintings following major volcanic eruptions, any structural differences seen in the paintings had not altered the results discussed above. However, a detailed study on all possible sources of uncertainties involved (such as the impact of clouds on R / G ratios) still needs to be studied. Because of the large number of paintings studied, we tentatively propose the conclusion that regardless of the school, red-to-green ratios from great masters can provide independent proxy AODs that correlate with widely accepted proxies and with independent measurements.

Vertical profiles of the optical (extinction and backscatter coefficients, lidar ratio and Ångström exponent), microphysical (mean effective radius, mean refractive index, mean number concentration) and geometrical properties as well as the mass concentration of volcanic particles from the Eyjafjallajökull eruption were retrieved at selected heights over Athens, Greece, using multi-wavelength Raman lidar measurements performed during the period 21–24 April 2010. Aerosol Robotic Network (AERONET) particulate columnar measurements along with inversion schemes were initialized together with lidar observations to deliver the aforementioned products. The well-known FLEXPART (FLEXible PARTicle dispersion model) model used for volcanic dispersion simulations is initiated as well in order to estimate the horizontal and vertical distribution of volcanic particles. Compared with the lidar measurements within the planetary boundary layer over Athens, FLEXPART proved to be a useful tool for determining the state of mixing of ash with other, locally emitted aerosol types. The major findings presented in our work concern the identification of volcanic particles layers in the form of filaments after 7-day transport from the volcanic source (approximately 4000 km away from our site) from the surface and up to 10 km according to the lidar measurements. Mean hourly averaged lidar signals indicated that the layer thickness of volcanic particles ranged between 1.5 and 2.2 km. The corresponding aerosol optical depth was found to vary from 0.01 to 0.18 at 355 nm and from 0.02 up to 0.17 at 532 nm. Furthermore, the corresponding lidar ratios (S) ranged between 60 and 80 sr at 355 nm and 44 and 88 sr at 532 nm. The mean effective radius of the volcanic particles estimated by applying inversion scheme to the lidar data found to vary within the range 0.13–0.38 μm and the refractive index ranged from 1.39+0.009i to 1.48+0.006i. This high variability is most probably attributed to the mixing of aged volcanic particles with other aerosol types of local origin. Finally, the LIRIC (LIdar/Radiometer Inversion Code) lidar/sunphotometric combined inversion algorithm has been applied in order to retrieve particle concentrations. These have been compared with FLEXPART simulations of the vertical distribution of ash showing good agreement concerning not only the geometrical properties of the volcanic particles layers but also the particles mass concentration.

The eruption of the Icelandic volcano Eyjafjallajökull in April–May 2010 represents a "natural experiment" to study the impact of volcanic emissions on a continental scale. For the first time, quantitative data about the presence, altitude, and layering of the volcanic cloud, in conjunction with optical information, are available for most parts of Europe derived from the observations by the European Aerosol Research Lidar NETwork (EARLINET). Based on multi-wavelength Raman lidar systems, EARLINET is the only instrument worldwide that is able to provide dense time series of high-quality optical data to be used for aerosol typing and for the retrieval of particle microphysical properties as a function of altitude. In this work we show the four-dimensional (4-D) distribution of the Eyjafjallajökull volcanic cloud in the troposphere over Europe as observed by EARLINET during the entire volcanic event (15 April–26 May 2010). All optical properties directly measured (backscatter, extinction, and particle linear depolarization ratio) are stored in the EARLINET database available at http://www.earlinet.org. A specific relational database providing the volcanic mask over Europe, realized ad hoc for this specific event, has been developed and is available on request at http://www.earlinet.org. During the first days after the eruption, volcanic particles were detected over Central Europe within a wide range of altitudes, from the upper troposphere down to the local planetary boundary layer (PBL). After 19 April 2010, volcanic particles were detected over southern and south-eastern Europe. During the first half of May (5–15 May), material emitted by the Eyjafjallajökull volcano was detected over Spain and Portugal and then over the Mediterranean and the Balkans. The last observations of the event were recorded until 25 May in Central Europe and in the Eastern Mediterranean area. The 4-D distribution of volcanic aerosol layering and optical properties on European scale reported here provides an unprecedented data set for evaluating satellite data and aerosol dispersion models for this kind of volcanic events.

A novel procedure has been developed to retrieve, simultaneously, the optical, microphysical and chemical properties of tropospheric aerosols with a multi-wavelength Raman lidar system in the troposphere over an urban site (Athens, Greece: 37.9° N, 23.6° E, 200 m a.s.l.) using data obtained during the European Space Agency (ESA) THERMOPOLIS project, which took place between 15–31 July 2009 over the Greater Athens Area (GAA). We selected to apply our procedure for a case study of intense aerosol layers that occurred on 20–21 July 2009. The National Technical University of Athens (NTUA) EOLE 6-wavelength Raman lidar system has been used to provide the vertical profiles of the optical properties of aerosols (extinction and backscatter coefficients, lidar ratio) and the water vapor mixing ratio. An inversion algorithm was used to derive the mean aerosol microphysical properties (mean effective radius (reff), single-scattering albedo ω) and mean complex refractive index (m)) at selected heights in the 2–3 km height region. We found that reff was 0.14–0.4 (±0.14) μm, ω was 0.63–0.88 (±0.08) (at 532 nm) and m ranged from 1.44 (±0.10) + 0.01 (±0.01)i to 1.55 (±0.12) + 0.06 (±0.02)i, in good agreement (only for the reff values) with in situ aircraft measurements. The water vapor and temperature profiles were incorporated into the ISORROPIA II model to propose a possible in situ aerosol composition consistent with the retrieved m and ω values. The retrieved aerosol chemical composition in the 2–3 km height region gave a variable range of sulfate (0–60%) and organic carbon (OC) content (0–50%), although the OC content increased (up to 50%) and the sulfate content dropped (up to 30%) around 3 km height; the retrieved low ω value (0.63), indicates the presence of absorbing biomass burning smoke mixed with urban haze. Finally, the retrieved aerosol microphysical properties were compared with column-integrated sun photometer CIMEL data.

A strong Saharan dust event that occurred over the city of Athens, Greece (37.9° N, 23.6° E) between 27 March and 3 April 2009 was followed by a synergy of three instruments: a 6-wavelength Raman lidar, a CIMEL sun-sky radiometer and the MODIS sensor. The BSC-DREAM model was used to forecast the dust event and to simulate the vertical profiles of the aerosol concentration. Due to mixture of dust particles with low clouds during most of the reported period, the dust event could be followed by the lidar only during the cloud-free day of 2 April 2009. The lidar data obtained were used to retrieve the vertical profile of the optical (extinction and backscatter coefficients) properties of aerosols in the troposphere. The aerosol optical depth (AOD) values derived from the CIMEL ranged from 0.33–0.91 (355 nm) to 0.18–0.60 (532 nm), while the lidar ratio (LR) values retrieved from the Raman lidar ranged within 75–100 sr (355 nm) and 45–75 sr (532 nm). Inside a selected dust layer region, between 1.8 and 3.5 km height, mean LR values were 83 ± 7 and 54 ± 7 sr, at 355 and 532 nm, respectively, while the Ångström-backscatter-related (ABR355/532) and Ångström-extinction-related (AER355/532) were found larger than 1 (1.17 ± 0.08 and 1.11 ± 0.02, respectively), indicating mixing of dust with other particles. Additionally, a retrieval technique representing dust as a mixture of spheres and spheroids was used to derive the mean aerosol microphysical properties (mean and effective radius, number, surface and volume density, and mean refractive index) inside the selected atmospheric layers. Thus, the mean value of the retrieved refractive index was found to be 1.49( ± 0.10) + 0.007( ± 0.007)i, and that of the effective radiuses was 0.30 ± 0.18 μm. The final data set of the aerosol optical and microphysical properties along with the water vapor profiles obtained by Raman lidar were incorporated into the ISORROPIA II model to provide a possible aerosol composition consistent with the retrieved refractive index values. Thus, the inferred chemical properties showed 12–40% of dust content, sulfate composition of 16–60%, and organic carbon content of 15–64%, indicating a possible mixing of dust with haze and smoke. PM10 concentrations levels, PM10 composition results and SEM-EDX (Scanning Electron Microscope-Energy Dispersive X-ray) analysis results on sizes and mineralogy of particles from samples during the Saharan dust transport event were used to evaluate the retrieval.

This study presents an inter-comparison of two active remote sensors (lidar and ceilometer) to determine the mixing layer height and structure of the Planetary Boundary Layer (PBL) and to retrieve tropospheric aerosol vertical profiles over Athens, Greece. This inter-comparison was performed under various strongly different aerosol loads/types (urban air pollution, biomass burning and Saharan dust event), implementing two different lidar systems (one portable Raymetrics S.A. lidar system running at 355 nm and one multi-wavelength Raman lidar system running at 355 nm, 532 nm and 1064 nm) and one CL31 Vaisala S.A. ceilometer (running at 910 nm). Spectral conversions of the ceilometer's data were performed using the Ångström exponent estimated by ultraviolet multi-filter radiometer (UV-MFR) measurements. The inter-comparison was based on two parameters: the mixing layer height determined by the presence of the suspended aerosols and the attenuated backscatter coefficient. Additionally, radiosonde data were used to derive the PBL height. In general, a good agreement was found between the ceilometer and the lidar techniques in both inter-compared parameters in the height range from 500 m to 5000 m, while the limitations of each instrument are also examined.

In this paper we present a statistical analysis on the geometrical and optical properties of Saharan dust layers observed over Athens, Greece, in a three-year period from 1 January 2004 up to 31 December 2006. The observations of the vertical aerosol profile were performed by the multi-wavelength (355-532-1064-387-607 nm) Raman lidar system of the National Technical University of Athens (NTUA) operated in the city of Athens (37°98' N, 23°77' E), Greece, in the frame of the European Aerosol Research Lidar Network (EARLINET-ASOS) project. The number of dust events was greatest in late spring, summer, and early autumn periods. This was evident also by aerosol observations during dust outbreaks obtained from the Moderate Resolution Imaging Spectroradiometer (MODIS). In our lidar measurements, multiple aerosol dust layers of variable thickness (680–4800 m) were observed. The center of mass of these layers was located in altitudes between 1600 and 5800 m. However, the mean thickness of the dust layer typically stayed around 2700 m and the corresponding mean center of mass was of the order of 2900 m. The top of the dust layer ranged from 2000 to 8000 m, with a mean value of the order of 4700 m. MODIS observations during dust outbreaks showed that the AOD values at 550 nm ranged between 0.3–0.6, while the corresponding Angström exponent (AE) values were of the order of 0.5–0.65, indicating the presence of rather large particles.

We present initial aerosol validation results of the space-borne lidar CALIOP -onboard the CALIPSO satellite- Level 1 attenuated backscatter coefficient profiles, using coincident observations performed with a ground-based lidar in Athens, Greece (37.9° N, 23.6° E). A multi-wavelength ground-based backscatter/Raman lidar system is operating since 2000 at the National Technical University of Athens (NTUA) in the framework of the European Aerosol Research LIdar NETwork (EARLINET), the first lidar network for tropospheric aerosol studies on a continental scale. Since July 2006, a total of 40 coincidental aerosol ground-based lidar measurements were performed over Athens during CALIPSO overpasses. The ground-based measurements were performed each time CALIPSO overpasses the station location within a maximum distance of 100 km. The duration of the ground–based lidar measurements was approximately two hours, centred on the satellite overpass time. From the analysis of the ground-based/satellite correlative lidar measurements, a mean bias of the order of 22% for daytime measurements and of 8% for nighttime measurements with respect to the CALIPSO profiles was found for altitudes between 3 and 10 km. The mean bias becomes much larger for altitudes lower that 3 km (of the order of 60%) which is attributed to the increase of aerosol horizontal inhomogeneity within the Planetary Boundary Layer, resulting to the observation of possibly different air masses by the two instruments. In cases of aerosol layers underlying Cirrus clouds, comparison results for aerosol tropospheric profiles become worse. This is attributed to the significant multiple scattering effects in Cirrus clouds experienced by CALIPSO which result in an attenuation which is less than that measured by the ground-based lidar.

A long-lasting Saharan dust event affected Europe on 18–23 May 2008. Dust was present in the free troposphere over Greece, in height ranges between the surface and approximately 4–5 km above sea level. The event was monitored by ground-based CIMEL sunphotometric and multi-wavelength combined backscatter/Raman lidar measurements over Athens, Greece. The dust event had the maximum of its intensity on 20 May. Three-dimensional dust spatial distribution over Greece on that day is presented through satellite synergy of passive and active remote sensing using MODIS and CALIPSO data, respectively. For the period under study, the ground-based measurements are used to characterize the dust event and evaluate the latest version of the BSC Dust Regional Atmospheric Modeling (BSC-DREAM) system. Comparisons of modeled and measured aerosol optical depths over Athens show that the Saharan dust outbreak is fairly well captured by BSC-DREAM simulations. Evaluation of BSC-DREAM using Raman lidar measurements on 20 May shows that the model consistently reproduces the dust vertical distribution over Athens.

An investigation of the Planetary Boundary Layer (PBL) height evolution over Greece, during the solar eclipse of 29 March 2006, is presented. Ground based observations were carried out using lidar detection and ranging devices and ground meteorological instruments, to estimate the height of the mixing layer (ML) before, during and after the solar eclipse in northern and southern parts of Greece exhibiting different sun obscuration. Data demonstrate that the solar eclipse has induced a decrease of the PBL height, indicating a suppression of turbulence activity similar to that during the sunset hours. The changes in PBL height were associated with a very shallow entrainment zone, indicating a significant weakening of the penetrative convection. Heat transfer was confined to a thinner layer above the ground. The thickness of the entrainment zone exhibited its minimum during the maximum of the eclipse, demonstrative of turbulence mechanisms suppression at that time. Model estimations of the PBL evolution were additionally conducted using the Comprehensive Air Quality Model with extensions (CAMx) coupled with the Weather Research and Forecasting model (WRF). Model-diagnosed PBL height decrease during the solar eclipse due to vertical transport decay, in agreement with the experimental findings; vertical profiles of atmospheric particles and gaseous species showed an important vertical mixing attenuation.

High aerosol optical depth (AOD) values, larger than 0.6, are systematically observed in the Ultraviolet (UV) region both by sunphotometers and lidar systems over Greece during summertime. To study in more detail the characteristics and the origin of these high AOD values, a campaign took place in Greece in the frame of the PHOENICS (Particles of Human Origin Extinguishing Natural solar radiation In Climate Systems) and EARLINET (European Aerosol Lidar Network) projects during August–September of 2003, which included simultaneous sunphotometric and lidar measurements at three sites covering the north-south axis of Greece: Thessaloniki, Athens and Finokalia, Crete. Several events with high AOD values have been observed over the measuring sites during the campaign period, many of them corresponding to Saharan dust. In this paper we focused on the event of 30 and 31 August 2003, when a dust layer in the height range of 2000-5000 m, progressively affected all three stations. This layer showed a complex behavior concerning its spatial evolution and allowed us to study the changes in the optical properties of the desert dust particles along their transport due to aging and mixing with other types of aerosol. The extinction-to-backscatter ratio determined on the 30 August 2003 at Thessaloniki was approximately 50 sr, characteristic for rather spherical mineral particles, and the measured color index of 0.4 was within the typical range of values for desert dust. Mixing of the desert dust with other sources of aerosols resulted the next day in overall smaller and less absorbing population of particles with a lidar ratio of 20 sr. Mixing of polluted air-masses originating from Northern Greece and Crete and Saharan dust result in very high aerosol backscatter values reaching 7 Mm-1 sr-1 over Finokalia. The Saharan dust observed over Athens followed a different spatial evolution and was not mixed with the boundary layer aerosols mainly originating from local pollution.

A two-wavelength ultraviolet (289–316nm) ozone Differential Absorption Lidar (DIAL) system is used to perform ozone measurements in the free troposphere in the Eastern Mediterranean (Northern Greece). The ozone DIAL profiles obtained during a Stratosphere-to-Troposphere Transport (STT) event are compared to that acquired by an electrochemical ozonesonde, in the altitude range between 2 and 10 km. The measurement accuracy of these two instruments is also discussed. The mean difference between the ozone profiles obtained by the two techniques is of the order of 1.11 ppbv (1.86%), while the corresponding standard deviation is 4.69 ppbv (8.16%). A case study of an STT event which occurred on 29 November 2000 is presented and analyzed, using ozone lidar, satellite and meteorological data, as well as air mass back-trajectory analysis. During this STT event ozone mixing ratios of 55–65 ppbv were observed between 5 and 7 km height above sea level (a.s.l.). Stratospheric air was mixed with tropospheric air masses, leading to potential vorticity (PV) losses due to diabatic processes. The ozone DIAL system can be used for following STT events and small-scale mixing phenomena in the free troposphere, and for providing sequences of vertical ozone profiles in the free troposphere. Keywords. Atmospheric composition and structure (Evolution of the atmosphere; Instruments and techniques) – Meteorology and atmospheric dynamics (Middle atmosphere dynamics; Turbulence)

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.

Routine lidar measurements of the vertical distribution of the aerosol extinction coefficient and the extinction-to-backscatter ratio have been performed at Thessaloniki, Greece using a Raman lidar system in the frame of the EARLINET project since 2000. Co-located spectral and broadband solar UV-B irradiance measurements, as well as total ozone observations, were available whenever lidar measurements were obtained. From the available measurements several cases could be identified that allowed the study of the effect of different types of aerosol on the levels of the UV-B solar irradiance at the Earth's surface. The TUV radiative transfer model has been used to simulate the irradiance measurements, using total ozone and the lidar aerosol data as input. From the comparison of the model results with the measured spectra the effective single scattering albedo was determined using an iterative procedure, which has been verified against results from the 1998 Lindenberg Aerosol Characterization Experiment. It is shown that for the same aerosol optical depth and for the same total ozone values the UV-B irradiances at the Earth's surface can show differences up to 10%, which can be attributed to differences in the aerosol type. It is shown that the combined use of the estimated single scattering albedo and of the measured extinction-to-backscatter ratio leads to a better characterization of the aerosol type probed.

A wide range of measurements was carried out in central and southeastern Europe within the framework of the EU project STACCATO (Influence of Stratosphere-Troposphere Exchange in a Changing Climate on Atmospheric Transport and Oxidation Capacity) with the principle goal to create a comprehensive data set on stratospheric air intrusions into the troposphere along a rather frequently observed pathway over central Europe from the North Sea to the Mediterranean Sea. The measurements were based on predictions by suitable quasi-operational trajectory calculations using ECMWF forecast data. A predicted deep Stratosphere to Troposphere Transport (STT) event, encountered during the STACCATO period on 20-21 June 2001, was followed by the measurements network almost from its inception. Observations provide evidence that the intrusion affected large parts of central and southeastern Europe. Especially, the ozone lidar observations on 20-21 June 2001 at Garmisch-Partenkirchen, Germany captured the evolution of two marked tongues of high ozone with the first one descending to nearly 2 km, thus providing an excellent data set for model intercomparisons and validation. In addition, for the first time to our knowledge concurrent surface measurements of the cosmogenic radionuclides 10Be and 7Be and their ratio 10Be/7Be are presented together as stratospheric tracers in a case study of a stratospheric intrusion. The ozone tracer columns calculated with the FLEXPART model were found to be in good agreement with water vapour satellite images, capturing the evolution of the observed dry streamers of stratospheric origin. Furthermore, the time-height cross section of ozone tracer simulated with FLEXPART over Garmisch-Partenkirchen captures many details of the evolution of the two observed high-ozone filaments measured with the IFU lidar, thus demonstrating the considerable progress in model simulations. Finally, the modelled ozone (operationally available since October 1999) from the ECMWF (European Centre for Medium-Range Weather Forecasts) atmospheric model is shown to be in very good agreement with the observations during this case study, which provides the first successful validation of a chemical tracer that is derived operationally from a weather forecast model. This suggests that coupling chemistry and weather forecast models may significantly improve both weather and chemical forecasts in the future.

A backscattering LIDAR system, the first of this kind in Brazil, has been set-up in a suburban area in the city of São Paulo (23º33' S, 46º44' W) to provide the vertical profile of the aerosol backscatter coefficient at 532 nm up to an altitude of 4-6 km above sea level (asl). The measurements have been carried out during the second half of the so-called Brazilian dry season, September and October 2001 and during the first half of the dry season in August and September 2002. The LIDAR data are presented and analysed in synergy with aerosol optical thickness (AOT) measurements obtained by a CIMEL sun-tracking photometer in the visible spectral region and with satellite measurements obtained by the MODIS sensor. This synergetic approach has been used, not only to validate the LIDAR data, but also to derive a typical value (45 sr) of the so-called extinction-to-backscatter ratio (LIDAR ratio) during the dry season. The satellite data analysis offers additional information on the spatial distribution of aerosols over Brazil including the determination of aerosol source regions over the country. The LIDAR data were also used to retrieve the Planetary Boundary Layer (PBL) height, aerosol layering and the structure of the lower troposphere over the city of São Paulo. These first LIDAR measurements over the city of São Paulo during the dry season showed a significant variability of the AOT in the lower troposphere (0.5-5 km) at 532 nm. It was also found that the aerosol load is maximized in the 1-3 km height region, although up to 3 km thick aerosol layers were also detected in the 2.5-5.5 km region in certain cases. Three-dimensional 96-hours air mass back-trajectory analysis was also performed in selected cases to determine the source regions of aerosols around São Paulo during the dry season.

We provide an overview of light detection and ranging (lidar) capability for describing and characterizing desert dust. This paper summarizes lidar techniques, observations, and fallouts of desert dust lidar measurements. The main objective is to provide the scientific community, including nonpractitioners of lidar observations with a reference paper on dust lidar measurements. In particular, it will fill the current gap of communication between research-oriented lidar community and potential desert dust data users, such as air quality monitoring agencies and aviation advisory centers. The current capability of the different lidar techniques for the characterization of aerosol in general and desert dust in particular is presented. Technical aspects and required assumptions of these techniques are discussed, providing readers with the pros and cons of each technique. Information about desert dust collected up to date using lidar techniques is reviewed. Lidar techniques for aerosol characterization have a maturity level appropriate for addressing air quality and transportation issues, as demonstrated by some first results reported in this paper.