Overview of the Chemistry-Aerosol Mediterranean Experiment/Aerosol Direct Radiative Forcing on the Mediterranean Climate (ChArMEx/ADRIMED) summer 2013 campaign

Atmos. Chem. Phys., 16, 455–504, 2016 www.atmos-chem-phys.net/16/455/2016/ doi:10.5194/acp-16-455-2016 © Author(s) 2016. CC Attribution 3.0 License. Overview of the Chemistry-Aerosol Mediterranean Experiment/Aerosol Direct Radiative Forcing on the Mediterranean Climate (ChArMEx/ADRIMED) summer 2013 campaign M. Mallet1, F. Dulac2, P. Formenti3, P. Nabat4, J. Sciare2,5, G. Roberts4, J. Pelon6, G. Ancellet6, D. Tanré7, F. Parol7, C. Denjean3,4, G. Brogniez7, A. di Sarra8, L. Alados-Arboledas9, J. Arndt10, F. Auriol7, L. Blarel7, T. Bourrianne4, P. Chazette2, S. Chevaillier3, M. Claeys4, B. D’Anna11, Y. Derimian7, K. Desboeufs3, T. Di Iorio8, J.-F. Doussin3, P. Durand1, A. Féron3, E. Freney12, C. Gaimoz3, P. Goloub7, J. L. Gómez-Amo8, M. J. Granados-Muñoz9, N. Grand3, E. Hamonou2, I. Jankowiak7, M. Jeannot13, J.-F. Léon1, M. Maillé3, S. Mailler14, D. Meloni8, L. Menut14, G. Momboisse4, J. Nicolas4,12, T. Podvin7, V. Pont1, G. Rea14, J.-B. Renard13, L. Roblou1, K. Schepanski15, A. Schwarzenboeck12, K. Sellegri12, M. Sicard16, F. Solmon17, S. Somot4, B Torres7, J. Totems2, S. Triquet3, N. Verdier18, C. Verwaerde7, F. Waquet7, J. Wenger10, and P. Zapf3 1Laboratoire d’Aérologie, Observatoire Midi-Pyrénées, 14 Avenue Edouard Belin, Toulouse, France 2Laboratoire des Sciences du Climat et de l’Environnement, LSCE/IPSL, CEA-CNRSUVSQ, Université Paris-Saclay, Gif-sur-Yvette, France 3Laboratoire Inter-Universitaire des Systèmes Atmosphériques (LISA), UMR CNRS 7583, Université Paris Est Créteil et Université Paris Diderot, Institut Pierre Simon Laplace, Créteil, France 4Météo-France, CNRM-GAME, Centre National de Recherches Météorologiques, UMR3589, Toulouse, France 5The Cyprus Institute, Energy, Environment and Water Research Center, Nicosia, Cyprus 6LATMOS-ISPL, UPMC Univ. Paris 06 et Université Versailles St-Quentin; CNRS/INSU, Paris, France 7LOA, Université Lille 1, Villeneuve d’Ascq, France 8ENEA, Laboratory for Earth Observations and Analyses, Via Anguillarese 301, Rome, Italy 9Department of Applied Physics, University of Granada, Granada, Spain 10Department of Chemistry and Environmental Research Institute, University College Cork, Cork, Ireland 11Institut de Recherches sur la Catalyse et l’Environnement de Lyon (IRCELYON), University of Lyon, Villeurbanne, France 12Laboratoire de Météorologie Physique, CNRS UMR6016, Observatoire de Physique du Globe de Clermont-Ferrand, Université Blaise Pascal, Aubière, France 13LPC2E-CNRS/Université d’Orléans, 3A Avenue de la recherche scientifique, Orléans, France 14LMD, IPSL, CNRS, Ecole Polytechnique, École Normale Supérieure, Université Paris 6, UMR8539, Palaiseau CEDEX, France 15Leibniz Institute for Tropospheric Research, Permoserstr. 15, 04318 Leipzig, Germany 16RSLab/CTE-CRAE-IEEC, Universitat Politècnica de Catalunya, Barcelona, Spain 17The Abdus Salam International Center for Theoretical Physics, Strada Costiera 11, Trieste, Italy 18Centre National d’Etudes Spatiales (CNES), DCT/BL/NB, 18 Avenue Edouard Belin, Toulouse, France Correspondence to: M. Mallet ([email protected]) Received: 19 June 2015 – Published in Atmos. Chem. Phys. Discuss.: 17 July 2015 Revised: 11 December 2015 – Accepted: 14 December 2015 – Published: 19 January 2016 Published by Copernicus Publications on behalf of the European Geosciences Union. 456 M. Mallet et al.: Overview of the ChArMEx/ADRIMED summer 2013 campaign Abstract. The Chemistry-Aerosol Mediterranean Experiment (ChArMEx; http://charmex.lsce.ipsl.fr) is a collaborative research program federating international activities to investigate Mediterranean regional chemistry-climate interactions. A special observing period (SOP-1a) including intensive airborne measurements was performed in the framework of the Aerosol Direct Radiative Impact on the regional climate in the MEDiterranean region (ADRIMED) project during the Mediterranean dry season over the western and central Mediterranean basins, with a focus on aerosol-radiation measurements and their modeling. The SOP-1a took place from 11 June to 5 July 2013. Airborne measurements were made by both the ATR-42 and F-20 French research aircraft operated from Sardinia (Italy) and instrumented for in situ and remote-sensing measurements, respectively, and by sounding and drifting balloons, launched in Minorca. The experimental setup also involved several ground-based measurement sites on islands including two ground-based reference stations in Corsica and Lampedusa and secondary monitoring sites in Minorca and Sicily. Additional measurements including lidar profiling were also performed on alert during aircraft operations at EARLINET/ACTRIS stations at Granada and Barcelona in Spain, and in southern Italy. Remote-sensing aerosol products from satellites (MSG/SEVIRI, MODIS) and from the AERONET/PHOTONS network were also used. Dedicated meso-scale and regional modeling experiments were performed in relation to this observational effort. We provide here an overview of the different surface and aircraft observations deployed during the ChArMEx/ADRIMED period and of associated modeling studies together with an analysis of the synoptic conditions that determined the aerosol emission and transport. Meteorological conditions observed during this campaign (moderate temperatures and southern flows) were not favorable to producing high levels of atmospheric pollutants or intense biomass burning events in the region. However, numerous mineral dust plumes were observed during the campaign, with the main sources located in Morocco, Algeria and Tunisia, leading to aerosol optical depth (AOD) values ranging between 0.2 and 0.6 (at 440nm) over the western and central Mediterranean basins. One important point of this experiment concerns the direct observations of aerosol extinction onboard the ATR-42, using the CAPS system, showing local maxima reaching up to 150Mm−1within the dust plume. Non-negligible aerosol extinction (about 50Mm−1) has also been observed within the marine boundary layer (MBL). By combining the ATR42 extinction coefficient observations with absorption and scattering measurements, we performed a complete optical closure revealing excellent agreement with estimated optical properties. This additional information on extinction properties has allowed calculation of the dust single scattering albedo (SSA) with a high level of confidence over the western Mediterranean. Our results show a moderate variability from 0.90 to 1.00 (at 530nm) for all flights studied compared to that reported in the literature on this optical parameter. Our results underline also a relatively low difference in SSA with values derived near dust sources. In parallel, active remote-sensing observations from the surface and onboard the F-20 aircraft suggest a complex vertical structure of particles and distinct aerosol layers with sea spray and pollution located within the MBL, and mineral dust and/or aged North American smoke particles located above (up to 6–7km in altitude). Aircraft and balloon-borne observations allow one to investigate the vertical structure of the aerosol size distribution showing particles characterized by a large size (>10µm in diameter) within dust plumes. In most of cases, a coarse mode characterized by an effective diameter ranging between 5 and 10µm, has been detected above the MBL. In terms of shortwave (SW) direct forcing, in situ surface and aircraft observations have been merged and used as inputs in 1-D radiative transfer codes for calculating the aerosol direct radiative forcing (DRF). Results show significant surface SW instantaneous forcing (up to −90Wm−2at noon). Aircraft observations provide also original estimates of the vertical structure of SW and LW radiative heating revealingsignificant instantaneousvalues of about5◦Kper day in the solar spectrum (for a solar angle of 30◦) within the dust layer. Associated 3-D modeling studies from regional climate (RCM) and chemistry transport (CTM) models indicate a relatively good agreement for simulated AOD compared with observations from the AERONET/PHOTONS network and satellite data, especially for long-range dust transport. Calculations of the 3-D SW (clear-sky) surface DRF indicate an average of about −10 to −20Wm−2(for the whole period) over the Mediterranean Sea together with maxima (−50Wm−2) over northern Africa. The top of the atmosphere (TOA) DRF is shown to be highly variable within the domain, due to moderate absorbing properties of dust and changes in the surface albedo. Indeed, 3-D simulations indicate negative forcing over the Mediterranean Sea and Europe and positive forcing over northern Africa. Finally, a multiyear simulation, performed for the 2003 to 2009 period and including an ocean–atmosphere (O–A) coupling, underlines the impact of the aerosol direct radiative forcing on the sea surface temperature, O–A fluxes and the hydrological cycle over the Mediterranean. 1 Introduction The Mediterranean region has been identified as one of the most prominent “Hot-Spots” in future climate change projections (Giorgi and Lionello, 2008). It is characterized by its vulnerability to changes in the water cycle (e.g., Chenoweth et al., 2011; García-Ruiz et al., 2011). General circulation model (GCM) and regional climate model (RCM) simulations show a substantial precipitation decrease and a warming of the region, especially in the long warm and dry Atmos. Chem. Phys., 16, 455–504, 2016 www.atmos-chem-phys.net/16/455/2016/ M. Mallet et al.: Overview of the ChArMEx/ADRIMED summer 2013 campaign 457 Mediterranean season. At the end of the 21st century, the average of the model outputs predicts a significant loss of freshwater (40% for the period 2070–2090 compared to 1950–1999; Sanchez-Gomez et al., 2009) over the Mediterranean region. More recently, Mariotti et al. (2015) have used the newly available Coupled Model Intercomparison ProjectPhase 5 (CMIP5) experiments and show a significant increase in the projected surface air temperature (by ∼ +2– 3◦C) for the 2071–2098 period compared to 1980–2005. These results need to be put in the context of an increasing anthropogenic pressure on the Mediterranean region, with an expected doubling of the population in countries around the Mediterranean basin in the next decades, with a contrast between a small decrease in European countries and a strong increase in African and Middle Eastern countries (Brauch, 2003). However, as highlighted by Mariotti et al. (2008), despite the high degree of model consistency, the results regarding the future climate projections for the Mediterranean Sea water budget from the global coupled models are still uncertain due to their horizontal spatial resolutions that are not capable of resolving the local to regional Mediterraneanspecific processes (air–sea exchanges, coastline, topography, north–south gradient of surface albedo). Indeed, the Mediterranean climate is affected by local processes induced by the complex physiography of the region and the presence of a large body of water (the Mediterranean Sea). For example, the Alpine chain is a strong factor in modifying traveling synoptic and mesoscale systems and the Mediterranean Sea is an important source of moisture and precipitation in the region (Gimeno et al., 2010; Schicker et al., 2010) and of energy for storms (Lionello et al., 2006). The complex topography, coastline and vegetation cover of the region are well known to modulate the regional climate signal at small spatial scales (e.g., Millán et al., 1997; Gangoiti et al., 2001; Lionello et al., 2006). So far, most global and regional climate simulations have investigated the impact of global warming on the Mediterranean climate without detailed considerations of the possible radiative influence and climatic feedback from the different Mediterranean aerosols (anthropogenic, marine, biomass burning, secondary biogenic and mineral dust particles). The Mediterranean region is rich in a variety of particles (natural and anthropogenic) from both continental and marine sources (Lelieveld et al., 2002). In Fig. 1, we illustrate the significant differences in aerosol loading between the eastern, central, and western sub-basins and between the north and the south of the Mediterranean shown by long-term aerosol satellite products. The aerosol optical depth (AOD), which represents the integration of the extinction by particles along the whole atmospheric column, displays annual mean values (Fig. 1) from 0.2 to 0.5 (in the visible wavelengths), depending on the aerosol types observed over the Euro-Mediterranean region (Nabat et al., 2013). Numerous studies have documented the AOD for polluted anthropogenic Mediterranean aerosols at a local scale over southeastern France (Mallet et al., 2006; Roger et al., 2006), Spain (Horvath et al., 2002; Alados-Arboledas et al., 2003, 2008), the western Mediterranean (Lyamani et al., 2015), Greece (Chazette and Liousse, 2001; Gerasopoulos et al., 2003), the Greek island of Crete (Fotiadi et al., 2006), and Italy (Tafuro et al., 2007; Ciardini et al., 2012). Under polluted conditions, they report low to moderate AOD values ranging between 0.1 and 0.5 (at visible wavelengths). In parallel, multi-year TOMS and MODIS observations over the eastern Mediterranean (Hatzianastassiou et al., 2009) or the Po Valley (Royer et al., 2010) indicate the occurrence of high AOD values (up to more than 0.8 at 500nm) over large urban areas surrounding megacities. Numerous studies (Markowicz et al., 2002; Ravetta et al., 2007; Liu et al., 2009; Kaskaoutis et al., 2011; Barnaba et al., 2011; Amiridis et al., 2012; Baldassarre et al., 2015) have also been dedicated to biomass burning aerosols over the Mediterranean, which are mainly observed in July and August (driest months of the year), when the development of forest fires is favored (Pace et al., 2005). Long-term observations of absorbing aerosols have clearly shown the major role of long-range transport of biomass (agriculture waste) burning in the eastern Mediterranean (Sciare et al., 2008). AOD data available for smoke particles show “intermediate” values between those observed for dust and anthropogenic particles. For example, AOD values ranging between 0.3 and 0.8 (Pace et al., 2005) have been observed at Lampedusa from 5 to 22 August 2003, in relation to intense fires developed in southern Europe and transported over the Mediterranean basin during a regional heat wave. In addition, the STAAARTE-MED experiment (August 1998 in the eastern Mediterranean) has also documented a mean AOD of 0.39 (at 550nm) for aged smoke plume from Canadian fires (Formenti et al., 2002). This kind of long-range transport has also been observed over the western Mediterranean (OrtizAmezcua et al., 2014). Concerning natural aerosols, different cases of Saharan mineral dust have been regularly documented with local optical measurements on the island of Lampedusa by Meloni et al. (2003, 2004), who indicate moderate AOD (at 415.6nm) of about 0.23–0.26 and one significantly larger event with AOD values of 0.51. Meloni et al. (2008) also report AOD (at 500nm) measurements ranging between 0.29 and 1.18 for the 1999 to 2006 period. For some extreme cases, dust AOD peaks may be even larger, reaching values up to 2 as observed by di Sarra et al. (2011). In parallel to Lampedusa observations, Kubilay et al. (2003) have also documented three dust intrusion events at Erdemli (Turkish coast), occurring in spring from the central Sahara, in summer from the eastern Sahara, and in the fall from the Middle East/Arabian Peninsula. In each case, the presence of dust particles significantly increased the AOD, up to 1.8. Over the western Mediterranean, different studies also reveal the impact of Saharan dust that occasionally can lead to extreme events with AOD (at 500nm) above 1 (Guerrero-Rascado et al., 2009). www.atmos-chem-phys.net/16/455/2016/ Atmos. Chem. Phys., 16, 455–504, 2016 458 M. Mallet et al.: Overview of the ChArMEx/ADRIMED summer 2013 campaign Figure 1. Aerosol optical depth (at 550nm) derived from the MODIS and MISR satellites for the 2003 to 2012 period. The AERONET/PHOTONS AOD is also indicated (in small circles). For sea-spray particles, which are the second main natural species observed over the Mediterranean, Nabat et al. (2013) report a relatively low monthly mean AOD derived from satellites and modeling data, with values lower than 0.05 (in the visible wavelengths). By using recent improvements in the sea-spray emission scheme, Spada et al. (2013) show an averaged sea-spray AOD of around 0.04 for the month of January (5-year period, 2002–2006), which is the favorable period for generating primary sea spray due to strong sea surface winds. Finally, and in the case of extreme wind episodes occurring over the western basin, Salameh et al. (2007) show that the amount of aerosol loading, solely due to the Mistral, Tramontane and Ligurian outflows, is as large as 3–4 times the background aerosol amount. They indicate that the contribution of sea-spray particles to the total aerosol loading and optical depth ranges from 1 to 10%. Salameh et al. (2007) report AOD around 0.15–0.20 (at 865nm) within the sea-spray aerosol plume during such strong wind events. In addition, Mulcahy et al. (2008) reported a high correlation between AOD and wind speed, with AOD values of 0.3–0.4 at moderately high wind speed. In addition to AOD, the knowledge of SSA is essential for estimating the aerosol direct and semi-direct radiative forcing. Concerning mineral dust particles observed over the Mediterranean, it should be noted that significant variations in SSA are reported, with values near 1 for purely scattering aerosols, and quite remarkable low values (0.74, 0.77 or 0.81) at Lampedusa (Pace et al., 2006; Meloni et al., 2003). At the high-altitude Alpine Jungfraujoch station, SSA values are generally higher than 0.9 in the case of African dust, but occasional SSA as low as 0.75–0.80 is reported by Collaud Coen et al. (2004). Intermediate values (0.85–0.92) have also been reported over the Mediterranean basin (Kubilay et al., 2003; Meloni et al., 2004; Saha et al., 2008). These estimates clearly indicate that significantly different SSA values are obtained following the dust particle origins and/or possible mixing of mineral dust with other species. For example, Kubilay et al. (2003) underlined the importance of mixing, showing SSA values clearly lower (0.85–0.90) in the case of mineral dust transport coincident with urban–industrial aerosols, as compared to pure dust (0.96–0.97). In addition, SSA observed in the case of urban/industrial regimes has also been well documented over the Mediterranean Sea and coastal regions. In most cases, moderate or low SSA (0.78–0.94) is observed due to emissions containing absorbing black carbon aerosols. Over southeastern France, optical computations performed by Mallet et al. (2004) and Saha et al. (2008) indicate SSA values of 0.85 and 0.83 (at 550nm) near the cities of Toulon and Marseilles, respectively. Aircraft observations performed over the Marseilles/Etang de Berre area during the ESCOMPTE campaign show values ranging between 0.88 and 0.93 (at 550nm) in the planetary boundary layer (PBL) (Mallet et al., 2005). These SSA values are close to those observed in southern Spain (0.86–0.90) by Horvath et al. (2002). Over southeastern Italy, Tafuro et al. (2007) reported a value of 0.94 during summertime corresponding to anthropogenic particles. Finally, polluted particles transported over the Mediterranean basin also have relatively low values as reported by Markowicz et al. (2002) over Crete (0.87) and by Di Iorio et al. (2003) (0.79–0.83) over Lampedusa for two cases (25 and 27 May 1999) of “aged” anthropogenic aerosols originating from Europe. As opposed to dust and polluted aerosols, few studies have derived the biomass burning SSA over the Mediterranean Sea. One estimate has been obtained during STAAARTEMED by Formenti et al. (2002), who reported a mean dry SSA of 0.89 (at 500nm) for aged smoke from North America. Meloni et al. (2006) report estimations at Lampedusa with values of 0.82±0.04 (at 415nm) for smoke aerosols over the Mediterranean region. The observed differences between SSA values could be due to the fact that the smoke events described by Meloni et al. (2006) are more “local” and not (or somewhat less) mixed with other secondary Atmos. Chem. Phys., 16, 455–504, 2016 www.atmos-chem-phys.net/16/455/2016/ M. Mallet et al.: Overview of the ChArMEx/ADRIMED summer 2013 campaign 459 species, as compared to biomass burning particles documented by Formenti et al. (2002), which were issued from very distant Canadian fires. Finally, at Palencia (Spain), Cachorro et al. (2008) reported a column-integrated SSA of 0.88 (at 440nm) for a biomass burning event occurring in 28 July 2004. One should remember that most estimations of SSA over the Mediterranean have been obtained from surface in situ or remote-sensing techniques. In that sense, the ChArMEx/ADRIMED project provides innovative observations of 3-D aerosol SSA, allowing investigating changes in its optical property during the transport of aerosols over the Mediterranean. Concerning the aerosol vertical profiles and apart from a few airborne in situ measurements (Formenti et al., 2002), most of the available information in the Mediterranean region comes from lidar observations, which provide highly resolved vertical profiles of aerosol backscattering at one or more wavelengths and, depending on the complexity of the instrumental setup, particles depolarization and extinction. Several sites are equipped with aerosol lidar systems and carry out regular observations in a coordinated way within EARLINET (European Aerosol Research Lidar Network; Pappalardo et al., 2014; Wang et al., 2014). Numerous studies have been specifically dedicated to the vertical distribution of Saharan dust during extended time periods and/or selected events from various Mediterranean regions, mainly from ground-based systems: (i) the eastern basin in Thessaloniki (Hamonou et al., 1999; Balis et al., 2004), Crete (Balis et al., 2006), the Aegean Sea (Dulac and Chazette, 2003), and Athens plus Thessaloniki (Papayannis et al., 2005; Balis et al., 2006); (ii) the central basin in Lampedusa (Di Iorio et al., 2003; Meloni et al., 2004), Lecce (Tafuro et al., 2006), and at Etna (Tafuro et al., 2006); and (iii) across the western basin with the first spaceborne lidar (Berthier et al., 2006) and at Observatoire de Haute Provence (Hamonou et al., 1999), and Barcelona (Pérez et al., 2006; Sicard et al., 2011). Finally, using data from 20 EARLINET lidar stations, Papayannis et al. (2008) indicate that African dust transport over the Mediterranean basin is layered. Their analysis confirms early observations by Hamonou et al. (1999) that not only different dust layers are superimposed at different altitudes, but that these layers have different source regions. The dust layers were generally detected between 1.8 and 9km altitude. Not only desert dust, however, can be transported above the marine atmospheric boundary layer. Balis et al. (2004) report non-dust aerosols within elevated layers over Thessaloniki, and Formenti et al. (2002) report a forest fire haze layer from Canada observed from airborne measurements between approximately 1 and 3.5km above the northeastern Mediterranean in August 1998. Pérez et al. (2004) describe the complex interaction among orography, sea breeze and pollution that causes the recirculation of pollutants and produces a strong layering with pollution aerosol layers above the boundary layer in the region of Barcelona. In addition, aerosol plumes are emitted sporadically in the Mediterranean free troposphere by Etna volcano. Such plumes have been observed to travel at altitudes between 4 and 5km (Pappalardo et al., 2004) or above (Sellitto et al., 2015) at a relatively short distance from Etna. To summarize, the lidar observations clearly show that only part of the aerosol transport occurs in the MBL demonstrating the need of using aircraft observations within the aerosol plume to determine the aerosol microphysical–chemical and optical properties of particles transported in altitude and so not detectable at the surface. Indeed, although lidar observations provide obviously crucial information on the aerosol vertical profiles, most of lidar systems cannot derive information on the aerosol size distribution, optical properties and chemical composition along the vertical. Such observations can only be obtained using in situ aircraft vertical profiles as proposed in this ChArMEx/ADRIMED experiment. As an example, this project provides interesting and unique observations of 3-D aerosol size distribution during the transport over the Mediterranean basin, allowing one to investigate changes in size distribution between mixed and pure mineral dust. In terms of radiative effects, such atmospheric aerosol characteristics (loadings, absorbing properties, vertical layering) are known (Nabat et al., 2012; Papadimas et al., 2012; Zanis et al., 2012) to significantly change the radiative budget of the Mediterranean region by (1) decreasing the sea surface incoming shortwave radiations, (2) increasing/decreasing outgoing shortwave fluxes depending on the surface albedo and (3) possibly heating turbid atmospheric layers when particles absorb solar light. This is the so-called aerosol direct radiative forcing (DRF). As for the AOD, many aerosol DRF calculations are now referenced over the Mediterranean, clearly showing that the DRF is significantly larger at daily timescales than the one exerted by the additional anthropogenic greenhouse gases. Concerning polluted aerosols, shortwave DRF has been estimated by many authors (Horvath et al., 2002; Markowicz et al., 2002; Meloni et al., 2003; Mallet et al., 2006; Roger et al., 2006; Saha et al., 2008; di Sarra et al., 2008; Di Biagio et al., 2009, 2010). Studies show significant decreases of surface solar fluxes of about 20–30Wm−2(daily mean) for different locations such as Almeria (southern Mediterranean coast of Spain), Finokalia (Crete), Lampedusa, Marseilles and Toulon (southeastern France). In parallel, the combination of surface and satellite remote-sensing observations performed at Lampedusa has been used to perform calculations of the DRF, both in the shortwave (SW; Di Biagio et al., 2010) and longwave (LW; di Sarra et al., 2011; Meloni et al., 2015) spectral regions for different cases of Saharan dust intrusions. These studies emphasize that the radiative effect of desert dust in the LW spectral range is significant, and offsets a large fraction of the SW forcing (di Sarra et al., 2011; Meloni et al., 2015). More recently, Sicard et al. (2014a, b) have also produced estimations of the dust LW radiative effect, www.atmos-chem-phys.net/16/455/2016/ Atmos. Chem. Phys., 16, 455–504, 2016 460 M. Mallet et al.: Overview of the ChArMEx/ADRIMED summer 2013 campaign based on remote-sensing observations in Barcelona and 1-D radiative transfer calculations. Concerning the smoke DRF, some calculations have been conducted over the Mediterranean region by Markowicz et al. (2002), di Sarra et al. (2008), Kaskaoutis et al. (2011) or Formenti et al. (2002). One estimate, proposed by Formenti et al. (2002) for an aged Canadian biomass-burning plume, reveals a significant SW surface dimming of about ∼ 60Wm−2. In addition, the DRF induced by smoke aerosols at Lampedusa between 3 and 23 August 2003, during the exceptionally hot and dry season, was derived by Pace et al. (2005) for the 300–800nm spectral range. The smoke atmospheric forcing was estimated to be between +22 and +26Wm−2, with a corresponding SW heating rate possibly exceeding 2Kd−1at the smoke plume altitude. At the regional scale, Papadimas et al. (2012) have proposed a recent estimation of the aerosol DRF using MODIS data from 2000 to 2007 for both all-sky and clear-sky conditions. They derived a multi-year regional mean surface of −19Wm−2, associated with a TOA DRF of −4.5Wm−2. Regional modeling studies have also been recently proposed by Nabat et al. (2012, 2014, 2015a) using the coupledchemistry RegCM and CNRM-Regional Climate System Model (RCSM) for multi-year simulations. These works reported a mean regional surface (TOA) forcing of about −12Wm−2(−2.4Wm−2) and −16Wm−2(−5.7Wm−2) for the RegCM and CNRM-RCSM models, respectively. RegCM has also been used to investigate direct and semidirect radiative effects of mineral dust over the Sahara and Europe in a test case of July 2003 (Santese et al., 2010). In this work, Santese et al. (2010) computed a daily mean SW DRF of −24Wm−2(−3.4Wm−2) on 17 July and −25Wm−2(−3.5Wm−2) on 24 July at the surface (TOA) on average over the simulation domain. Zanis et al. (2012) also proposed a regional estimate of the DRF of anthropogenic particles over the 1996–2007 period using RegCM and showed a significant forcing of up to −23Wm−2at TOA over eastern Europe. In addition, Péré et al. (2011) have used the CTM-CHIMERE model coupled to the WRF model, for estimating the DRF of anthropogenic particles during the heat wave of summer 2003 and showed significant effects with implications on the planetary boundary layer height (decrease up to 30% in the presence of anthropogenic aerosols) and local air quality. In addition to their important effects on the surface and TOA DRF, most of the Mediterranean aerosols are also able to absorb more or less effectively the solar radiations leading to a significant atmospheric forcing and associated SW heating rate. Local studies previously mentioned (Roger et al., 2006; Saha et al., 2008; Pace et al., 2005; Péré et al., 2011; Meloni et al., 2015) clearly report significant SW heating rate due to absorbing particles with values reaching up to 2–3K per day, depending on the aerosol types. Finally, aerosols also have a significant effect on photolysis rates that may affect tropospheric chemistry and ozone production over the basin (Casasanta et al., 2011; Mailler et al., 2015). With regards to such surfaces, TOA and atmospheric forcings, there is a need to investigate how the change in the radiative budget due to natural/anthropogenic aerosols influences surface temperature (both over land and sea), relative humidity profiles, exchanges (latent heat fluxes) between ocean and atmosphere, cloud cover (semi-direct effect of absorbing particles), precipitation and finally the whole Mediterranean hydrological cycle. The induced perturbations in the sea surface–atmosphere fluxes is expected to be important despite the relatively small size of the Mediterranean Sea, since this basin plays an important role at much larger scale by providing moisture for precipitation to its surroundings land region extending to northern Europe and northern Africa (Gimeno et al., 2010; Schicker et al., 2010). Indeed, and as shown by Ramanathan et al. (2001) for the Indian region or Foltz and McPhaden (2008) and Yue et al. (2011) for the Atlantic Ocean, a modification of the sea surface evaporative fluxes, due to the dimming radiative effect of aerosols at the sea surface could significantly influence the lower troposphere moisture content and the associated precipitation distribution around the Mediterranean. In parallel, the absorbing particles over the Mediterranean (Mallet et al., 2013) could exert a semi-direct effect that could modify the vertical profiles of relative humidity and cloud cover, which has to be quantified. To our knowledge, there is no regional climate simulation over the Mediterranean basin at this time that includes an ocean–atmosphere (O–A) coupled system model for investigating this specific question. In that context of the referenced modeling and observation research over the Mediterranean basin, the main objectives of the ChArMEx/ADRIMED project were the following: –to conduct an experimental campaign, based on surface and aircraft observations, for creating a rich 3-D database of physical, chemical and optical properties of the main Mediterranean aerosols, including (i) original in situ aircraft observations of extinction coefficients, size distribution, black carbon concentrations as well as (SW and LW) radiative fluxes and associated heating rates, (ii) balloons observations of aerosol size distribution and (iii) surface measurements including original characterization of chemical properties –to investigate how the aerosol size distribution and optical (especially SSA) properties evolve along the vertical, between the MBL and elevated layers, and during the transport over the Mediterranean –to use experimental surface and aircraft observations to estimate the 1-D local DRF and forcing efficiency of different aerosols at the surface, TOA and within the atmospheric layer –to investigate how the modifications of the radiative budget due to aerosols affect the sea surface evaporation Atmos. Chem. Phys., 16, 455–504, 2016 www.atmos-chem-phys.net/16/455/2016/ M. Mallet et al.: Overview of the ChArMEx/ADRIMED summer 2013 campaign 461 fluxes, relative humidity profiles, cloud cover, precipitation and more largely the Mediterranean hydrological cycle The present article describes the experimental setup of the campaign and the meteorological context and illustrates important results detailed in a series of companion papers. The rest of this article is divided into six different parts. In the first and second parts (Sects. 2 and 3), we describe the in situ and remote-sensing instrumentation deployed at the two super-sites (Ersa and Lampedusa) and secondary sites (Minorca, Capo Granitola and the Barcelona and Granada EARLINET/ACTRIS stations), the additional AERONET/PHOTONS (AErosol RObotic NETwork/PHOtométrie pour le Traitement Opérationnel de Normalisation Satellitaire, http://aeronet.gsfc.nasa.gov/; Holben et al., 1998) and EARLINET/ACTRIS (European Aerosol Research Lidar Network/Aerosols, Clouds, and Trace gases Research InfraStructure Network, http://www. actris.net/; Pappalardo et al., 2014) network stations that we used, and the airborne observations obtained onboard the two French research aircraft (ATR-42 and F-20) and with sounding and drifting balloons. Sect. 4 is dedicated to presenting the main meteorological conditions, cloud cover and precipitation, which controlled the aerosol emission and transport during the period of observations. Sect. 5 presents some examples of results concerning the in situ and remote-sensing observations, in terms of aerosol physical, chemical, optical properties, and vertical profiles, as well as 1-D DRF SW and LW calculations. In the last part (Sect. 6), the modeling effort is presented. Different models are involved in this project, from high-resolution meteorological and chemistry transport models to regional climate models. The modeling results are used to describe the anthropogenic (carbonaceous, secondary inorganic and organic species) and natural (dust and sea-spray) loading and the estimated DRF at the regional scale for the period of experiment. An example of results of longer (inter-seasonal and inter-annual) aerosol-climate simulations is presented in Sect. 6, based on the work of Nabat et al. (2015a). 2 Overview of the surface observation network The regional experimental setup deployed in the western and central Mediterranean during the ChArMEx SOP-1a campaign is shown in Fig. 2. Two super-sites (Cape Corsica and Lampedusa) and 10 secondary sites (Fig. 2) have been used in this project. 2.1 The Cape Corsica and Lampedusa surface super-sites Two super-sites were fully equipped for documenting the aerosol chemical, physical and optical properties as well as their possible mixing and their vertical structure at local scale Figure 2. The regional experimental setup deployed in the western and central Mediterranean during the ChArMEx SOP-1a campaign. The two aircraft were based at Cagliari. (Table 1). The main characteristics of these two surface stations are presented here. The first station was located in Ersa on Cape Corsica (42◦5801000 N, 9◦2204900 E), near the northern tip of Corsica. This station was primarily instrumented for investigating polluted air masses transported over the Mediterranean basin from the highly industrialized regions of the Po Valley (Royer et al., 2010) and/or the Marseille– Fos–Berre (Cachier et al.,2005) zone and Rhone Valley. This ground-based remote station is located at an altitude of about 530m above mean sea level (a.m.s.l.) on a ridge equipped with windmills and benefits from a direct view to the sea over a northern sector of ∼270◦extending from the SW to the SE. The Cape Corsica peninsula is a remote site ensuring that the in situ measurements are not contaminated by local anthropogenic pollution. The Lampedusa super-site (35◦310500 N, 12◦3705100 E) was established at the Roberto Sarao station permanently operated by ENEA on the small island of Lampedusa (∼20km2), and it was augmented during the field campaign by the observations of the PortablE Gas and Aerosol Sampling UnitS (PEGASUS) mobile station operated by LISA. This surface station was mainly used for documenting very aged air masses in southwesterly flow from Europe, southern air masses from northern Africa (Tunisia, Algeria and Libya) possibly laden with mineral dust, as well as marine aerosols. It is situated on a cliff at about 45ma.m.s.l. on the northeastern tip of the island. The complete instrumentation deployed during the SOP1a experiment for both super-sites is detailed in Table 1. Briefly, it served to determine the complete aerosol physical, chemical and optical properties as well as vertical profiles, and to measure radiative fluxes (broadband SW and LW, and spectral SW). www.atmos-chem-phys.net/16/455/2016/ Atmos. Chem. Phys., 16, 455–504, 2016 462 M. Mallet et al.: Overview of the ChArMEx/ADRIMED summer 2013 campaign Table 1. List of the instrumentations deployed over the two super-sites (Ersa and Lampedusa) during the SOP-1a experiment for the characterization of physical, chemical and optical properties of aerosols, vertical profiles, columnar-averaged properties and radiation measurements. Meteorological parameters and gas concentrations are not included in this table. Ersa Lampedusa Instruments Frequency Instruments Frequency Number concentration 1 CPC (0.01–3µm) Continuous (10) 1 W-CPC (0.01–3µm) Continuous (20) CCN concentration 1 CCN counter Continuous 1 CCN counter Continuous Mass concentration 1 PM2.5Continuous 1 PM40 (TEOM) Continuous 1 PM10 Continuous Number size distribution 1 OPC (0.3–5µm) Continuous 2 GRIMM (0.25–32µm) Continuous 1 APS (TSI) Continuous 1 APS (TSI) (0.5–20µm) Continuous 1 SMPS (3–300nm) Continuous 2 (dry/ambient) SMPS Continuous Mass size distribution 2 Impactor DEKATI (13 stages) 48h 2 Impactor DEKATI (13 stages) 48h 1 Impactor Nano-MOUDI 24h PM1composition 1 PILS Continuous AMS (Aerodyne) Continuous 1 PILS Continuous PM10 composition 1 FAI Hydra Sampler 12h Mass BC concentration 1 (7-λ) aethalometer Continuous 1 PSAP Continuous (1h) 1 MAAP Continuous Vertical profiles 1 (1-λ355nm) Leosphere Continuous 1 (1-λ) Leosphere ALS 300 Continuous (200) 2 (3-L) ENEA/Univ. of Rome lidar Continuous (10) microwave radiometer (p,T, RH) Continuous (150) radiosondes on event Scattering coefficient 1 (3-λ) TSI nephelometer Continuous (10) 1 (3-λ) TSI nephelometer Continuous (10) (450-550-700nm) (450-550-700nm) Absorbing coefficient 1 (7-λ) aethalometer Continuous 1 (7-λ) aethalometer Continuous (370-420-490-520-660-880-950nm) (370-420-490-520-660-880-95nm) Extinction coefficient 1 (1-λ) (860nm) PAX Continuous (10) Column optical properties 1 (9-λ) AERONET/PHOTONS Continuous (150for AOD) 1 (9-λ) AERONET/PHOTONS Continuous (150for AOD) 2 (12-L) MFRSRs Continuous (15s) Mineral aerosol deposition 1 CARAGA Continuous (7-day) 1 CARAGA Continuous (7-day) Downward shortwave irradiance 1 pyranometer Continuous (30s) 1 (CMP 21) pyranometer Continuous (30s) Downward longwave irradiance 1 pyrgeometer Continuous (30s) 1 (CGR4) pyrgeometer Continuous (30s) Downward window (8–14µm) irradiance 1 modified CG3 pyrgeometer Continuous (60s) Direct solar radiance 1 CHP1 Pyrheliometer Continuous (30s) Direct spectral solar radiation 1 PMOD Precision SpectroRad. Continuous (30s) Spectral downward global solar irradiance 1 HyperOCR spectrometer Continuous (30s) Spectral downward diffuse solar irradiance 1 HyperOCR spectrometer Continuous (30s) Spectral direct solar irradiance 1 spectroradiometer Continuous (60s) Downward spectral actinic flux 1 Diode array spectrometer Continuous (60s) Atmos. Chem. Phys., 16, 455–504, 2016 www.atmos-chem-phys.net/16/455/2016/ M. Mallet et al.: Overview of the ChArMEx/ADRIMED summer 2013 campaign 463 2.1.1 In situ measurements at super-sites Both super-sites measured the mass concentration online using tapered element oscillating microbalance (TEOM) analyzers. The number size distribution of particles are also measured, including fine and coarse fractions (radius ranges and corresponding instruments are reported in Table 1). The aerosol composition was derived from chemical analyses of filters and cascade impactors (DEKATI and MOUDI) with time resolution varying from 12 to 48h (depending on the aerosol load), but also from high time resolution online measurements by an ACSM (Aerosol Chemical Speciation Monitor) at Ersa, a C-TOF-AMS (Time of Flight Aerosol Mass Spectrometer) at Lampedusa, and two PILS (Particle Into Liquid Sampler) systems at both sites (Table 1). The original observations of aerosol chemical properties obtained from PM10-PILS instrument at Ersa are detailed in Claeys et al. (2016). Concerning aerosol optical properties, scattering and absorption coefficients (at wavelengths listed in Table 1) have been estimated for both super-sites using a 3-λ nephelometer and a 7-λaethalometer, respectively. At Ersa station, the extinction coefficient (at 870nm) was also estimated using a Photoacoustic Extinctiometer (PAX) instrument, while it has been estimated at 2-λ(450 and 630nm) at Lampedusa using two Cavity Attenuated Phase Shift Spectroscopy (CAPS) systems. Additional in situ measurements were performed at the Ersa station. The mixing state of fine particles (at the two selected diameters of 50 and 110nm in dry conditions) has also been estimated from their hygroscopic behavior using a VHTDMA (volatilization and humidification tandem differential mobility analyzer) system (Johnson et al., 2004). In parallel, a TSI (model 3800) aerosol time of flight mass spectrometer (ATOFMS) (Gard et al., 1997) was used to measure the size-resolved chemical composition of single particles in the vacuum aerodynamic diameter (dva)size range 100– 3000nm. 2.1.2 Remote-sensing and radiation measurements at super-sites A Leosphere Raman lidar model RMAN510 was setup at low altitude (∼11m above sea level) in the small village of Macinaggio (42◦5704400 N, 9◦2603500 E) located on the eastern coast of Cape Corsica. The lidar was operated at about 6km east of the Ersa station and less than 700m from the shoreline. The RMAN510 uses a laser emitting at 355nm. It measures the total and polarized backscatter at 355nm and the Raman nitrogen signal at 387nm at night-time. A second ALS300 510 lidar system has been deployed in Lampedusa (Formenti et al., 2016) as well as a more powerful University of Rome–ENEA homemade lidar measuring backscatter at 532 and 1064nm (Di Iorio et al., 2016). The main characteristics of lidar systems are provided and detailed in Table 1. At each station, a multi-wavelength sun-photometer from the AERONET/PHOTONS network was operated, allowing the operational retrieval of column-integrated AOD at 340, 380, 440, 500, 675, 870, and 1020nm (and also at 1650nm at Ersa) and aerosol optical and microphysical properties such as the single scattering albedo, refractive index and particle size volume distribution (Dubovik and King, 2000; Dubovik et al., 2000, 2002, 2006). The Ersa sun-photometer has been positioned since June 2008 near the navy semaphore on the northwestern tip of Cape Corsica (43◦0001300 N, 9◦2103300 E, alt. ∼75ma.m.s.l.) at about 4.2km NNW of the Ersa surface station. Both super-sites were complemented by a pyrgeometer and a pyranometer for monitoring longwave and shortwave downward fluxes measurements, respectively. Additional radiation measurements were performed at Lampedusa (Table 1). Spectral measurements of global, diffuse, and direct radiation were carried out with other instruments deployed by ENEA and the Physikalisch-Meteorologisches Observatorium Davos, World Radiation Center (PMOD/WRC, Switzerland). Multi-filter rotating shadowband radiometer observations were carried out jointly with AERONET sunphotometer (di Sarra et al., 2015) and allowed the derivation of the AOD at several wavelengths. By combining these two measurements, a long-term series of AOD, started in 2001, was obtained. Measurements of the spectral actinic flux, allowing the determination of the photolysis rates (Mailler et al., 2015), were carried out with a diode array spectrometer. Measurements of broadband irradiance included a CG3 pyrgeometer sensitive to radiation in the atmospheric infrared window. Finally, the total ozone and spectral UV irradiance were obtained with a Brewer spectrophotometer. Several radiosondes were also launched from Lampedusa during the SOP-1a, and vertical profiles of temperature and humidity were continuously measured by a microwave radiometer. 2.2 The secondary sites 2.2.1 Montesoro station The Cape Corsica station was complemented by an additional remote-sensing setup at the peri-urban air quality station of Montesoro, southward of Bastia at about 45ma.m.s.l. (Léon et al., 2015), including a Leosphere model EZ lidar operating at 355nm (42◦4001700 N, 9◦2600500 E) and a Cimel AERONET/PHOTONS sunphotometer (42◦4001900 N, 9◦2600600 E). In addition, some airquality parameters were monitored by Qualitair Corse, including PM2.5and PM10. This station is less than 1km from the shore on the northeastern coast of Corsica, about 32km south of Macinaggio. www.atmos-chem-phys.net/16/455/2016/ Atmos. Chem. Phys., 16, 455–504, 2016 470 M. Mallet et al.: Overview of the ChArMEx/ADRIMED summer 2013 campaign Figure 5. Total AOD (550nm) obtained from the MODIS, MISR and SEVIRI (sea only) sensors for the June–July 2013 period. The AERONET/PHOTONS AOD is also indicated. Biscay induced a strong southerly flow at 700hPa between the Balearic Islands and Corsica associated with large dust optical depth concentrated in this zone as shown by SEVIRI AOD (AOD of 0.3–0.4) for 19 June (Fig. 6). This period of the SOP-1a corresponds to the two ATR-42 flights 33 and 34 (Fig. 3). After 20 June, this low-pressure system moved eastward, generating a trough located between France and Italy, and inducing a waving westerly flow over the northwestern Mediterranean. As a result, the aerosol loading over the western basin decreased between 21 and 24 June, but the westerly (northerly) winds observed at 700hPa in Minorca (Ersa) (Fig. 8) reinforced the transport of dust aerosols over the central basin and the Lampedusa station (where winds were from the northwesterly direction at 3km height). These meteorological conditions lead to an increase of the dust optical depth over the central Mediterranean as shown by the SEVERI instrument and AERONET/PHOTONS data. Between 25 and 29 June, a northwesterly flow set up between the Gulf of Lion and Sicily. The vertical profiles of the wind direction reveal a remarkable transition on 29 June with significant changes in direction from westerlies to north-northwesterlies, notably over the Minorca and Ersa stations above 850hPa. The 700hPa geopotential field on 29 June at 12:00UTC from the ERA-Interim reanalysis (Fig. 7) shows a maximum over the Atlantic Ocean, whereas a deep low-pressure system was located over southern Algeria. This strong geopotential gradient leads to intense northerly to northwesterly winds at 700hPa over the western basin, leading to significant AOD over Libya (AOD of 0.4–0.5) and the Alboran Sea (AOD of 0.5–0.6) as shown in Fig. 6. These meteorological conditions lead to low dust optical thickness over the central Mediterranean as observed by AERONET/PHOTONS data. Finally, during the last period of the SOP-1a experiment (30 June– 5 July), weather conditions became more anticyclonic over the region, while low systems were confined to northern Europe. Figure 8 shows northwesterly winds in the whole troposphere in Lampedusa and Minorca, limiting the presence of dust aerosols to the southern part of the northwestern Mediterranean. 4.2 Surface temperature, cloud cover and precipitation In terms of surface temperature, which is one of the most important meteorological variables that control biogenic or biomass burning aerosol emissions over the EuroMediterranean region, the summer 2013 was mostly characterized by moderate values as shown in Fig. 9. Indeed, during the SOP-1a period, surface temperatures (in ◦C and at 12:00UTC) derived from NCEP reanalysis (Kalnay et al., 1996) for different days reveal moderate values, especially over the western Mediterranean region (southwestern France and Spain). One can observe temperatures of about 15–20◦C (at 12:00UTC) over Spain and Portugal, which are one of the main regions of the Mediterranean where large fire events occur. In addition, part of France was also characterized by moderate surface temperature but slightly higher than over Spain especially over northeastern regions. A strong west to east gradient is observed over Europe with strongest values over the eastern regions (around 30◦C over Greece and the Balkans) compared to the western basin. A similar conclusion is obtained over the Mediterranean Sea with differences of about 5◦C between the eastern (around 25◦C for the SOP-1a period) and the western (around 20◦C) basins. Among other factors (such as cloud fraction and shortwave radiations), such moderate surface temperatures do not create favorable meteorological conditions to produce intense Mediterranean biomass burning events and/or significant production of secondary organic and inorganic aerosols. Concerning smoke aerosols, GAFS-V1 emission data, analyzed for the SOP-1a period, do not reveal important primary BC and OC flux emissions (not shown). This is consistent with the APIFLAME biomass burning emission estimates (Turquety et al., 2014) data as reported by Menut et al. (2015). During the SOP-1a, the cloud cover retrieved over the Euro-Mediterranean region (excluding the Mediterranean Sea) from CRU (Climate Research Unit) data (Harris et al., 2013) (Fig. 10) indicates the largest values (between 75 and 95%) over the French, Benelux and eastern Europe regions. In parallel, southern France, as well as western Spain and the Atmos. Chem. Phys., 16, 455–504, 2016 www.atmos-chem-phys.net/16/455/2016/ M. Mallet et al.: Overview of the ChArMEx/ADRIMED summer 2013 campaign 471 Figure 6. AOD MSG/SEVIRI observations for 5 different days during the SOP-1a experiment (16, 19, 22, and 29 June and 3 July). The AERONET/PHOTONS AOD is also indicated. Balkans are characterized by moderate cloud cover with values around 50–60% for June 2013. Over the Mediterranean coast, the cloud cover strongly decreases for most of countries, with values lower than 40%. Such spatial cloud cover (observed during the SOP-1a) over the Euro-Mediterranean could limit the photochemical processes over the main anthropogenic sources (such as the Benelux and Po Valley) and the associated production of secondary aerosols. This could explain for a part the low to moderate contribution of fine anthropogenic particles to the total atmospheric loading during the SOP-1a. In parallel, the mean precipitation (averaged for June 2013), obtained from the TRMM (Tropical Rainfall Measuring Mission) instrument over land and sea (CRU observations are only available over land; see Fig. 10), are found to be very heterogeneous over the EuroMediterranean continental region, with some important values over the Balkans, Alps and eastern Europe (from 100 to 250mm for the month of June 2013) and moderate values over Italy, Croatia, western France and Benelux (80 to 100mm, as shown in Fig. 11). Over the Mediterranean Sea, southern Spain and northern Africa, the precipitation was smaller, with most of the values lower than 20mm during the SOP-1a. To summarize, this global view of the synoptic situation, cloud cover and regional precipitation patterns indicates that the meteorological conditions during the experimental campaign were favorable to moderate mineral dust emissions, associated with a weak contribution of anthropogenic aerosols over the western basin. This important characteristic of the SOP-1a is well observed in Fig. 12, which indicates the AOD anomalies (calculated for the period 2000–2013) of summer 2013 compared to all AOD summer derived from MODIS and MISR data. Indeed, negative AOD anomalies of about −0.05 are found over the western Mediterranean basin for the summer of 2013, both from MODIS and MISR observations. To conclude, it appears that the period of observations during the SOP-1a was characterized by aerosol concentration slightly lower but in the same range of magnitude that usually observed during summer over the western Mediterranean. The level of aerosol concentration was found to be moderate but allows investigating several dust and seaspray events as well as an interesting intense biomass burning plume advected from North America. 4.3 An aged smoke plume advected over Europe During the SOP-1a, several large forest fires occurred in North America (Colorado, Alaska, Canada) from 17 to 24 June 2013, as identified by the MODIS instrument. The absorbing aerosol index produced from GOME-2 by KNMI (http://www.temis.nl/aviation/aai-pmd-gome2b.php? year=2013) shows that a large smoke plume crossed the North Atlantic and reached western European coasts on 25 June. Main fire areas, with fire radiative power higher than 50MW (Shroeder et al., 2010), have been detected over Canada (Ancellet et al., 2015). Average MODIS AOD during the same period (23 to 28 June 2013) indicate values as high as 1 over the Atlantic Ocean, suggesting that a significant fraction of the aerosol produced by the fires was transported to western Europe during the ChArMEx/ADRIMED field campaign. To investigate whether the western Mediterranean has been impacted by these fires, a forward simulation of the FLEXPART (Ancellet et al., 2015) Lagrangian plume dispersion model has been conducted to quantify the spatial extent of the fire plume transport for 11 days. Fire emissions www.atmos-chem-phys.net/16/455/2016/ Atmos. Chem. Phys., 16, 455–504, 2016 472 M. Mallet et al.: Overview of the ChArMEx/ADRIMED summer 2013 campaign Figure 7. Geopotential and wind intensity at 700hPa (ERA-Interim reanalysis), and mass dust concentration (in mgm−3) simulated from the ALADIN-Dust model for 6, 19, 22, 29 June and 2 July at 12:00UTC. areas were identified by MODIS observations over several locations in Canada and Colorado. The aerosol mass is emitted in the transport model from 17 to 28 June in a 3km layer as suggested by the CALIOP lidar observations over Canada. The biomass burning plume reaches much lower latitudes over Europe, down to the western Mediterranean 4–10 days after the emission in Canada. During the SOP-1a, the plume was mainly present in the altitude range of 2.5–4.5km and has been sampled by many remote-sensing and in situ instruments on 27 and 28 June: at the Minorca and Cagliari surface stations, and between Sardinia and Lampedusa onboard the ATR-42 aircraft. 5 Overview of aerosol physical-chemical-optical properties, vertical profiles and local direct radiative forcing 5.1 Aerosol physical and chemical properties 5.1.1 Aerosol mass and number concentration at the two super-sites First, PM concentrations between the two different stations are reported in Fig. 13, which reports the daily time series of PM1and PM10 at Ersa, as well as PM10 and PM40 at Lampedusa. The results indicate a significantly higher mass concentration at Lampedusa compared to Ersa. Indeed, the mass concentration observed at Lampedusa is comprised between 10 and 30µgm−3, with a mean of 21µgm−3, which is two times higher than the averaged PM10 (∼9µgm−3) measured at Ersa. One can note the significant peak of PM40 (maxima of 75µgm−3) at Lampedusa during the 24 to 26 June period that corresponds to a significant production of primary Atmos. Chem. Phys., 16, 455–504, 2016 www.atmos-chem-phys.net/16/455/2016/ M. Mallet et al.: Overview of the ChArMEx/ADRIMED summer 2013 campaign 473 Figure 8. Wind profiles between 1000 and 200hPa during the SOP1a experiment for three different sites (Ersa, Lampedusa and Minorca) simulated from the ALADIN-Dust model. The wind intensity (in ms−1) is also reported at the different stations. marine aerosols. Finally, the PM1concentration at Ersa is found to be almost constant during the period of the campaign, with a mean value of 6µgm−3. In order to take into account the difference of altitudes between the two sites of Lampedusa and Ersa, we have applied a correction factor to PM10 observed at Ersa (530m) for estimating a new PM10 concentration corresponding to the altitude of Lampedusa. In that sense, we have applied the logarithmic law reported by Piazzola et al. (2015) using a value of 0.75 for the factor s to correct the mass concentration of sea-spray aerosols only. The calculated mean value of PM10 is about 12µgm−3(Fig. 13), closer to the mean value observed at Lampedusa (21µgm−3). In addition, the background aerosol number concentrations (for Dp>0.01µm) observed within the boundary layer in Corsica averaged ∼2000cm−3(not shown). The lowest concentrations (∼200cm−3) resulted from aerosol activation to cloud droplets, and scavenging from cloud droplets and rain drops, while high concentrations as high as 10000cm−3were observed during pollution events from continental European air masses. The number concentrations showed a diurnal cycle suggesting that the site was situated within the marine boundary layer during daytime and within the free troposphere during night-time. The analysis of the diurnal variation of the particle number size distribution further indicates that nucleation events also increased the particle number concentration during daytime, about one-third of the time (Sellegri et al., 2016). The periods of high aerosol number concentrations detected between 12 and 25 June were also dominated by a single mode with diameters between 30 and 150nm. The small Aitken mode (dg<50nm) associated with pollution events suggests a relatively fresh aerosol that has been formed during transport from the European continent. The largest mode (dg∼150nm) occurred during the dust event on 18 June. 5.1.2 Columnar particle volume size distribution We have used the column-integrated particle size volume distributions derived from AERONET/PHOTONS sky radiance measurements (Dubovik et al., 2000). These size distributions allow investigation of the changes in aerosol size distribution between different stations during the SOP-1a and over the western basin. Four different stations have been studied, which include the two super-sites of Lampedusa and Ersa, as well as the aircraft and balloon base stations, Cagliari and Cap d’En Font, respectively. Daily volume size distributions for both sites are represented in Fig. 14, as well as the averaged (red curve) size distribution for the whole period (1 June to 5 July) and the number of observations. In addition, the mean values of the volume radius, concentration of fine and coarse mode and the standard deviations of the volume size distribution are reported in Table 6. It should be noted that the scales of the yaxis are different for each figure. One can note the bimodal size distribution for both stations with large spread of radius values, especially for the coarse mode. The most important concentrations are obviouslyobserved in Lampedusa, near the mineral dust sources, with maxima of ∼0.12µm3µm−2for the coarse mode. In parallel, the lowest concentrations are observed at the Ersa station due to the absence of intense polluted photochemical or smoke aerosol events over southern France and Italy during the SOP-1a. In that sense, the mean contribution (red curve) of the coarse mode to the aerosol volume size distribution appears to be predominant at most sites, except at the Ersa station. However, the inclusion of the corrected factor (Piazzola et al., 2015) for taking into account the altitude of the Ersa site reduces slightly the differences in the concentration of the coarse mode with the Lampedusa station (see Table 6). This point is well noted for the Cap d’En Font station, where the concentration of each modes appear as equivalent, due to the absence of pollution from the Iberian Peninsula during the period of observations. For this site, it is interesting to note the intense peak of concentration (0.08µm3µm−2) for 27 June, which is due to the transport of an important smoke plume over the Mediterranean (see Ancellet et al., 2015; Chazette et al., 2015). Finally, the contribution of the coarse www.atmos-chem-phys.net/16/455/2016/ Atmos. Chem. Phys., 16, 455–504, 2016 474 M. Mallet et al.: Overview of the ChArMEx/ADRIMED summer 2013 campaign Figure 9. Surface temperature (at 12:00UTC) obtained from NCEP re-analysis for 16, 17, 19, 20 and 22 June. mode clearly increases for the two other, more southern Italian sites of Cagliari and Lampedusa, which are more affected by the mineral dust compared to Ersa and Cap d’En Font. The variability of AERONET products collected over a period of 4 years at Ersa and Palma de Mallorca, near Cap d’En Font, is reported in Sicard et al. (2015b). It is interesting to note the variability (±0.05) in the derived size of the coarse mode at Lampedusa (see Table 6), which will be analyzed with regards to dust sources in a future study. The derived volume concentrations over these two stations highlight the moderate dust activity occurring during the SOP-1a experiment, when compared to stations under high dust conditions. As an example of comparisons, Dubovik et al. (2002) reported a large range of concentration for the coarse mode for dusty sites (such as Cape Verde or Solar Village), which are characterized by larger concentrations, close to 0.30µm3µm−2. In parallel, the Bahrain (Persian Gulf) AERONET station is characterized by a concentration of 0.14–0.15µm3µm−2. 5.1.3 Particle size distribution during transport Figure 15 presents an example of the evolution of the aerosol particle number concentrations in the 19 particle size classes of the LOAC instrument as measured along the northward trajectory of BPCL balloon B74 from Minorca to the French coast (see Fig. 4). The balloon was launched at 09:46UTC on 16 June 2013 during a moderate desert dust event shown Atmos. Chem. Phys., 16, 455–504, 2016 www.atmos-chem-phys.net/16/455/2016/ M. Mallet et al.: Overview of the ChArMEx/ADRIMED summer 2013 campaign 475 Figure 10. Monthly cloud cover and precipitation (over land only) derived from the Climate Research Unit (CRU) data for June 2013. Table 6. Main aerosol volume size distribution characteristics: rvf (µm), σf,rvc (µm), σc,Cvf, and Cvc, for the four different AERONET/PHOTONS stations: Ersa, Lampedusa, Cagliari and Cap d’En Font. Cvi denotes the particle volume concentration, rvi is the median radius, and σiis the standard deviation. Each average value in the table is accompanied by its standard deviation (this is not an accuracy of the retrieval). As mentioned in the text, the concentration of the coarse mode at Ersa has been corrected to be comparable to results at other stations closer to the sea surface, using the logarithmic law proposed by Piazzola et al. (2015). Ersa Ersa corrected Lampedusa Cagliari Cap d’En Font Number of observations 25 18 20 17 rvf (µm) 0.16±0.02 n/a 0.14±0.01 0.15±0.03 0.17±0.03 σf0.43±0.03 n/a 0.50±0.06 0.46±0.04 0.45±0.04 rvc (µm) 2.49±0.43 n/a 2.36±0.48 2.52±0.28 2.48±0.30 σc0.69±0.03 n/a 0.68±0.05 0.71±0.04 0.71±0.04 Cvf (µm3/µm2)0.02±0.01 n/a 0.02±0.01 0.02±0.01 0.02±0.01 Cvc (µm3/µm2)0.03±0.01 0.04 0.08±0.05 0.05±0.03 0.04±0.03 Figure 11. Same as Fig. 10 but for the Tropical Rainfall Measuring Mission (TRMM) precipitation observations. at the top of Fig. 6 (AERONET-derived AOD at 550nm of 0.15). It drifted at a constant altitude of ∼2.1km at the bottom of the African dust layer observed with the WALI lidar at Minorca (not shown; see Chazette et al., 2015), and was automatically forced to land on the sea before reaching the coast south of Marseilles, after a 12h flight of 368km. The dominant mineral dust nature of the particles was confirmed by the LOAC particle typology measurements (Renard et al., 2015b). The figure illustrates that LOAC has detected large particles of up to 50µm in diameter, although the plume originated from northern Africa a few days before (Renard et al., 2015b). The concentrations of particles remained relatively constant during the flight, suggesting either no significant sedimentation of the largest particles during the flight or compensation by particles coming from above. BPCL balloon B70 launched a few minutes later drifted at an upper altitude of ∼3.1km and followed a different trajectory towards the east (Fig. 4), but showed a quite similar extended particle size range with larger concentrations in almost all channels except the extremes (not shown). The four other drifting balloons launched in the dust layer during this event on 17 and 19 June (Table 5) did confirm the presence of very large particles (>20µm), which cannot be reported by the AERONET particle size distribution retrieval algorithm (Hashimoto et al., 2012). In addition, observations of large particles (>15µm) were systematically found during www.atmos-chem-phys.net/16/455/2016/ Atmos. Chem. Phys., 16, 455–504, 2016 476 M. Mallet et al.: Overview of the ChArMEx/ADRIMED summer 2013 campaign MISR observations MODIS observations Figure 12. AOD (550nm) anomaly for summer 2013 estimated from the MODIS and MISR sensor data. Figure 13. Time series of daily PM mass concentrations estimated at the Lampedusa (PM40 and PM10) and Ersa (PM1and PM10) super-stations. Problems in PM10 data acquisition that occurred at Ersa explain the gaps. PM10 Ersa corrected curves correspond to PM10 estimated at an altitude of 45m to be comparable with Lampedusa results, following the logarithmic law provided by Piazzola et al. (2015) (see text in Sect. 5.1.1 for details). all other LOAC balloon flights drifting in African dust layers, which will need further analysis to better understand the process that can maintain such large particles in suspension during several days. Concerning the aerosol microphysical properties, aircraft observations have allowed us to investigate the vertical structure of aerosol size distribution showing particles characterized by large sizes (>10µm in diameter) within dust plumes. In addition, in most of the cases, a coarse mode of mineral dust particles, characterized by an effective diameter Deff,c ranging between 5 and 10µm, has been detected within the dust layer located above the MBL. Such values are found to be larger than those referenced in dust source regions during FENNEC, SAMUM1 and AMMA, as well as measurements in the Atlantic Ocean in the Cape Verde region during SAMUM-2 and at Puerto Rico during PRIDE. The complete analysis of aerosol size distribution is detailed in Denjean et al. (2015). 5.1.4 Aerosol chemical composition In terms of aerosol chemical properties, an example of averaged mass size distributions for carbonaceous (elemental and organic carbon, EC and OC) species (mass size distributions of inorganic and mineral dust aerosols are not shown) obtained at Ersa from a 12-stage cascade impactor (DEKATI system; see Table 1) is reported in Fig. 16. The aerosol chemical properties obtained from the PILS instrument at Ersa are detailed in Claeys et al. (2016). As mentioned in Table 1, the measurements were obtained by using a 2-day collection period in order to obtain a sufficient aerosol mass on filters for chemical analyses. This system provides the speciation of the mass size distribution, including fine and coarse fractions. Such information is very useful for deriving optical properties using Mie calculations (Mallet et al., 2011) for the main particle types (sulfates, ammonium, nitrates, sea spray, dust, black and organic carbon). This provides crucial information on key radiative properties that are classically used in regional climate models (mass extinction efficiencies, SSA and asymmetry parameter). Furthermore, it allows one to assess the spectral dependence of radiative properties, which cannot always be estimated from in situ instrumentation. Concerning OC (blue curves), observations clearly report a bi-modal mass size distribution with two different peaks for the majority of cases. The first (almost constant) peak is found in the 0.4–0.5µm size range in diameter and more occasionally a second one occurs in the coarse fraction around 3µm. Compared to the few available data over the western Mediterranean, these mass size distributions are found to be different from those obtained over southern France, especially for the accumulation mode. Indeed, during the ESCOMPTE experiment in southern France, Mallet Atmos. Chem. Phys., 16, 455–504, 2016 www.atmos-chem-phys.net/16/455/2016/ M. Mallet et al.: Overview of the ChArMEx/ADRIMED summer 2013 campaign 477 Figure 14. AERONET/PHOTONS volume size distribution derived at four different stations: Ersa, Lampedusa, Cagliari and Cap d’En Font (the red curve represents the mean of observations). The characteristics of the volume size distribution are provided in Table 6. et al. (2003) also observed a bi-modal size distribution for OC aerosols but with a finer accumulation mode observed in the 0.1–0.2µm size range. Differences between the two observations is likely due to the proximity of anthropogenic sources during the ESCOMPTE experiment compared to the Ersa station, where the possible ageing of carbonaceous particles could affect the size of aerosols. On the contrary, the coarse mode of OC appears in the same range of size, around 3µm, for both experiments. Compared to data obtained in the eastern Mediterranean basin, the OC mass size distributions are in good agreement with those estimated by Sciare et al. (2003) in Crete during the MINOS campaign, with two modes around 0.4 and 3µm. The EC (green curves in Fig. 16) mass size distribution is also characterized by a bi-modal size distribution, with two modes well correlated with the mass size distribution of OC, except for the 16–19 June period (dust episode), where the size of EC fine mode is higher (∼0.5–0.6µm) than OC aerosols, the EC coarse mode remaining similar at ∼3µm. This reveals a possible external mixing of carbonaceous aerosols for this event. It should also be noted that the EC concentrations observed at the Ersa station are logically (due at least to the altitude of the station and the absence of intense pollution during the SOP-1a; see Sect. 4) lower (0.39µgm−3) than EC concentrations (PM2.1) reported by Eleftheriadis et al. (2006) from the eastern Mediterranean during the summer season (0.6µgm−3) in July 2000. The same ascertainment is obtained on OC concentrations with higher values (4.2µgm−3) reported by Eleftheriadis et al. (2006) compared to observations at Ersa (1.5µgm−3). Concerning the modes of the OC and EC particle mass size distributions, the two identified modes detected in Ersa are consistent with those reported by Mallet el al. (2011) at the Porquerolles coastal island (southeastern France), who also detected two (fine and coarse) different modes of the mass size distributions for EC (0.3–0.4 and 4–6 µm) and OC (0.3 and 5–6µm) aerosol particles. In most cases, we observed at Ersa lower concentrations of EC particles for both modes compared to OC aerosols. The mass of OC and BC observed during the SOP-1a, for both modes, are found to be equivalent with those observed by Sciare et al. (2003) in Crete in summer 2001. They report mean values of 0.30 and 0.15µgm−3for fine OC and BC, respectively. During the MINOS experiment, the mean concentrations for OC and BC coarse modes were about 0.1 and 0.02–0.03µgm−3, what is also consistent with the observations at Ersa. Finally, the mass concentrations obtained for each mode at Ersa are logically lower than those obtained during the ESCOMPTE experiment, located much closer to pollution sources. For example, EC and OC fine-mode concentrations were, respectively, between 0.8 and 2.8µgm−3 and between 3.1 and 6.9µgm−3during ESCOMPTE (Mallet et al., 2003). In addition, and as discussed in Sects. 4.1 and www.atmos-chem-phys.net/16/455/2016/ Atmos. Chem. Phys., 16, 455–504, 2016 478 M. Mallet et al.: Overview of the ChArMEx/ADRIMED summer 2013 campaign Figure 15. Particle size distribution measured with a LOAC during the ∼12h flight of BPCL balloon B74 drifting from Minorca towards Marseilles (see the trajectory in Fig. 4). The first and last 20min correspond to the ascending and descending phases of the quasi-Lagrangian flight that occurred at a constant altitude of 2091±10m. 4.2, the meteorological conditions (surface temperature, meteorological synoptic situations) observed during the SOP1a campaign were not favorable to producing large concentrations of polluted or smoke aerosols, compared to the ESCOMPTE campaign, where AOD as large as 0.3–0.5 (in the visible range) has been observed due to important concentrations of anthropogenic-polluted particles. It should be noted that, in parallel to filter analyses, higher time-resolved observations from the PILS systems have been deployed at the two stations of Lampedusa and Ersa (Claeys et al., 2016) during the SOP-1a. In parallel to filter chemical analysis, over 700000 single particle mass spectra were generated by the A-TOFMS instrument during the sampling period (not shown). A Kmeans algorithm (K=80), as described in detail by Healy et al. (2010) and Gross et al. (2010), was used to classify aerosol mass spectra into different particle classes. More than 40 distinct ATOFMS particle classes were identified and subsequently grouped into 8 general categories for clarFigure 16. EC and OC (48h mean) aerosol mass size distributions obtained at Ersa from the DEKATI impactor instrument for the entire SOP-1a period. ity. Elemental carbon-containing particles dominated the data set (55% of total spectra), followed by K-rich particles (30%) and sea spray (7%). The remaining particle categories include organic carbon (OC)-containing (3%), Atmos. Chem. Phys., 16, 455–504, 2016 www.atmos-chem-phys.net/16/455/2016/ M. Mallet et al.: Overview of the ChArMEx/ADRIMED summer 2013 campaign 479 Figure 17. Vertical profiles of rBC concentrations estimated from the SP2 instrument for five different zones (Granada, Minorca, Lampedusa, southern France and Ersa). trimethylamine (TMA)-containing (3%), shipping (2%), Fecontaining (0.5%) and Ca-containing (0.3%). EC particles dominated the first third of the sampling period, decreased noticeably for approx. 6 days and then dominated the rest of the sampling period again. In contrast, K-rich particle (associated with biomass burning and dust) numbers were high only for the latter half of the campaign, with a peak on 27– 28 June. The profiles of these two particle categories suggest transport from regional sources. Sea-spray particle numbers were at their highest during the period where EC particles were at their lowest, and were generally low when EC particle numbers were high. OC-containing particles were present during the same period K-rich numbers peaked, suggesting an association with the transport of biomass burning particles. TMA particles were present in low numbers throughout the sampling period, suggesting a less regional source, independent of the air masses influencing EC and sea-spray particle occurrence. The same can be said of Feand Cacontaining particles, likely to be local dust, while shipping particle numbers were slightly higher during the first half of the sampling period. Finally and concerning the aerosol chemical properties, an interesting aspect of the obervations deployed during the SOP-1a concerns the rBC concentrations obtained from the SP2 instrument onboard the ATR-42. Despite its importance, studies on rBC were until now limited to surface-based measurements in the Mediterranean region. Measurements of vertical distribution of rBC concentrations provide crucial information for assessing the rBC radiative effects in the region. Figure 17 shows the vertical distributions of rBC mass concentrations measured by the SP2 in the five areas (Granda, Minorca, Lampedusa, southern France and Ersa). For the different vertical soundings, rBC mass concentrations ranged between 20 and 690ngm−3close to the surface. The surface rBC concentrations were generally less than 200ngm−3, typical for continental and regional background sites in the western Mediterranean basin (Ripoll et al., 2015). The lowest surface concentrations of rBC (∼20ngm−3) were found in southern France over the open sea, with almost no local contribution of anthropogenic aerosols. Maxima surface concentrations (∼690ngm−3) were recorded over Granada where frequently heavy traffic emissions are occurring. These observations were obtained between 07:15 and 07:45UTC, when the convection was not fully developed, which probably did not favor the vertical transport of local emissions over Granada. A prominent feature in vertical profiles is the presence of significant concentrations of rBC up to 5–6km altitude. Therefore the regional transport of rBC particles was not only limited to the MBL but occurred also at higher altitude. In most of the observed cases, the rBC vertical distribution in the free troposphere reveals a strongly stratified structure characterized by either single isolated plumes or more uniform layers. It is worth noting the presence of rBC layers above the MBL in the open sea that could be attributed to convective transport from distant sources. Only in few observed cases, rBC mass concentration decreased monotonically with increasing altitude, most likely due to vertical transport of air masses from surface to higher heights. 5.2 Aerosol optical properties 5.2.1 In situ optical properties at the surface Figure 18 reports the (daily mean) time series of nephelometer observations obtained at the surface for the Ersa and Lampedusa stations. Daily scattering coefficients (at the three nephelometer wavelengths of 450, 550 and 700nm) are reported, as well as the scattering Angström exponent (AE) calculated between 450 and 700nm. At 550nm and at Ersa, the scattering coefficient presents a significant variability during the SOP-1a, with peaks of about 35–40Mm−1 during the dust event (19–20 June) transported over Corsica, associated with low values (15Mm−1) for certain periods of time, as for 21–22 June. The mean scattering coefficient (at 550nm) is 24Mm−1. Such scattering coefficient values are comparable to observations reported by Vaishya et al. (2012) at the Mace Head station for Atlantic marine air, with scattering coefficient (at 550nm) ranged between 10 and 25Mm−1 during the summer period. In terms of scattering spectral dependence, the calculated scattering AE is found to be almost constant, with AE∼1.5–2 and a mean value of 1.71 (indicating that scattering is mostly dominated by fine aerosols) during the SOP-1a, except for 23–24 June. The lowest values (AE∼0.3–0.5) observed during this period are the result of a large contribution of coarse sea-spray aerosols (Claeys et al., 2016) due to moderate (5ms−1) westerly winds (see Fig. 8) at the Ersa station, which is also observed from the filter chemical size-resolved analyses and detected on the ATOFMS and VHTDMA data. In parallel, we observe that the dust event occurring in Ersa on 18–20 June is not correlated with low scattering AE, revealing a possible contribution of fine dust particles only to scattering, the result of a possible deposition of the coarse dust fraction during transport. The www.atmos-chem-phys.net/16/455/2016/ Atmos. Chem. Phys., 16, 455–504, 2016 486 M. Mallet et al.: Overview of the ChArMEx/ADRIMED summer 2013 campaign Figure 25. Particle concentrations as a function of size and altitude in the troposphere and lower stratosphere from the LOAC flight under meteorological balloon BLD9 launched from Minorca at the end of a dust event on 19 June 2013, 10:12UT (Table 4; see the daytimeaveraged aerosol optical depth over the sea in Fig. 6). al. (2016) and Nicolas et al. (2016). Figure 26 indicated the average of all instantaneous (clear-sky) DRFs (in Wm−2) estimated during a day for both AERONET/PHOTONS stations. Estimates are reported at the surface (bottom left), at TOA (bottom right) and within the total atmosphere (down). Averaged values of the DRF are also indicated in Fig. 26. As mentioned above, sun-photometer retrievals demonstrate a significant DRF during the SOP-1a experiment. As an example and at the surface, the mean forcing is comprised between −15Wm−2(Barcelona, not affected by dust transport) and −35Wm−2in Burjassot. Such values are consistent with independent 1-D estimates reported by Nicolas et al. (2016) and Meloni et al. (2016). AERONET/PHOTONS data also reveal a negative DRF at TOA over most of sites, meaning that aerosols exert in majority a cooling effect at TOA, with values around ∼ −6 to −12Wm−2. These negative values are also due to the fact that most of the AERONET/PHOTONS stations are located over islands, which are characterized by low surface albedo. Logically and due to the moderate values of aerosol absorption observed during the SOP-1a (Denjean et al., 2015), a positive atmospheric forcing is observed with mean values from +7 to +30Wm−2(with maxima in Burjassot), which could affect the vertical profiles of temperature and relative humidity as shown recently by Nabat et al. (2015a). 5.4.3 Estimates using in situ radiative flux observations As shown by di Sarra et al. (2011), an estimate of the aerosol radiative forcing can be obtained by comparing irradiance measurements made during days characterized by different aerosol loads. In particular, the identification of a cloud-free day with low aerosol amounts is important to provide a reference for pristine conditions. During the SOP-1a, 17 June at Lampedusa displayed a very low aerosol optical depth (daily average of 0.064 at 500nm) and cloud-free conditions throughout the day, and was identified as the reference day forpristineconditions.3 July, conversely, wasoneof the days characterized by the presence of desert dust, with moderate values of the AOD (0.28). As shown in Fig. 22a, dust was present above 2km altitude and there were no major changes in the aerosol vertical distribution during the day, as it also appears from the limited daily variability of the AOD (daily standard deviation of the AOD at 500nm of 0.015). Cloudfree conditions were present throughout the day. Figure 27 displays the downward solar irradiance measured on 3 July, compared with the one measured on the pristine reference day (17 June). The irradiance measurements were corrected for the radiometer thermal offset as discussed by Di Biagio et al. (2009). The sharp narrow peak occurring on 17 June around 06:30 was related to a small isolated cloud, and these data were discarded from the analysis. The differences between the downward irradiances measured on these 2 days were calculated as a function of the solar zenith angle; these differences are due to the effect of aerosol and, to a smaller extent, column water vapor. The effect of water vapor was estimated by means of a radiative transfer model (see, e.g., di Sarra et al., 2011), and the remaining difference was integrated over 24h to obtain the daily average effect, 1I, on the downward solar irradiance. The daily aerosol radiative forcing RF can be derived as RF =1I (1−A), where 1I is the difference between the two curves of Fig. 27 integrated over 24h, and Ais the surface albedo. For a surface albedo of 0.07 (di Sarra et al., 2011), the estimated surface RF is −14.8Wm−2. The radiative forcing efficiency (RFE), which is the radiative forcing produced by a unit AOD, was calculated as RFE =RF/(AOD2−AOD1), where AOD1and AOD1are the measured daily average aerosol optical depth on 17 June and 3 July, respectively. The estimated RFE is −67.4Wm−2. Di Biagio et al. (2010), Atmos. Chem. Phys., 16, 455–504, 2016 www.atmos-chem-phys.net/16/455/2016/ M. Mallet et al.: Overview of the ChArMEx/ADRIMED summer 2013 campaign 487 Figure 26. One-dimensional (clear-sky) instantaneous (shortwave only) DRF calculations (in Wm−2) based on the AERONET/PHOTONS data set for the different stations listed in Table 2 (BOA, TOA and ATM refer to the bottom of the atmosphere, top of atmosphere and atmospheric forcings). Figure 27. Time evolution of the downward solar irradiance observed at Lampedusa on 17 June and 3 July 2013. based on a multi-year data set at Lampedusa, derived a similar value for desert dust (−68.9Wm−2) at the equinox; di Sarra et al. (2011), for an intense desert dust event occurring in March 2010 found values between −70 and −85Wm−2. For a desert dust event associated with the propagation of a gravity wave, with values of AOD similar to those of 3 July, di Sarra et al. (2013) derived an RFE equal to −79Wm−2. Valenzuela et al. (2012) determined RFE for Saharan dust episodes over the western Mediterranean with different origins, showing values in the range from −74Wm−2(for air masses coming from northern Morocco) to −65Wm−2(for air masses coming from Algeria and Tunisia). Values of the dust RFE at the surface in the same range were obtained by Derimian et al. (2006), although they were derived in different conditions for which the influence of surface albedo should be taken into account. The downward LW irradiance measured on 3 July was higher than on 17 June by 23Wm−2. Most of this effect is due to differences in the water vapor column amount (about 1cm difference between the 2 days, with larger values on 3 July). Once the water vapor contribution was subtracted by means of radiative transfer calculations, we found a net positive effect induced by the aerosol of about +5.5Wm−2. This is, on the daily timescale, about 35% of the SW effect. The resulting aerosol RFE in the LW spectral range is +25.5Wm−2, in agreement with previous results by di Sarra et al. (2011), who found values between +25.9 and www.atmos-chem-phys.net/16/455/2016/ Atmos. Chem. Phys., 16, 455–504, 2016 488 M. Mallet et al.: Overview of the ChArMEx/ADRIMED summer 2013 campaign 0 1000 2000 3000 4000 5000 6000 Altitude (m) 0 200 400 600 800 1000 1200 Irradiance E SW (Wm -2 ) E SW Dwn E SW Up E SW Net < θ s >=29.7° < θ s >=29.7° < θ s >=29.7° 0 1000 2000 3000 4000 5000 6000 Altitude (m) -200 -100 0 100 200 300 400 Irradiance E LW (Wm -2 ) E LW Dwn E LW Up E LW Net dT/dt-ATR35 0 1000 2000 3000 4000 5000 6000 Altitude (m) -10 -8 -6 -4 -2 0 2 4 6 8 10 Cooling/heating rate dT/dt (K/day) LW SW Total < θ s > =29.7° (a) (b) (c) Figure 28. SW (a) and LW (b) upward and downward radiative fluxes observed over the Lampedusa station for 22 June and estimated SW and LW heating rate (c) in the two spectral regions (see Sect. 5.4.4 for details). +27.9Wm−2, or Antón et al. (2014), who reported RFE values around +20Wm−2(in reference to AOD at 675nm). 5.4.4 Estimations of the SW and LW radiative heating rate along the vertical One important original aspects of this study concerns the estimates of the vertical profiles of SW and LW radiative heating rate. To our knowledge, all the referenced estimates of this important parameter, which controls for a part the semidirect radiative effect of aerosols, have been conducted using remote-sensing techniques or in situ observations of aerosol optical properties, coupled with radiative transfer modeling. Here, we propose a first estimate of the SW and LW heating rate derived directly from upward and downward (SW and LW) radiative fluxes obtained onboard the ATR-42 aircraft. Because of the nature mainly diffuse of longwave upward and downward irradiances (irradiances in thermal infrared), and of the upward shortwave irradiance (irradiance in solar domain), in first approximation, no correction due to the altitude of the aircraft will be applied to these measurements. Only shortwave downward irradiances will be corrected. Three kinds of corrections are applied: –Correction of the aircraft attitude (unavoidable movements due to the aircraft pitch and roll) –Correction of cosine response of the pyranometer –Correction due to the non-horizontal position of the sensor when a leg stabilizes (i.e., determination of offsets on roll and pitch). Let θmthe angle between the sun direction and the normal to the pyranometer sensor (depending on pitch, roll and aircraft heading given by the inertial navigation system), and θsthe solar zenith angle, the attitude correction coefficient is Xn d=cosθm cosθs. Finally, we obtain the global (direct plus diffuse) downward irradiance, for the solar zenith angle θs:E↓ SW (θs)= Em↓ SW(θm) (Xn d[1−c(θs)]−D)f(θs)+D. In this equation, Em↓ SW (θm)is the measured global irradiance, c(θs)is the cosine response of the pyranometer, D=21 R 0[1−C(θm)]µdµ and f(θs)is the part of direct downward irradiance in the global (estimation obtained from radiative transfer code). Taking into account these corrections, Fig. 28a shows downward EDwn SW , upward EUp SW, and net ENet SWshortwave irradiances obtained from measurements performed onboard ATR-42 aircraft on 22 June between 10.35 and 11.30UT. Irradiances are reduced to the mean solar zenith angle θs=29.7◦. Similarly, Fig. 28b shows corresponding measurements of downward EDwn LW , upward EUp LW, and net ENet LWlongwave irradiances. Total Atmos. Chem. Phys., 16, 455–504, 2016 www.atmos-chem-phys.net/16/455/2016/ M. Mallet et al.: Overview of the ChArMEx/ADRIMED summer 2013 campaign 489 net irradiances are then determined versus the aircraft altitude for the mean air mass factor of the considered studied flight phase. The radiative cooling/heating rate is finally derived and shown in Fig. 28c, in which the longwave (LW) and shortwave (SW) parts are distinguished. Concerning the SW heating rate vertical profiles (Fig. 28c), one can observe the significant increase of the calculated instantaneous SW heating rate in the two different aerosol layers detected for this case (Fig. 21), especially above 4km, which corresponds to the maximum of extinction coefficient (up to 100Mm−1) due to the presence of mineral dust. For this specific layer, the values of SW heating rate peak at 4–5◦K per day for a solar angle of 29.7◦. We can also observe a similar tendency in the second aerosol layer, located between 1.5 and 3km (see Fig. 21). Concerning the LW heating rate, Fig. 28c indicates instantaneous values ranging between −2 and −4◦K per day, which is also consistent with the well-known cooling effect of mineral dust in the longwave spectrum (Mallet et al., 2006; Zhu et al., 2007). As shown in Fig. 28c, the net heating rate is dominated by the SW heating (the maximum LW cooling is less than 60% of the SW heating), which leads to net SW radiative heating ranging between +0.5 and +2K per day inside the dust layer above the MBL. Such unique and original database of SW and LW radiative heating obtained over the western Mediterranean should be now used to evaluate the ability of the different models involved in the ChArMEx/ADRIMED project (see the following Sect. 6) to simulate this important radiative property for the different identified dust cases. 6 Overview of modeling activities Several models are used to analyze the SOP-1a period: the meso-scale meteorological COSMO-MUSCAT model, the chemistry transport model (CTM) CHIMERE model, and two regional climate (RegCM and CNRM-RCSM) models. These models differ in terms of horizontal and vertical resolutions, physical parameterizations, aerosolchemical schemes and are able to deliver complementary information to address key scientific questions of the ChArMEx/ADRIMED experiment. Their main characteristics are summarized in Table 7. 6.1 COSMO-MUSCAT model The COSMO-MUSCAT (Wolke et al., 2012) parallelized multi-scale regional model system consists of the COSMO (Consortium for Small-scale Modelling) non-hydrostatic atmosphere model that is on-line coupled to the MUSCAT (MUltiScale Chemistry Aerosol Transport Model) 3-D chemistry tracer transport model. The atmospheric dust cycle consisting of the emission, transport and deposition of dust particles is simulated within MUSCAT using meteorological and hydrological fields from COSMO. Dust emission is calculated using the emission scheme by Tegen et al. (2002) and depends on local surface wind friction velocities, surface roughness length, soil texture and soil moisture. Calculated dust emission fluxes depend on particle diameter for individual size classes that are assumed to be log-normally distributed. Following Marticorena and Bergametti (1995), dust emission is considered as threshold function of local friction velocities and thus initial dust emission is computed as a function of soil particle size distribution. Dust emission is limited to regions where active dust sources have been identified during 2006–2009 from MSG SEVIRI observations (Schepanski et al., 2007). The advection of dust particles is described by a third-order upstream scheme; dust particles are transported as a passive tracer in five independent size classes with limiting radius at 0.1, 0.3, 0.9, 2.6, 8, and 24µm. The removal of dust particles from the atmosphere is described by dry and wet deposition taking particle size, particle density, and atmospheric conditions into account. Here, the simulations of the atmospheric dust cycle are performed at a 28km horizontal grid and 40 vertical layers covering northern African dust sources, the eastern North Atlantic, the Mediterranean basin and Europe. 6.2 The CHIMERE chemistry transport model CHIMERE is a chemistry transport model able to simulate concentrations fields of gaseous and aerosol species at a regional scale. The model is off-line and thus needs precalculated meteorological fields to run. In this study, we used the version fully described in Menut et al. (2013), forced by the WRF meso-scale model. The horizontal domain is the same as the one of WRF, and, for the vertical grid, the 28 vertical levels of WRF are projected on the 20 levels of the CHIMERE mesh. The gaseous species are calculated using the MELCHIOR 2 scheme and the aerosols using the scheme developed by Bessagnet et al. (2004). This module takes into account species such as sulfate, nitrate, ammonium, primary organic (OC) and black carbon (BC), secondary organic aerosols (SOA), sea spray, mineral dust, and water. These aerosols are represented using ten bins, from 40nm to 20µm, in diameter. The life cycle of these aerosols is completely represented with nucleation of sulfuric acid, coagulation, adsorption/desorption, wet and dry deposition and scavenging. This scavenging is both represented by coagulation with cloud droplets and precipitation. The formation of SOA is also taken into account. The anthropogenic emissions are estimated using the same methodology as the one described in Menut et al. (2013) but with the HTAP masses as input data. These masses were prepared by the EDGAR Team, using inventories based on MICS-Asia, EPA-US/Canada and TNO databases (http://edgar.jrc.ec.europa.eu/htap_v2). Biogenic emissions are calculated using the MEGAN emissions scheme (Guenther et al., 2006), which provides fluxes of isoprene, terpene and pinenes. In addition to this 2013 version, several processes were improved and added in the framework www.atmos-chem-phys.net/16/455/2016/ Atmos. Chem. Phys., 16, 455–504, 2016 490 M. Mallet et al.: Overview of the ChArMEx/ADRIMED summer 2013 campaign of this study. First, mineral dust emissions are now calculated using new soil and surface databases, as described in Menut et al. (2013). Second, chemical species emissions fluxes produced by vegetation fires are estimated using the new highresolution fire model presented in Turquety et al. (2014). Finally, the photolysis rates are explicitly calculated using the FastJ radiation module (Mailler et al., 2015). 6.3 The RegCM regional climate model The RegCM system is a community model designed for use by a varied community composed of scientists in industrialized countries as well as developing nations. It is supported through the Regional Climate Network, or RegCNET, a widespread network of scientists coordinated by the Earth System Physics section of the Abdus Salam International Centre for the Theoretical Physics (ICTP, Giorgi et al., 2012). RegCM is a hydrostatic, compressible, sigma-p vertical coordinate model. As a limited area model, RegCM requires initial and boundary conditions that can be provided both by NCEP or ECMWF analyses. The horizontal resolution used need to be higher than 10km, due to the hydrostatic dynamic core of the model, associated with 23 vertical levels. A simplified aerosol scheme specifically designed for application to long-term climate simulations has been incrementally developed within the RegCM system. Solmon et al. (2006, 2008) first implemented a first-generation aerosol model including sulfates, organic carbon, and black carbon. Zakey et al. (2006) then added a four-bin desert dust module, and Zakey et al. (2008) implemented a two-bin sea-spray scheme. In RegCM, the dust emission scheme accounts for sub-grid emissions by different types of soil. The dust emission size distribution can now also be treated according to Kok (2011). When all aerosols are simulated, 12 additional prognostic equations are solved in RegCM, including transport by resolvable scale winds, turbulence and deep convection, sources, and wet and dry removal processes. In RegCM, the natural/anthropogenic aerosols are radiatively interactive both in the solar and infrared regions and so are able to feedback on the meteorological fields. 6.4 The CNRM-RCSM regional climate model The fully coupled RCSM (Regional Climate System Model), which is developed at CNRM, has also been used within the ChArMEx/ADRIMED project. This model includes the ALADIN-Climate (Déqué and Somot, 2008) regional climate atmospheric model, the NEMOMED8 (Beuvier et al., 2010) regional ocean model and the ISBA (Noilhan and Mahfouf, 1996) land-surface model. We used here the version described in Nabat et al. (2015b) with a 50km horizontal resolution. ALADIN-Climate includes the Fouquart and Morcrette radiation scheme based on the ECMWF model incorporating effects of greenhouse gases as well as direct effects of aerosols. The NEMOMED8 ocean model is the regional eddy-permitting version of the NEMOV2.3 ocean model that covers the Mediterranean Sea. Concerning the aerosol phase, the ALADIN-Climate model incorporates a radiative scheme to take into account the direct and semidirect effects of five aerosol types (sea spray, desert dust, sulfates, black and organic carbon aerosols) through either AOD climatologies or a prognostic aerosol scheme (Nabat et al., 2013, 2015b). On the one hand, Nabat et al. (2013) have proposed a new AOD monthly climatology over the period 2003–2009, based on a combination of satellite-derived and model-simulated products. The objective is to have the best estimation of the atmospheric aerosol content for these five most relevant aerosol species. On the other hand, a prognostic aerosol scheme has been recently implemented in ALADIN-Climate, and has shown its ability to reproduce the main patterns of the aerosol variability over the Mediterranean (Nabat et al., 2015b). Using CNRM-RCSM with the new AOD monthly climatology over the period 2003–2009 (Nabat et al., 2013), Nabat et al. (2015a) have notably highlighted the response of the Mediterranean sea surface temperature (SST) to the aerosol direct and semi-direct radiative forcing. Figure 29a presents the annual average difference in SST over the period 2003– 2009 between a simulation ensemble including aerosols and a second one without any aerosol. Aerosols are found to induce an average decrease in SST by 0.5◦C, because of the scattering and absorption of incident radiation. As a consequence, the latent heat loss is also reduced by aerosols (Fig. 29b), as well as precipitation (Fig. 29c). This result also underlines the importance of taking into account the ocean– atmosphere coupling in regional aerosol-climate studies over the Mediterranean. 6.5 SOP-1a multi-model aerosol simulations 6.5.1 Aerosol optical depth Figure 30 reports the AOD (in the visible range) simulated for the SOP-1a period and for the COSMO-M (550nm), RegCM (between 440 and 670nm), CNRM-RCSM (550nm) and CHIMERE (500nm) models. Except for the CTMCHIMERE model, which includes all the secondary species (SOA and inorganic), the others have different aerosol schemes and take into account both natural (COSMO-M) or natural plus a part of anthropogenic aerosols as described in Table 7. The configurations used for each models are listed in Table 7. One can observe the large variability of AOD simulated by models over the Mediterranean region, with highest values clearly simulated by the COSMO-M (AOD∼1– 1.5 in the visible wavelengths) over the northern African region. The CHIMERE model indicates two different regions where AOD peaks around 1, over Algeria–Tunisia and southern Morocco. For COSMO-M and CHIMERE, no intense dust AOD is simulated over northeastern Africa (Libya and Egypt), and values are below 0.25, contrary to RegCM and Atmos. Chem. Phys., 16, 455–504, 2016 www.atmos-chem-phys.net/16/455/2016/ M. Mallet et al.: Overview of the ChArMEx/ADRIMED summer 2013 campaign 491 (a) (b) (c) Figure 29. Annual average difference in (a) sea surface temperature (SST), latent heat loss (b) and precipitation (c) over the period 2003– 2009 between a simulation ensemble including aerosols and a second one without any aerosol. Figure 30. AOD averaged for the 15 to 25 June 2013 period from the models RegCM, COSMO-MUSCAT, CHIMERE and CNRM-RCSM. Details about the model configurations are provided in Table 7. CNRM-RCSM that simulate moderate AOD over this region, with more intense peaks (∼0.7 for CNRM-RCSM simulations). Some identified regions with important AOD over Tunisia, Algeria, and southern Morocco are well captured by all models, except COSMO-M, which shows more intense AOD south of Algeria. It should be noted that this regional pattern of AOD is found to be consistent with MODIS observations as shown by Menut et al. (2015) for the CHIMERE model. Averaged over the SOP-1a period, all models simulate low to moderate AOD over the EURO-Mediterranean region, which is consistent with AERONET/PHOTONS observations (Fig. 14). Once again and as noted by Menut et al. (2015), this modeling exercise clearly shows that the summer of 2013 was not characterized by intense dust plumes or intense anthropogenic or forest fire emissions. However, modeling results indicate regular dust intrusions during the SOP-1a characterized by moderate atmospheric loads. Over Europe, the CHIMERE CTM obviously simulates anthropogenic aerosol AOD (AOD∼0.3), especially over Benelux and the Po Valley that are not simulated by www.atmos-chem-phys.net/16/455/2016/ Atmos. Chem. Phys., 16, 455–504, 2016 492 M. Mallet et al.: Overview of the ChArMEx/ADRIMED summer 2013 campaign Figure 31. Time series of AOD comparisons between AERONET/PHOTONS observations and COSMO-MUSCAT, CHIMERE, CNRMRCSM and RegCM model outputs over the two stations of Ersa and Lampedusa. The MODIS retrievals are also indicated. the two other regional models. Indeed, CNRM-RCSM simulations reveal a more diffuse AOD about 0.2 over Europe with maximum over Western France certainly due to the advection of primary marine particles generated over the Atlantic Ocean. RegCM simulations indicate a plume of anthropogenic aerosols over the Balkan region mainly due to secondary inorganic species. As RegCM does not use the spectral nudging technique in this simulation and are only forced at the boundaries during the period of simulation, some biases in meteorological fields could appear (as for the precipitation location and intensity), which need to be evaluated. Finally and in addition to analysis of the AOD regional pattern, a specific comparison with in situ observations and remotesensing (AERONET/PHOTONS and satellite) data has been made for the CTM-CHIMERE model (Menut et al., 2015) and is planned in accompanied studies for the COSMO-M, RegCM and CNRM-RCSM models, associated with an intercomparison exercise for evaluating the dust emissions, vertical distribution, size distribution and dry/wet deposition using all data collected in the framework of the SOP-1a. In parallel to time-averaged AOD simulated at the regional scale, we report comparisons of simulated AOD with AERONET/PHOTONS data for the two reference stations (Lampedusa and Ersa). As reported in Table 7, it should be reminded here that all models did not take into account aerosol species in a similar way. As an example, COSMOMUSCAT includes mineral dust only in this simulation, while CNRM-RCSM and the RegCM model include natural (sea spray and dust) and sulfates as well as secondary ammonium and nitrate particles (treated as bulk aerosols), but for RegCM only. The most complete regional model is the CTM-CHIMERE, which takes into account natural and all anthropogenic particles (including secondary organics and inorganic) resolved in size by using a number of bins (Menut et al., 2013) higher than used in RegCM, CNRM-RCSM or COSMO-MUSCAT (number of dust bins between 3 and 4 bins) models. Figure 31 reports the time evolution of simulated and observed AOD at 550nm for the two sites (Ersa and Lampedusa) during the SOP-1a. Time correlation, as well as bias, is calculated after removing AERONET/PHOTONS data for 27 June, strongly affected by smoke aerosols transported from North America biomass burning sources that are not included in the different domains. Figure 31 indicates that all models are able to simulate AOD in the range of magnitude of observations. For the dusty Lampedusa site, CNRM-RCSM and CHIMERE reveal high temporal correlations (0.82 and 0.85, respectively), with standard deviations close to AERONET/PHOTONS data, especially for CHIMERE. For this station, COSMO-M and RegCM display moderate temporal correlation (0.55 and 0.49, respectively) compared to CNRM-RCSM and CHIMERE. As already mentioned, one reason for the lowest time correlation for these models is related to the fact that they are only forced at the boundaries, and the synoptic conditions inside the domain can be derived during the simulation. This effect is limited for CNRM-RCSM that used the spectral nudging technique and for CHIMERE forced by WRF meteorological field (Menut et al., 2015). For each models, biases are shown to be low, both positive (for CNRM-RCSM and CHIMERE) and negative (for COSMO-M and RegCM). For the Ersa station, less influenced by long-range transport of mineral dust during this period, temporal correlations are lowest and found to be moderate (0.40) for CHIMERE and COSMO-M and low for RegCM and CNRM-RCSM. In terms of bias, values are positive and low (0.02 to 0.04) for all models, except for COSMO-M (−0.07), which does not include anthropogenic aerosols or sea spray in the present simulation (Table 7). For each model, calculated standard deviations are in the same range of magnitude but slightly higher than observations, especially for RegCM (bias of 0.08) that simulated a large AOD for the 19–20 June period. By comparison with the values obtained in Lampedusa, these low correlations at Ersa reveal the limitations of these models in terms of horizontal resolution with respect to the Atmos. Chem. Phys., 16, 455–504, 2016 www.atmos-chem-phys.net/16/455/2016/ M. Mallet et al.: Overview of the ChArMEx/ADRIMED summer 2013 campaign 493 Figure 32. Averaged surface and TOA SW DRF simulated in clear-sky conditions and over the SOP-1a period by the CNRM-RCSM and RegCM models. representativeness of the site. Lampedusa being isolated in the middle of the Mediterranean and under the main pathways of African mineral dust, AOD is mostly related to longrange transport. On the other hand, the site of Ersa in Corsica may under several types of aerosol contributions (anthropogenic, biogenic) be more intense and more spatially variable than in Lampedusa. Ersa being closer to large industrial areas, the models with a horizontal resolution of tens of kilometers are probably not highly enough resolved to catch small-scale aerosol plumes from the continent. 6.5.2 Regional SW 3-D direct radiative forcing The SW (clear-sky) DRF, averaged for the SOP-1a period, has been estimated from the RegCM and CNRM-RCSM models, both at the surface and TOA, as shown in Fig. 32. For this discussion, we only consider these two models as they estimate the clear-sky SW DRF by taking into account natural and anthropogenic aerosols, contrary to the COSMOMUSCAT model in this study. At the surface first, one can observe the large regional dimming due to anthropogenic (especially over Europe) and natural (northern Africa and Mediterranean) particles over the Euro-Mediterranean. Concerning the northern African region, both models simulate large surface forcing ∼ −20Wm−2(with local maxima of −50Wm−2associated with higher AOD). CNRM-RCSM is shown to simulate higher surface radiative forcing for the whole domain, especially over Algeria. Although such RCM climate models are not designed to simulate finely the size distribution and the chemical composition of aerosols as an A-Q (Air-Quality) system (Menut et al., 2013), a first estimate of the radiative effect of polluted particles over Europe is provided. Figure 32 displays a negative forcing, obviously lower than for mineral dust, of about −10 to −15Wm−2for RegCM, especially over the Balkans, and no significant radiative effect over the Benelux region for this period. Over the continental region, CNRM-RCSM simulated a more diffuse surface forcing with values around −10Wm−2, including a large part of Europe (France, Benelux and eastern Europe). As shown recently by Nabat et al. (2015a), this decrease in SW radiations due to aerosols could perturb the surface continental temperature, SST and latent heat fluxes over the Mediterranean Sea and more largely on meteorological fields. At TOA, the dipole of the direct forcing between the north and the south of the domain is well reproduced by the two RCM systems, with more intense values for CNRMRCSM. One can clearly observe positive forcing at TOA (heating) over northern Africa and negative forcing (cooling) over the Mediterranean and Europe. This represents one of www.atmos-chem-phys.net/16/455/2016/ Atmos. Chem. Phys., 16, 455–504, 2016 494 M. Mallet et al.: Overview of the ChArMEx/ADRIMED summer 2013 campaign Table 7. Main characteristics (period of simulations, horizontal resolution, number of vertical layers, main aerosol (primary and/or secondary) species, radiative transfer codes) of the four different 3-D models used during the SOP-1a experiment (see Sect. 6) (GME is for the global model of the German Weather Service). Models Time of Horizontal Number of Aerosol species Boundary layer Radiative transfer code simulation resolution vertical layers forcing CHIMERE 1 Jun–31 Jul 50km 20 Dust, sea salt, secondary organic and inorganic, primary OC–BC WRF FastJX CNRM-RCSM 1 Jun–31 Jul 50km 31 Dust, sea salt, sulfates, primary OC–BC ERA-Interim SW: FMR (six bands, Morcrette et al., 1989) LW: RRTM (Mlawer et al., 1997) RegCM 13 Jun–5 Jul 25km 23 Dust, sea salt, secondary inorganic, primary OC–BC NCEP reanalysis CCM3 or RRTM COSMO-MUSCAT 15 May–31 Jul 28km 40 Dust GME Ritter and Geleyn (1992) the characteristics of the Euro-Mediterranean region, with a large variability of surface albedo from the south (with higher values) to the north (low to moderate albedo). Due to this gradient in the surface albedo, moderate absorbing dust aerosols emitted over northern Africa (characterized by high surface albedo) decrease the shortwave radiations reflected at TOA, compared to a non-turbid atmosphere. When advected above low surface reflectance as marine or dense forest over Europe, dust aerosols increase the upward SW radiations at TOA, leading to a cooling effect. One can see the transition from positive to negative TOA forcing that occurs over northern Algeria and Morocco as soon as dust particles are transported over darker surfaces. This TOA radiative forcing gradient is well captured by such RCM models that use a finer resolution than GCM. Over Europe and Mediterranean, the TOA forcing is simulated to be negative for both RCM with lower values around −5 to −10Wm−2. Such results are consistent with the study of Nicolas et al. (2016), who performed two different simulations using different surface albedo (from marine to continental), based on the ATR-42 observations above the Balearic Islands and the Granada station. The inclusion of high surface albedo (0.27 at 870nm) in the 1-D radiative transfer model compared to low sea surface albedo (0.02 at 870nm) contributes to decrease the TOA radiative effect at Granada. The last important point to mention here concerns the fact that most SW radiation losses at the surface are not completely compensated for by fluxes reflected back to space. Hence, this gain in solar energy within dusty layers (due to moderate dust SW absorption; see Denjean et al., 2015) has been shown to result in significant feedbacks on the temperature and relative humidity profiles over the Mediterranean region, with some important implications for its climate (Nabat et al., 2015a). 7 Conclusions The special observing period (SOP-1a) performed during the Mediterranean dry season (11 June to 5 July 2013) over the western and central Mediterranean basins has been described in detail, as well as the 1-D to 3-D modeling effort, involved in the ChArMEx/ADRIMED project focused on aerosol-radiation-climate interactions. Details of the in situ and remote-sensing instrumentation deployed at the different sites and the main meteorological conditions that occurred during the campaign have been provided. Some results from the in situ and remote-sensing observations, vertical profiles, 1-D and 3-D aerosol direct radiative forcing (DRF) computations have also been presented. Concerning the aerosol loading during the SOP-1a, our results indicate that numerous but moderate mineral dust plumes were observed during the campaign with main sources located in Morocco, Algeria and Tunisia, leading to AOD between 0.1 and 0.6 (at 440nm) over the western and central Atmos. Chem. Phys., 16, 455–504, 2016 www.atmos-chem-phys.net/16/455/2016/ M. Mallet et al.: Overview of the ChArMEx/ADRIMED summer 2013 campaign 495 Mediterranean. Analysis of synoptic situations demonstrates unfavorable conditions to produce large concentrations of polluted smoke particles during the SOP-1a, but interesting sea-spray events have been observed. Aerosol extinctions measured onboard the ATR-42 show local maxima reaching up to 150Mm−1within the dust plume, associated with extinctions of about 50Mm−1within the marine boundary layer (MBL), possibly due to the presence of sea-spray aerosols. By combining ATR-42 extinction, absorption and scattering measurements, complete optical closures have been made, revealing an excellent agreement in estimated optical properties. This additional information on extinction properties has allowed calculation of the dust single scattering albedo (SSA) with a high level of confidence over the western Mediterranean. Our results show a surprising moderate variability from 0.90 to 1.00 (at 530nm) for all flights studied, corroborated by AERONET/PHOTONS SSA retrievals. The SSA derived during the ChArMEx/ADRIMED project has also been compared with referenced values obtained near dust sources, showing a relatively low difference in this optical parameter at 530nm. Concerning the aerosol vertical structure, active remotesensing observations, at the surface and onboard the F-20, indicate complex vertical profiles of particles with sea spray and pollution located in the MBL, and mineral dust and/or even aged North American smoke particles located above (up to 6–7km in altitude). Microphysical properties of aerosols measured onboard the ATR-42 and balloon-borne observations for transported/aged mineral dust reveal particle volume size distributions with diameters greater than 10µm. In most cases, a coarse mode of mineral dust particles, characterized by an effective diameter Deff,cranging between 5 and 10µm, has been detected within the dust layer located above the MBL. Such values are found to be larger than those referenced in dust source regions during FENNEC, SAMUM1 and AMMA, as well as measurements in the Atlantic Ocean in the Cape Verde region during SAMUM-2 and at Puerto Rico during PRIDE. In terms of shortwave (SW) and longwave (LW) DRF, in situ surface and aircraft observations have been merged and used as inputs in different radiative transfer codes for calculating the 1-D DRF. Modeling results show significant surface (instantaneous) SW radiative forcing down to as much as −90Wm−2over super-sites. In parallel, AOD together with surface radiative flux observations have also been used to directly estimate the local daily surface forcing in SW (and LW) spectral regions, showing a significant effect with values of −15Wm−2(+5.5Wm−2) over Lampedusa. Such DRFs are consistent with those previously referenced over the Mediterranean basin. In parallel, aircraft observations also provide original and new estimates of SW and LW radiative heating vertical profiles with significant values of SW heating of about 5◦K per day within the dust layer (for a solar angle of 30◦). Associated 3-D modeling studies, using regional climate (RCM) and chemistry transport (CTM) models, indicate a relatively good agreement between simulated AOD and that determined from AERONET/PHOTONS data. Such models allow 3-D calculations of the daily SW DRF revealing a regional DRF of −10 to −20Wm−2(at the surface and in clear-sky conditions), when averaged over the SOP-1a period. At TOA, a significant dipole in the DRF is estimated between the north and the south of the domain, with positive (heating) over northern Africa and negative (cooling) DRF over the Mediterranean basin and Europe, reflecting changes in surface albedo associated with moderately absorbing aerosols. A first multi-year simulation (conducted for the 2003 to 2009 period) that takes into account the ocean– atmosphere coupling has demonstrated that the significant aerosol radiative forcing is responsible for a decrease in sea surface temperature (on average −0.5◦C for the Mediterranean). In addition, the latent heat loss is shown to be weaker in the presence of aerosols, resulting in a decrease in specific humidity in the lower troposphere, and a reduction in cloud cover and precipitation. This unprecedented data set of aerosol microphysical, chemical, optical properties and vertical profiles obtained over the western Mediterranean will now be used for evaluating regional models to reproduce such properties. In addition to classical model evaluations based generally on the AOD, new comparisons between models and in situ observations on aerosol absorbing (SSA and AAOD) properties and SW and LW heating rates, which control the semi-direct effect of aerosols, should be conducted. Comparisons will also be performed on the aerosol size distribution for investigating the ability of regional models to simulate the observed large dust particle size during the transport over the Mediterranean, which could be helpful for improving the representation of deposition in such models. In parallel, in situ observations of sea-spray particles obtained at the surface and from ATR-42 measurements will also be used to evaluate the different primary sea-spray generation schemes, in terms of concentration and size distribution. The objective is to improve the representation of microphysical and optical properties of aerosols in regional climate models that will be used in multi-year simulations to assess the impact of natural and anthropogenic aerosols on the climate in this region. Acknowledgements. This research has received funding from the French National Research Agency (ANR) project ADRIMED (contract ANR-11-BS56-0006). This work is part of the ChArMEx project supported by ADEME, CEA, CNRS-INSU and MétéoFrance through the multidisciplinary programme MISTRALS (Mediterranean Integrated Studies aT Regional And Local Scales). The station at Ersa was partly supported by the CORSiCA project funded by the Collectivité Territoriale de Corse through the Fonds Européen de Développement Régional of the European Operational Program 2007–2013 and the Contrat de Plan Etat-Région. We acknowledge the AERONET/PHOTONS sun-photometer networks www.atmos-chem-phys.net/16/455/2016/ Atmos. Chem. Phys., 16, 455–504, 2016 502 M. Mallet et al.: Overview of the ChArMEx/ADRIMED summer 2013 campaign Pérez, C., Nickovic, S., Baldasano, J. M., Sicard, M., Rocadenbosch, F., and Cachorro, V. E.: A long Saharan dust event over the western Mediterranean: Lidar, sun photometer observations, and regional dust modeling. J, Geophys. Res., 111, D15214, doi:10.1029/2005JD006579, 2006. Petzold, A., Onasch, T., Kebabian, P., and Freedman, A.: Intercomparisonof a CavityAttenuatedPhase Shift-based extinctionmonitor (CAPS PMex) with an integrating nephelometer and a filterbased absorption monitor, Atmos. Meas. Tech., 6, 1141–1151, doi:10.5194/amt-6-1141-2013, 2013. Piazzola, J., Tedeschi, G., and Demoisson, A.: A model for a transport of sea-spray Aerosols in the coastal zone, Bound. Lay. Meteorol., 155, 329–350, 2015. Ramanathan, V., Crutzen, P. J., Lelieveld, J., Mitra, A. P., Althausen, D., Anderson, J., Andreae, M. O., Cantrell, W., Cass, G., Chung, C. E., Clarke, A. D., Coakley, J. A., Collins, W. D., Conant, W. C., Dulac, F., Heintzenberg, J., Heymsfield, A. J., Holben, B., Howell, S., Hudson, J., Jayaraman, A., Kiehl, J. T., Krishnamurti, T. N., Lubin, D., MacFarquhar, G., Novakov, T., Ogren, J. A., Podgorny, I. A., Prather, K., Priestley, K., Prospero, J. M., Quinn, P. K., Rajeev, K., Rasch, P., Rupert, S., Sadourny, R., Satheesh, S. K., Shaw, G. E., Sheridan, P., and Valero, F. P. J.: Indian Ocean experiment: An integrated analysis of the climate forcing and effects of the great Indo-Asian haze, J. Geophys. Res., 106, 28371– 28398, doi:10.1029/2001JD900133, 2001. Ravetta, F., Ancellet, G., Colette, A., and Schlager, H.: Longrange transport and tropospheric ozone variability in the western Mediterranean region during the Intercontinental Transport of Ozone and Precursors (ITOP-2004) campaign, J. Geophys. Res., 112, 1–12, doi:10.1029/2006JD007724, 2007. Rea, G., Turquety, S., Menut, L., Briant, R., Mailler, S., and Siour, G.: Source contributions to 2012 summertime aerosols in the Euro-Mediterranean region, Atmos. Chem. Phys. Discuss., 15, 8191–8242, doi:10.5194/acpd-15-8191-2015, 2015. Renard, J.-B., Dulac, F., Berthet, G., Lurton, T., Vignelles, D., Jégou, F., Tonnelier, T., Thaury, C., Jeannot, M., Couté, B., Akiki, R., Mineau, J.-L., Verdier, N., Mallet, M., Gensdarmes, F., Charpentier, P., Mesmin, S., Duverger, V., Dupont, J.-C., Elias, T., Crenn, V., Sciare, J., Giacomoni, J., Gobbi, M., Hamonou, E., Olafsson, H., Dagsson-Waldhauserova, P., Camy-Peyret, C., Mazel, C., Décamps, T., Piringer, M., Surcin, J., and Daugeron, D.: LOAC: a small aerosol optical counter/sizer for ground-based and balloon measurements of the size distribution and nature of atmospheric particles – Part 1: Principle of measurements and instrument evaluation, Atmos. Meas. Tech. Discuss., 8, 1203– 1259, doi:10.5194/amtd-8-1203-2015, 2015a. Renard, J.-B., Dulac, F., Berthet, G., Lurton, T., Vignelle, D., Jégou, F., Tonnelier, T., Thaury, C., Jeannot, M., Couté, B., Akiki, R., Mineau, J.-L., Verdier, N., Mallet, M., Gensdarmes, F., Charpentier, P., Mesmin, S., Duverger, V., Dupont, J.-C., Elias, T., Crenn, V., Sciare, J., Giacomoni, J., Gobbi, M., Hamonou, E., Olafsson, H., Dagsson-Waldhauserova, P., Camy-Peyret, C., Mazel, C., Décamps, T., Piringer, M., Surcin, J., and Daugeron, D.: LOAC: a small aerosol optical counter/sizer for ground-based and balloon measurements of the size distribution and nature of atmospheric particles – Part 2: First results from balloon and unmanned aerial vehicle flights, Atmos. Meas. Tech. Discuss., 8, 1261–1299, doi:10.5194/amtd-8-1261-2015, 2015b. Ripoll, A., Minguillón, M. C., Pey, J., Pérez, N., Querol, X., and Alastuey, A.: Joint analysis of continental and regional background environments in the western Mediterranean: PM1and PM10 concentrations and composition, Atmos. Chem. Phys., 15, 1129–1145, doi:10.5194/acp-15-1129-2015, 2015. Roger, J. C., Mallet, M., Dubuisson, P., Cachier, H., Vermote, E., Dubovik, O., and Despiau, S.: A synergetic approach for estimating the local direct aerosol forcing: application to an urban zone during the Experience sur Site pour Contraindre les Modèles de Pollution et de Transport d’Emission (ESCOMPTE) experiment, J. Geophys. Res., 111, D13208, doi:10.1029/2005JD006361, 2006. Royer, P., Raut, J.-C., Ajello, G., Berthier, S., and Chazette, P.: Synergy between CALIOP and MODIS instruments for aerosol monitoring: application to the Po Valley, Atmos. Meas. Tech., 3, 893– 907, doi:10.5194/amt-3-893-2010, 2010. Ryder, C. L., Highwood, E. J., Rosenberg, P. D., Trembath, J., Brooke, J. K., Bart, M., Dean, A., Crosier, J., Dorsey, J., Brindley,H.,Banks,J., Marsham, J. H.,McQuaid, J.B., Sodemann, H., and Washington, R.: Optical properties of Saharan dust aerosol and contribution from the coarse mode as measured during the Fennec 2011 aircraft campaign, Atmos. Chem. Phys., 13, 303– 325, doi:10.5194/acp-13-303-2013, 2013. Saha, A., Mallet, M., Roger, J. C., Dubuisson, P., Piazzola, J., and Despiau, S.: One year measurements of aerosol optical properties over an urban coastal site: Effect on local direct radiative forcing, Atmos. Res., 90, 195–202, doi:10.1016/j.atmosres.2008.02.003, 2008. Sanchez-Gomez, E., Somot, S., and Mariotti, A.: Future changes in the Mediterranean water budget projected by an ensemble of regional climate models, Geophys. Res. Lett., 36, L21401, doi:10.1029/2009GL040120, 2009. Salameh, T., Drobinski, P., Menut, L., Bessagnet, B., Flamant, C., Hodzic, A., and Vautard, R.: Aerosol distribution over the western Mediterranean basin during a Tramontane/Mistral event, Ann. Geophys., 25, 2271–2291, doi:10.5194/angeo-25-22712007, 2007. Santese, M., Perrone, M. R., Zakey, A. S., De Tomasi, F., and Giorgi, F.: Modeling of Saharan dust outbreaks over the Mediterranean by RegCM3: case studies, Atmos. Chem. Phys., 10, 133– 156, doi:10.5194/acp-10-133-2010, 2010. Sassen, K.: Lidar backscatter depolarization technique for cloud and aerosol research, in: Light Scattering by Nonspherical Particles: Theory, Measurements, and Applications, edited by: Mishchenko, M., Hovenier, J. W., and Travis, L. D., Academic Press, San Diego, 393–417, 1999. Schepanski, K., Tegen, I., Laurent, B., Heinold, B., and Macke, A.: A new Saharan dust source activation frequency map derived from MSG-SEVIRI IR-channels, Geophys. Res. Lett., 34, 18803, doi:10.1029/2007GL030168, 2007. Sciare, J., Cachier, H., Oikonomou, K., Ausset, P., Sarda-Estève, R., and Mihalopoulos, N.: Characterization of carbonaceous aerosols during the MINOS campaign in Crete, July–August 2001: a multi-analytical approach, Atmos. Chem. Phys., 3, 1743– 1757, doi:10.5194/acp-3-1743-2003, 2003. Sciare, J., Oikonomou, K., Favez, O., Liakakou, E., Markaki, Z., Cachier, H., and Mihalopoulos, N.: Long-term measurements of carbonaceous aerosols in the Eastern Mediterranean: evidence of Atmos. Chem. Phys., 16, 455–504, 2016 www.atmos-chem-phys.net/16/455/2016/ M. Mallet et al.: Overview of the ChArMEx/ADRIMED summer 2013 campaign 503 long-range transport of biomass burning, Atmos. Chem. Phys., 8, 5551–5563, doi:10.5194/acp-8-5551-2008, 2008. Schicker, I., Radanovics, S., and Seibert, P.: Origin and transport of Mediterranean moisture and air, Atmos. Chem. Phys., 10, 5089– 5105, doi:10.5194/acp-10-5089-2010, 2010. Schroeder, W., Csiszar, I., Giglio, L., and Schmidt, C. C.: On the use of fire radiative power, area, and temperature estimates to characterize biomass burning via moderate to coarse spatial resolution remote sensing data in the Brazilian Amazon, J. Geophys. Res., 115, D21121, doi:10.1029/2009JD013769, 2010. Sellegri, K., Rose, C., Culot, A., Sauvage, S., Roberts, G., Marchand, N., Pey, J., Sciare, J., Bourriane, T., Mallet, M., and Dulac, F.: Spatial extent, occurrence and precursors of nucleation events over the western Mediterranean basin, in preparation, 2016. Sellitto, P., di Sarra, A., Corradini, S., Boichu, M., Herbin, H., Dubuisson, P., Sèze, G., Meloni, D., Monteleone, F., Merucci, L., Rusalem, J., Salerno, G., Briole, P., and Legras, B.: Synergistic use of Lagrangian dispersion modelling, satellite and surface remote sensing measurements for the investigation of volcanic plumes: the Mount Etna eruption of 25–27 October 2013, Atmos. Chem. Phys. Discuss., 15, 31335–31383, doi:10.5194/acpd-1531335-2015, 2015. Sicard, M., Rocadenbosch, F., Reba, M. N. M., Comerón, A., Tomás, S., García-Vízcaino, D., Batet, O., Barrios, R., Kumar, D., and Baldasano, J. M.: Seasonal variability of aerosol optical properties observed by means of a Raman lidar at an EARLINET site over Northeastern Spain, Atmos. Chem. Phys., 11, 175–190, doi:10.5194/acp-11-175-2011, 2011. Sicard, M., Bertolín, S., Mallet, M., Dubuisson, P., and Comerón, A.: Estimation of mineral dust long-wave radiative forcing: sensitivity study to particle properties and application to real cases in the region of Barcelona, Atmos. Chem. Phys., 14, 9213–9231, doi:10.5194/acp-14-9213-2014, 2014a. Sicard, M., Bertolín, S., Muñoz, C., Rodríguez, A., Rocadenbosch, F., and Comerón, A.: Separation of aerosol fineand coarsemode radiative properties: Effect on the mineral dust longwave, direct radiative forcing, Geophys. Res. Lett., 41, 6978–6985, doi:10.1002/2014GL060946, 2014b. Sicard, M., Barragan, R., Muñoz-Porcar, C., Comerón, A., Mallet, M., Dulac, F., Pelon, J., Alados-Arboledas, L., Amodeo, A., Boselli, A., Bravo-Aranda, J. A., D’Amico, G., GranadosMuñoz, M. J., Leto, G., Guerrero-Rascado, J. L., Madonna, F., Mona, L., Pappalardo, G., Perrone, M. R., Burlizzi, P., Rocadenbosch, F., Rodríguez-Gómez, A., Scollo, Spinelli, N., Titos, G., Wang, X., and Zanmar Sanchez, R.: Contribution of EARLINET/ACTRIS to the summer 2013 Special Observing Period of the ChArMEx project, Óptica Pura y Aplicada, submitted, 2015a. Sicard, M., Barragan, Dulac, F., Alados-Arboledas, L., and Mallet, M.: Aerosol optical, microphysical and ratiative properties at three regional background insular sites in the western Mediterranean Basin, Atmos. Chem. Phys. Discuss., in press, 2015b. Solmon, F., Giorgi., F., and Liousse, C.: Aerosol modelling for regional climate studies: application to anthropogenic particles and evaluation over a European/African domain, Tellus, 58B, 51–72, doi:10.1111/j.1600-0889.2005.00155.x, 2006. Solmon, F., Mallet, M., Elguindi, N., Giorgi, F., Zakey, A., and Konaré, A.: Dust aerosol impact on regional precipitation over western Africa: mechanisms and sensitivity to absorption properties, Geophys. Res. Lett., 35, L24705, doi:10.1029/2008GL035900, 2008. Spada, M., Jorba, O., Pérez García-Pando, C., Janjic, Z., and Baldasano, J. M.: Modeling and evaluation of the global sea-salt aerosol distribution: sensitivity to size-resolved and sea-surface temperature dependent emission schemes, Atmos. Chem. Phys., 13, 11735–11755, doi:10.5194/acp-13-11735-2013, 2013. Tafuro, A. M., Barnaba, F., De Tomasi, F., Perrone, M. R., and Gobbi, G. P.: Saharan dust particle properties over the central Mediterranean, Atmos. Res., 81, 67–93, 2006. Tafuro, A. M., Kinne, S., De Tomasi, F., and Perrone, M. R.: Annual cycle of aerosol direct radiative effect over southeast Italy and sensitivity studies, J. Geophys. Res., 112, D20202, doi:10.1029/2006JD008265, 2007. Tanré, D., Kaufman, Y. J., Herman, M., and Mattoo, S.: Remote sensing of aerosol properties over oceans using the MODIS/EOS spectral radiances, J. Geophys. Res., 102, 16971–16988, 1997. Tanré, D., Bréon, F. M., Deuzé, J. L., Dubovik, O., Ducos, F., François, P., Goloub, P., Herman, M., Lifermann, A., and Waquet, F.: Remote sensing of aerosols by using polarized, directional and spectral measurements within the A-Train: the PARASOL mission, Atmos. Meas. Tech., 4, 1383–1395, doi:10.5194/amt-4-1383-2011, 2011. Tegen, I., Harrison, S. P., Kohfeld, K. E., and Prentice, I. C.: Impact of vegetation and preferential source areas on global dust aerosol: Results from a model study, J. Geophys. Res., 107, 4576, doi:10.1029/2001JD000963, 2002. Thieuleux, F., Moulin, C., Bréon, F. M., Maignan, F., Poitou, J., and Tanré, D.: Remote sensing of aerosols over the oceans using MSG/SEVIRI imagery, Ann. Geophys., 23, 3561–3568, doi:10.5194/angeo-23-3561-2005, 2005. Torres, B., Dubovik, O., Fuertes, D., Lapyonok, T., Toledano, C., Schuster, G. L., Goloub, P., Blarel, L., Barreto, A., Mallet, M., and Tanré, D.: Advanced characterization of aerosol properties from measurements of spectral optical thickness of the atmosphere, in preparation, 2016. Turquety, S., Menut, L., Bessagnet, B., Anav, A., Viovy, N., Maignan, F., and Wooster, M.: APIFLAME v1.0: high-resolution fire emission model and application to the Euro-Mediterranean region, Geosci. Model Dev., 7, 587–612, doi:10.5194/gmd-7-5872014, 2014. Vaishya, A., Jennings, S. G., and O’Dowd, C.: Wind-driven influences on aerosol light scattering in north-east Atlantic air, Geophys. Res. Lett., 39, L05805, doi:10.1029/2011GL050556, 2012. Valenzuela, A., Olmo, F. J., Lyamani, H., Antón, M., Quirantes, A., and Alados-Arboledas, L.: Aerosol radiative forcing during African desert dust events (2005–2010) over Southeastern Spain, Atmos. Chem. Phys., 12, 10331–10351, doi:10.5194/acp12-10331-2012, 2012. Vialard, J., Duvel, J.-P., McPhaden, M., Bouruet-Aubertot, P., Ward, B., Key, E., Bourras, D., Weller, R., Minnett, P., Weill, A., Cassou, C., Eymard, L., Fristedt, T., Basdevant, C., Dandoneau, Y., Duteil, O., Izumo, T., de Boyer Montegut, C., Masson, S., Marsac, F., Menkes, C., and Kennan, S.: Cirene: Air-sea interactions in the Seychelles-Chagos thermocline ridge region, B. Am. Meteorol. Soc., 90, 45–61, doi:10.1175/2008BAMS2499.1, 2009. Villani, P., Picard, D., Marchand, N., and Laj, P.: Design and validation of a 6-volatility tandem differential mobilwww.atmos-chem-phys.net/16/455/2016/ Atmos. Chem. Phys., 16, 455–504, 2016 504 M. Mallet et al.: Overview of the ChArMEx/ADRIMED summer 2013 campaign ity analyzer (VTDMA), Aerosol Sci. Technol., 41, 898–906, doi:10.1080/02786820701534593, 2007. Wang, Y., Sartelet, K. N., Bocquet, M., Chazette, P., Sicard, M., D’Amico, G., Léon, J. F., Alados-Arboledas, L., Amodeo, A., Augustin, P., Bach, J., Belegante, L., Binietoglou, I., Bush, X., Comerón, A., Delbarre, H., García-Vízcaino, D., GuerreroRascado, J. L., Hervo, M., Iarlori, M., Kokkalis, P., Lange, D., Molero, F., Montoux, N., Muñoz, A., Muñoz, C., Nicolae, D., Papayannis, A., Pappalardo, G., Preissler, J., Rizi, V., Rocadenbosch, F., Sellegri, K., Wagner, F., and Dulac, F.: Assimilation of lidar signals: application to aerosol forecasting in the western Mediterranean basin, Atmos. Chem. Phys., 14, 12031–12053, doi:10.5194/acp-14-12031-2014, 2014. Waquet, F., Cornet, C., Deuzé, J.-L., Dubovik, O., Ducos, F., Goloub, P., Herman, M., Lapyonok, T., Labonnote, L. C., Riedi, J., Tanré, D., Thieuleux, F., and Vanbauce, C.: Retrieval of aerosol microphysical and optical properties above liquid clouds from POLDER/PARASOL polarization measurements, Atmos. Meas. Tech., 6, 991–1016, doi:10.5194/amt-6-991-2013, 2013. Wolke, R., Schroeder, W., Schroedner, R., and Renner, E.: Influence of grid resolution and meteorological forcing on simulated European air quality: A sensitivity study with the modeling system COSMO-MUSCAT, Atmos. Environ., 53, 110–130, 2012. Yue, X., Liao, H., Wang, H. J., Li, S. L., and Tang, J. P.: Role of sea surface temperature responses in simulation of the climatic effect of mineral dust aerosol, Atmos. Chem. Phys., 11, 6049– 6062, doi:10.5194/acp-11-6049-2011, 2011. Zakey, A. S., Solmon, F., and Giorgi, F.: Implementation and testing of a desert dust module in a regional climate model, Atmos. Chem. Phys., 6, 4687–4704, doi:10.5194/acp-6-4687-2006, 2006. Zakey, A. S., Giorgi, F., and Bi, X.: Modeling of sea salt in a regional climate model: Fluxes and radiative forcing, J. Geophys. Res., 113, D14221, doi:10.1029/2007JD009209, 2008. Zanis, P., Ntogras, C., Zakey, A., Pytharoulis, I., and Karacostas, T.: Regional climate feedback of anthropogenic aerosols over Europe using RegCM3, Clim. Res., 52, 267–278, doi:10.3354/cr01070, 2012. Zhu, A., Ramanathan, V., Li, F., and Kim, D.: Dust plumes over the Pacific, Indian, and Atlantic oceans: Climatology and radiative impact, J. Geophys. Res., 112, D16208, doi:10.1029/2007JD008427, 2007. Atmos. Chem. Phys., 16, 455–504, 2016 www.atmos-chem-phys.net/16/455/2016/