Abstract

Dust intrusions from African desert regions have an impact on the Mediterranean Basin (MB), as they cause an anomalous increase of aerosol concentrations in the tropospheric column and often an increase of particulate matter at the ground level. To estimate the Saharan dust contribution to PM10, a significant dust intrusion event that occurred in June 2006 is investigated, joining numerical simulations and specific measurements. As a first step, a synoptic analysis of this episode is performed. Such analysis, based only on meteorological and aerosol optical thickness observations, does not allow the assessment of exhaustive informations. In fact, it is not possible to distinguish dust outbreaks transported above the boundary layer without any impact at the ground level from those causing deposition. The approach proposed in this work applies an ad hoc model chain to describe emission, transport and deposition dynamics. Furthermore, physical and chemical analyses (PIXE analysis and ion chromatography) were used to measure the concentration of all soil-related elements to quantify the contribution of dust particles to PM10. The comparison between simulation results and in-situ measurements show a satisfying agreement, and supports the effectiveness of the model chain to estimate the Saharan dust contribution at ground level.

1. Introduction

Aerosols have direct and indirect effects on global climate [1], altering the radiative balance of the Earth-atmosphere system [24] and changing the microphysical and the radiative properties of clouds. In fact, aerosol particles can act as cloud condensation nuclei, modifying cloud lifetime and amount [35]. The mineral dust from desert areas, which represent a great source of aerosol injected into the atmosphere, suppresses precipitation in thin low-altitude clouds [68]. In addition, dust deposition can modify the ocean biogeochemical cycle, providing an important source of micronutrients [9]. It has also an impact on terrestrial ecosystems, providing nutrients as phosphorus to the soil [10]. Saharan desert is one of the most important sources of mineral dust [11], which has a considerable impact both at the global [12] and at the regional scale, as the Mediterranean Basin (MB) [13].

