Chemically aged and mixed aerosols over the Central Atlantic Ocean – potential impacts

Chemically aged and mixed aerosols over the Central Atlantic Ocean – potential impacts M. Astitha, G. Kallos, C. Spyrou, W. O’Hirok, J. Lelieveld, and H. A. C. Denier van der Gon University of Athens, School of Physics, Division of Environmental Physics-Meteorology, Bldg PHYS-V, 15784, Athens, Greece Energy, Environment and Water Research Centre, The Cyprus Institute, Nicosia, Cyprus Max Planck Institute for Chemistry, Becherweg 27, 55128 Mainz, Germany TNO Built Environment and Geosciences, Princetonlaan 6, P.O. Box 80015, 3584 CB Utrecht, The Netherlands


Introduction
Aerosols are complex mixtures of atmospheric constituents that play a role in many environmental and climate issues, ranging from air quality and human health impacts to potential effects on the hydrological cycle (Ramanathan et al., 2001;Lau and Kim, 2007).Aerosols originate from natural (desert dust, sea salt particles) and anthropogenic sources (urban, industry, biomass burning), and can travel long distances from Europe and Africa (Lelieveld et al., 2002) reaching the Central Atlantic Ocean with different age and composition (Prospero et al., 2001;Kaufman et al., 2002b;Formenti et al., 2003;Kallos et al., 2006Kallos et al., , 2007)).Chemical aging of dust aerosols occurs through heterogeneous interactions between particles and gases of different origin.The aged dust particles can have distinctly different physico-chemical properties and hygroscopicity than the emitted primary particles.For example, desert dust particles being initially not very soluble and ineffective cloud condensation nuclei (CCN) can become coated with soluble material turning them into effective CCN (Levin and Cotton, 2009).On the other hand, non-coated and less hygroscopic dust particles can act as ice nuclei (IN) (Levin and Cotton, 2009).Anthropogenically increased aerosol concentrations cause a direct radiative forcing of climate (Forster et al., 2007 and references therein) and indirect climate effects through CCN interactions with clouds also indirect climate effects (Ramanathan et al., 2001;Denman et al., 2007).The links between aerosols, radiation, clouds and precipitation are a main subject of theoretical, experimental and modelling studies during recent years.Many of these studies were summarized in the IPCC 4th Assessment Report (Forster et al., 2007); yet major uncertainties remain.
During the warm season, air masses from Europe and the Mediterranean region can be transported towards the Inter-Tropical Convergence Zone (ITCZ, located over N-Africa) and subsequently towards the Central Atlantic Ocean by the trade wind easterlies (Millan et al., 1997;Kallos et al., 1998;Gangoiti et al., 2006;Kallos et al., 2007 and references therein).The aerosols entering the Central Atlantic region, where intense convective activity can develop, may influence cloud formation through their CCN Introduction

Conclusions References
Tables Figures

Back Close
Full and IN activity depending on their number, size and chemical composition.Furthermore, particulate air pollution, including desert dust, has been indicated as potentially reducing the sea surface temperature (SST) in the Central and Northern Atlantic Ocean (Lau and Kim, 2007;Evan et al., 2009) and affecting tropical cyclone activity (Evan et al., 2006).Aerosol and CCN measurements near the Azores (Garrett and Hobbs, 1995) have shown that polluted continental air can alter cloud structures through enhancements in the CCN number concentrations and shifts in the droplet size spectra.More recent field campaigns, such as the African Monsoon Multidisciplinary Analysis (AMMA), were dedicated to understanding the role of the African monsoon (Redelsperger et al., 2006) in the West Africa region.A continuation of the AMMA project has been the field campaign in the East Atlantic (NASA-AMMA) to investigate the linkage between the African easterly waves, the Saharan Air Layer and tropical cyclogenesis (Zipser et al., 2009).Although this is a complex and controversial topic, encompassing a wide range of scales of atmospheric (micro)physics, chemistry and dynamics, there is general agreement that changes in aerosols can affect clouds and precipitation.The sign and extent of the aerosol influences is subject of ongoing field and modelling studies.The subject of this work is the characterization of the aerosols over the Central Atlantic Ocean in terms of number and size distribution and chemical composition.The methods used involve chemistry-transport modelling, the use of aerosol measurements and data retrieved from satellite observations, focusing on August 2005, a month with intense hurricane and tropical storm activity over the Atlantic Ocean.Implementing state-of-the-art regional modelling systems and satellite retrievals may be considered as a first step towards the characterization of aerosols and their possible climate effects.It should be noted that the modelling tools used here are not suited for studying the intricate mechanisms of aerosol/radiation/cloud processes, which therefore are not directly addressed.The focus of this work is on the chemical composition and number concentration of the aerosols entering the eastern Central Atlantic Ocean tropical storm genesis region.The models selected for this study, the simulation period chosen Introduction

Conclusions References
Tables Figures

