Five Years of Dust Episodes at the Southern Italy GAW Regional Coastal Mediterranean Observatory: Multisensors and Modeling Analysis

The Mediterranean area is a climate-change hotspot because of the natural and anthropogenic pollution pressure. The presence of natural aerosols, such as dust, influences solar radiation and contributes to the detection, in storm episodes, of significant concentrations of PM10 in Southern Italy, where generally fresh and clean air is due to local circulation, and particulate matter concentrations are very low. We present the results of medium-term observations (2015–2019) at Lamezia Terme GAW (Global Atmospheric Watch) Regional Observatory, with the purpose of identifying the dust incursion events by studying the aerosol properties in the site. To achieve this goal, the experimental data, collected by several instruments, have been also correlated with the large-scale atmospheric patterns derived by the ERA5 reanalysis dataset, in order to study the meteorological conditions that strongly influence dust outbreaks and their spatio-temporal behavior. An intense dust-outbreak episode, which occurred on 23–27 April 2019, was chosen as a case study; a detailed analysis was carried out considering surface and column optical properties, chemical properties, large-scale pattern circulation, air-quality modeling/satellite products, and back-trajectory analysis, to confirm the capability of the modeled large-scale atmospheric fields to correctly simulate the conditions mainly related to the desert dust-outbreak events.


