Analysis of the ETNA 2015 Eruption Using WRF–Chem Model and Satellite Observations

The aim of the present work is to utilize a new functionality within the Weather Research and Forecasting model coupled with Chemistry (WRF–Chem) that allows simulating emission, transport, and settling of pollutants released during the Etna 2015 volcanic activities. This study constitutes the first systematic application of the WRF–Chem online-based approach to a specific Etna volcanic eruption, with possible effects involving the whole Mediterranean area. In this context, the attention has been focused on the eruption event, recorded from 3–7 December 2015, which led to the closure of the nearby Catania International Airport. Quantitative meteorological forecasts, analyses of Etna volcanic ash transport, and estimates of the ash ground deposition have been performed. In order to test the performance of the proposed approach, the model outputs have been compared with data provided by satellite sensors and Doppler radars. As a result, it emerges that, as far as the selected eruption event is concerned, the WRF–Chem model reasonably reproduces the distribution of SO2 and of volcanic ash. In addition, this modeling system may provide valuable support both to airport management and to local stakeholders including public administrations.


Introduction
Airport security is constantly under threat due to extreme weather events, such as strong wind-shear conditions and heavy rainfall, as well as due to natural hazards, like dust intrusions from desert and volcanic fallout. These latter may cause interruptions to flight operations and damage to ground infrastructures [1]. Indeed, the particulate expelled by volcanoes, consisting of small particles (typically ranging from a few millimeters to a few micrometers), can remain suspended in the air for usually long times, during which the particles retain their identity even if they are involved in physical-chemical processes in the atmosphere. For instance, suspended particles act as condensation nuclei for water droplets, favoring the occurrence of extreme weather phenomena, and acid rains [2,3]. In addition, the ash plume can cause serious damage to aircrafts during their flights. For example, it can be mentioned the emblematic case occurred on June 24th, 1982. The British Airways Boeing 747, with 263 people on board, crossed a cloud of volcanic ash expelled during the eruption of the volcano Galunggung (Indonesia). This event, never recorded before, caused the shutdown of the four airplane's engines [4]. The siliceous ashes sandblasted on the windshield and on the landing lights. Finally, The WRF-Chem package has been already applied to model emissions and transport of ash and SO 2 from volcanic eruptions all over the globe. In particular, Stuefer et al., [27] made a significant development on the WRF-Chem software architecture exploring the functionality that permits simulating emissions, transport, and settling of volcanic particles and gases. Webley et al. [28] analyzed the Eyjafjallajökull 2010 eruption simulating with the WRF-Chem model the ash cloud dispersing toward mainland Europe. Beside WRF-Chem, the 2010 Eyjafjallajökull eruption has also been analyzed through LIDAR observations [29], dispersion models [30,31] including FALL3D [32], and inverse modeling [33]. The WRF-Chem simulations performed by Webely et al. [28] used the Eruption Source Parameters (ESP) introduced by the fundamental work of Mastin et al., 2009 (hereinafter M09) [6].
A coupling of aerosol clouds and radiative processes in the WRF-Chem model has already been investigated [34,35], as well as the effect of biomass burning [36]. Further significant contributions to the development of the WRF-Chem package have been performed in Refs. [37][38][39][40].
In the present study, in order to perform a preliminary validation of the above-mentioned modeling package, we compare the meteorological and aerosols model predictions with experimental data obtained from reanalysis and satellite retrievals.
The paper is organized as follows: in Section 2, we describe in detail the considered case study, while the WRF-Chem model is introduced in Section 3. A semi-empirical description of eruption processes and input parameters of the WRF-Chem model is presented in Section 4, with results and discussion in Section 5. Finally, conclusions and perspectives are given in Section 6.

