The volcanic eruption of Cumbre Vieja (La Palma Island, Spain), started on 19 September 2021 and was declared terminated on 25 December 2021. A complete set of aerosol measurements were deployed around the volcano within the first month of the eruptive activity. This paper describes the results of the observations made at Tazacorte on the west bank of the island where a polarized micro-pulse lidar was deployed. The analyzed two-and-a-half months (16 October–31 December) reveal that the peak height of the lowermost and strongest volcanic plume did not exceed 3 km (the mean of the hourly values is 1.43 ± 0.45 km over the whole period) and was highly variable. The peak height of the lowermost volcanic plume steadily increased until week 11 after the eruption started (and 3 weeks before its end) and started decreasing afterward. The ash mass concentration was assessed with a method based on the polarization capability of the instrument. Two days with a high ash load were selected: The ash backscatter coefficient, aerosol optical depth, and the volume and particle depolarization ratios were, respectively, 3.6 (2.4) Mm⁻¹sr⁻¹, 0.52 (0.19), 0.13 (0.07) and 0.23 (0.13) on 18 October (15 November). Considering the limitation of current remote sensing techniques to detect large-to-giant particles, the ash mass concentration on the day with the highest ash load (18 October) was estimated to have peaked in the range of 800–3200 μg m⁻³ in the lowermost layer below 2.5 km.

In this paper a statistical study of cirrus geometrical and optical properties based on 4 years of continuous ground-based lidar measurements with the Barcelona (Spain) Micro Pulse Lidar (MPL) is analysed. First, a review of the literature on the two-way transmittance method is presented. This method is a well-known lidar inversion method used to retrieve the optical properties of an aerosol–cloud layer between two molecular (i.e. aerosol and cloud-free) regions below and above, without the need to make any a priori assumptions about their optical and/or microphysical properties. Second, a simple mathematical expression of the two-way transmittance method is proposed for both ground-based and spaceborne lidar systems. This approach of the method allows the retrieval of the cloud optical depth, the cloud column lidar ratio and the vertical profile of the cloud backscatter coefficient. The method is illustrated for a cirrus cloud using measurements from the ground-based MPL and from the spaceborne Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP). Third, the database is then filtered with a cirrus identification criterion based on (and compared to) the literature using only lidar and radiosonde data. During the period from November 2018 to September 2022, 367 high-altitude cirrus clouds were identified at 00:00 and 12:00 UTC, of which 203 were successfully inverted with the two-way transmittance method. The statistical results of these 203 high-altitude cirrus clouds show that the cloud thickness is 1.8 ± 1.1 km, the mid-cloud temperature is −51 ± 8 ^∘C and the linear cloud depolarization ratio is 0.32 ± 0.13. The application of the transmittance method yields an average cloud optical depth (COD) of 0.36 ± 0.45 and a mean effective column lidar ratio of 30 ± 19 sr. Statistical results of the errors associated with the two-way transmittance method retrievals are also provided. The highest occurrence of cirrus is observed in spring and the majority of cirrus clouds (48 %) are visible (0.03 < COD < 0.3), followed by opaque (COD > 0.3) with a percentage of 38 %. Together with results from other sites, possible latitudinal dependencies have been analysed together with correlations between cirrus cloud properties. For example, we noted that in Barcelona the COD correlates positively with the cloud base temperature, effective column lidar ratio and linear cloud depolarization ratio and negatively with the cloud base height.

Cloud–radiation interaction still drives large uncertainties in climate models, and its estimation is key to make more accurate predictions. In this context, high-altitude cirrus clouds play a fundamental role because (1) they have a high occurrence frequency globally and (2) they are the only cloud that can readily cool or warm the atmosphere during daytime, depending on their properties. This study presents a comprehensive analysis of optical scattering properties and direct radiative effect of cirrus clouds based on 4 years of continuous ground-based lidar measurements with the Barcelona (Spain) Micro Pulse Lidar. First, we introduce a novel approach of a self-consistent scattering model for cirrus clouds to determine their optical scattering properties at different wavelengths using only the extinction coefficient and cloud temperature. Second, we calculate the direct radiative effects of cirrus clouds with the discrete ordinates method, and we validate our results with SolRad-Net pyranometers and NOAA-20 measurements. Third, we present a case study analysing the direct radiative effect of a cirrus cloud along its back-trajectory using data from the Chemical LAgrangian Model of the Stratosphere with microphysics scheme for Ice clouds formation (CLaMS-Ice). The results show that the cirrus clouds with an average ice water content of 4.97 ± 5.53 mg m−3, in the nighttime, have a positive direct radiative effect at the top of the atmosphere (TOA; +40.4 W m−2) almost twice than at the bottom of the atmosphere (BOA; +22.1 W m−2); in the daytime, they have generally a negative direct radiative effect at BOA (−11.5 W m−2, 82 % of the cases) and always a positive effect at TOA (+14.2 W m−2). In these simulations, the influence of the lower layer aerosols is negligible in the cirrus direct radiative effects, with a bias (denoted BIAS) of −1.2 %. For the case study, the net direct radiative effects produced by the cirrus cloud, are from 0 to +40 W m−2 at TOA and from −51 to +20 W m−2 at BOA. This study reveals that the complexity of the cirrus cloud direct radiative effect calculation lies in the fact that it is highly sensitive to the cirrus scene properties.