Introduction
The contribution of desert dust emission in the troposphere to global atmospheric aerosols (1.5 × 10 9 tons/year) is due to arid regions [1,2]. The principal dust-producing regions on Earth are the Sahara Desert, the Middle East, Northwest China, Southwest Asia, central Australia, basins in southern Africa (in Namibia and Botswana), some areas in Bolivia, and the Great Basin (USA) [3,4].
In Europe, the main regions responsible for desert dust transport events are the northern African regions, in particular the northwestern part (Sahara Desert) and the northeast part (i.e., Egypt, Libya, Sudan, etc.) [5,6].
The most frequent trajectory patterns related to the transport of Saharan dust are four [7,8]. More than half (about 60%) of the dust loaded from Saharan sources remains in Africa (deposited in the

Data and Methods
In this section, we describe a comprehensive dataset, involving several instruments, used for our study. Satellite and model used tools are also described together with the method adopted to identify dust events occurring at our experimental site.

The Experimental Site
The CNR-ISAC observatory of Lamezia Terme is a Regional Global Atmospheric Watch of World Meteorological Organization (GAW-WMO) site (LMT, 38.88 N 16.23 E; 6 m.a.s.l.) is managed by the Institute of Atmospheric Science and Climate (ISAC) of the National Research Council (CNR). LMT is a coastal site located about 600 m inland from the Tyrrhenian coastline (west side; Figure 1). The area is characterized by anthropogenic pollution emissions related to transportation and agriculture due, in particular, to the presence of Lamezia Terme International Airport (north direction) (14,699 flights in 2015) and the city of Lamezia Terme (north-east direction). The highway (north-east direction) runs around the observatory location clockwise from N to S, and is located 7 km (northward) to 3.5 km (southward) from the observatory. The area is characterized by moderate wind breezes, converging on Marcellinara gap between the Jonian and Thyrrenian Sea from the sea (NW-SW) (see Figure 1), which mainly develop during daytime, while NE gentle wind breezes from land mainly affecting the night-time period [39].

The Instrumental Setup
Several instruments within different monitoring programs have been active at LMT Observatory since 2015.
A unique pump is used for aerosol sampling, maintaining a uniform flow rate of 200 L/min pumped in the main inlet tube, at 4 m height, equipped with a PM 10 head (50%) transmission efficiency at 10 µm aerodynamic diameter. The aerosol stream is split isokinetically into several instruments such as nephelometer, Multi-Angle Absorption Photometer (MAAP), Optical Particle Counter (OPC), Beta Attenuation Sampling and Measurement System (SWAM).
Regarding quality data checking, the following procedures were guaranteed: • raw data with relative humidity greater than 60% condition were rejected according to [31,32]; • a periodical full nephelometer calibration was performed using CO2 as span gas; • a zero calibration was carried out once an hour using internally filtered particle-free air.
The MAAP (5012 Model Thermo Scientific™ (Franklin, MA, USA) operates by measuring black carbon and aerosol short-wave radiation absorption properties. The aerosol absorption (short-wave radiation) coefficient (sa) and continuous measurements of equivalent black carbon (eBC) concentrations are measured at 637 nm at 1 min resolution time [41].
The MAAP determines aerosol absorption by illuminating a particle-loaded filter and simultaneously measuring the radiation passing through the filter and the backscattered light at three different angles [42]. The sample flow provided by the internal pump was set to 16.7 L/min. A mass closure, for the black carbon mass concentration measured by MAAP, was verified by ensuring that the eBC mass concentration was lower than 50% of the PM 2.5 mass concentration.
Nephelometer and MAAP operate according to SOP GAW-WMO specification [43] The OPC (Multichannel Monitor, FAI Instruments, Rome, Italy) [44] measures the number size distribution of particles in 22 size bins with diameters ranging from 0.28-10 µm. OPC operates with a 1 L/min flow rate at 1 min temporal resolution and the relative humidity of its air samples is kept lower than 50% using dilution (1:2) (a clean dried air and a smart heater placed in the diluter along the mixing chamber is automatically activated only when needed). It ensures proper conditioning of the aerosol and it minimizes of the effects of the water vapor. PM 10 and PM 2.5 mass concentration measurements are performed by β-ray attenuation method using a low volume (2.3 m 3 /h) FAI SWAM 5a-Dual Channel Monitor [45].
An automatic weather station (Vaisala WXT520, Finland) measures, at 10 m.a.s.l., the following meteorological parameters: temperature, relative humidity, wind speed and direction, pressure, and rain (10 averaged minutes).
A sun photometer (CIMEL CE-318, Paris, France) [46], operative within AErosol RObotic Network (AERONET) [47], is used for ground-based aerosol measurements and in particular atmospheric aerosol columnar properties. It is a multichannel, automatic sun-and-sky scanning radiometer measuring the direct solar irradiance and sky radiance at the Earth's surface. The sun photometer allows the obtaining of, with direct sun measurements with a 1.2 • full, the Aerosol Optical Depth (AOD) at 340, 380, 440, 500, 675, 870, 1020 and 1640 nm, along with the water vapor column content and the estimation of several optical and microphysical aerosol properties, such as the refractive index and the size distribution. The system is fully automatic and powered by a solar panel. The solar irradiances received by the photometer are inverted into the columnar aerosol optical and microphysical properties through the implementation of retrieval algorithms [48]. The used standardized methods concerning instrument maintenance, calibration, cloud screening and data processing allow aerosol monitoring and comparison between different study periods and areas [47]. The uncertainty in the retrieval of AOD under cloud-free conditions is ±0.01 for wavelengths larger than 440 nm and ±0.02 for shorter wavelengths [49]. Errors in aerosol size distribution retrievals depend on particle size, aerosol type, and actual values of the size distribution. For particles in the size range 0.1 < r < 7 µm, retrieval is around 10%-35%, while for diameters/particle sizes lower than 1 µm and higher than 7 µm retrieval errors rise to 80%-100% [48,50]. For measurements under low aerosol load conditions (i.e., AOD (440 nm) < 0.2), the retrieval accuracy decreases significantly [48,50]. In this work, the AERONET level 2.0 data are used, in relation to the dust-outbreak case, including the AOD and inversion products derived, SAE at 440-675 nm. Moreover, the accuracy of the retrieved AOD by the CIMEL radiometer is spectrally dependent, being better (<0.01) for wavelengths longer than 440 nm and lower (<0.02) for the UV wavelengths [49].
A Scanning Electron Microscope (SEM) Zeiss EVO MA15, equipped with Energy Dispersive X-ray Analysis (EDX) Microanalysis, BRUKER Detector Xflash 6-10, located at Research Center of the National Institute for Insurance against Accidents at Work in Lamezia Terme IT, was used to analyze the filter samples collected with the SWAM. SEM is managed by "Smart SEM" software and the EDX microanalysis is managed by "Bruker Quantax" software. A lanthanum filament is used, and the instrument calibration is repeatedly checked by the operator using a copper standard. The filament current is adjusted to an energy voltage of 20 keV.
The analysis was carried out for dust particles collected on quartz filters by means of a SWAM 5a Dual Channel Beta Attenuation Sampling and Measurement Systems, using a PM 10 cut-off.
The portions of each filter are mounted on (12.5 mm diameter) aluminum stubs and metalized using gold. To proceed with recognition of individual particles, a portion of the blank filter acts as a reference. To minimize the recognition of false particles, the acquired images, of backscattered electrons (BSE), are compared with secondary electron (SE) ones. The SEM analysis technique, combined with the EDX Energy Dispersive X-ray Analysis, is frequently used to obtain geochemical and morphology information on single dust particles laid on filters. It is also possible to perform the morphology of dust particles and related chemical elemental analysis.

Modeling Products
With the aim of studying the main large-scale atmospheric patterns linked to the revealed dust events, we analyzed the European Centre for Medium-Range Weather Forecasts (ECMWF) ERA5 (Copernicus Climate-Change Service (C3S) (2017)) [51]; ERA5 fields have a horizontal resolution of about 31 km and 137 vertical levels from the surface to 0.01 hPa.
The considered fields, in this work, are geopotential height at 500 hPa (Z500), mean sea level pressure (MSLP), and U and V wind component(s) at 925 hPa (in order to calculate wind speed and direction at this level; WIND925). These parameters are extracted at 00 UTC, for the whole 5-year period.
To confirm the occurrence of the dust events, to evaluate the location of the dust sources for the investigated cases and their simulated spatio-temporal propagation, we also analyzed (and plotted) some CAMS (Copernicus Atmosphere Monitoring Service) model products [52,53]. In particular, we used the CAMS Global archived analysis (6 hourly analysis, approximately 40km horizontal resolution, on 60 vertical levels from the surface up to 0.1 hPa). The main parameter taken into account is the Dust Aerosol Optical Depth, obtained from CAMS simulations vertically integrating the aerosol extinction coefficient at 550 nm.
For a single-event discussion, the back-trajectories computed with the HYSPLIT (Hybrid Single-Particle Lagrangian Integrated Trajectory) model [54] are considered and commented upon; the trajectories are obtained using the READY (Real-time Environmental Applications and Display sYstem) platform (https://www.ready.noaa.gov/index.php) [55].

Satellite Products
To search for further confirmation of the characteristics of the dust events, images from the Moderate Resolution Imaging Spectrometer (MODIS) on NASA's Aqua and Terra satellite are also considered.
The parameter taken into account is the MODIS Near Real-Time NRT Level 3 value-added Aerosol Optical Depth (MCDAODHD) [56][57][58][59][60]; the sensor resolution is 0.5 degrees, the imagery resolution is 2 km, and the temporal resolution is daily.

Optical Properties
On the entire dataset, we analyzed the total scattering and backscattering coefficients data from nephelometer, for the entire period. The non-Lambertian and truncation error was applied at raw signal of σsc and σbsc by using the method in [38,[61][62][63].
We averaged raw data at 5-minute intervals. For σsc, we used detection limits values, signal-to-noise ratio of 2, of 0.44, 0.17, and 0.26 Mm −1 and for σbsc the detection limits were 0.29, 0.11, and 0.21 Mm −1 both respectively at 450, 550, and 700 nm [40]. Based on σsc and σbsc coefficients, three intensive (independent of the amount of particles and dimensionless) properties were calculated: the hemispheric backscattering fraction (b), asymmetry parameter (g) and Angstrom coefficient SAE. Using Equation (1), b λ was calculated at the wavelength of λ 1 = 450 nm, λ 2 = 550 nm and λ 3 = 700 nm for the whole period: The b λ parameter is weakly dependent on concentration and it could be very useful for describing the aerosol-scattered radiation reaching the ground and it is very useful in radiative transfer models [64].
An analytic relationship by [65] is adopted in this manuscript (Equation (2)) to calculate g using b λ .
The asymmetry parameter ranges from −1 (entirely backscattered light) to 1 (entirely forward-scattered light) [64] and g is equal to zero value for symmetric (e.g., Rayleigh) scattering light. This parameter is one of the optical properties widely used to estimate aerosol forcing of climate. In this study, in the discussion of results, we used g at λ = 550nm. SAE was calculated for each pair of wavelengths: SAE 450,700 , SAE 550,700 , and SAE 450,550 by using Equation (3) from the corrected scattering data.
A prevalence of coarse particles (dust, sea salt) occurs when the SAE values are small (SAE 450,700 < 0.5), while SAE 450,700 > 1.5 indicate a preponderance of small particles such as secondary sulphates, nitrate aerosols, and biomass burning particles [50]. Although this interpretation could be unclear if a mixture of fine and coarse mode is present in the sampled air [66,67]. In this work, SAE 450,700 will be used to discuss our results. The absorption coefficient (sa) was calculated at 637 nm [48] from MAAP measurements according to Equation (3).
where eBC is the black carbon mass (µg/m 3 ) measured by MAAP and MAC = 6.6m 2 /g is the Mass Absorption Coefficient (defined in [68]).

The Identification of Dust Events: The Method
African dust events were identified considering the aerosol optical properties at the surface, using an approach adopted in past studies [69,70]. The results were matched with available ground measurements of PM coarse fraction and PM coarse and fine fraction to find the number of days when dust events occurred. Also, advection of African dust and air mass back-trajectory analysis [71,72] was used to identify the specific period of dust events in the Mediterranean site. In our case study, we used such an approach to reconstruct the event, already selected by in situ observations, and identify the possible direction sources of sampled dust at LMT. Additional data based on the concentration of mineral elements (mainly Fe, Al, Si, Ti), typically contained in dust particles, are available for the chosen case study, using the SEM instruments. A characterization of the aerosol types, present at surface, can be based on measurements of σsc and σbsc , Ångström exponent (SAE).
SAE coefficient mainly depends on the particle size distribution (low values of SAE indicate a prevailing role of coarse particles). The combined use of SAE, with other optical parameters, allows the identification of different aerosol types at ground, including dust; usually, low values of SAE (<0.5) are typical of Saharan dust [71]. In [69] different kinds of particles, using surface optical properties, were identified and studied. The authors show that during dust events, the mean values of SAE 450,700 , g and SSA are respectively less than 0.7 ± 0.34, 0.62 ± 0.05 and 0.86 are thresholds to characterize dust events at the Lamezia Terme site. Therefore, in this study, we selected dust events largely dominant at the surface, using threshold values for SAE 450,700 less than 0.7, and, with the application of this threshold, the number of total cases found is 143 days (with probably dust events). These 143 days matched the results of PM coarse fraction dataset (SWAM instruments). Regarding PM coarse fraction, a weighting daily mean was performed. To identify, with a high probability, a dust event occurrence, we adopted a threshold value. Such a threshold, for our site, according to the coarse fraction mean value, is the 95th percentile (19.9 µg/m 3 ). This value is likely associated with a dust event considering PM 10 concentration. The initial 143 events that matched the 19.9 µg/m 3 mass were then reduced to 62 days.
Continuous mixing with cleaner (marine) air masses, typical of the marine-coastal environments, activates a frequent lowering of eBC concentration [73] at LMT, with eBC values much lower (0.64 µg/m 3 ) than the PM 10 (10.69 µg/m 3 ), on the whole measured period. During the dust events, the concentration of PM 10 usually rose drastically with peaks of about 50-60 µg/m 3 while the eBC concentration remained close to the mean values with peaks of about 1.6 µg/m 3 . On the 62 dust events, we performed an additional check of the eBC concentration with the aim of excluding the events with high mean daily eBC mass concentration values, which contribute to the coarse fraction of the PM. This will allow us to study only cases of dust not affected by the presence of local pollutants due to the products of combustion, such as biomass burning and road traffic. From the 62 events, all those with a mean daily of eBC concentration values above the 90 percentile (1.01 µg/m 3 ), as shown in Table 1, were therefore excluded from the analysis. With this further discrimination from 62 dust episodes, we pass to 52 in total, as in Table 2 within the initial identified 143 days with dust. Using such an approach, we identified the events as reported in Table 2, their duration and the identification of the season in which they occurred. We divided cases occurring in the cold from the warm season by considering the respective following grouping of months: respectively for October, November, December, January, February and March (ONDJFM), and April, May, June, July, August, and September (AMJLAS). In Table 3 the PM coarse fraction statistics in warm and cold season are shown. Regarding the mean values from warm and cold season, they are very similar, with differences of 2.51 µg/m 3 , the maximum value is reached in cold season with a peak value of 60.04 µg/m 3 . To verify and confirm the dust events identified by the optical surface characteristics and PM deposition, in Tables 2 and 3, we discussed above the synoptic conditions over the Mediterranean area, associating each of the 15 "dust events" with the large-scale ERA5 fields at 00 UTC on the first day of each case.
In summary, to use our approach, two main kinds of measure are necessary: • dust collection, through impactors or dry deposition with the opportunity to recognize by SEM the presence of dust particles; and • optical measurements to derive related properties regarding AOD and SAE.
The integrated use of available model output and satellite measurement allow completion the method here proposed.

Results and Discussion
In this section, we show results of statistics applied to selected cases by adopting the methodology defined in Section 3 according to seasonality and prevailing wind direction. Furthermore, analysis of large-scale patterns and the selected case of dust outbreaks is described in a detailed comprehensive analysis of measures and tools described in Section 2.

General Features
To apply our methodology to dust-outbreak case selection, we analyzed the statistical results on the entire 2015-2019 dataset considering SWAM OPC, MAAP and nephelometer and quality check of their measures as shown in Table 4. Regarding nephelometer data, availability of the initial 87.48% was reduced to 49.35% as quality check of data about relative humidity caused a high percentage loss of data, as described in [32].
Statistics of aerosol optical properties, eBC mass concentrations and PM coarse fraction (mean, median, minimum and maximum values, 10th, 90th percentiles) over the whole measurement period (2015-2019), is presented in Table 5 for LMT Observatory. Table 5   Over these five years of measurement, the highest median (1.94) with respect to the mean (1.84) values of SAE 450,700 were observed and the asymmetry factor g median values ranged between 0.49 g 700 and 0.54 g 550 . The observatory site shows high mean scattering values >40 Mm −1 and SAE values greater than one. The site, with a strategic position in the central Mediterranean Basin, depending on the local wind circulation and breeze system (west-east), is influenced by multiple sources with emissions of both coarse and fine fraction of aerosols. African dust events, sea spray intrusion, and episodic dust coming from Etna and Stromboli volcanic eruptions influence the area of study [74]. Furthermore, a large contribute of fine particles is evident (SAE 450,700 > 1) due to the near city and airport area [69,74,75]. Similarly, we analyzed eBC parameters obtaining a mean value of 0.643 µg/m 3 and a maximum value of 11.266 µg/m 3 resulted as shown in Table 5. To quantify the actual amount of dust, and other solid particles, we used the SWAM instruments collecting PM coarse fraction at the surface. For sample collection by SWAM, a quartz filter is used for 24h sampling in spring and summer seasons (when events are most frequent) and a 72h sampled in other seasons. The material collected can also be used for further analysis including chemical ones. Regarding the PM coarse and fine fraction, to homogenize the dataset due to different sampling time intervals, we calculate the daily weighted mean of the samples. For the whole period, the mean values of PM coarse fraction results are 10.69 µg/m 3 as in Table 5.

Climatology of Dust Events and Seasonality
In Table 6 we report the statistics of all dust events that occurred during the considered period (2015-2019). We note that the mean value of SAE coefficient is 0.614 ± 0.012, lower than the average SAE value calculated for whole period, while the g values is about 0.596 ± 0.001, and does not change much compared to the average calculated on 5 years according [28,50,64,69]. Regarding the σsc mean value, it is higher than the average over 5 years, confirming a greater presence of coarse particles in the atmosphere during the dust event. The contribution to the measured concentrations and the duration of the events are highly variable. Selected events present difference considering a seasonal analysis whose statistics are presented in Tables 7 and 8. During the cold season, the mean values of σsc, SAE and g in Table 7 are lower than the average on all dust events in Table 6.
Furthermore, analyzing the warm season in Table 8, an increase in SAE 450,700 , g and σsc compared to the average values described in Table 6 is observed. These results suggest that during dust event occurrences, we should analyze deeper optical parameters that are influenced by the circulation of meteorological fields due to different possible sources hence to the presence of mixed aerosol. This result explains the highest SAE values, compared to the averaged values present in Table 6, when the dust is transported from the S-E. In this case, dust particles are mixed with continental source aerosols.
Often, the synoptic configurations drive the air masses to flow together from SW and SE (see next section), generating a single convergent trajectory over LMT. During such an event, transported dust particles are mixed from the two opposite directions and the information regarding their source is lost.
In Table 9 we show the mean values of aerosol optical properties at surface. Here, 6 events with dust in the SW sector and 9 events in the SE sector are evident. The mean values of SAE (0.4) and g (0.570), considering only the SW sector, are in agreement with several other works [64,76,77].
Analyzing the dust events from the SE sectors, the SAE and g values increase to confirm change in the shape and size of dust particles according to [69].  Table 9. Statistics (σsc, σbsc, SAE, g), at LMT Observatory, all selected dust events with South-West and South-Est wind prevalent direction.

Large-Scale Atmospheric Patterns
A detailed analysis was carried out on the single identified event, producing and analyzing the upper-level and surface maps of the main atmospheric fields derived from the ERA5 reanalysis.
From an atmospheric point of view, wanting to analyze the synoptic forcing, we started considering 15 "dust events". As previously seen, these cases led to multiple and contiguous daily and 3-day dust events. To confirm that the 15 "cases" are, really, representative of all the 52 days (over 5 years), all daily maps were subjectively analyzed (not shown for brevity). After this first check, we confirmed the validity of our methodology and the mean values of the large-scale patterns were calculated and analyzed for the 15 "dust events".
We associated each of the 15 "dust events" with the large-scale ERA5 fields at 00 UTC on the first day of the considered period (see Table 2 for the dates). The analyzed parameters are geopotential height at 500 hPa (Z500), mean sea level pressure (MSLP), and U and V wind component ad 925 hPa (in order to calculate wind speed and direction at this level; WIND925). Once the single parameters were considered, we performed a spatial average for the all "dust events" (15 cases, 00 UTC; see Table 2 for the dates considered for the average), over a sufficiently large domain (see Figures 2-4) that can be considered representative of dust phenomena in southern Italy.
After this first check, we confirmed the validity of our methodology and the mean values of the large-scale patterns were calculated and analyzed for the 15 "dust events"; to do this, we performed a spatial average of the atmospheric parameters for all cases (fields at 00 UTC; see tab 2 for the dates considered for the average), over a sufficiently large domain (-20W 50E 0N 60N) that can be considered representative of dust phenomena in southern Italy.
In Figures 2-4 we show the averaged Z500, MSLP, and WIND925 (speed and direction), respectively. The averaged geopotential height at 500 hPa is shown in Figure 2. The domain appears sharply divided into two distinct parts, with opposite baric values. The demarcation is denoted by the zone of latitude between 20 • and 35 • . Only in this area are the contour lines are relatively dense, demonstrating the occurrence of intense upper-level winds. Elsewhere, upper-level leveled pressure fields and consequent weak winds prevail.
Upper-level winds in southern Italy and, thus, over Calabria, are guided by a pronounced high-pressure ridge over Greece/southeastern Europe. Following this configuration, the currents come from the Sahara Desert area towards our site of interest (from south-west to north-east). The averaged mean sea level pressure is shown in Figure 3. During the period taken into account, the mean transport pathway of the individuated African dust events towards southern Italy is mainly determined by a low-pressure system over the central Mediterranean Basin, which directly affects the Italian peninsula, and by two high-pressure systems over southeastern Europe (most importantly) and in the Atlantic. The surface circulation driven by these configurations facilitates the coming of southerly winds from North Africa toward southern Italy. In particular, the counter-clockwise circulation of winds around shows how the currents mostly flow from the main source of African dust, the Sahara Desert. At the same time, the high pressure over southeastern Europe facilitates the flowing of winds from north-eastern African regions.
An obvious condition for the existence of a dust source is the availability of fine-grained material, which can be lifted from the ground when the surface wind speed exceeds a certain threshold. In the literature it is found to vary from about 6 to 9 m/s; (e.g., [78]). In our case, the winds reached values higher than this threshold (Figure 4), further confirming that these were clear dust events. The averaged wind field at 925 hPa (Figure 4), which is directly connected to the dust emission and transport, reveals an important occurrence for our study cases and a specific peculiarity of the dust events in southern Italy. The currents follow, on average, two different (but confluent) prevalent directions. The first one is the prevalent mechanism responsible for the detected (at Lamezia Terme site) Saharan dust events, with currents directed from SW to NE. Another important contribution is represented by currents directed from south-east to north/north-west towards the Lamezia Terme site. These individuated pathways are in clear agreement with the mean sea level pressure fields shown in Figure 3. It is interesting to note the merging of these two pathways in the southern Mediterranean, and the subsequent flow together in southern Italy. This, considering the averaged values, makes the dust coming from the Sahara Desert indistinguishable from the dust coming from north-eastern African regions.
The desert dust events individuated in this work are mainly characterized by the synoptic patterns described above. Following these forcings, dust plumes were transported northward, resulting in intense dust events over the Mediterranean and Southern Italy.
Analyzing air mass circulation, it arises that most of the dust transported to our site, during the hot season, is associated with wind fields from the south-east direction, while in the cold season the prevalence of dust intrusions occurs in S-W air masses incursion. To confirm our results, we identified all the events according to which the wind fields at 925hPa flows exclusively from the south-west sector. In such a case, Saharan contribution we observe at LMT is evident, and, when the wind flows from the SE sector mixing particle occurrence is more evident.

Case Study: The Dust Outbreak 23-27 April 2019
In this section, we present a selected dust-outbreak case that is the most significant that occurred at LMT. The case selection and the related methodology was presented in Section 2.3. Here we show results of a detailed analysis in which all described instruments were operating simultaneously and filtered SEM analysis was available for one of the most high values of PM 10 registered at the site. Model and satellite output and images were also available for the period.

Derived Intensive Parameters
In addition to the direct measurements in situ, already analyzed, during the selected case studied we calculated their intensive aerosol parameters and we used sun photometer data, in particular SAE values available from the AERONET site. For the nephelometer measurements, SAE is calculated using Equation (3). Similarly, for column-integrated measurements, we used SAE in the spectral range 440-870 nm.
During intense dust outbreaks, a reduction of SAE 450,700 and a simultaneous increasing of PM coarse fraction is evident; this depends on the degree of local pollution and intensity of the dust event [79].  We observed a delay from the column-integrated values obtained from the sunphotometer with respect to in situ surface SAE 450,700 values measured from nephelometer. This delay corresponds to the permanence at the surface of a dust plume over LMT station; in fact, on 26 and 27 April, we still observed high daily mean concentration of PM 10 with values 46.1 µg/m 3 (from SWAM instruments) and 26.6µg/m 3 respectively, as shown in Figure 6. Both used methodologies (in situ and sunphotometer) calculate values close to 0.5 (even close to zero for columnar measures) and this is an indication of coarse particle predominance. The difference between columnar and in situ SAE measurements, during the dust event, suggests a stronger predominance of coarse particles in the atmospheric column (SAE 440,870 close to 0) [47] than on the ground (SAE 450,700 < 0.7) due to the influence of continental particles mixed with dust, as analyzed later with SEM in Figure 7.
To determine the starting, ending, and peak times for PM 10 concentration levels, we analyzed the temporal behavior of the PM 10 values derived from daily mean OPC and SWAM instruments at LMT. On 25th and 26th April 2019, the highest daily PM 10 values were recorded at LMT station ( Figure 6).
During these days SWAM measured one of the most relevant concentrations, as depicted in Figure 6, with respect to the mean PM 10 values during the five years. The statistics is shown in Table 10. Also, the maximum value of eBC concentration mass in this period is 0.05 µg/m 3 , which excludes a possible contribution of eBC due to a biomass burning. When compared with measurements from the previous day, daily mean PM 10 values turned from 7.71 µg/m 3 to maximum peaks values during 25 and 26 on April (50.4 µg/m 3 ). This variation, according to the behavior of SAE 450,700 coefficients, shows a minimum peak around 25 April. The highest values of SAE 450,700 (>0.65) during the days of the intrusion of dust show an important component of fine particles due to the mixing aerosols that we measured at surface (Figure 6). In Figure 6 the OPC fine and coarse contribution whose peaks correspond respectively to 31.97 µg/m 3 and 61.15 g/m 3 is depicted. However, the dust contribution at the surface aerosol sites is remarkable.
The comparison between the observed parameters mainly conductive of dust events collected at the Lamezia Terme experimental site, with different types of model outputs, were confirmed by other analysis of a sample collected on filters. This allows an investigation of the model's ability to reproduce both the horizontal and the vertical displacement of the dust plume, as well as its evolution in time, taking into account the large-scale atmospheric forcing in driving the dust event. Regarding the sampled filters the ratio Al/Si, for three selected dust episodes (see Table 11), was higher than 0.3 and this condition is always indicative of a desert origin of the particles [80]. In the literature, common mineralogical markers are often used to identify the origin of particulates. The prevalence of Si and Al, during the in-dust days, ascribes a prevailing Saharan origin to the particles [81]. Moreover, elevated values of the Al/Si ratio was used to discriminate the aerosol origins [82,83]. Regarding sampled filters for the Al/Si ratio, for three selected dust episodes (see Table 11), the result were higher than 0.3 and this condition is always indicative of a desert origin of the particles [80]. SEM-EDX analysis on selected samples reveals the presence of Si, Al, Mg, K, Ca and Fe elements. Significant amounts of Na and Cl elements were also detected in the two filters sampled in the spring season for the following two events: (14)(15)(16)(17) and (25 April 19), where the presence of NaCl salt crystals reveals the contribution of marine sea spray due to local circulations (see Figure 8). The SEM-EDX spectra of the single dust particles shows a mixed elemental composition, not homogeneous and spherical shape as in Figure 7 in agreement with the studies by [84][85][86], in which dust particles transported from arid regions are often aggregated.

Synoptic Features
In Figures 9-11 we show Z500, MSLP, and WIND925 (speed and direction), at the 00 UTC on 23 April 2019.
From a meteorological point of view this day, and the selected time, are indicative of the start of the dust event.
The geopotential height map at 500 hPa ( Figure 9) shows a strong similarity with the averaged field mentioned above (Figures 10 and 11); this confirms that the event is representative of the all dust events recorded at the experimental site.
A pronounced high-pressure ridge over the Aegean drives the upper-level currents to flow from the Sahara Desert area towards southern Italy. At the surface, the mean sea level pressure ( Figure 10) slightly differs from the average, but the main patterns are confirmed. In this case, the low-pressure area on the Mediterranean (gulf of Genoa) is less deep than the average (see Figure 10), while a wide area of high pressure is very pronounced over Eastern Europe. The trajectory of the surface winds, as confirmed by the denser isobars along the western side of the high, is more shifted to east with respect to the averaged fields, and the associated wind speeds are more intense. This fact is confirmed by the wind field map (Figure 11), at 925 hPa that shows the existence of moderate-to-high winds that flow across a pathway oriented from south to north, with values up to 25 m/s. Despite the orientation of the pathway, we can see that in northern Africa there are two different pathways coming, respectively, from the Sahara Desert and from northwestern regions; also, in this case the merging of these two pathways in the Ionian Sea and the consequent flow together in southern Italy is apparent, although the trajectory is slightly shifted to the east with respect to the average fields.

Satellite and Model Products
The Moderate Resolution Imaging Spectrometer (MODIS) on NASA's Aqua and Terra satellite permits confirmation of the characteristics of the event.
The MODIS value-added AOD is shown in Figure 12, superimposed on the MODIS Corrected Reflectance, for 25 April 2019, in a central day of the considered dust event period.
It is important to highlight that this product is a result of a specific retrieval method and is not a direct measurement. Furthermore, the aerosol forecasting community largely uses the product to quantitatively describe some specific events with high-resolution satellite products. As we can see, southern Italy is affected by a high value of AOD, in particular on 25 April. Observing the positions of the higher AOD values (reddish colors), and observing its evolution in the whole period (all daily maps, not reported for brevity), we deduce that the storm comes from the Sahara and moves east-northeast, according to synoptic analysis and with CAMS-modeled AOD as discussed below.
To evaluate the location of the dust sources that are directly connected to the investigated dust intrusion, we show the CAMS-modeled dust AOD at 550 nm. These fields, despite being modeling products, could be useful also to confirm the meteorological analysis discussed above and the selection of the events to study in more detail.
CAMS analysis confirmed the severity of this dust event. In particular, in Figure 13 the 24-h time interval simulated AOD (from 00 UTC on 22 April to 00 UTC on 26 April) is reported. Although in the previous section (ERA5 analysis) we analyzed the fields starting from 00 UTC on 23 April (the first day of the recognized event), in this part we show the CAMS outputs starting from the day before, in order to highlight the source of the dust. Values higher than 1.6 were reached on 23 April over southern Italy, but larger values (about 2) were simulated in Northern Africa and in the Mediterranean the day before.
Following the spatio-temporal evolution of the simulated dust storm, we can recognize the Saharan nature of the case, mainly on the first phase of the event. In a second phase, large values of AOD are also measured at Lamezia Terme (according to CAMS simulations, Figure 13), but in this case, simultaneously to the eastward translation of the synoptic system, even the dust source appears to be more shifted to the central/eastern African regions.
In Figure 14 the Total AOD at 00 UTC on 23 April, derived by CAMS, is also reported; this map is directly comparable with the Figure 13 referring to the same time. From this map we can see the high values of the Total AOD (>1.8) also in southern Italy, proving the existence of a small anthropogenic contribution to the aerosol detection at the site of interest. Comparing the MODIS Combined Value-Added AOD (Figure 12, referred to on 25 April) and the CAMS Dust AOD at 00 UTC on 25 April, we note a slight difference in the values and position of the larger AOD values. It is important to highlight that none of the two products is a direct observation and, moreover, that the satellite layer is referred to for the whole day, while the CAMS analysis is indicative of a single time (00 UTC in this case).  One of the principal results obtained in this preliminary analysis, regarding the main directions of the dust intrusions, regards the recognized existence of two pathways in the southern Mediterranean, and the subsequent fusion of them, that makes the dust sources (Sahara Desert or north-eastern African regions) often indistinguishable. This peculiarity is unique in the context of similar studies on Mediterranean regions, and characterizes the experimental site of Lamezia Terme.
This dust-outbreak case is a typical example of this occurrence; to confirm this, we show the HYSPLIT [54] back-trajectories, computed for 48 h run, for 23-25 April 2019, respectively (Figure 15a-c); the trajectories are forced with the archived NCEP Global Data Assimilation System (GDAS), at 0.5 resolution. The trajectories are obtained using the READY (Real-time Environmental Applications and Display sYstem) platform (https://www.ready.noaa.gov/index.php) [55] and are forced with the archived NCEP Global Data Assimilation System (GDAS), at 0.5 • resolution. In Figure 15a relating to 23 April, the dust is coming from the Sahara Desert, as confirmed in the previous discussions done with CAMS and MODIS outputs. As anticipated, in the following days, synergistically with the eastward movement of the storm, the dust becomes advected from south-southeast, thus highlighting the predominance of the second pathway in forcing the dust from central/north-eastern Africa regions.

Conclusive Remarks and Future Works
Here, we presented a comprehensive study of aerosol optical properties, performing statistical analysis integrated with modeled large-scale atmospheric fields to characterize dust events at our LMT GAW Regional site.
The paper focused on an incremental method to select dust-outbreak events that occurred in 2015-2019 period. We adopted thresholds based on a basic statistical analysis on the entire dataset, for aerosol optical properties, and collected particles in order to help us to efficiently identify dust outbreaks. The dust events were identified by the intersection between two different datasets using, respectively, SWAM and nephelometer data. SWAM data allowed us to select days and events (grouping at least three consecutive days) by choosing collected filters that present the double of the mean value (i.e., the 95th percentile) of particle concentration corresponding to 10.69 µg/cm 3 , this represents the aerosol mass threshold in which desert dust transport has a high probability to occur, and in our case this is confirmed by SEM observation. By optical parameters, obtained from nephelometer, we choose within the already selected days, through SWAM-considered mass threshold, those days in which SAE results were lower than 0.65. This selection technique plausibly reduces the total number of events that occurred, but allows us to keep track of only those occurrences objectively recorded by the sensors, according to the technique adopted. A total of 62 days of possible dust events were first identified by SWAM analysis. In 10 of these days, the values of eBC concentration are greater than a threshold choice of the 90th percentile of the daily mean (1.01 µg/cm 3 ), and this indicates the presence of local pollution; therefore, these days were discarded (52 days are retained). SWAM sampling was sometimes done with daily (one-day) collection times and, other times, with three-day collection times. By unifying all the three-day collection times, the sample was reduced to 26 sampling periods, now consisting of both one-day and three-day collection times. These sampling periods are, in some cases, contiguous. By further merging the all periods with contiguous dates, thus associating the same plausible large-scale atmospheric configuration to the single dust event, the sample is further reduced to 15 dust periods.
From analysis of data results that the LMT site is subject to frequent dust incursions and the most significant events occur during warm and cold seasons.
The collected dust particles generally present a coarse grain characterization with optical mean values that do not present relevant difference.
Due to the peculiarity of our site, dust outbreaks are related to two different pathways coming, respectively, from the Sahara Desert and from the northwestern African regions.
A difference exists for events coming from the two main directions in which western storms transport coarse dust particles, while eastern storms present mixed resulting aerosols that cross anthropized surroundings.
The analyzed case event, using all available instruments at LMT and by comparing the observed parameters with different types of model outputs and confirmed by other analysis on sample collected on filters, allows us to investigate the model's ability to reproduce both the horizontal and the vertical displacement of the dust plume, as well as its evolution in time, taking in a major account the large-scale atmospheric forcing in driving the dust event.
In this last case, the convergence of two air masses coming from both privileged directions supports our results: • SAE parameter with value close to 0.65, shows that the PM coarse fraction at the surface is mixed with fine/coarse due to the influence of the transport of continental aerosol; • g parameter does not show a significant change, but there is a small increase when we considering the S-E event sector; • filter images by SEM demonstrate the presence of mixed particles.
The ERA5 dataset was found to reproduce the synoptic meteorological conditions driving the dust outbreak well. Analysis of future and further dust outbreaks will give us a new chance to test further our incremental method.
Future work is needed, mainly related to mineralogical studies, to understand if this circulation pattern, that sometimes occurs, transports mass from both main dust sources, the Sahara Desert and Middle East.

Conflicts of Interest:
The authors declare no conflict of interest.