Case Study: December 2015
From December 3rd-7th, 2015, at Etna volcano, an explosive event occurred; it lasted about a week and was characterized by various episodes of different intensity. Moreover, these events had repercussions on the air traffic; in fact, the Catania International Airport reported numerous closures due to the volcanic ashes accumulated on the runway. In the weeks before the explosive event, constant strombolian activity (explosions of varying intensity and frequency, with the casting of incandescent material sometimes accompanied by small amounts of volcanic ash) was observed inside the Voragine crater (often called "central"), one of the four summit craters of Etna. The strombolian activity intensified on the evening of December 2th, culminating in a brief, but very violent, paroxysm on the night of December 3th. Around 02:20 UTC, the monitoring tools recorded the beginning of the eruption, which remained active until around 03:10 UTC. The eruption column exceeded 3 km in height above the vent, accompanied by volcanic explosions and lightning inside the ash cloud formed above the Voragine crater. The recorded data, obtained through the Voldorad Doppler radar, detected a signal/noise ratio of +18 dB and a longitudinal speed at the beam axis equal on average to 50 m s −1 , with peaks of 60 m s −1 . The eruption of December 3th was followed by three other similar paroxysmal episodes: two occurred on December 4th, respectively between 09:00 and 10:15 UTC and between 20:27 and 21:17 UTC, while the other one was detected in the afternoon of December 5th between 14:40 and 16.20 UTC. Figure 2 shows the intensity measurements of the volcanic activity (power), expressed in decibel-milliwatt (dBm) and taken from the database of the Physique du Globe de Clermont-Ferrand University [41]. Such data, captured by the OPGC Doppler radar (VOLDORAD 2B), provides immediate information on the eruption intensity by means of near-source measurements of echo power.
These three events were the most significant of the whole series of eruptions that affected Etna in this period of intense explosive activity. The eruption that occurred on the morning of December 4th was characterized by eruptive columns that reached about 7 km in height above the top of the volcano in about 15 min and at the same time frequent emissions of brown-greyish ash from a new crater were observed, along with high amplitude of volcanic tremor and strombolian explosions. The coarser pyroclastic material was deposited on the high south-western side of the volcano, at an elevation of about 2 km, while the ash fallout occurred in its eastern sector, covering the Giarre-Zafferana Etnea area. A representative map showing the position of Giarre and Zafferana Etnea with respect to Etna volcano is given in Figure 3.  The paroxysms that occurred on the evening of December 4th and in the afternoon of December 5th had almost the same characteristics as the previous ones, but the radial speed of the particles ejected by the volcano and the intensity of the volcanic activity (still measured by means of echo power) substantially decreased. In the late evening of December 5th, only degassing phenomena have been observed. At the dawn of the following day, the strombolian activity moved to the New Southeast Crater (NSEC), followed by landslides of unstable material (mixed with fragments of hot lava), which in a short time turned into pyroclastic flows. At the same time, a seismic swarm occurred in the north-eastern slope of Etna, probably connected with the volcano deflation due to the removal of a significant amount of magma: this deflation can be connected to the response of the tectonic structures in the eastern flank of the volcano. After about a week, the situation became normal again, and the explosive event finished.