Back Close
Full and the input data implemented in the model are detailed in Sect. 2. In the subsequent section, the results are presented and discussed focusing on the number distribution of natural and anthropogenic aerosols near the cloud formation areas, and the final section presents the conclusions.
2 Models and data used The modelling system used in the present study is the combination of the latest versions of the SKIRON/Dust model (Kallos et al., 1997(Kallos et al., , 2006(Kallos et al., , 2009;;Nickovic et al., 2001) and the Comprehensive Air Quality Model with Extensions-CAMx (Environ 2006), with additions/improvements to the aerosol schemes (Astitha and Kallos, 2009).A recently updated version of the SKIRON/Dust model includes the treatment of dust particles using a lognormal size distribution partitioned into eight size bins from 0.1 µm to 10 µm diameter (Kallos et al., 2009).It also includes the online calculation of aerosol optical depth and its influence on radiative transfer by utilizing look-up tables (Kaufman et al., 2002a).Since the dust module is radiatively-dynamically coupled to the atmospheric model, prognostic meteorological calculations are expected to realistically treat the injection of dust into the lower atmospheric levels.The performance of the SKIRON/Dust model has been tested in several recent studies (Kallos et al., 2006(Kallos et al., , 2007;;Astitha et al., 2008).
The Comprehensive Air Quality Model with Extensions-CAMx (Environ, 2006) is a Eulerian chemistry-transport model that includes gas and aerosol chemistry, aqueous phase processes, thermodynamic partitioning of the inorganic aerosols (Nenes et al., 1998)  model, impacts of dust optical depth on photochemical reactions, online production of sea salt and online heterogeneous reactions on the surface of dust particles.
An important new feature is the use of the desert dust fluxes from the SKIRON/Dust model.These fluxes are converted to emissions of crustal material and they are imported into the CAMx model for the simulation set up.The mentioned new features are described in detail in Astitha and Kallos (2009) and are not re-analyzed in this work.In Astitha and Kallos (2009) there is also a statistical analysis of the model performance for the on-line sea salt production.It should be noted that organic aerosols (primary and secondary) are not accounted for in this study.Even though we acknowledge their significance in the atmospheric composition (Kanakidou et al., 2000), in this study we emphasize inorganic aerosols and desert dust and their interactions.As this system is already complex to analyze, the inclusion of organics in the aerosol mixture is associated with large uncertainties and the analysis of their contribution is a subject for future work.
Information on the input data and configuration of the modelling system appear in The chemistry-transport simulations with the CAMx model used the 3-D meteorological fields from the SKIRON/Dust model as input (wind, temperature, pressure, water vapour, cloud, precipitation, diffusivities).The initial and boundary conditions are temporally and spatially constant concentrations based on climatological data.For the gasphase chemical reactions the CBIV (Carbon Bond IV) mechanism is chosen and for aerosol chemistry the multi-section scheme (CMU) is used that includes ISORROPIA Introduction

Conclusions References
Tables Figures

Back Close
Full and RADM for the thermodynamics and the aqueous phase processes.The aerosol particles are partitioned into three size sections: 0.03-0.1,0.1-2.5 and 2.5-10 µm.These sections were chosen in accordance with the sea salt production scheme, even though dust is produced from the SKIRON/Dust model using eight size bins.The emission fields necessary for the model simulations are obtained from various sources and are described below.
Because of the large spatial area required to perform the simulations, ranging from Northern Europe to Southern Africa and from the Arabian Peninsula to the Western Atlantic Ocean, it was necessary to derive the emission inventory from a combination of available sources.Specifically, the emissions for the European region are from TNO (Netherlands Organisation for Applied Scientific Research) in 0.25 • ×0.125 • horizontal analysis for the year 2000 (Visschedijk and Denier van der Gon, 2005); the emissions for the shipping sector are from the Convention on Long-range Transboundary Air Pollution-EMEP (Vestreng, 2003) and for the remaining areas they come from the Global Emissions Inventory Activity-GEIA database.To accomplish a uniform gridding of the different inventories, an interpolation scheme is used to re-grid the different resolutions of the databases.The emission input data include SO 2 , NH 3 , CO, NO x , Non-Methane Volatile Organic Compounds (NMVOC, speciated in seven species based on the requirements of the CBIV chemical mechanism) and isoprene.
The time-period selected for simulation is August 2005.This was a period of intense convective activity in the Central Atlantic Ocean, accompanied with tropical storms, which in several occasions developed into hurricanes (Lau and Kim, 2007).During August 2005, two severe hurricanes (Irene, Katrina) and three tropical storms (Harvey, Jose, Lee) developed according to the National Hurricane Center, with more than 20 being identified for the period of June to November 2005 (http://www.nhc.noaa.gov).Introduction

Conclusions References
Tables Figures

Back Close
Full Screen / Esc Printer-friendly Version Interactive Discussion (Ramanthan et al., 2001;Lau and Kim 2007;Evan et al., 2009) and the chemical composition of dust entering the Central Atlantic (Husar et al., 1997;Formenti et al., 2003;Zipser et al., 2009).In fact, this is only a subset of the studies performed for this region, emphasizing the scientific interest and continued research interest in the area.The transported amount of aerosols near cloud formation regime is an important characteristic of this period in the selected region, as addressed in this study.

Results and discussion
In the absence of precipitation, fine mode aerosol particles can travel thousands of kilometres from their origin.Furthermore, aerosols in this size range can potentially act as CCN.Andreae et al. (2009) suggest that the typical number size distribution for the aerosols activated as CCN at 0.3% supersaturation (a typical median supersaturation in clouds) is in the range of the Aitken to accumulation size mode.Dusek et al. (2006) show that the critical size range for CCN activation is 0.04-0.12µm at a supersaturation between 0.25 and 1.5%.Ramanathan et al. (2001) state that aerosols such as sulphates, larger than 0.05 µm in diameter, can provide the nuclei for most cloud drops and ice crystals.Therefore, as size and chemical composition of the aerosols are important aspects of the particle efficiency to serve as CCN, the fine mode is selected for the investigation of the aerosol number density.In this work we emphasize the fine mode which includes the first two size sections of the chemistry-transport model, with defined diameter ranges of 0.03-0.1 µm (Aitken mode) and 0.1-2.5 µm (accumulation mode).

Meteorological conditions
During August 2005, the quasi-stationary anticyclone over the Azores exhibited high variability with respect to its position and intensity.As indicated in Fig. 1 keeping in mind that the model does not account for organic aerosols which can also play an important role in the fine mode aerosol mass concentration.Finally, sulphates within the coarse mode are evaluated as they are an essential part of the current study due to their ability to travel long distances and their hygroscopic nature.
The results of the model evaluation are presented in Table 2 and Fig. 2. The performed statistical analysis includes calculation of the mean predicted, mean observed, root mean square error (RMSE), mean bias (MB), mean absolute gross error (MAGE), normalized mean error (NME), normalized mean bias (NMB), fractional error (FE), fractional bias (FB) and the fraction of data within a factor of 2 (FAC2) for a series of modelobservation pairs (Table 2).For the comparison with the measurements, values of the species concentrations at the first model layer were used (0-50 m above ground).The definition of each statistical metric is provided in the Appendix (EPA 2007).The above metrics followed by the scatter diagrams in Fig. 2 provide a general overview of the model performance for the particulate matter concentrations.The PM 10 model concentration is approximated by summing the modelled concentrations of desert dust, sulphates and nitrates (anthropogenic and formed on dust), ammonium, sodium and chloride aerosols of all size sections.PM 2.5 model concentration is approximated by summing the mentioned species of the first two size sections (Table 1).PM 10 and PM 2.5 are underestimated, as indicated by the negative mean, fractional and normalized mean biases (MB, FB and NMB) which are metrics sensitive to systematic errors (MB=−7.3µg/m 3 , NMB=−30.6% for PM 10 and MB=−3.8 µg/m 3 , NMB=−35.67% for PM 2.5 ).The underestimation of the measured concentrations is expected, since the model simulation did not include organic aerosols.This is confirmed by the scatter diagrams in Fig. 2a-b .This statistical metric is considered more robust than the traditional correlation coefficient, since it is not sensitive to outlier data pairs (high or low) as the correlation coefficient.Nevertheless, both metrics are calculated and stated in Fig. 2 and Table 2.
The statistical analysis for the coarse sulphate aerosols shows that the model also performs reasonably well for the sulphates.The mean observed and modelled values are in good agreement (Table 2), the RMSE (1.54µg/m 3 ), MAGE (1.18 µg/m 3 ) and MB (0.64 µg/m 3 ) are quite low, as well as all the other statistical metrics.81.75% of the modelled values are within 2 and 0.5 times the observations (FAC2) and the linear regression of the modelled versus observed values (Fig. 2c) shows a good performance with a correlation coefficient of 0.71.The metrics discussed above provide reasonable confidence in the model performance.
When addressing issues of transport and transformation of pollutants, it is common to use the mass concentration of species or the integrated mass load or both.When the aim is to address the complex interactions between aerosols and clouds, it is preferable to use the number concentration.The conventional use of mass does not account for the chemical and physical characteristics of the aerosols which are important in assessing their effects on clouds and/or precipitation (Andreae et al., 2009).The number distribution for each size section is calculated with the use of the modelled concentration fields, the mass of a single particle (assumed to be spherical) and the geometric Introduction

Conclusions References
Tables Figures

Back Close
Full , where i is the lower size cut point of each section bin).Consequently, in addition to the aerosol concentrations, a comparison has been performed between the calculated aerosol number concentrations and measurements publicly available for August 2005 (Fig. 3).The number of stations providing such data for 2005 is limited and we acquired hourly data from four stations in France, Switzerland and Germany from the database maintained by the Norwegian Institute for Air Research (NILU-EBAS database) and through personal communication with the data providers.Since the model produced 2 h average output fields, the comparison between modelled and observed values was realized on a 2 h average basis.The station in France is Puy de Dome (FR0030R, 45.77 N-2.97 E) and the one in Switzerland is Jungfraujoch (CH0001G, 46.55 N-7.985 E).They both used a condensation particle counter instrument (CPC TSI 3010) for the aerosol number distribution, covering aerosol sizes larger than 0.01 µm for FR0030 and in the range 0.005-10 µm for CH0001 (diameter).The comparison with the total number concentration calculated from the modeled mass gave a reasonable agreement with a correlation coefficient approximately 0.6 (within a total number of 677 pair values) as shown in Fig. 3a.Average and standard deviation values were 1175±735 cm −3 and 1402±761 cm −3 for measurements and model values respectively.For the data pairs of the two stations the statistical analysis has shown that 72% of the modeled data is within 2 and 0.5 times the observations (FAC2).
The small overestimation of the measured average number concentration originates mostly from the comparison with the Jungfraujoch station (CH0001G).For some days in August the model overestimates the number concentration measured in this station, probably due to overestimation of the fine particulate sulfate (the largest portion of the number concentration comes from the smallest size section).Jungfraujoch is an alpine research station located at 3580 m above sea level, surrounded by highly industrialized regions at lower altitudes.Cozic et al. (2008) analyzed time series of measurements in Jungfraujoch station and have shown that the PM 1 mode mass concentration is mainly Introduction

Conclusions References
Tables Figures

Back Close
Full Screen / Esc

Printer-friendly Version
Interactive Discussion composed of organic matter and sulphates with a lower contribution from nitrates, ammonium and black carbon for August 2005.The available number concentration measurements span a diameter range of 0.005-10 µm, which does not enable evaluating the model results based on the model's diameter size bins.As organic aerosols are not part of the current simulation, it was expected that the comparison with the measurements would exhibit deviations from a 1:1 relationship, especially near the species source locations.Overall, the number concentration calculated from the model for the first two stations which are located in the heart of anthropogenic emission sources, can be considered reasonable without large discrepancies.
The other two stations used for the comparison of the aerosol number concentration are located in Germany.The aerosol measurements from the Melpitz station (DE44, 51.53 N-12.90 E) are obtained with the use of a twin differential mobility particle sizer (TDMPS) (Engler et al., 2007;Putaud et al., 2004).A Single Scanning Mobility Particle Sizer (SMPS) was used for the measurements in the Zugspitze station (47.42 N-10.98 E) (Birmili et al., 2009a, b andpersonal communication, 2009).These two techniques are similar but the size cut-offs of the aerosols differ which makes it appropriate to provide different comparison plots and metrics for these stations (Fig. 3b and c).For the alpine station in Zugspitze the comparison is also considered reasonable and the modeled number concentrations agree quite well with the observations (Fig. 3c).The site is located at 2600 m a.s.l. on mount Zuspitze and the aerosol number concentration is given for the diameter size range of 0.0115-0.585µm in 40 size bins.The first size section of the model is used (0.03-0.1 µm, Aitken mode) against the summation of the first 23 size sections of the measurements (0.0115-0.105 µm).The average and standard deviation values for measured and modeled values are 2020±1167 cm −3 and 1899±1366 cm −3 respectively (total number of 326 pair values).The 71.78% of the modeled values are within 0.5 and 2 times the observations (90.8% of the modeled values are lower than 2 times the observed values).For all the stations used in the evaluation, the results from the comparison between calculated and measured aerosol number concentrations show an acceptable agreement.This result should be evaluated keeping in mind that the calculated aerosol number concentration from the modeled concentration is an assumption and it can not produce values as detailed and accurate as a lognormal size distribution.Also, an important issue for the underestimation of the number concentration is the exclusion of the organic aerosols (primary and secondary) which are part of the total aerosol number concentration in the observations.Overall, the calculated aerosol number concentration compared to the observed values give confidence in the model's ability to reproduce the number concentration.This is important for the main goal of this work which concerns the characterization of coefficients of the aerosol models used in the Collection 5 retrieval algorithm (Remer et al., 2006).Whereas the aerosol optical depths are considered to be "validated", the aerosol mass loads are derived products since they do not represent solutions to the retrieval inversion procedure.Hence, their usefulness for a direct pixel-to-gridbox comparison is rather limited, but for spatial structures and relative magnitudes they do provide a useful means for assessing the model performance.
In Figs. 4 and 5 presented here (left column panels), the aerosol mass loads represent a composite of retrievals based on both the land and ocean algorithms.Gaps in the data are primarily caused by clouds or the inability of the algorithm to perform retrievals over bright surfaces.To maximize the number of observations for a given day, data from both Terra during the morning and Aqua in the afternoon are used.Where valid retrievals are available for the same pixel for the two different observational times, the value provided represents an average of the two retrievals.The satellite data are available for the whole period of August 2005; however specific days are selected for discussion in August 2005: 14th-16th in Fig. 4; 27th and 28th in Fig. 5.It should be noted that the aerosol optical depth is not a standard model output field; hence it could not be used in the evaluation procedure.
The modelled daily PM 10 values are determined by the sum of sulphates, nitrates, dust, sodium chloride, and the chemically aged and mixed aerosols (sulphates and nitrates on dust).The modelled mass load represents the vertical integration of the simulated atmospheric column from 0-10 km whereas the satellites cover the entire atmospheric column.As can be seen from Figs. 4 and 5 the model calculated aerosol plumes over the Central Atlantic region follow similar spatial patterns compared to the satellite-based.In particular, the model agrees very well with the observations for the aerosol plume structure as it expands out into the Atlantic.The weaker match for the 14th along the coast and the spatial dislocation on the 28th may be related to temporal mismatches since the satellite observations occur in the middle of the day and the model output represents daily averages.Also there is a strong maximum in the satellite data on 15 and16 August, close to the Canary Islands (Fig. 4), which is underestimated Introduction

Conclusions References
Tables Figures

Back Close
Full Screen / Esc

Printer-friendly Version
Interactive Discussion by the model.The maximum in the model is inland and rather close to the coast while the satellite shows the maximum above the sea.There may be two reasons for this difference; one is spatial dislocation of the simulated aerosol mass load and the second could be the underestimation of aerosols like dust and sulphates (together with the absent organics and DMS fluxes from the sea).While the peak values are of the same magnitude, the model generally produces higher mass concentrations than those inferred from the satellite observations.This result is not unexpected since the satellite based estimates rely on a mass concentration coefficient that is set to a fixed value for the aerosol models used in the retrievals and cannot vary with humidity.Overall, this comparison together with the previous statistical analysis gives reasonable confidence of the ability of the modelling system to simulate the production, transformation and transport of aerosols.Two days were selected for the extended discussion about the atmospheric conditions during the second half of August, the 16 and 27 August 2005.On the 16 and the 27 August 2005, cloud formation was evident, beginning with low clouds in the area of interest, and expanding in the vertical direction as they were travelling west towards the Central Atlantic Ocean.Concentrating on the CWC and the aerosol load in two vertical layers of the model (300-450 m and 450-600 m), two interesting areas were found, one for each day (Fig. 6).The CWC reached 0.8 to 0.9 g/m 3 in these areas (denoted with a black circle in Fig. 6) and all the studied aerosols enter the area of interest (from the daily average load of the species under investigation in Fig. 7).The region illustrated in Fig. 6 and 7 is a zoomed area of the simulation domain, indicated with the rectangular box in Fig. 8.The maximum daily average loads for each aerosol type in the zoomed area of Fig. 7 are quite similar for the two days (16 and 27 August).The desert dust maximum is 255 µg/cm 2 on the 16th and 301 µg/cm 2 on the 27 August.The sea salt maximum is 6.15 µg/cm 2 and 7.48 µg/cm 2 for the two days, respectively.The sulphate on dust maximum load is 0.09 µg/cm 2 on both days and anthropogenic sulphate is 2.81 µg/cm 2 and 2.21 µg/cm 2 , respectively.Large differences occur among the different aerosol types, with desert dust having the largest daily average load as mentioned above.This result, in conjunction with the total daily values in Figs. 4 and 5, shows that the maximum values are related to desert dust and the rest of the species contribute with lower loads.

Case selection for the typical aerosol number concentrations
During all simulated days, daily values of desert dust mass load reach 200-300 µg/cm 2 along the West African coastline.Over the Atlantic Ocean the maximum daily sodium chloride mass load is approximately 9 µg/cm 2 .For the daily anthropogenic aerosols in this region (sulphates and nitrates) the mass load reaches 8-9 µg/cm 2 and for the chemically aged and mixed dust aerosols 3-4 µg/cm 2 are computed.The contributions by the various aerosol types appear rather different when we consider number distribution, as analyzed below.Introduction

Conclusions References
Tables Figures

Back Close
Full Screen / Esc Printer-friendly Version Interactive Discussion

Aerosol vertical profiles -aerosol number distribution
The vertical profiles are shown for the lowest 5 km of the atmosphere, as significant aerosol transport towards the cloud layers takes place.During the simulation days, no deep convective activity was evident in the areas of high aerosol concentrations.Our study mainly associates the aerosols with stratus, stratocumulus and cumulus cloud types, according to the CWC vertical profiles.As expected, the aerosols in the smaller size section -Aitken mode (0.03-0.1 µm) -exhibit higher number concentrations than the ones in the range of 0.1-2.5 µm (accumulation mode).To examine the time evolution of the aerosol/cloud coexistence, vertical profiles were derived every 6 h for all the selected locations in the domain of interest.These locations are shown in Fig. 8 with the green points indicating profiles for 16 August (P1-P3) and the red showing locations for 27 August (P4-P6).The locations were selected using the criterion of simultaneous presence of cloud, all of the studied aerosol types and the absence of rain.The main interest was that all the different aerosol types enter the region of cloud formation, and we selected areas where aerosol/cloud interaction processes commence.For each 0.6 particles/cm3 respectively.For both modes, it appears that the aerosols are mostly affecting the cloud layers with high concentrations from above.A most interesting vertical distribution is of the anthropogenic sulphates (PSO4 1), being stratified above the cloud top, having a maximum of 1260 particles/cm 3 at 1.3 km altitude.
The second location P2 (lat=21.0 • N, lon=18.2 • W) is closer to the coast and the time evolution of the vertical profile indicates larger amounts of Aitken mode aerosols towards the clouded area (Fig. 10, left panels).CRST 1 reaches 1400 particles/cm 3 comparable to PSO4 1, which exhibits a maximum at about 1460 particles/cm 3 but not simultaneously with the non-sulphate-coated dust.The NA 1 maximum is approximately 12 particles/cm 3 near the cloud base and the DSO4 1 particles have only limited vertical extent (max=24 particles/cm 3 ).Both the desert dust and the anthropogenic sulphate aerosols in the Aitken mode can influence the cloud tops with high concentrations, although the sulphate quickly decreases with time.In contrast, the accumulation mode sulphate decreases at slower pace throughout the day.The sulphate aerosol in both modes exhibits different vertical structure which can be explained from the different processes that govern each mode.The accumulation mode has a contribution from the sea-salt sulphate (not present in the smaller size section) that can explain the higher concentration up to 1km from the surface.In the Aitken mode, the sulphate aerosol exhibits a vertical layering throughout the previous day (15 August) which diappears at the beginning of the 16 August.Since the species did not undergo wet deposition (absence of rain at that specific period), the main reasons for the different vertical structure can be the advection/diffusion process and consequently the long-range transport pattern and/or the nucleation process that assigns the nucleated mass to the first section of the distribution (the Aitken mode for this simulation).For the accumulation mode, the number density (Fig. 10, right panels) of CRST 2 reaches Introduction

Conclusions References
Tables Figures

Back Close
Full Screen / Esc Printer-friendly Version

Interactive Discussion
Location three P3 (lat=23.5 • N, lon=17.5 • W) is also close to the coast and the time evolution of the vertical profile is quite similar to that at the previously discussed location for the accumulation mode (Fig. 11).Vertical profiles can be discerned for sulphates and desert dust, which are exhibiting a maximum mostly above the cloud top.The maximum values reached are 466, 9, 1400 and 24 particles/cm 3 for CRST 1, NA 1, PSO4 1 and DSO4 1, respectively (Fig. 11, left panels).Likewise, for the accumulation mode, the maximum values are 77, 2.7, 51 and 0.46 particles/cm 3 for CRST 2, NA 2, PSO4 2 and DSO4 2, respectively.
During the second selected day, 27 August 2005, cloud formation events were much less frequent near the Cape Verde region and along the West African coast.In general, the dense cloud formation during this day occurred mostly over the open ocean towards the west (western part of the domain in Fig. 6).The area investigated with the locations P4-P6 represents the entrance of the aerosols (mainly dust and sulphate on dust) into the oceanic region and in all locations the clouds are shallow.In location P4 (lat=17.0 • N, lon=17.0 • W), CRST 1 reaches 1430 particles/cm 3 (Fig. 12, left panels) and the CRST 2 maximum is 203 particles/cm 3 (Fig. 12, right panels).All other aerosol types have a lower number density compared to the conditions discussed above, with PSO4 1 being almost absent in the Aitken mode.Another important feature that differs from the previously discussed locations is the transport of dust.It typically occurs throughout the lower kilometres with a maximum towards the surface and a strong influence throughout the cloud layer, as the P4 location is very close to the continental outflow region.
The final two locations P5 (lat=15.3 • N, lon=22.5 • W) and P6 (lat=17.0 • N, lon=23.0 • W) are very close to the Cape Verde Islands.For both sites, we compute transport of desert dust towards the cloud layers near the cloud top, with similar concentrations for the accumulation and Aitken modes as those discussed previously (Figs. 13 and 14).DSO4 has the same vertical structure as mentioned earlier for dust, and is present at rather low number densities.The Saharan dust layer extends from approximately 500 m to 3 km, in agreement with that reported by Formenti et al. (2003).Introduction

Conclusions References
Tables Figures

Back Close
Full The flight measurements by Formenti et al. (2003) took place in the same area during September 2000 and showed that the Saharan Air Layer (SAL) base was between 0.5 and 1 km.The top of the layer reached 4-4.5 km in some locations.These findings and the ones in the present work are also in accord with Karyampudi et al. (1999) and Kallos et al. (2006).Anthropogenic sulphate (PSO 4 ) and sodium (NA) can influence the cloud formation from below the cloud base and their number density is lower than before for the sulphates and similar for the sodium aerosol.In terms of number distribution, the predominance of the Aitken mode over the accumulation mode is also a common finding with Formenti et al. (2003) and with Garrett and Hobbs (1995).
To summarize these findings, typical ranges of number concentrations for the different aerosol species for the six locations in the case studies are tabulated in Table 3.These ranges refer to the values occurring at the different vertical layers up to 5 km, except for the last row.When the minimum number concentration was below one, the range of values was substituted by the less-equal sign.The common characteristic of the six locations is the appearance of clouds in the area.The last row in this table (named after "other" in the location column) indicates the maximum values for each aerosol type found in the area under investigation (small box in Fig. 8), regardless of the CWC amount.The values are indicative of the aerosols that can be found in this region throughout the period addressed by this modelling study.Typically, the Aitken mode (0.03-0.1 µm) particle numbers dominate over the accumulation mode (0.1-2.5 µm).Aitken number concentration ranges from 6 to 1400 particles/cm 3 , depending on the aerosol type.Accumulation number concentration ranges from lower than 1 to approximately 200 particles/cm 3 .These values refer to the selected locations P1-P6.When looking at the generic local maxima, the Aitken mode again dominates with a range of 42-5000 particles/cm 3 and the accumulation mode is in the range of 1-600 particles/cm 3 approximately.The most abundant species are desert dust and anthropogenic sulphate with comparable number densities in several cases.The third most important is sulphate on dust and the fourth sodium containing sea salt aerosol.Introduction

Conclusions References
Tables Figures

Back Close
Full

Conclusions
In this study we focused on aerosol number distributions of desert dust, sodium (representing sea salt), sulphate from anthropogenic sources and sulphate produced on the surface of dust particles in the Eastern-Central Atlantic Ocean area.The modelling results suggest that rather small amounts of accumulation mode desert dust, sodium and chemically aged mixed aerosols entered the region where tropical storms and/or hurricane genesis usually takes place (Central Atlantic Ocean).Aerosols of smaller size (Aitken mode) are more abundant in the area and in some occasions sulphates of anthropogenic origin and desert dust are of the same order in terms of number concentration.Of all the considered aerosols, desert dust and anthropogenic sulphate exhibit the highest mass and number concentrations although the other aerosol types are non-negligible as well.Additional findings include: -Desert dust is present in the cloudy area with a maximum number density around 1400 particles/cm 3 in the Aitken mode (0.03-0.1 µm) and 200 particles/cm 3 in the accumulation mode.These maximum values can be higher in the extended area, reaching 3000 particles/cm 3 and 600 particles/cm 3 in the Aitken and accumulation mode, respectively.The results are in good agreement with measurements by Formenti et al. (2003) for the Cape Verde region and Karyampudi et al., (1999) for the Saharan Air Layer.In addition, dust reaches the cloud layers mostly towards the cloud top, and sometimes throughout the lower atmosphere, especially close to the African coast.It is thus conceivable that in the case of storm development more remote from the coast, the dust enters the convective clouds in the active growth phase by entrainment above the boundary layer.Since this is a region where air parcels are likely to enter strong updrafts, the dust particles may serve as ice nuclei in the mixed cloud phase and can potentially affect the cloud water phase transitions and the latent energy budget.
-Sulphates of anthropogenic origin arrive at the Eastern and Central Atlantic Ocean region in significant amounts, reaching approximately 1400 particles/cm Introduction

Conclusions References
Tables Figures

Back Close
Full in the Aitken mode and 60 particles/cm 3 in the accumulation mode.As with dust, these maximum values can be higher in the extended area, reaching 5500 particles/cm 3 and 70 particles/cm 3 in the Aitken and accumulation mode, respectively.Transport towards the stratiform cloud formation area can occur either near the top or the base of the clouds.This depends on the large scale transport component that leads the aerosols to travel either through or above the boundary layer.For the period of August, the Azores anticyclone governs the circulation in the area establishing a strong N-NE flow from Western Europe along the Western Africa coastline.These air masses carry anthropogenic sulphate aerosols above the boundary layer (as can be seen in the locations P1-P3 during the 16 August).
On the contrary, when the anticyclonic system weakens and shifts slightly to the west (as occurs on the 27 August), the sulphate aerosols show a different vertical pattern with the main characteristic the transport through the lower atmospheric layers up to 1 km.In both cases the transported particles strongly enhance the number concentration of CCN as compared to non-polluted and non-dust influenced conditions.
-Sulphates produced on dust particles exhibit relatively low number densities in the accumulation mode as they enter the Central Atlantic Ocean region, with no significant variability in their vertical distribution.In particular, the maximum percentage of sulphate produced on dust is 0.035% of the dust number distribution in the accumulation mode for the six studied locations.This is a rather small percentage of dust being coated with sulphates.On the other hand, this aerosol type is represented by about 20-40 particles/cm 3 in the Aitken mode, with maximum values reaching ∼50 particles/cm 3 .Sulphate on dust in the Aitken mode represents a percentage of the uncoated dust number distribution ranging from 1% to 55% for the six studied locations and for number concentrations higher than 1 particle/cm −3 (average maximum varies between 18% to 38%).This result is significant because the sulphate cover of the dust turns the particles into hygroscopic so that they can also serve as efficient CCN whereas their ability to serve 5207 Introduction

Conclusions References
Tables Figures

Back Close
Full as IN may diminish.Transport occurs mostly through the upper atmospheric layers and rarely through the lower layers, following the dust movement in the area.
and secondary organic aerosol formation/partitioning.The meteorological fields that drive the air quality simulations are taken from the SKIRON/Dust model in an offline mode, with a time resolution of one hour.There is no feedback from the air quality simulation (CAMx) to the meteorological model (SKIRON/Dust).The version of CAMx implemented in the current study is an updated version of the original model.New features include processes that adopt the desert dust fluxes from the SKIRON/Dust Introduction size section (d g = d i * d i +1 The fine fraction of the aerosol number distribution was available (0.028-0.86 µm) in 40 size bins in the Melpitz station.To be in accordance with the model, we added the size bins in the range 0.032-0.1 µm (8 bins) and compared with the first size section of the modeled number concentration (0.03-0.1 µm, Aitken mode).As shown in Fig.3b, the calculated number distribution correlates reasonably well with the measured values (r=0.52).The average and standard deviation values for measured and modeled values are 1979±1138 cm −3 and 1898±1165 cm −3 respectively (total number of 348 pair values).66.66% of the modeled values are within 0.5 and 2 times the observations (88% of the modeled values are lower than 2 times the observed values).The Melpitz sampling site in eastern Germany is surrounded by populated areas with anthropogenic aerosol sources (Engler et al., 2007) which can be indicative of the good correlation with the modelled values.
the aerosols over the Central Atlantic Ocean in terms of number and size distribution and chemical composition.Continuing with the ability of the model to simulate the columnar aerosol mass load, a comparison is made between predicted values and those inferred from observations by the Moderate Resolution Imaging Spectroradiometer (MODIS) aboard the National Aeronautics and Space Administration's (NASA) Aqua and Terra satellite platforms.Tropospheric columnar aerosol mass loads are obtained from the MODIS Collection 5 Level-3 global 1×1 degree dataset.They are based on the weighted combination of retrieved fine and coarse mode aerosol optical depths and the mass concentration Introduction For a few selected days of August 2005 we perform a detailed analysis based on the triple criterion that in the storm genesis region (the area of the Central Atlantic Ocean where cloud formation and convective activity are initiated), clouds, aerosols and small amounts of rain coexist.These days are identified by examining the integrated aerosol column load and the cloud water content in the model vertical layers.Cloud water content is used instead of the cloud cover, since it is explicitly calculated by the atmospheric model SKIRON/Dust.Cloud cover in the CAMx model is diagnosed during the simulation.Sub-grid scale clouds are not included in the model.Clouds are either occupying the entire cell or are absent, depending on the condensed cloud water content (CWC) values.When the CWC is above the threshold of 0.05 g/m 3 the grid cell is designated as cloudy.As indicated in Seinfeld and Pandis (2006) the liquid water content of a cloud varies approximately from 0.05 to 3 g/m 3 , with the most frequent observed values in the range 0.1-0.3g/m 3 .Also, the CWC has a characteristic range of values for different cloud types.Hence, the CWC is selected as the marker of cloud formation in the region under investigation.
location and time two vertical distributions are presented; one for the Aitken and one for the accumulation mode (Fig.9-14).The numerical scale in the Figs.9-14 for each mode is different because the values can differ up to 1 or 2 orders of magnitude for the same species.At the first location P1 (lat 23.5• N, lon 18.5• W) the vertical profile for the Aitken mode (Fig.9, left panels) shows a clear stratification above the cloud top mostly for sulphates, whereas the other species exhibit a much smoother vertical distribution.Dust (CRST 1), sodium containing sea salt (NA 1), anthropogenic sulphates (PSO4 1) and sulphate on dust (DSO4 1) have maximum number concentrations of 188, 9.4, 1260 and 30 particles/cm 3 , respectively.The accumulation mode does not indicate significantly enhanced accumulation mode aerosols in the vicinity of the cloud formation (Fig.9, right panels).Desert dust (CRST 2), sodium dominated sea salt (NA 2), sulphate (PSO4 2) and sulphate on dust (DSO4 2), reach 26, 2.5, 45 and Notes: (1) The notation 1 refers to diameters 0.03-0.1 µm and the notation 2 refers to diameters in the range 0.1-2.5 µm.(2) The "other" in the location column refers to generic maximum values of aerosol number density regardless of the location and time.(3) CRST=crustal material (desert dust), DSO4=sulphate produced on dust surface, PSO4=anthropogenic sulphate with the contribution of ss-sulphate, NA=sodium aerosol.(4) Maximum values for each species are denoted with bold font.12h accumulated precipitation and geopotential heights for August 9 and 16, 2005 t 12:00UTC.

Fig. 4 .Figure 5 :
Fig. 4. Daily average mass load of aerosols from MODIS satellite (left) and from the model simulation (right) (µg/cm 2 ) for 14-16 August 2005.The grey areas in the satellite images are areas with no satellite data.

Figure 6 :
Figure 6: Hourly condensed cloud water content (CWC) during 16 (top) and 27 (bottom) August 2005 (g/m 3 ) at 12UTC.The left to right panels indicate two vertical layers, 300-450m and 450-600m respectively.The domain shown is the one indicated with a rectangular box in Figure 7.

Fig. 6 .Figure 7 :
Fig. 6.Hourly condensed cloud water content (CWC) during 16 (top) and 27 (bottom) August 2005 (g/m 3 ) at 12:00 UTC.The left to right panels indicate two vertical layers, 300-450 m and 450-600 m respectively.The domain shown is the one indicated with a rectangular box in Fig. 8.

Fig. 7 .
Fig. 7. Integrated column mass load (µg/cm 2 ) for desert dust (crustal), sea salt aerosol, sulphate on dust and anthropogenic sulphate particles.The left column is for 16 August 2005 and the right for 27 August 2005.

Figure 8 :
Figure 8: Location of the points used for the vertical profile of the aerosol number distribution (the green points are for August 16, 2005 and the red points for August 27, 2005) (image from Google Earth).The rectangular box denotes the area under investigation for the chemical composition and vertical number distribution of the aerosols.

Fig. 8 .Figure 9 :Fig. 9 .Figure 10 :Fig. 10 .Figure 11 :Fig. 11 .Figure 12 :Figure 13 :Fig. 13 .Figure 14 :Fig. 14 .
Fig. 8. Location of the points used for the vertical profile of the aerosol number distribution (the green points are for 16 August 2005 and the red points for 27 August 2005) (image from Google Earth).The rectangular box denotes the area under investigation for the chemical composition and vertical number distribution of the aerosols.

Table 1
and are briefly discussed below.For the performed simulations with the SK-IRON/Dust model, we used initial and lateral boundary conditions (0.5• ×0.5• ) from the European Centre for Medium-Range Weather Forecasts (ECMWF) at 16 isobaric levels.Also, daily SSTs are adopted from the ECMWF.The topography and vegetation input fields were obtained from the US Geological Survey (USGS) 30 ×30 horizontal analysis dataset.The soil texture was taken from ZOBLER and FAO/UNESCO at 30 ×30 resolution.The model horizontal resolution applied is 0.24 • ×0.24• , with 38 vertical levels up to 22 km altitude.

2 Evaluation of the modelled aerosol concentrations
, during the ninth day of the simulation (9 August 2005) the anticyclone was not present in the area, Introduction Atlantic Ocean, these stations are important for identifying the dust transport towards the west.Unfortunately, from the 37 available stations in this region, only 11 provided PM 10 data for August 2005 and from those only 4 could be used because the coarse resolution of the model included more than one station in a single grid cell.The selection of the four stations was based on their type, preferring the use of Introduction , where the correlation coefficients are 0.47 Introduction general, aerosols with different chemical composition and size reach the Eastern and Central Atlantic cloudy region in different vertical layers.The size and composition is related to the origin of the air masses (anthropogenic or natural sources) and the transformations occurring during transport.The aerosol mixture, hygroscopicity (CCN and IN properties) and the altitude at which they arrive may crucially influence the microphysics and the life cycle of different cloud types.Introduction Frank, G. P., Hildebrandt, L., Curtius, J., Schneider, J., Walter, S., Chand, D.,  Drewnick, F., Hings, S., Jung, D., Borrmann, S., and Andreae, M. O.: Size Matters More Than Chemistry for Cloud-Nucleating Ability of Aerosol Particles, V., and Anderson, B.: The Saharan Air Layer and the Fate of African Easterly Waves -NASA's AMMA Field Study of Tropical Cyclogenesis, B. Am.Meteor.Soc., 90(8), Introduction

Table 2 .
Comparison of modelled versus measured daily concentrations of PM 10 , PM 2.5 , PM 10 Sulfates for the period 2-30 August 2005.Data adopted from the EMEP database (www.emep.

Table 3 .
Typical ranges of number concentrations for desert dust, sulphates and sodium aerosols (cm −3 ) in the vertical column of 0-5 km.Model results for the 16 and 27 August 2005 (2 h average values).