More specifically, Saharan dust outbreaks can cause an anomalous rise of PM10 concentrations over large parts of the MB, with impacts on Spain, Greece, and Italy [1416]. The exceedance of daily thresholds established by the European Union (Directive 2008/50/EC) may depend not only on human emissions but also on the mineral dust contribution. The same directive includes the possibility of removing exceedances due to natural sources, such as atmospheric resuspension or the transport of natural particles from dry regions. Although the assessment of a recognised procedure to quantify and subtract the effective Saharan dust contribution to PM10 concentration levels is still under discussion, preliminary guidelines are proposed in a working paper of the European Commission (http://www.ec.europa.eu/environment/air/quality/legislation/pdf/sec_2011_0208.pdf), mainly following the approach reported in Escudero et al. [17].

The method proposed in the European guidelines suggests to combine the use of satellite retrievals, model systems, and PM10 measurements at the ground level. The meteorological analysis (ECMWF http://www.ecmwf.int/) checks the evolution of synoptic systems prone to the generation of strong gradient winds able to transport dust plumes from North Africa towards Europe. The satellite observations (http://oceancolor.gsfc.nasa.gov/SeaWiFS/HTML/dust.html) and the aerosol index (TOMS http://toms.gsfc.nasa.gov/ozone/ozone_v8.html) can be useful in the detection of the dust episode and in the identification of the area affected by the plumes. Daily results of numerical models (e.g., SKIRON http://forecast.uoa.gr/, BSC-DREAM http://www.bsc.es/projects/earthscience/BSC-DREAM/ and NAAPs http://www.nrlmry.navy.mil/aerosol/) can be checked to identify the occurrence and the duration of the dust episode. Using only satellite observations and model simulations of aerosol optical thickness, it is not possible to clearly distinguish episodes that involve dust transport above the boundary layer without any impact at the ground level from those that cause dust deposition. For this reason, the European guidelines also require an analysis of PM10 mass concentration, which is routinely measured at the ground level by air quality station networks. Distinguishing between the Saharan dust contribution is not trivial since the PM10 concentration may be contributed by many different sources [16]. In a near future, the results of this method will be provided by GMES-Atmosphere service (http://www.gmes-atmosphere.eu/).

Within this framework, a study of Saharan dust intrusions over the MB is proposed, focusing on the deposition mechanism to assess the dust impact at the ground level. For this reason, a comprehensive model chain is configured to reconstruct the dynamic evolution of the Saharan desert dust. Numerical models play an important role in the description of the process of dust emission, transport, and deposition from desert zones. It is possible to reconstruct the spatial distributions of the dust with a good description of the vertical concentration profile. To characterise the PM10 mass concentration and composition at the ground level, the model-based analysis is complemented by physical and chemical analyses of the PM10 sampled in different sites. In particular, the concentration of all soil-related elements is measured by PIXE (Particle-Induced X-ray Emission) to assess the mineral dust concentration. The ionic composition is also used as ancillary data to subtract the sea-salt contributions to Na and Mg. Because the impact of desert dust is well characterised by an increase of all soil-related elements (and by changes in elemental ratios), this approach is useful for assessing the real impact of dust episodes on PM10. The measurement of all crustal elements may allow a quantitative assessment of the desert dust contribution, which is more accurate with respect to estimations that may be obtained simply by the analysis of PM10 mass concentration data. In addition to the specific quantitative information obtained by the PM10 measurements in the selected sampling sites, the integrated use of model simulations and in situ experimental data may provide an overall picture of the mineral dust distribution.

This approach, which joins numerical simulations and specific measurements, is employed to study the impact of Saharan dust event of June 2006, a month characterised by significant dust outbreaks over the MB [18, 19]. Since at the moment, the GMS-Atmosphere service is not operational, a preliminary analysis of this long spell is conducted using the available data as the synoptic meteorological data from NCEP, the satellite observation images (MODIS/AQUA, MODIS/TERRA http://modis.gsfc.nasa.gov/), and the aerosol index (TOMS). The global model GOCART (http://disc.sci.gsfc.nasa.gov/gesNews/gocart_data_V006), developed by the Georgia Institute of Technology-Goddard to simulate aerosol optical thickness [20, 21], is used to integrate the satellite data. The HYSPLIT model (Hybrid Single Particle Lagrangian Integrated Trajectory model, http://ready.arl.noaa.gov/HYSPLIT.php) is used to generate backward trajectories to trace back the sources of the air masses at different levels and for different hours of the day.

The results of model chain, configured to describe the spatiotemporal evolution of this specific dust outbreak, are compared with GOCART data. Physical and chemical analyses are conducted on PM10 samples collected from three sites in central Italy (Tuscany). Finally, the numerical simulation results are compared with these specific measurements to evaluate the effectiveness of the model chain in the quantitative estimation of the Saharan dust contribution at the ground level.

2. Instruments and Methods

2.1. The Model Chain

The natural phenomena involved in the dust cycle in the atmosphere consist of two major physical processes: a wind stress lifting mechanism enabling the dust particles to rise up from bare soil surfaces and then the transport and deposition of this mineral dust [2224]. To provide regional characteristics of Saharan dust intrusion over the MB, an atmospheric emission and dispersion model chain is developed. The model chain is based on three different modules: the atmospheric, dust emission, and transport/deposition modules. The Regional Atmospheric Modeling System (RAMS) [19], used by CNR-IBIMET in operational mode [2527], provides the input data for the other modules. The DUST Emission Model (DUSTEM), specifically developed for this aim, simulates the dust emission from the desert. The Comprehensive Air quality Model with extensions (CAMx) (http://www.camx.com/home.aspx) takes the meteorological inputs from RAMS and the emission rate from DUSTEM [28], providing the dynamical transport and deposition of the dust particles.

This section describes the model chain configuration used for this case study from 1 June to 5 July 2006.

2.2. The Regional Atmospheric Modeling System (RAMS)

The RAMS 6.02 version is run over a domain including a large part of the Northern Hemisphere and performed by parallel computing. The initial and boundary atmospheric conditions need to be set as well as the forcing data during the period of simulation. To address these necessities, the Reanalysis2 dataset [29], with a 2.5-degree horizontal resolution, is employed. In particular, the geopotential height, temperature, relative humidity, and zonal and meridional wind component fields are used and forced as boundary conditions every 6 hours throughout the simulation period. The Kain-Fritsch convective scheme is also adopted. The configuration is set with the following features: the domain is centered on 40°N–5°E, with 200 × 80 grid points, 32 sigma vertical levels (with a stretching factor to obtain a greater resolution near the soil and a lesser one above 2000 m) and 11 soil levels. The horizontal resolution is 0.54 degrees, and the time step is 120 seconds, with a temporal output of 1 hour.

2.3. The Dust Emission Model (DUSTEM)

Dust emissions are estimated by developing an ad hoc model called the DUST Emission Model (DUSTEM) [28]. This model estimates dust emission rates using empirical relationships based on soil texture and friction velocity [30]. DUSTEM can take into account four different soil types: clay, small silt, large silt, and sand [30, 31]. The main features of these classes are shown in Table 1. For the selected case study, only the first two dimensional classes, clay and small silt, are considered because they are the only ones involved in long-range transport [31]. DUSTEM takes the soil information from the GLC2000 land cover [32] and FAO Textural Map [33] as input data to obtaining the bare soil map. The hourly meteorological fields (soil moisture and friction velocity) are provided by the RAMS model. The computational domain, in Polar Stereographic coordinates (with the pole at 40 degrees north and 5 degrees east), is formed by 380 × 340 cells with a 30 km resolution to provide the input data to CAMx with a temporal resolution of 1 hour. The maps of emission rates cumulated on the whole simulation period, relative to the first and second soil types, are shown in Figure 1.

2.4. The Comprehensive Air Quality Model (CAMx)

The Comprehensive Air quality Model (CAMx), developed by ENVIRON International Corporation California, is an Eulerian photochemical dispersion model that allows for an integrated “one-atmosphere” assessment of gaseous and particulate air pollution over many scales, ranging from urban to superregional (http://www.camx.com/home.aspx). CAMx simulates the dispersion, chemical reactions, and removal of pollutants in the troposphere by solving the pollutant continuity equation for each species. The aerosol deposition is handled by adopting the algorithm of S. A. Slinn and W. G. N. Slinn [34] and Seinfeld and Pandis [35] for the dry and wet deposition, respectively. To simulate the transport and deposition of mineral dust, the chemical module is switched off in this study. The meteorological input data are provided by RAMS, and the emission rates of clay and small silt types are provided by DUSTEM (Table 1). The initial and boundary conditions are set to zero. To consider the atmospheric dust loading, the simulation starts on 1 June, some days in advance of the outbreak over Europe. The horizontal computational domain is the same as that used for DUSTEM. There are 18 vertical levels, from 10 m to 10,500 m, with a finer resolution near the ground. The concentration outputs for clay and small silt soil types are provided with a temporal resolution of 1 hour.

2.5. PM10 Sampling and Analysis

To characterise the PM10 mass concentration and composition at the ground level, the concentration of all soil-related elements is measured by PIXE to assess the mineral dust concentration [3640]. The simultaneous detection, with high sensitivity, of all the elements that compose mineral dust makes PIXE highly effective for these investigations. A simultaneous increase in the concentration of all crustal elements is indeed a first indication of the occurrence of dust events; clearly, it is a stronger signature than the increase of the PM mass concentration, which may be due to many other sources, both natural and anthropogenic. However, an increase of crustal element concentrations may also be due to local soil dust resuspension: in this context, the analysis of changes in elemental ratios may help distinguish between these two kinds of events. Moreover, when the aerosol is sampled at more than one site, the observation of a simultaneous increase of soil-related elements at the different sites may also suggest the impact of long-range transported soil dust affecting a wide area, while the observation of different patterns from site to site clearly hints at the contribution of local sources.

2.6. PM10 Sampling

Aerosol samples are collected on 47 mm Teflon filters on a daily basis (from midnight to midnight) using low-volume (2.3 m3/h) sequential samplers (HYDRA Dual Sampler, FAI Instruments) equipped with a PM10 inlet in accordance with the European rule EN 12341.

The PM10 daily mass concentrations are obtained by weighing the filters with an analytical balance (sensitivity 1 μg) before and after the sampling, always after a storage period (48 hours) in a temperature and humidity controlled room with an ambient temperature 𝑇 = (20 ± 1)°C and a relative humidity RH = (50 ± 5)%. Electrostatic effects are avoided by using a deionising gun.

In particular, a comprehensive dataset of PM10 measured concentrations and composition at the ground level is obtained in the framework of the PATOS project, the first extensive field campaign for PM10 characterisation in Central Italy (Tuscany), which was supported by the Regional Government (http://servizi.regione.toscana.it/aria/index.php?idDocumento=18348). PM10 samples are collected from September 2005 to September 2006 at six sampling sites in Tuscany that are representative of different types of areas: Florence (urban background), Prato (urban traffic), Capannori-Lucca (urban background), Arezzo (urban traffic), Grosseto (urban background), and Livorno (suburban background). In this field campaign, the three available samplers are used to collect the PM10 data in the six sites previously cited. These samplers are moved every 15 days from three sampling sites (Arezzo, Prato, and Livorno) to the other three (Florence, Capannori-Lucca, and Grosseto). In this way, data relative to the long-range transport dust intrusion of June and the background level due to local dust resuspension is sampled.

2.7. PM10  Compositional Analysis

The PM10 samples are analysed by Particle-Induced X-ray Emission (PIXE) to determine the aerosol elemental composition and, in particular, the concentrations of all the main soil-related elements (Na, Mg, Al, Si, K, Ca, Ti, Mn, Fe, and Sr). The same samples are also analysed by ion chromatography (IC) to assess the soluble fraction of inorganic ions, specifically the Na+ and Mg2+ concentrations. These values are used to determine the sea-salt contributions to Na and Mg to be subtracted for the calculation of the soil dust fraction of these elements.

The PIXE measurements are performed at the 3 MV Tandetron accelerator of the LABEC laboratory of INFN in Florence using the external beam setup dedicated to environmental applications [39, 41]. For this case study, each sample is irradiated for 500 s with a 3.2 MeV proton beam with an intensity of 5 nA over a spot of 2 mm2. During irradiation, the filter is moved in front of the beam so that most of the area of the deposit is analysed. The PIXE spectra are fitted using the GUPIX code [42] and elemental concentrations are obtained via a calibration curve from a set of thin standards (Micromatter Inc.). The minimum detection limits (MDLs) are 10 ng/m3 for low-Z elements and 1 ng/m3 for medium-high-Z elements. The elemental concentration uncertainty is determined by taking into account the independent uncertainties on the standard sample thickness (±5%), the aerosol deposition area (±2%), the airflow (±2%), and the X-ray counting statistics (from ±2 to ±20% or higher when concentrations approach MDLs).

Ion analysis is performed on the water extract obtained from a quarter of each Teflon filter. Each sample is analysed for cations (Na+, NH4+, K+, Mg2+, and Ca2+), inorganic anions (F, Cl, NO3, and SO42), and some organic anions (methanesulphonate (MSA), acetate, formate, glycolate, and oxalate) by 3 Dionex ion chromatographic systems operating in parallel under the working conditions summarised by Becagli et al. [43]. The detection limits are several orders of magnitude lower than the concentrations found in these samples. The uncertainty is mainly determined by the ion chromatography accuracy, which is typically ±5%.

3. Results and Discussion

3.1. Synoptic Conditions Featuring the Dust Outbreak.

The synoptic atmospheric conditions favourable for dust raising are those associated with air mass advection coming from Northern Europe or the Balkans regions and heading towards Algeria, Libya, and Egypt up to Chad [44]. These air masses are characterised by strong and constant wind in the lower layers. The associated thermal gradient has two important effects for dust emissions: the first is the increase of the soil friction in correspondence with the air mass front due to the acceleration of the thermal wind; the latter is the reduction of the soil moisture content due to the positive surface, which favours the raising of dust.

Once the dust is raised above the boundary layer, a strong circulation between 850 and 700 hPa is necessary to transport the dust far from the emission area. Throughout the year, massive airborne plumes of desert dust from the Sahara and surrounding regions are exported to the tropical Atlantic Ocean and the Mediterranean Sea. The majority of the dust intrusions over the MB are usually associated with the passage of either a cold or a warm low-pressure system. Saharan depressions develop most readily when a Polar or Arctic air mass from the northwest (Maritime Polar or Maritime Arctic) or northeast (Continental Polar or Continental Arctic) flows behind the dry desert air [45].

To analyse the synoptic circulation over the Mediterranean area during the case study of June 2006, the sea level pressure, air temperature, and wind fields are investigated using a daily NCEP Reanalysis database provided by the NOAA/OAR/ESRL PSD (website at http://www.esrl.noaa.gov/psd/). The anomalies are calculated with respect to the long-term mean (LTM) period of 1981–2010. The screening of 6-hourly sea level pressures (SLPs), wind vectors, and air temperatures at different pressure levels reveals major dust outbreaks originating from the North African deserts towards Europe and provides insight into their temporal evolution. The means and anomalies of these 3 meteorological fields, compared with their LTMs, are calculated and discussed as follows.

From 10 to 21 June 2006, the average SLP distribution over the MB highlights a pronounced high-pressure system over the West Mediterranean Basin (4 mb in excess of the LTM for this period) and a low-pressure system over southern Algeria (Figure 2(a)). The West Mediterranean (WM), exposed to the influence of both systems, is subjected to the generation of SE gradient winds at 850 hPa, resulting in this episode being 6-7 m s−1 stronger than the normal wind flow over this region (Figure 2(b)). On 23 June, an Atlantic perturbation advances toward the Central Mediterranean area and weakens the SW flux and the dust transport over Italy. The following period, from 25 to 29 June 2006, is characterised by an extension of the Azorean ridge deep into northern Europe together with a weakening of the African Monsoon over southern Algeria (Figure 2(c)). This synoptic configuration yields a negative anomaly of 1-2 hPa over the whole MB (Figure 2(d)) with one of its two lowest pressure cores centred over the North Western African coast. This low-pressure system induces a cyclonic flow at shallow atmospheric layers over the WM, characterised by anomalous strong Westerlies over this region accompanied by stronger than usual SE winds over Sardinia and Corsica (Figure 2(e)). This counterclockwise flow led to a positive anomaly in air temperature over the WM and the Central Mediterranean (Figure 2(f)).

The satellite observations, as shown in Figure 3(a) for 15 June, and the TOMS aerosol index reveal air masses transporting significant amounts of dust over large parts of the WM originating from southern Algeria and SW Africa. The TOMS data and the AQUA satellite observations (Figure 3(b)) for 27 June 2006 show the plume of dust originating from SW Algeria migrating in a cyclonic direction path through Sicily and Central Italy, crossing Sardinia and Corsica on its way to the Genoa Bay and eastern Spain.

The GOCART dust aerosol optical depth (AOD) daily maps are also used to follow the evolution of the dust outbreak on a daily basis [20]. As an example, Figure 4 represents the AOD for 15, 21, 24, and 26 June, showing the dust transport over the MB.

Finally, the Lagrangian trajectory HYSPLIT (Hybrid Simple Particle Lagrangian Integrated Trajectory) model is used to generate backward trajectories to trace back the sources of the air masses at different levels and for different hours of the day. Figure 5 shows the trajectories ending in Florence on 29 June 2006 at three different levels of A.G.L. (500, 1000, and 1500 m) starting from Northern Africa.

3.2. Numerical Simulations of Mineral Dust Events.

The model chain is performed from 1 June to 5 July 2006. An analysis of the model results, based on the daily vertically integrated dust concentrations for the particle size 1–20 μm (as a sum of clay and small silt soil types), is presented. From 6 June, the Saharan desert dust is transported to the Iberian Peninsula, and in the following days, the dust reaches Northern Europe (France, UK, and Norway). From 16 to 22 June, a large area ranging from the West MB to the Scandinavian Peninsula and Russia is affected by a mineral dust intrusion that also reached Italy. In particular, during 17 and 18 June, Spain, France, Switzerland, Germany, Italy, and the Balkan Peninsula are affected by elevated dust concentrations. During 23 and 24 June, the Italian peninsula experienced a decrease of dust concentration due to a newly formed cyclone coming from the Atlantic Ocean, which interrupted the southwestern flux over Italy. Afterwards, the transport mechanism proceeds again until the first days of July, although the concentration values are lower than those in the previous period. Figure 6 shows the daily vertically integrated dust concentrations for 15, 21, 24, and 26 June 2006. The qualitative comparison between these simulation results, the satellite retrievals, and the TOMS index data show a good agreement with the location of the zones that encountered dust intrusions. Furthermore, the GOCART dust aerosol optical depth daily maps are compared with simulation’s results for the whole period, from 1 June to 5 July, to evaluate whether the model chain correctly reproduces the dust event. There is temporal and spatial agreement between the CAMx and GOCART results, as shown in Figures 4 and 6.

The RAMS/DUSTEM/CAMx model chain, configured with a finer resolution than GOCART, provides more detailed information on the vertical distribution of dust concentration and on the deposition at the ground level. In fact, the vertical sections of the mean daily dust concentration, for example at latitude 43.78°N (Florence), show dust advancing towards Spain and France, at both low and high levels up to 5000–7000 meters. The greatest dust concentration over Europe is reached during the period of 15–22 June between 1000 and 8000 metres, whereas under the boundary layer, the concentrations are lower (30–40 μg/m3 versus 200–500 μg/m3) (Figure 7). During 23-24 June, the Atlantic cyclonic flow temporarily interrupts the dust transport over Tuscany. On 25 June, the dust is transported only above the boundary layer without deposition. Observing the daily modeled dust concentration at the ground level, a strong decrease is evident, as also proved by ground measurements (next paragraphs). Such information cannot be obtained from the results based on the daily vertically integrated dust concentrations or from an analysis of the optical depth maps. Finally, in the last part of the episode (26 June–2 July), the vertical extent of the dust concentration is lower, reaching 5000–6000 meters, and the dust event is less intense (Figure 7).

3.3. Chemical and Physical Analysis

During the month of June 2006, the concentrations of all the main crustal elements (Na, Mg, Al, Si, K, Ca, Ti, Mn, Fe, and Sr) show a clear increase together with the WM Saharan intrusion simultaneously in all sampling sites in Tuscany (Figure 8). The concentrations start to increase from 18 June and return to background values on 2 July. The time series show a two-peak shape, with a central minimum (24-25 June). As mentioned previously, the samplers are moved from three sampling sites to the other three every 15 days: in June, these displacements occurred on the 14th and 29th. Consequently, most of the episode (18–28 June) is measured in Florence, Capannori-Lucca, and Grosseto, one day (28 June) is missed, and the last part of the event (30 June–2 July) is observed in Arezzo, Livorno, and Prato. The average concentration values of the crustal elements during the episode are 2–4 times higher than their background concentrations. Al, Si, and Ti show a higher increase: their maximum concentration values during the episode (2.5 μg/m3, 6.0 μg/m3, and 0.2 μg/m3) are approximately 6 times higher than their background values.

Elemental ratios among the crustal elements also show significant changes, thus further reinforcing the Saharan intrusion hypothesis. For example, the Si/Al ratio mean and standard deviation during the episode is 2.3 ± 0.1 in Florence, 2.4 ± 0.1 in Capannori-Lucca, and 2.2 ± 0.1 in Grosseto, while the same ratios calculated for the other days (i.e., excluding the days of the episode) are 15% higher (2.7 ± 0.2, 2.8 ± 0.2, and 2.5 ± 0.2, resp.); the Ti/Fe ratio during the episode is 0.08 ± 0.01 in Florence, 0.08 ± 0.01 in Capannori-Lucca, and 0.09 ± 0.01 in Grosseto, while it is 50–70% lower during the other days (0.04 ± 0.01, 0.05 ± 0.01, and 0.06 ± 0.01, resp.); the Al/Ca ratio during the episode is 0.41±0.07 in Florence, 0.49 ± 0.11 in Capannori-Lucca, and 0.57 ± 0.12 in Grosseto, while it is 40–80% lower during the other days (0.29 ± 0.06, 0.28 ± 0.07, and 0.32 ± 0.08, resp.). Similar results are found at two remote sites of Northern and Central Italy by Bonelli et al. [36]. In both studies, a decrease of the Si/Al ratio and an increase of the Ti/Fe and Al/Ca ratios are indicated as highly representative fingerprints of African dust transport events. Their reported values for the Si/Al ratio during these events (and during days not affected by Saharan intrusions), that is, 2.3-2.4 (2.7-2.8), match well with those found in this study. The Si/Al ratio decrease may be explained by a higher contribution during the Saharan intrusions of agglomerated clay minerals with respect to bulk crustal material, while the increase of the Ti/Fe and Al/Ca ratios may be ascribed to an iron and Ca enrichment in local soil dust.

An estimation of the soil dust component may be calculated as the sum of the contributions of the main oxides of all the crustal elements (Na2O2, MgO, SiO2, Al2O3, TiO2, K2O, CaO, and Fe2O3): [][][][][K][][][],soildust=1,7nssNa+1,67nssMg+1,89Al+2,14Si+1,2+1,4Ca+1,43Fe+1,67Ti(1) where [nssNa] and [nssMg] are the concentrations of “non-sea-salt” Na and Mg (i.e., excluding the sea-salt contribution).

This expression may overestimate soil dust due to the contribution of other sources, such as biomass burning (for K), traffic, and other anthropogenic sources for Fe, and Ca. During the studied period the enrichment factors of K, Fe, and Ca (with respect to Al using the crust composition reported in Mason [18]) turn out to be low, suggesting the absence of strong anthropogenic contributions. Consequently, no correction is applied to the calculation of their contribution to soil dust. The sea-salt fractions of Na and Mg are calculated using the IC data, as the measured Na+ and Mg2+ soluble ion concentrations may actually be considered good estimations for these contributions.

The daily time series of soil dust concentration, calculated as reported above, are shown in Figure 8. During the Saharan intrusion, values of 20–30 μg/m3 are reached, to be compared with background values of the order of a few μg/m3. The Saharan dust contribution, for all the days of the episode, is estimated by subtracting the background soil dust concentration, which is calculated, site by site, as an average from mid-May to the end of July, excluding the Saharan intrusion days. The results together with the PM10 concentration values are reported in Table 2. The Saharan dust contribution is quite high, with concentrations up to ~20 μg/m3, and it can be considered the main cause of the PM10 50 μg/m3 limit value exceedances (values shown in bold in the table) occurring in this period.

3.4. Comparison between Model and Measurements

The comparison with in situ measurements plays a fundamental role for an evaluation of the model chain effectiveness in reproducing the deposition mechanism.

As explained in Section 3.3, the in situ measurements identify the dust outbreak from 18 June to 1 July, with a small residual contribution on 2 July. For the comparison, only three sites are used (Firenze, Capannori-Lucca, and Grosseto), as the relative sampling periods cover almost the entire extent of the episode.

In Figure 9, the time series of the daily simulated dust concentration at the ground level are compared with the Saharan dust contribution to the PM10 obtained by the measurements. Concerning the period extension, in the inland cities (Florence and Capannori-Lucca), the model simulation identifies the dust episode from 16 June to 2 July. The onset of the dust outbreak is thus anticipated in 2 days with respect to the measured data, while in the coastal city (Grosseto), the simulation and the measurements identify the same temporal range.

At all the sites, both the simulated and observed time series show a good agreement in the representation of the concentration decrease on 24-25 June and the following increasing phase up to the end of the episode.

Because the DUSTEM emission model provides particle sizes of 1–20 μm, the comparison with PM10 measurements is limited to a qualitative analysis and may lead to simulated concentration values that are higher than the soil dust concentrations obtained from the measurements. This behaviour may actually be observed in the inland cities. However, in the coastal cities, the simulation gives lower values than the measurements, which could be attributed to the model representation of the boundary layer’s vertical extent, which could be underestimated near the coast, obstructing the dust intrusions from higher to lower levels. The meteorological model RAMS has been configured with a horizontal and vertical resolution useful to reproduce the raising and transport phenomena, but it is not fine enough to describe the boundary layer near the coastline. Perhaps it will be necessary to introduce a finer nested grid into the meteorological model (and, consequently, the CAMx model) to describe the boundary layer with a better resolution over the target area.

4. Conclusions

Within the framework of the European Commission’s guidelines to assess the natural contribution to PM10, the aim of this work is to study the Saharan dust intrusions over the MB to quantify the dust impact at the ground level. The selected case study considers the large dust outbreak over the MB of June 2006. After a preliminary synoptic analysis, an ad hoc model chain is configured and physical-chemical analyses are conducted on PM10 samples collected in Central Italy. The model chain properly reproduces the dust emission and transport dynamics, as proved by comparison with available model maps. Furthermore, it describes the dust distribution on vertical profiles and the deposition at the ground level, adding detail to the satellite observations and optical thickness data analysis. Regarding ground-based measurements, the specific assessment of all the crustal elements by PIXE gives a more accurate quantitative estimate of the desert dust contribution than that obtained solely by PM10 mass concentration. The comparison between simulation results and specific insitu measurements shows a satisfying overall agreement, especially for the temporal evolution of the studied episode. Nevertheless, in the sampling site near the coast, the model simulation gives lower values than the measured ones. Near the coast, the model representation of the boundary layer’s vertical extent could be underestimated, obstructing the dust intrusions to the ground level. To improve the performance of the model chain, it may be useful to introduce a nested grid in the meteorological model (and, consequently, in the CAMx model) with a higher resolution over the target area, where the measurement sites are located. However, the comparison with PM10 measurements is limited to a qualitative analysis because the DUSTEM model provides emission rates for particle sizes of 1–20 μm. As a further development, the DUSTEM will be improved, taking into account a better soil type description to provide PM10 emission rates.

Acknowledgments

The authors would like to acknowledge the Tuscany Regional Government (Italy), which supported the PATOS project and provided the experimental data used here, and the support received within the CNR-Italy Short Term Mobility Program 2010. The authors would also like to acknowledge the NOAA Air Resources Laboratory for the provision of the HYSPLIT model results and NASA for the TOMS maps, satellite images, and the Giovanni online data system, developed and maintained by the NASA GES DISC for the dust maps used in this study.