The WRF-Chem Model
The numerical simulation of the eruptive event was carried out by means of the WRF-Chem model [18] in version 3.6.1 [42]. It consists of a type of "online" model, as it allows the evaluation of the meteorological fields together with the emission, transport, dispersion, and sedimentation of all anthropogenic and natural pollutants. In such a framework, it is possible to consider both physical processes and chemical transformations of aerosol particles released during volcanic eruptions. This key characteristic has been shown to significantly improve the accuracy of model forecasts in comparison with offline models [43]. Indeed, it is worth noting that, although weather conditions are the main factor that determines air quality, they may also depend on the direct and indirect effects that aerosol compounds have on the solar radiation and on the cloud microphysics. Most of the Volcanic Ash Transport and Dispersion (VATD) approaches are "offline" models, since they separately describe the physics and chemistry characterizing the dispersion of volcanic emissions into the atmosphere, and generally, given their low computational cost, they are largely adopted by the operational forecasting centers. However, it is more realistic to privilege the "online" approach, since an inappropriate treatment of atmospheric processes and their feedback with aerosol components could lead to incorrect predictions concerning the deposition of the ash and the diffusion of aerosols [43,44].
The WRF-Chem is based on the WRF model [17,[45][46][47][48], which supports both research applications and operational meteorological forecasting applications [49], including various options for dynamic cores and physical parameterizations [50]: therefore, it can be used to simulate atmospheric conditions on a wide range of space and time scales. The chemical package that simulates tracer gases and aerosol interactively with the meteorological fields making use of different photochemical treatments and aerosols models has been integrated by using the WRF-ARW numerical core [51]. Among others, the WRF-Chem model can be used to forecast and simulate meteorological/climatic conditions on a regional and local scale, predicting the release and transport of pollutants and natural aerosols [52], estimating the air quality, and studying important processes for global climate change [35]. The capability of the model to predict the transport and concentration of ash clouds and SO 2 depends on the initial information regarding volcanic emissions, such as the scale of the eruption (which includes the expelled mass), the height of the plume above the vent, the eruption rate, the date, and duration of the event [27,53].
In our simulations, the spectrum of the grain size distribution of the ash particles has been taken from Stuefer et al. [27]. In particular in their Table 1, adapted from M09 [6], they classified Eruption Source Parameters (ESP) for Etna as "small mafic", or M1 [24].  Figure 4 shows the numerical domain adopted in the present work: this domain includes a part of North Africa, Italy, and a part of Balkan Europe, with a 180 × 180 grid centered at a latitude of 37.74 • and longitude of 15.18 • . The horizontal space of the grid is 10 km in both directions with 40 vertical levels up to 50 hPa. The initial and boundary conditions were acquired by NCAR/NCEP Final Analysis (FNL by GPS) (ds083.2), with resolution of 1 degree, every six hours [54].  Table 1 shows the physical schemes selected for the setting of the parameterizations, specifically Mellor-Yamada-Nakanishi and Niino (MYNN, level 2.5) for the planetary boundary layer [55], surface similarity [56], and the surface model RUC (rapid update cycle) [57]. These schemes have been chosen to represent the physics of the surface layer and the parametrization of the Earth's surface. The Rapid Radiative Transfer Model (RRTMG) for both short-wave (ra_sw_physics = 4) and long-wave (ra_lw_physics = 4) radiation is used for the direct radiative effect of aerosol [58]. The Purdue Lin scheme (mp_physics = 2) is used for microphysics processes. This configuration is compatible with short-wave radiative feedback [59].
The implementations of the other parameters follow the guidelines of Stuefer et al., 2009 [27] (such as the mass of volcanic ash, set to 1.7 × 10 9 kg), or are based on estimates derived from the literature (as for the mass of SO 2 set to 1.0 × 10 9 kg) [60]. These parameters are reported in Table 2: for the injection heights above the vent and duration, the observations taken from the cameras and radar Doppler VOLDORAD 2B [39] have been adopted. Such a radar is placed at approximately 3 km from the top of the volcano and has been made available by the Observatoire de Physique do Globe de Clermont-Ferrand (OPGC) in collaboration with the Istituto Nazionale di Geofisica e Vulcanologia (INGV) of Catania.

Semi-Empirical Description of Eruption Processes and Input Parameters of the WRF-Chem Model
The classification drawn up by M09 is based on the study of well-documented past eruptions and includes eleven types of eruptions. If new observations reveal data different from the predefined ones, one needs to find the type of eruption closest to the one under study and use its parameters. Stuefer et al., 2013 [27], divided the granulometric size spectrum of the ash particles into 10 bins (in their Table 2) and coded them within the WRF-Chem software. This classification considers particles' diameter in a range going from 2 mm to less than 3.9 µm and the corresponding percentage of the mass fraction depending on the type of Eruption Source Parameters (ESP). The latter are needed to predict the transport and dispersion of volcanic ash clouds during eruptions, with the aim of identifying risks for the air navigation. For each ash-bin (Table 3), it is possible to define the aerodynamic radius as half of the arithmetic mean between the limits of the diameters of each ash-bin. The erupted mass is distributed in the so-called umbrella cloud model, for which it is assumed that 75% of the erupted mass is detrained in the umbrella cloud and 25% beneath, with a linear distribution from the umbrella to the vent [27]. According to M09, all the possible eruption categories may be classified as M-type, or mafic types, which include basaltic and ultramafic magmas; and S-type, or silicic types, which include andesite, dacite, rhyolite, and others such as phonolite that can produce high ash columns (M09, Table 3, pag.16). As far as Etna is concerned, the M09 classification assigned the class M1 to its eruptions, which are mainly of an effusive type, i.e., they are characterized by basaltic dense-lava fountains, with relatively high temperatures (>1000 • C). In Table 3, it is reported the particle size distribution assigned to the emitted particles diameters in the range between 1000 µm and 31.25 µm.
It has been formulated considering for Etna a grain distribution that has been obtained by the semi-empirical methodology defined by Scollo et al. [24]. According to the classification reported in the table, 4% of particles have a diameter between 0.5 and 1 mm (vash_2), 10% between 0.25 and 0.5 mm (vash_3), 50% between 125 and 250 µm (vash_4), 34% between 62.50 and 125 µm (vash_5), and 2% between 31.25 and 62.50 µm (vash_6). No particles with diameter less than 31.25 µm are considered [27]. It is evident from Table 3 that smaller particles (bins 7-10) are not considered in M1 classification. In this context, it may be crucial for future studies to consider in-situ measurements, for example from an Optical Particle Counter, for the direct measurement of the size distribution of volcanic aerosols after the Etna eruption [61].
In Figure 5, the distributions of volcanic ash emitted by Etna (M1) are compared with those observed in different classes of eruptions, according to the database by M09. In this figure, each point represents the midpoint of each bin reported in Table 3, as stated above. It emerges that the distribution of Etna ash shows a well-defined peak at a particle size of 175 µm, with about 50% of emitted particles having this size. In addition, it is worth noting that this is the less polydisperse distribution, with only five values other than zero. These findings suggest that this particular class of eruptions may attain very high temperatures in order to generate a huge amount of particles of relatively large size with a low degree of polydispersity and this result should likely be verified by local measurements as stated above. It is worth noting that also different total grain size distributions could be adopted for the Etna volcano, such those reported in Ref. [62]. In future studies, we plan to compare these different distributions in order to assess how the ash dispersion and transportation depends on their spectrum in the WRF-Chem model. The input parameters of the WRF-Chem model are determined by using a software pre-processor called PREP-CHEM-SRC [63], developed to prepare the grid emission parameters for the most important anthropogenic pollutants. In this work, the Version 1.5 of the PREP-CHEM-SRC has been employed. Such a pre-processor is specifically designed to include input parameters concerning volcanic eruptions. Once the height of the plume above the vent and the eruption duration are defined, the total erupted mass can be estimated according to the following equation, holding for volcanic ash only (see M09): where M is the mass eruption rate (in kg), ρ is the ash density (in kgm −3 ), d the eruption duration (in s), and ρ the plume height (in km) above the vent. Moreover, the PREP-CHEM-SRC provides the position of the volcano in the numerical domain considering the nearest grid cell.
It is important to mention that in addition to the ash, during a volcanic eruption, a considerable amount of SO 2 is also emitted, which, once in the atmosphere, oxidizes and transforms into sulfuric acid (H 2 SO 4 ); the latter condenses into sulfate aerosols, characterized by a residence time in the atmosphere proportional to the gases containing sulfur present in the volcanic plume. Unlike the ash that settles within few days, SO 2 can last up to several months [12].
The data provided by the international program AEROCOM (started in 2002 with the target to reduce the uncertainty on the impact of aerosols on the climate system) contain the volcanic emissions of SO 2 and other variables for the time interval going from January 1979 to December 2009 for all the volcanoes listed in the Global Volcanism Program database provided by the Smithsonian Institution. In particular, there is a file for each year containing the number of events occurred, along with the name of the volcano, the date, the height above the average sea level, the height of the plume above the vent, the longitude, the latitude, and the daily emission rate of SO 2 , distinguishing between volcanic eruptive and non-eruptive emissions. In analogy with the preprocessing phase of ash emissions, the SO 2 emissions are placed in a grid point of the spatial domain, which includes the position of the investigated volcano. The total SO 2 emission corresponds to the amount of SO 2 emitted by the volcano inside the grid point, and it is expressed in kg m −2 dy −1 [27].
The pre-processing tool, as in the case of volcanic ash, places the SO 2 emissions in the grid point of the WRF-Chem domain that surrounds the geographical position of the investigated volcano. The total emission is calculated by adding the emissions of all the volcanoes inside the cell and is expressed in Kg km −2 hr −1 . However, it is important to point out that this value is often only approximated, in consideration of the difficulty usually met in accurately estimating the correct total amount of SO 2 emitted during a volcanic eruption.

Geopotential
In Figure 6a-d, the comparison between the geopotential height at an atmospheric pressure of 500 hPa predicted by the WRF-Chem (left panels) and data downloaded from ECMWF/ERA5-reanalysis (right panels) is reported. ERA5-reanalysis are available from the Copernicus Climate Change Service (C3S, 2017) at the following URL (https://www.ecmwf.int/en/forecasts/datasets/reanalysis-datasets/ era5) [64].
During December 4th (Figure 6a), ERA5-reanalysis (right panel) shows a pressure minimum localized between the southern coast of Sardinia and northern coast of Tunisia that fosters a south-east circulation near Etna. The following day, December 5th ( Figure 6b, right panel), ERA5-reanalysis depicts that the pressure minimum is further dislocated toward the south entering the African continent near the Tunisia/Libyan border. The consequent circulation at Etna is rotated in the north-west direction. On December 6th, again considering ERA5-reanalysis (Figure 6c, right panel), a deepening of the low pressure between Balearic Islands and Sardinia producing northward wind near Etna is evident. Finally, on December 7th, when the eruptive intensity started to decrease, ERA5-reanalysis (Figure 6d, right panel) shows that the cyclonic system moved in the north-west direction toward the southern coast of France, inducing a northern flow near Etna. The close agreement between the ERA5 (Figure 6, right panels) geopotential at 500 hPa and the WRF-Chem model is evident (Figure 6, left panels). This is quite important, as the dispersion of eruptive plumes of ashes and gases (SO 2 ) is strictly connected with wind speed and direction.

Transport of SO 2
Sulfur dioxide (SO 2 ) is a colorless, irritating, non-flammable gas, very soluble in water and with a pungent smell. This compound is one of the most widespread, dangerous, and studied atmospheric pollutants. It derives from the oxidation of sulfur during the combustion processes of substances that contain this element. Due to the characteristics of sulfur dioxide, it is interesting to analyze its distribution in the course of the studied event. At this purpose, the daily data relating to the columnar density of SO 2 , expressed in Dobson Units (DU), have been downloaded from the NASA/Earthdata program [66] through one of its dedicated web-portals [67]. The DU indicates the concentration of gas in a column of air, above a certain point on the Earth's surface. The daily SO 2 total column data have been provided by the EOS-Aura satellite [68] via the OMI sensor (Ozone Monitoring Instrument) with a spatial resolution of 0.25 • × 0.25 • . The SO 2 observational data have been compared with the output of the WRF-Chem model: the latter have been first averaged over a single day and then column-integrated in every surface grid point, finally providing values expressed in DU. The conversion factor is: 1 DU ≈ 28,000 µg m −2 .
The WRF-Chem simulations allow the analysis of the spatial distribution of daily averaged sulfur dioxide emitted during the eruptive event under investigation. The scale of the daily average SO 2 concentrations shown in Figure 7 identifies values ranging from 0 to 1 DU. In particular, on December 4th (Figure 7a), there is an intense emission of SO 2 that spreads out mainly to the south-east with a plume that crosses the Mediterranean until it almost reaches the coast of Africa. Between December 5th and 7th (Figure 7b-d) another intense flow of SO 2 from Etna is transported to the north, covering the entire Italian southern peninsula. In this context, it may be evidenced that the results of the simulations provide a spatial distribution of the daily average SO 2 that reflects in a fairly precise way the synoptic analysis and the description of the event discussed above (see Figure 6).
A comparison between the results obtained by using the WRF-Chem model and the daily data retrieved by the Aura satellite via the OMI sensor has been also performed. At 13:00 UTC on December 4th, immediately after the paroxysm, the OMI map reports approximately the same distribution of SO 2 , characterized by a shift of the flow towards the south-east (Figure 8a). According to the simulation results, the SO 2 concentration amounts to approximately 1.8 DU, close enough to the observed value of about 2 DU. Due to the lack of satellite data, it is not possible to make a comparison on December 5th. On December 6th (Figure 8b), there is a good qualitative correspondence between the spatial distribution predicted by the model and that provided by the satellite, albeit the amount of SO 2 appears slightly underestimated by the simulations. In particular, it is possible to observe from both maps the transport of the SO 2 to the north, up to the Gargano and the Adriatic regions. Finally, with reference to December 7th (Figure 8c), both maps show that sulfur dioxide, spread to central Italy, decreased in concentration, although it is not possible to explain the presence of an abundant quantity of SO 2 recorded by the OMI sensor between Sicily and Calabria. The origin of this discrepancy can be ascribed to the fact that in the WRF-Chem approach, we simulate only the three main events (see Table 2), whereas puffs of erupted material can be expected along the entire duration of the event (degassing). Despite the discrepancy observed on the last day of the eruption, the overall distribution of SO 2 along the whole duration of the process appears fairly well reproduced by the WRF-Chem model.

Distribution of Volcanic Ash
In addition to the SO 2 , in this study the transport of ash particles of different sizes has been analysed for the time interval from December 4th to 7th.
It is important to mention that in its default configurations, the WRF-Chem model uses the Stokes law corrected by the Cunningham slip factor [27] for calculating the terminal fall velocity. This is well suited for small particles only (typically with a diameter < 0.125 mm). Such an issue could be solved by implementing in the WRF-Chem model different formulations of the settling velocity, as detailed, for instance, in Refs. [69][70][71][72]. In this context, for the following analysis, we showed results only for vash_5 and vash_6 particles. As reported in Figure 9, the average values of the different vash sizes range between 10 −1 and 10 8 µg m −2 . These data, expressed as mass per unit area, have been obtained by integrating the column load. It is observed that the larger particles (vash_5), being heavier, have fallen back into the areas adjacent to the volcano, while the finer ones (vash_6), and therefore lighter, have been transported to greater distances, reaching central Italy and the African coast. In order to validate the model predictions, the simulation results have been compared with observed data characterized by the Aerosol Index (AI), explicitly developed to detect the presence of aerosols in the air. The aerosol measurements were carried out again by the OMI sensor on board the above-described NASA/Aura spacecraft. In Figure 10, the area circumscribed by the dashed line highlights the regions affected by the dispersion of aerosols due to the eruption of Etna, showing how the suspension of particles is spread mainly over Sicily, the Mediterranean Sea, and the Gargano, in agreement with the wind circulation derived from the synoptic analysis previously described. There are also high concentrations of aerosols (orange coloring) in the south-east part of Sicily, but since the wind mainly blew northwards during the paroxysmal episodes, this can be ascribed to dust from the Sahara. A qualitative agreement with the WRF-Chem results is clearly visible, albeit a more detailed comparison is not possible due to the different units adopted. A further analysis has been done considering only the volcanic ash with diameter between 15.62 and 31.25 µm (vash_7) that, according to M09 distribution reported in Table 3, should not be present for M1 class. Therefore, the starting amount of vash_7 has been set by choosing that corresponding to the M2 classification (see M09); according to this prescription, the percentage of vash_7 particles is 1.5% of the total ash. This supplementary study aims to investigate the transport of small particles in the numerical domain. In Figure 11a, the WRF-Chem predictions on the transport of vash_7 are shown: as observed, in this case, the finest volcanic particles can be found in the eastern regions of Europe. Finally, in Figure 11b, the amount of vash_7 that exceeds the detection limit of satellite-based ash retrieval algorithms of 0.1 g m −2 is reported. Even if we only consider vash_7, the large distribution of volcanic ash above the threshold limit that may be considered potentially dangerous for any flight operations is evident. This aspect needs to be further investigated and verified by future studies by considering the LIDAR measurements network in Europe [73,74] or the satellite retrievals from the Spinning Enhanced Visible and Infrared Imager (SEVIRI) onboard the Meteosat Second Generation spacecraft.

Conclusions
In this work, we have presented an application of the WRF-Chem model to systematically investigate the Etna volcanic eruption registered at the beginning of December 2015 and its role in affecting the security of the airports involved as well as the transportation of pollutants and volcanic ash in the Mediterranean area. As far as the case study considered is concerned, the WRF-Chem model has proved to be an excellent candidate for the prediction of the transport and dispersion of natural and anthropogenic aerosols, as well as for volcanic ash and sulfur dioxide emitted during the volcanic eruption. In this study, the output obtained from the numerical model agrees with the satellite retrievals and with the meteorological reanalysis. It has been possible to reconstruct the dynamics of Etna volcanic ash transport and to confirm the data concerning both the concentration and the direction of propagation of ash and gases emitted into the atmosphere.
More in detail, immediately after the first paroxysm, the ejected material was transported towards the north-east by winds. Then, the circulation underwent a south-east rotation and, at the end of the eruptive event, the volcanic cloud overlooked the skies of southern Italy, up to the Gargano, the Balkans, and the African coast. Concentrations of SO 2 , of the order of 1.8 DU, were recorded in the time bands of maximum explosive activity, which progressively decreased to a few tenths. The particles expelled during the eruption reached maximum concentrations 2.0 DU, and the finest ones (vash_6) went as far as the Balkans and Africa, while those with a larger diameter (vash_5) fell back in the areas surrounding Etna mountain.
Although the model has exhaustively reproduced the event, it is advisable to modify the WRF-Chem input settings, especially with regard to the particle size distribution to be adopted. The results, obtained by simulating the transport of volcanic ash including particles having dimensions between 15.62 and 31.25 µm (vash_7), suggest the possible spreading out of Etna volcanic ash until the east of Europe. Upon summarizing the model predictions and the relative comparison with the observed data, it emerges that the WRF-Chem model can be a valid support both to airport management companies affected by volcanic emission phenomena and to other local stakeholders such as public administrations. In addition, the capability of the model to describe the diffusion of Etna volcanic ash and SO 2 in the atmosphere is particularly useful in order to better understand the impact of this volcano on the meteorological phenomena involving the whole Mediterranean area. Additional developments of the WRF-Chem model aiming to further increase the accuracy of the predictions are currently underway, as well as the comparison with different models, including FALL3D [32], HYSPLIT [75], and Ash3d [76].
Author Contributions: All authors contributed to conceptualization, methodology, investigation, writing, reviewing, and editing of the present work. All authors have read and agreed to the published version of the manuscript.