Small-scale volcanic aerosols variability, processes and direct radiative impact at Mount Etna during the EPL-RADIO campaigns

Sep 16, 2020

Summit observations: emissions

A first target of EPL-RADIO is the proximal aerosol plume characterisation, i.e. at Etna summit. Figure 1a,b shows (hbox {PM}_1) and (hbox {PM}_{10}) (particulate matter with diameter less than 1 and (10,upmu hbox {m})) concentrations measured with the OPC at Etna summit (BNC and VOC craters) on 20/07/2017. Etna’s plume was directed towards the south-east direction. Two plume traverses were realised on that particular day. The (hbox {PM}_1) concentrations were around (100,upmu hbox {g/m}^3) for most of the traverse and reached up to (600,upmu hbox {g/m}^3) in what can be identified as the approximate center of the summit plume. The (hbox {PM}_{10}) concentration reached up to 8,000(,upmu hbox {g/m}^3) in a limited-sized location. Results from chamber testing have shown that these low-cost sensors can underestimate (hbox {PM}_1) mass concentration and overestimate (hbox {PM}_{10}) at these extreme values (see “Methods” section). However, the observed concentrations are of the same order of magnitude as for previous measurements at Etna’s summit20,21. At similar passive degassing conditions, Allen et al.20 observed mass concentrations between (sim 😉 1,400 and (200,upmu hbox {g/m}^3), in the coarse mode (diameter (hbox {D}_{p} > 3.5,upmu hbox {m})), and between (sim 150) and 4,500(,upmu hbox {g/m}^3), in the fine mode ((hbox {D}_{p} < 3.5,upmu hbox {m})), at varying distances from the near-summit to (sim ;10,hbox {km}) distance. Similarly, very high concentrations have been observed for other volcanoes, like the Popocatepetl volcano (total suspended particles (sim ;1,440,upmu hbox {g/m}^3))22 or the Soufrière Hills ((hbox {PM}_{10} = 100)(500,upmu hbox {g/m}^3))23. Our observations confirm (and exceed) extremely large PM values for Etna’s summit and for similar volcanoes. These extremely big PM sources (including the most dangerous small-sized particle matter) have been associated in the past to degraded air quality and connected pulmonary diseases of the neighbouring people22,23,24. These health risks are still to be systematically assessed for Etna area, including for the persistent passive degassing activity.

One can note that the (hbox {PM}_{1}) and (hbox {PM}_{10}) hot-spots are not localized exactly in the same plume section. Therefore, the sources of fine and coarse particles may be different. Moreover, the column integrated (hbox {SO}_{2}) concentrations were simultaneously measured along the same path; time series of (hbox {PM}_{1}) and (hbox {SO}_{2}) concentrations are shown in Fig. 1c,d. The (hbox {SO}_2) measurements show very limited correlations with both (hbox {PM}_1) and (hbox {PM}_{10}) concentrations (see a scatterplot of (hbox {PM}_1) and (hbox {SO}_2) concentrations in Fig. 1e). This evidence contrasts with previous evidences of a strong correlation of total particle column and (hbox {SO}_2) et Etna near-summit15. Even if during this measurement session the volcanic plume was mostly flattened to the ground due to the moderate wind at the summit, it has to be mentioned that part of the plume can be missed from the PM observations because it is located at higher altitudes. While the wind was not extreme, during the session, it is still possible that some of the coarser aerosol observations are moderately affected by surface particles re-mobilisation. Due to these limitations and the lack of complementary observations of atmospheric oxidants or ultra-fine particles, the precise origin of the observed aerosol is not possible. Nevertheless, these results do not show clear evidence of formation of secondary aerosols. Primary aerosol emissions were found dominating, with respect to secondary aerosol formation, at proximal locations, for other passive degassing plumes like for Masaya volcano in the past25.

Distal and itinerant observations: small-scale variability

Getting new insights into the volcanic aerosol plume at larger proximal scale, i.e. under (sim 20,hbox {km}) from the degassing craters, is a further target of EPL-RADIO. Figure 2a shows a summary map of the MIIOM observations of the volcanic aerosols plume’s optical properties during C2 campaign. Aerosol optical depths (AOD) for Etna’s plume were derived at 320.5 nm (hereafter referred to as UV AOD) and 1,020.0 nm (hereafter referred to as NIR AOD), with the method described by Sellitto et al.19. Using UV and NIR AODs, the Ångström exponent (AE) and the Ångström (beta) coefficient are subsequently derived19. In general, shorter-wavelength AODs carry more information on the presence of smaller particles than longer-wavelength AODs26. The spectral variability of the AOD (and then of the extinction of radiation due to the presence of the aerosol) can be represented using the empirical Ångström law, as a function of AE and (beta):

begin{aligned} mathrm {AOD}(lambda )=beta lambda ^{-mathrm {AE}} end{aligned}

(1)

Starting from the fundamental theory of particles-radiation interaction (Mie theory26), AE and (beta) can be considered as optical proxies for particles mean size and burden/composition, in the aerosol layer. Bigger AEs are linked to smaller particles, on average, in the sampled plume. While this relationship of AE and average particle size is generally true, precise studies on the average particle size evolution by means of AE observations may be complicated by the possible evolution of the overall particle size distribution, and in particular of the number of its modes. Log-normal size distributions with two20 or three modes15 have been found at near-crater locations for passive degassing volcanoes. On the contrary, mono-modal size distributions, in the accumulation mode (typical mean particles size around 0.1–(0.5,upmu hbox {m})), have been found at distal locations27 . The (beta) parameter is representative of the AOD value at wavelength of (1.0,upmu hbox {m}). For the particle types in volcanic plumes and the typical size distributions discussed above, (beta) is representative of the coarse particles (mainly ash) in the plume. Thus, while its interpretation is more complex than for the AE, in our present study (beta) can be linked to the burden of the coarse aerosol in the layer, bigger (beta)s being representative of bigger burdens. Despite the intrinsic complexity and uncertainty in their interpretation, the two Ångström parameters are widely recognised as diagnostic tools for volcanic plumes morphological and microphysical properties and their inherent evolution processes28. The average values of AE and (beta) are estimated for each measurement session at each location, during C2, as listed in Table 2 and are shown on the map. Thus, Fig. 2a represents the small-scale variability and evolution with distance of the average size and coarse aerosol burden of a volcanic plume. During C2 campaign a clear image of progressively more dominant smaller particles, associated with a smaller burden of coarse particles, is apparent from crater to increasing distances. At the same time, the UV AOD has a different behavior during both campaigns C2 and C3, with slightly increasing values in the distal field, thus indicating a possible slight increase of the finer aerosol content in the plume. This result is consistent with previous airborne remote observations of Spinetti and Buongiorno16, that have shown a similar increasing trend, for increasing distance from the crater, of the AE and a near-constant to slightly increasing trend of the distal (starting from a few km downwind the craters) short-wavelength AOD (at about 400 nm), during passive degassing activity of Mount Etna volcano. A substantial variability of the plume properties is found at spatial scales that are smaller than typical grid-points of chemistry/transport models and pixels of satellite observations. This reveals that processes are at play which are not represented with current regional scale modelling and observations. To get further insights into the small-scale aerosol variability, in Fig. 2b,c we show the MIIOM AE (C2 campaign) and MIISP AE (C3 campaign), as a function of the distance from craters and of simultaneous measurements of (beta). The error bars in Fig. 2b,c represent the standard deviation of the average AE and (beta), thus representing their variability for each measurement session/location. Results shown in Fig. 2b,c confirm the picture drawn before for C2 campaign: progressively smaller particles and decreasing coarse aerosol burden are found along plume dispersion in the first (sim)10 km around Etna’s degassing craters. The correlation of the occurrence of smaller particles with increasing distance and with thinner plumes is significant (R(^2) correlation coefficients 0.90 and 0.71, respectively). This can be attributed to relatively quick sedimentation of coarser ash particles (bigger than a few micrometers) emitted by discontinuous ash puffs (observed during C3 campaign and on 20/07/2016 during C2 campaign) and, possibly, by sustained gas-to-particle conversion of (hbox {SO}_2) volcanic emissions to secondary SA and subsequent grow by condensation (up to a few hundreds nanometers). These two processes are expected to concur to the progressive modification of the size distribution towards smaller average sizes. During C2 campaign, volcanic activity was based on a prevalent passive degassing and (hbox {SO}_2) emission rates were significant yet very variable (average values of 2,240.5 ± 882.5 t/day, observed by the FLAME network). Nevertheless, the same picture is not found for C3 campaign. Even if smaller particles sizes and coarse aerosol burdens are found with increasing distance from the emitting craters, the regression lines have less pronounced slopes and much smaller correlations (0.21 and 0.33) than for C2 campaign. This may be linked to the larger amount of emitted coarse ash particles and/or a smaller signal of gas-to-particles conversion to SA. For C3 campaign, sporadic explosive activity is found, with significantly larger ash emissions and with smaller (hbox {SO}_2) emission rates (average values of 1,315.9 ± 254.2 t/day), which can be limiting factors to secondary SA in-plume production. A similar negative correlation between the volcanic aerosol amount and average size, stronger for ash-free that ash-bearing plumes, as shown in Fig. 2c, was found during past field campaigns at different volcanoes, including Mount Etna28,29, and attributed to ash sedimentation and in-plume SA formation.

Simultaneous side-by-side measurements of MIIOM and MIISP volcanic AODs, and UVS (hbox {SO}_2) column amounts, were carried out during both C2 and C3 campaigns, and are compared in Fig. 3. The UV AOD correlates much better than the NIR AOD to UVS (hbox {SO}_2) estimations (Fig. 3a,b). In general, the correlation decreases with wavelength and is only significant in the UV-to-shorter-visible range (Fig. 3c). As shorter-wavelengths AOD values are more strongly correlated with the presence of small-sized particles in the aerosol layer than longer-wavelengths AOD, and (hbox {SO}_2) is the main precursor of the tiny secondary SA, this is a strong indication of (hbox {SO}_2)-to-SA conversion processes at play. In addition, UV and NIR AODs correlate both better to (hbox {SO}_2) during C2 than C3 campaign. This gives a further indication that possible (hbox {SO}_2)-to-SA conversion is more marked during volcanic passive degassing conditions, with stronger (hbox {SO}_2) emissions, than in mildly explosive conditions, with weaker (hbox {SO}_2) emissions. New SA particle formation has been recently identified for Etna’s degassing plume, at larger spatial scales, by Sahyoun et al.30. It is important to mention that the presence of in-plume and background atmospheric aerosols in the line-of-sight of the UVS spectrometer might have an impact leading to either over- or under-estimating (hbox {SO}_2) column amounts measurements, by multiple scattering and effective path lenght reduction due to light dilution, respectively31. These effects are expected to be negligible for the observations carried out in this study. Multiple scattering by in-plume aerosols is expected to be important for aerosol extinctions significantly larger than those observed in our study. The light dilution effect depends on the background atmosphere and is not related to the plume AOD. In addition, for vertical observation geometry, this effect is largely spectrally independent. As we took care to carry out observations at as much as possible vertical geometry, this has a negligible effect on the results discussed above.

Besides a general short-term depiction of the geographical near-source variability of volcanic AOD for Etna, in Fig. 2a different experiments are identified: two plume transects (TS1 and TS2, see details in Table 2), i.e. rapid longitudinal plume scans, and a plume half-traverse, i.e. a rapid perpendicular plume scan from the plume centre to one periphery. For the plume half-traverse, additional indications of the measured UV AODs are given in the map. The dramatic change in average size and coarse aerosol burden of the plume during the transects (details in Table 2), in such a short time and small distance, can be readily attributed to the sedimentation of possible coarser ash particles and the rapid formation of secondary SA. Coarser particles are found during TS2 than TS1, due to the visible ash puffs during this measurement session. Some residual relatively-coarse ash particles might be present at distal TS2 location due to these more pronounced ash emissions, which are unlikely present during TS1-distal observation ((hbox {AE} = 1.90pm 0.67), typical of extremely small particles). The larger relative uncertainties at proximal ((sim 95%), TS1, and (sim 60%), TS2) than distal UV AOD ((sim 45%), TS1, and (sim 35%), TS2), see Table 2, point at more inhomogeneous plumes near the sources. This plumes tend to get more and more temporally homogeneous with increasing distance, as different particles types (like fine ash) are removed from the plume along dispersion. Larger variabilities of volcanic plumes’ optical properties have been associated with puffiness of the plume and to the presence of ash29. The short-term atmospheric processes smooth the plume, in terms of its properties. The plume half-traverse revealed the perpendicular structure of Etna’s plume, during a typical passive degassing situation. The burden is maximum at the plume core (UV AOD (sim ;0.2)) and steadily decreases towards the plume periphery (UV AOD (sim ;0.1), i.e. half the burden than at the plume core). Correspondingly, the particle mean size decreases from the core to the periphery. New SA particles formation is expected to be more effective at the plume periphery, i.e. in presence of lower concentrations of pre-existing aerosols that can act as condensational sink30.

Distal and itinerant observations: vertical distribution

The vertical distributions of volcanic aerosols is further studied using LiDAR observations from the fixed SLN station. Figure 4 shows typical aerosol vertical structures at passive degassing conditions (19/07/2016). Back-trajectories analyses (not shown here) show that this profile is mostly affected by local air masses, at all altitudes. No larger scale features, like desert dust transport events, are observed at all altitude ranges. Zenith- and crater-pointing observations allow a detailed three-dimensional characterisation of the passive degassing plume. These observations reveal an aerosol layer at low altitude over the station (from near-ground to (sim) 2 km, aerosol signal peaking at (sim) 1 km, zenith-pointing observations), possibly covering the whole line of sight, from station to crater (crater-pointing observations, aerosol signal over the whole line-of-sight and peaking at (sim) 5–7 km, i.e. at and near the crater itself). The core of the volcanic degassing plume is characterised by weakly depolarising aerosols ((sim 3) and 1%, for zenith- and crater-pointing observations, respectively), thus indicating the predominant presence of spherical droplets, e.g. SA. The weak signal-to-noise ratio after plume crossing, i.e. visible in the average crater-pointing observations, may point at partially absorbing particles near the crater area, possibly SA embedding a more absorbing fine ash core.

On 19/07/2016, quasi-simultaneous MIIOM and LiDAR zenith observations have been realised at SLN station. This allows the direct comparison of the AODs obtained with the two instruments. The MIIOM UV AOD (at 320.5 nm) is (0.14pm 0.06), for the volcanic plume. Using a LiDAR ratio (LR) of 48 sr, as done previously, for volcanic plumes, at the same station32, we obtain a LiDAR VIS AOD (at 532 nm) of 0.06. This value, scaled using the quasi-simultaneous AE estimations made with MIIOM at SLN, produces a full-column LiDAR UV AOD (at 320.5 nm) of 0.16. While a limited impact of other local aerosol sources at the very lowest altitudes cannot be excluded, we attribute to the volcanic plume the aerosols observed up to 2 km over SLN station and obtain a LiDAR UV AOD of 0.12, for the volcanic plume. The consistency of quasi-simultaneous MIIOM photometric and LiDAR UV AODs, for the volcanic plume, validates our choice of LR and the LiDAR observations themselves.

Radiative transfer modelling: direct impact on the regional climate

These estimations of the radiative forcing can be compared with more distal estimations of the radiative forcing for Etna’s emissions during explosive eruptions. The radiative forcing efficiency (radiative forcing per AOD unit) has been estimated, at TOA and surface: (1) for a SA-dominated plume, at (sim 300,hbox {km}) downwind Etna for a moderate eruption (25/10/2013), at (sim) 40–50 and 50–(70,hbox {W}/(hbox {m}^2hbox {AOD}))11, and (2) for an ash-containing plume, at (sim 400,hbox {km}) downwind Etna for a relatively strong eruption (03/12/2015), at (sim ;112) and (145,hbox {W}/(hbox {m}^2hbox {AOD}))36. These values, scaled at the AOD observed in this case, lead to the forcings, at TOA and surface, of (sim ;-5)–6 and (-6)(8,hbox {W/m}^2), for the SA-dominated, and (sim) (-13) and (-17,hbox {W/m}^2), for the ash-dominated plume. Our estimations, for a proximal passive degassing plume, are more in the magnitude range of purely SA plumes. The proximal radiative forcing efficiency of Mount Etna’s passive degassing plume is comparable to the very distal (some hundreds km) efficiency of explosive SA-dominated plumes, though with a factor ten smaller AOD37.