The versatile Generalized Retrieval of Aerosol and Surface Properties (GRASP) algorithm exploits the advantages of synergic ground-based aerosol observations such as radiometric (sensitive to columnar aerosol optical and microphysical properties) and lidar (sensitive to vertical distribution of the optical properties) observations. The synergy is possible when the complementary data is mutually constrained by GRASP parametrization that includes, for the first time ever, the degree of linear polarization (DoLP) parameter measured by a polarized sun-sky-lunar AERONET photometer (380, 440, 500, 675, 870, 1020, and 1640 nm) in synergy with the vertical profiles from an elastic lidar (355, 532, and 1064 nm). First, a series of numerical tests is performed using simulated data generated using a climatology of data and ground-based measurements. The inversions are performed with and without random noise for five different combinations of input data, starting from the AERONET-like dataset and increasing to the complex one by adding more information for three aerosol scenarios: I—high aerosol optical depth (AOD) with dominant coarse mode; II—low AOD with dominant coarse mode; III—high AOD with dominant fine mode. The inclusion of DoLP improves (i) the retrieval accuracy of the fine-mode properties when it is not dominant; (ii) the retrieval accuracy of the coarse-mode properties at longer wavelengths and that of the fine-mode properties at shorter wavelengths; (iii) the retrieval accuracy of the coarse-mode real part of the refractive index (up to 36% reduction), but has no effect on the retrieval of the imaginary part; (iv) reduces up to 83% the bias of the sphere fraction (SF) retrieval in coarse-mode dominated regimes; and (v) the root mean square error (RMSE) of the retrieval for most of the parameters in all scenarios. In addition, the addition of more photometer channels in synergy with a three-wavelength elastic lidar reduces the RMSE for the real part (67% in the coarse mode) and the imaginary part (35% in the fine mode) of the refractive index, the single scattering albedo (38% in the fine mode), the lidar ratio (20% in the coarse mode), and the SF (43%).

This research aims to estimate long-term aerosol radiative effects by combining radiation and Aerosol Optical Depth (AOD) observations in Barcelona, Spain. Aerosol Radiative Forcing and Aerosol Forcing Efficiency (ARF and AFE) were estimated by combining shortwave radiation measurements from a SolRad-Net CM-21 pyranometer (level 1.5) and AERONET AOD (level 2), using the direct method. The shortwave AFE was derived from the slope between net solar radiation and AOD at 440, 675, 879, and 1020 nm, and the ARF was computed by multiplying the AFE by AOD at six solar zenith angles (20°, 30°, 40°, 50°, 60°, and 70°). Clear-sky conditions were selected from all-skies days by a quadratic fitting. The aerosol was classified to investigate the forcing contributions from each aerosol type. The aerosol classification was based on Pace and Toledano’s thresholds from AOD vs. Ångström Exponent (AE). The GRASP inversions were performed by combined AOD, radiation, Degree of Linear Polarization (DoLP) by zenith angles from the polarized sun–sky–lunar photometer and the elastic signal from the UPC-ACTRIS lidar system. The long-term AFE and ARF are both negative, with an increasing tendency (in absolute value) of +24% (AFE) and +40% (ARF) in 14 years. The yearly AFE varied from −331 to −10 Wm⁻²τ⁻¹, and the ARF varied from −64 to −2 Wm⁻², associated with an AOD (440 nm) from 0.016 to 0.690. The three types of aerosols on clear-sky days are mixed aerosols (61%), desert dust (10%), and urban/industrial-biomass burning aerosols (29%). Combined with Gobbi’s method, this classification clustered the aerosols into four groups by AE analysis (two coarse- and two fine-mode aerosols). Then, the contribution of the aerosol types to the ARF showed that the desert dust forcing had the largest cooling effect in Barcelona (−61.5 to −37.4 Wm⁻²), followed by urban/industrial-biomass burning aerosols (−40.4 to −20.4 Wm⁻²) and mixed aerosols (−31.8 and −24.0 Wm⁻²). Regarding the comparison among Generalized Retrieval of Atmosphere and Surface Properties (GRASP) inversions, AERONET inversions, and direct method estimations, the AFE and ARF had some differences owing to their definitions in the algorithms. The DoLP, used as GRASP input, decreased the ARF overestimation for high AOD.