Impact of mineral dust on shortwave and longwave radiation: evaluation of different vertically resolved parameterizations in 1-D radiative transfer computations

. Aerosol radiative properties are investigated in southeastern Spain during a dust event on 16–17 2013 in the framework of the ChArMEx/ADRIMED (Chemistry-Aerosol Mediterranean Experiment/Aerosol Direct Radiative Impact on the regional climate in the MEDiterranean region) optical and microphysical properties from create a of input parameterizations, feed the 1-D transfer model (RTM) GAME (Global Atmospheric consider three datasets: (1) a ﬁrst parameterization based on the retrievals by an advanced aerosol inversion code (GRASP; Generalized Retrieval of Aerosol and Surface Properties) applied to combined photometer and lidar data, (2) a parameterization based on the photometer columnar optical properties and vertically resolved lidar retrievals with the two-component Klett–Fernald algorithm, and (3) a parameterization based on vertically resolved optical and microphysical aerosol properties measured in situ by the aircraft instrumentation. Once retrieved, the outputs of the RTM in terms of both shortwave and longwave radiative ﬂuxes are compared against ground and in situ airborne measurements. In addition, the outputs of the model in terms of the aerosol direct radiative effect are with respect to the different input parameterizations. Results show that calculated radiative ﬂuxes differ no more than 7 % from the measured ones. The three parameterization datasets produce a effect due to mineral both at the and the Aerosol radiative effects with differences of up to 10 W m − 2 in the shortwave spectral range and W m − Abstract. Aerosol radiative properties are investigated in southeastern Spain during a dust event on 16–17 2013 in the framework of the ChArMEx/ADRIMED (Chemistry-Aerosol Mediterranean Experiment/Aerosol Direct Radiative Impact on the regional climate in the MEDiterranean region) optical and microphysical properties from create a of input parameterizations, feed the 1-D transfer model (RTM) GAME (Global Atmospheric consider three datasets: (1) a ﬁrst parameterization based on the retrievals by an advanced aerosol inversion code (GRASP; Generalized Retrieval of Aerosol and Surface Properties) applied to combined photometer and lidar data, (2) a parameterization based on the photometer columnar optical properties and vertically resolved lidar retrievals with the two-component Klett–Fernald algorithm, and (3) a parameterization based on vertically resolved optical and microphysical aerosol properties measured in situ by the aircraft instrumentation. Once retrieved, the outputs of the RTM in terms of both shortwave and longwave radiative ﬂuxes are compared against ground and in situ airborne measurements. In addition, the outputs of the model in terms of the aerosol direct radiative effect are with respect to the different input parameterizations. Results show that calculated radiative ﬂuxes differ no more than 7 % from the measured ones. The three parameterization datasets produce a effect due to mineral both at the and the Aerosol radiative effects with differences of up to 10 W m − 2 in the shortwave spectral range and W m −

Abstract. Aerosol radiative properties are investigated in southeastern Spain during a dust event on 16-17 June 2013 in the framework of the ChArMEx/ADRIMED (Chemistry-Aerosol Mediterranean Experiment/Aerosol Direct Radiative Impact on the regional climate in the MEDiterranean region) campaign. Particle optical and microphysical properties from ground-based sun/sky photometer and lidar measurements, as well as in situ measurements on board the SAFIRE ATR 42 French research aircraft, are used to create a set of different levels of input parameterizations, which feed the 1-D radiative transfer model (RTM) GAME (Global Atmospheric ModEl). We consider three datasets: (1) a first parameterization based on the retrievals by an advanced aerosol inversion code (GRASP; Generalized Retrieval of Aerosol and Surface Properties) applied to combined photometer and lidar data, (2) a parameterization based on the photometer columnar optical properties and vertically resolved lidar retrievals with the two-component Klett-Fernald algorithm, and (3) a parameterization based on vertically resolved optical and mi-crophysical aerosol properties measured in situ by the aircraft instrumentation. Once retrieved, the outputs of the RTM in terms of both shortwave and longwave radiative fluxes are compared against ground and in situ airborne measurements. In addition, the outputs of the model in terms of the aerosol direct radiative effect are discussed with respect to the different input parameterizations. Results show that calculated atmospheric radiative fluxes differ no more than 7 % from the measured ones. The three parameterization datasets produce a cooling effect due to mineral dust both at the surface and the top of the atmosphere. Aerosol radiative effects with differences of up to 10 W m −2 in the shortwave spectral range (mostly due to differences in the aerosol optical depth) and 2 W m −2 for the longwave spectral range (mainly due to differences in the aerosol optical depth but also to the coarse mode radius used to calculate the radiative properties) are obtained when comparing the three parameterizations. The study reveals the complexity of parameterizing 1-D RTMs as sizing and characterizing the optical properties of min-eral dust is challenging. The use of advanced remote sensing data and processing, in combination with closure studies on the optical and microphysical properties from in situ aircraft measurements when available, is recommended.

Introduction
The radiative effect by atmospheric aerosol is estimated to produce a net cooling effect of the Earth's climate. However, an accurate quantification of this cooling is extremely difficult. In fact, the aerosol radiative effect (ARE) is affected by large uncertainties. Due to the direct aerosol-radiation interaction, the ARE is estimated to be −0.27 W m −2 on average at the global scale, with an uncertainty range of −0.77 to −0.23 W m −2 ; whereas the radiative effect related to cloud adjustments due to aerosols is −0.55 W m −2 (−1.33 to −0.06 W m −2 ; , which is the largest unknown in the radiative forcing of the atmosphere. The extent to which the ARE uncertainty range reported is due to physical processes or due to the measurement uncertainty itself is still hard to quantify. In previous studies, the AREs in longwave (LW) spectral range were commonly neglected due to the complexity of an accurate quantification of the optical properties in this spectral range Mallet et al., 2008;. However, the contribution of the LW component to the ARE is nonnegligible for large aerosol particles, i.e., marine aerosol or mineral dust (e.g., Markowicz et al., 2003;Otto et al., 2007;Sicard et al., 2014a, b;Meloni et al., 2018).
The contribution of mineral dust to the ARE in the infrared spectral range is especially relevant because of its large size and abundance (Meloni et al., 2018). Mineral dust is estimated to be the most abundant aerosol type in the atmosphere by mass (e.g., , with global emission between 1000 and 3000 Mt yr −1 . The high temporal and spatial variability in dust concentrations and the variability in their microphysical and optical properties present a significant challenge to our understanding of how these particles impact the environment . Many measurements worldwide have been made using different approaches, including satellites, which can provide global coverage of mineral dust properties. However, the retrievals of particle properties are still affected by large uncertainties (Levy et al., 2013) and the information on mineral dust properties is quite scarce .
One of the areas frequently influenced by mineral dust is the Mediterranean Sea region, affected by dust intrusions from the close by Sahara or the Middle East region (Moulin et al., 1998; producing significant perturbations to the shortwave (SW) and the LW radiation balance Papadimas et al., 2012; as well as the regional climate (Nabat et al., 2015). The ARE in the Mediterranean region can be responsible for a strong cooling effect both at the surface (or bottom of the atmosphere, BOA) and at the top of the atmosphere (TOA). The so-called forcing efficiency (FE), which is defined as the ratio between the ARE and the aerosol optical depth (AOD) for the SW spectral component ranges between −150 and −160 W m −2 for solar zenith angles (SZAs) in the range 50-60 • (di , being able to reach values larger than 200 W m −2 at the BOA during strong dust events in the Mediterranean region . The LW component accounts for an effect of up to 53 % of the SW component and with an opposite sign . The Aerosol Direct Radiative Impact on the regional climate in the MEDiterranean region (ADRIMED) field campaign within the Chemistry-Aerosol Mediterranean Experiment (ChArMEx, http://charmex.lsce.ipsl.fr, last access: 25 January 2018) took place in the Mediterranean region from 11 June to 5 July 2013 . It aimed at characterizing the different aerosol particles and their radiative effects using airborne and ground-based measurements collected in the Mediterranean Basin, with special focus on the western region. In particular, two ChArMEx/ADRIMED flights, F30 and F31 from the French ATR 42 environmental research aircraft of SAFIRE (http://www.SAFIRE.fr, last access: 25 January 2018), took place above southeastern Spain during a Saharan dust episode on 16 and 17 June 2013.
In this paper, we present an analysis of the mineral dust radiative properties during this particular episode and take advantage of the thorough database that is available. Multiple datasets are used as input in a radiative transfer model (RTM) to evaluate the influence of the different measurements and data processing in the retrieved direct ARE. The model used here is the Global Atmospheric ModEl (GAME; Dubuisson et al., 1996, which allows for calculating both the solar and thermal infrared fluxes. An evaluation against aircraft in situ measurements of radiative fluxes is also presented. Two main goals are pursued: (i) the quantification of the direct ARE for two case studies within a dust transport episode and (ii) the evaluation of the model estimate sensitivity to the aerosol input used.
The paper is structured as follows: Sect. 2 includes a description of both the ground-based and in situ aircraft instrumentation and a short description of the retrieval algorithms used for the present study, Sect. 3 is devoted to the description of GAME and the input datasets used here, and results are presented in Sect. 4; finally, a short summary and concluding remarks are included in Sect. 5.

Ground-based measurements
Ground-based measurements used in this work were carried out at the Andalusian Institute for Earth System Research (IISTA-CEAMA) of the University of Granada, Spain (lat 37.16, long −3.61; 680 m a.s.l.) by the Atmospheric Physics Group of the University of Granada (GFAT-UGR). This experiment site is located in the western Mediterranean basin, near the African continent (∼ 200 km). Therefore, long-range transport of mineral dust particles from north Africa is a main source of natural atmospheric aerosol in the region (e.g., Lyamani et al., 2005;. The station is also affected by long-range transported smoke  and fresh smoke from nearby biomass burning . Anthropogenic sources such as pollution from Europe, the Iberian Peninsula, and the Mediterranean Sea (Pérez-Ramírez et al., 2016) also affect the station. Local sources are mainly road traffic and central heating systems . IISTA-CEAMA station is equipped with a CE-318-4 (Cimel Electronique) sun/sky photometer, which belongs to the AERONET network . This instrument performs direct solar irradiance measurements, used to derive AOD, and sky radiance measurements both measured at least at the following nominal wavelengths (λ): 440, 670, 870, and 1020 nm. The AOD product provided by AERONET has uncertainties of ±0.01 for λ > 440 nm and of ±0.02 for λ < 440 nm . AERONET also provides aerosol optical and microphysical properties such as columnar particle size distribution (PSD), real and imaginary parts of the refractive indices (RRI and IRI, respectively), asymmetry factor (g), and single scattering albedo (SSA) using the AOD and sky radiance values in an inversion algorithm . For the present study, AERONET Version 2 Level 1.5 (Level 2.0 when available) data are used. The uncertainty in the retrieval of SSA is ±0.03 for high aerosol load (AOD 440 > 0.4) and SZAs > 50 • ; while for measurements with low aerosol load (AOD 440 < 0.2), the retrieval accuracy of SSA drops down to 0.02-0.07 . For high aerosol load and SZAs > 50 • , errors are about 30 %-50 % for the IRI. For particles in the size range 0.1 < r < 7 µm (r being the aerosol radius), errors in PSD retrievals are around 10 %-35 %, while for sizes lower than 1 µm and higher than 7 µm retrieval errors rise up to 80 %-100 %. The inversion code provides additional variables such as the volume concentration; effective radius, r eff ; and geometric standard deviation of the equivalent lognormal distribution, σ , for fine and coarse modes of the retrieved PSD that will be used in the current study.
The multiwavelength aerosol Raman lidar MULHACEN, based on a customized version of LR331D400 (Raymetrics S.A.) is operated at Granada station as part of EAR-LINET/ACTRIS (European Aerosol Research Lidar Network/Aerosols, Clouds, and Trace Gases Research Infrastructure Network; https://www.actris.eu/default.aspx, last access: 15 June 2018; Pappalardo et al., 2014) since April 2005. The system has a monostatic biaxial configuration, which usually requires an overlap correction to minimize the incomplete overlap effect . The system emits vertically to the zenith by means of a pulsed Nd:YAG laser, with second-and third-harmonic generators, that emits simultaneously at 1064, 532, and 355 nm. The receiving system consists of several detectors, which can split the radiation according to the three elastic channels at 355, 532 (parallel-and perpendicular-polarized;, and at 1064 nm; two nitrogen Raman channels at 387 and 607 nm; and a water vapor Raman channel at 408 nm (Navas-Guzmán et al., 2014). The aerosol backscatter coefficient profiles (β aer (z, λ), z being the vertical height) obtained from the multiwavelength lidar were calculated with the Klett-Fernald method . For the retrieval of the aerosol extinction coefficient profiles (α aer (z, λ)), a height-independent lidar ratio (LR) obtained by forcing the vertical integration of α aer (z, λ) to the AOD from the AERONET photometer  was assumed. The assumption of a constant LRs introduces uncertainty in α aer (z, λ) retrievals, especially when different types of aerosol appear at different layers. In our case, the LR used for the Klett-Fernald retrieval are very similar to those provided by GRASP (see . Considering the different uncertainty sources, total uncertainty in the profiles obtained with the Klett-Fernald method is usually 20 % for β aer (z, λ) and 25 %-30 % for α aer (z, λ) profiles .
Additionally, surface temperature and pressure are continuously monitored at IISTA-CEAMA by a meteorological station located 2 m above the ground. At the same location, the global and diffuse downward radiative fluxes for the SW component are continuously measured with a CM11 pyranometer (Kipp & Zonen) and diffuse downward radiative fluxes for the LW component are measured with a precision infrared radiometer pyrgeometer (Eppley), both being instruments regularly calibrated at the site .

Airborne measurements
The SAFIRE ATR 42 aircraft performed two overpasses above Granada on 16 (flight F30) and 17 June (flight F31) in 2013 during the ChArMEx/ADRIMED campaign. During F30, the SAFIRE ATR 42 descended performing a spiral trajectory from 14:15 to 14:45 UTC; whereas during flight F31, the aircraft ascended in the early morning (from 07:15 to 07:45 UTC) at around 20 km from Granada station (see Fig. 1 from . Additional flight details can be found in previous studies Mallet et al., 2016;. The airborne instrumentation includes a scanning mobility particle sizer (SMPS) and an ultra-high sensitivity aerosol spectrometer (UHSAS) for measuring aerosol number size distribution in the submicron range. The forward-scattering spectrometer probe model 300 (FSSP-300) and the GRIMM optical particle counter (sky-OPC 1.129) were used to measure the optical size distributions in the diameter nominal size range between 0.28 and 20 µm and between 0.3 and 32 µm, respectively. A nephelometer (TSI Inc, model 3563) was used to measure the particle scattering coefficient at 450, 550, and 700 nm, and a cavity attenuated phase shift extinction monitor (CAPS-PMex, Aerodyne Inc.) was employed to obtain the aerosol extinction coefficient (α aer ) at 530 nm. For more details on the aircraft instrumentation see  and references therein. The PLASMA (Photomètre Léger Aéroporté pour la Surveillance des Masses d'Air) system, which is an airborne sun-tracking photometer, was additionally used to obtain AOD with wide spectral coverage (15 channels between 0.34 and 2.25 µm) with an accuracy of approximately 0.01, as well as the vertical profiles of the aerosol extinction coefficient .
Airborne radiative fluxes (F ) were measured with Kipp & Zonen CMP22 pyranometers and CGR4 pyrgeometers. Upward and downward SW fluxes ( ↑ F SW and ↓ F SW ) were measured in the spectral range 297-3100 nm by two instruments located above and below the aircraft fuselage. The same setup was used for the pyrgeometers, which provided the LW upward and downward radiative fluxes ( ↑ F LW and ↓ F LW ) for wavelengths larger than 4 µm. Both pyranometers and pyrgeometers were calibrated in January 2013 and data were corrected for the temperature dependence of the radiometer's sensitivity following .
Radiation measurement data from the aircraft were filtered out for large pitch and roll angles and corrected from the rapid variations in the solar incidence angle around the SZA due to the aircraft attitude (pitch and roll). This correction also depends on aircraft heading angle and solar position. It should be noted that, beforehand, roll and pitch offsets must be determined (the axis sensor is not necessarily vertical on average during a horizontal leg). Cosine errors were taken into account. Finally, data were corrected from variations in the SZA during the flight to ease the comparison with GAME retrievals. After these various corrections, an estimated uncertainty of ±5 W m −2 is considered to affect the data, taking into account the accuracy of the calibration and the acquisition system together with the consistency of airborne measurements (Meloni et al., 2018).

The GRASP code
The GRASP (Generalized Retrieval of Aerosol and Surface Properties) code  provides aerosol optical and microphysical properties in the atmosphere by combining the information from a variety of remote sensors (e.g., . In our case, GRASP was used to invert simultaneously coincident lidar data (range-corrected signal, RCS, at 355, 532, and 1064 nm) and sun/sky photometer measurements (AOD and sky radiances both from AERONET at 440, 675, 870, and 1020 nm) providing a detailed characterization of the aerosol properties, both column-integrated and vertically resolved. It is worthy to note that this GRASP scheme, based on Lopatin et al. (2013), presents the main advantage that it allows for retrieving aerosol optical and microphysical properties for two distinct aerosol modes, namely fine and coarse. The α aer , β aer , SSA (all at 355, 440, 532, 675, 870, 1020, and 1064 nm), and aerosol volume concentration (VC) profiles obtained as output from GRASP will be used as input to GAME in the present study, together with the column-integrated PSD properties (namely r eff and σ for fine and coarse modes). A more in-depth analysis of GRASP output data retrieved using the lidar and sun/sky photometer data at Granada station for the two inversions coinciding with the aircraft overpasses during flights F30 and F31 during ChArMEx/ADRIMED campaign can be found in .
3 GAME radiative transfer model

GAME description
The GAME code is widely described by Dubuisson et al. (2004 and . It is a modular RTM that allows for calculating upward and downward radiative fluxes at different vertical levels from the ground up to 20 km (100 km) in the SW (LW) spectral range. The solar and thermal infrared fluxes are calculated in two adjustable spectral ranges, which in this study were fixed to match those of the aircraft radiation measurements (namely 297-3100 nm for the SW and 4.5-40 µm for the LW) by using the discrete ordinates method . Note that the GAME code has a variable spectral sampling in the SW (depending on the spectral range considered) and a fixed spectral sampling (115 values) in the LW spectral range (Table 1).

GAME input data parameterization
The two considered SAFIRE ATR 42 flights, F30 and F31, took place on 16 and 17 June 2013, respectively, coinciding with ground-based lidar and sun/sky photometer measurements performed at the station. On these days, mineral dust with origin in the Sahara region (southern Morocco near the border with Algeria) reached Granada after ∼ 4 days of traveling, according to back-trajectory analysis (see Supplement Fig. S1) and the results presented in . A homogenous dust layer reaching up to 5 km a.g.l. was ob- served on 16 June, whereas on 17 June the dust layer was decoupled from the boundary layer and located between 2 and 4.5 km a.g.l. . A very similar vertical structure was observed for the same dust event above Minorca . Daily maps of Meteosat Second Generation-derived AOD over the Mediterranean from 15 to 18 June during the dust event shown in Fig. 4  , and simultaneous lidar and sun/sky photometer data were available. Unfortunately, the airborne vertical profile of extinction by the CAPS measurements was not available during this second flight. Clouds were detected by the lidar on 17 June after 15:00 UTC. Furthermore, a sky camera and the ground-based pyranometer and pyrgeometer data indicate cloud contamination (but not in the zenith) in the radiation data much earlier (around 09:00 UTC), also preventing satellite retrievals in the region. A summary of the experimental data used as input for GAME calculations during these two case studies is presented in Table 2. This input includes surface parameters and atmospheric profiles of meteorological variables, main gas concentrations, and aerosol properties. The aerosol properties used in the present study are parameterized using three different datasets, based on the different instrumentation and retrievals available, i.e., Dataset 1 (DS1), Dataset 2 (DS2) and Dataset 3 (DS3). A more detailed description of the different parameters is provided next.

Surface parameters and profiles of meteorological variables
The surface parameters required for GAME are the surface albedo (alb(λ)) and land-surface temperature (LST). The alb(λ) for the SW range is obtained from the sun/sky photometer data using the AERONET retrieval at 440, 675, 880, and 1020 nm, and for the LW from the integrated emissivity between 4 and 100 µm provided by the Single Scanner Footprint (SSF) Level2 products of the CERES (Clouds and the Earth's Radiant Energy System (http://ceres.larc.nasa. gov/, last access: 10 December 2018) instrument (Table 3). LST values are obtained from MODIS (Moderate Resolution Imaging Spectroradiometer) 1 km daily level-3 data (Wan, 2014) on 16 June. Unfortunately, on 17 June MODIS data were not available due to the presence of clouds and the local surface temperature was obtained from temperature measurements at the Granada site, where the meteorological station is located at 2 m above the ground. LST and alb(λ) values used for the two analyzed cases are included in Table 3. Table 2. Summary of the data sources used to obtain the input data parameterizations for GAME computations both in the SW and LW spectral ranges, including the surface parameters (albedo, alb, and land-surface temperature, LST), profiles of meteorological variables, and main gases and the aerosol parameters. For the aerosol parameters (aerosol extinction, α aer ; single scattering albedo, SSA; and asymmetry parameter, g) three different datasets are used (DS1, DS2, and DS3) based on different instrumentation and retrievals. The indications below the sources of the aerosol parameters indicate whether the parameter is column-integrated (col) or if it is vertically resolved (z) and the number of wavelengths at which it is given (nλ   Figure 1 shows the pressure (P ), temperature (T ), and relative humidity (RH) profiles obtained from the SAFIRE ATR 42 measurements. Data from the meteorological station located at IISTA-CEAMA are used to complete these profiles at the surface level; whereas at altitudes above the aircraft flight, a scaled US standard atmosphere is used for completion. The concentration profiles of the main absorbing gases (O 3 , CH 4 , N 2 O, CO, and CO 2 ) are also taken from the US standard atmosphere, while for the gaseous absorption coefficients the HITRAN database is used (as in Sicard et al., 2014a, b). Variations in the concentration profiles of the main absorbing gases have low impact of the radiative fluxes and the ARE, thus small uncertainty is introduced by this approach. A sensitivity test performed in the present study, varying the O 3 profiles up to double concentrations, indicates maximum differences of 4 W m −2 in the F SW and 3.6 W m −2 in the case of the F LW . For the ARE, differences are negligible (below 0.2 W m −2 ).

Aerosol parameterization
As for the aerosol parameterization, α aer (λ, z), SSA(λ, z), and g(λ, z) are required as GAME input data (Table 2). For the SW wavelengths, these properties can be obtained from the measurements performed with the instrumentation available during the campaign; namely the lidar, the sun/sky pho-tometer, and the in situ instrumentation on board the aircraft. On the other hand, direct measurements of the aerosol properties in the LW are not so straightforward and thus scarce. Hence, the aerosol LW radiative properties are calculated by a Mie code included as a module in GAME. According to , the dust particles nonsphericity effect at the thermal infrared wavelengths is not significant on the LW direct ARE, thus the shape of the mineral dust can be assumed as spherical for the Mie code retrievals introducing negligible uncertainties.
For the SW simulations, we run GAME using three different aerosol input datasets, i.e., DS1, DS2, and DS3 (Table 2), in order to evaluate their influence on the ARE calculations. DS1 relies on a parameterization based on the advanced postprocessing GRASP code, which combines lidar and sun/sky photometer data to retrieve aerosol optical and microphysical property profiles; DS2 relies on Klett-Fernald lidar inversions and AERONET products and corresponds to a reference parameterization (easily reproducible at any station equipped with a single-or multiwavelength lidar and an AERONET sun/sky photometer and without the need of an advanced postprocessing algorithm); and DS3 relies on in situ airborne measurements and corresponds to an alternative parameterization to DS1 and DS2.   For DS1 ( Fig. 2a and d), α aer profiles at seven different wavelengths obtained with GRASP are used as input data in GAME. In DS2 ( Fig. 2b and e), the α aer profiles are obtained from the lidar data using Klett-Fernald retrievals and adjusting the lidar ratio to the AERONET retrieved AODs, as mentioned in Sect. 2.1. Finally, for DS3 ( Fig. 2c and f) the α aer values are obtained from the aircraft in situ measurements (CAPS and PLASMA data on 16 June and PLASMA on 17 June). A detailed analysis and discussion on the comparison between α aer profiles provided by the aircraft measurements, GRASP, and the lidar system at Granada is already included in . In general, the lidar, GRASP, and the CAPS data are in accordance, observing the same aerosol layers and similar values, with differences within 20 %. GRASP slightly overestimates CAPS data by 3 Mm −1 on average, whereas the differences with PLASMA are larger, reaching 30 % (or 11 Mm −1 ). In the case of the Klett-Fernald retrieval, values are lower than those retrieved with GRASP by up to 19 %. Considering that the uncertainty in α aer is around 30 % for both GRASP and the Klett-Fernald retrieval and 3 % for the CAPS data, this discrepancy is well below the combined uncertainty in the different datasets. Dif-ferences in the α aer profiles translate into differences in the integrated extinction and, hence, into differences in the AOD values used as input in the radiative flux retrievals. The AOD values presented here (included in Table 4) are obtained by integrating the α aer profiles at 550 nm from the surface up to the considered top of the aerosol layer (4.3 km on 16 June and 4.7 km on 17 June). In GRASP retrieved α aer profiles, values above this top of the aerosol layer are slightly larger than zero since GRASP takes into account stratospheric aerosols by an exponential decay (Lopatin et al., 2013), thus the approach used here to calculate the AOD leads to lower values compared to the column-integrated AOD provided by the sun photometer. Differences among the three datasets are more noticeable on 16 June, when the AOD for DS1 is 0.05 lower than for DS2 and DS3; whereas on 17 June the maximum difference is 0.03, obtained between DS1 and DS2. The AOD values at 550 nm reveal that GRASP input data (DS1) and to a lesser extent the aircraft in situ data (DS3) underestimate the aerosol load in the analyzed dust layer compared to AERONET (DS2) due to the differences in the retrieval techniques, e.g., although AERONET provides integrated AOD for the whole column, low α aer values above the aerosol layer are neglected for the AOD calculations in DS1 and DS3.     Fig. 3b). The mean SSA at 440 nm is equal to 0.92 on 15 June, whereas on 17 June it is 0.85. On 17 June the SSA profiles present lower values and more variation with height than on 16 June; the lower SSA values indicate the presence of more absorbing particles on 17 June. The vertical variation on 17 June is associated with the presence of two different layers, whereas a more homogeneous dust layer is observed on 16 June. For DS2, the SSA are taken from AERONET columnar values and assumed to be constant with height (Fig. 4a). The SSA at 440 nm was 0.89 and 0.83 on 16 and 17 June, respectively; as already observed in Fig. 3, SSA values are lower on 17 June due to the intrusion of more absorbing particles. For DS3, SSA values at 530 nm are obtained from the nephelometer and the CAPS or PLASMA on board the ATR. In order to reduce the uncertainty in the measured data, only averaged values for the column will be considered, being 0.88 and 0.83 on 16 and 17 June, respectively (Fig. 4). Therefore, differences of up to 0.04 and 0.02 are observed on 16 and 17 June, respectively, among the SSA values obtained with the three datasets. Despite these differences, the retrieved SSA values obtained here are within the range of typical values for dust aerosols Lopatin et al., 2013) and differences are still within the uncertainty limits, which range between 0.02 and 0.07 depending on the aerosol load for AERONET data  and is 0.04 for the aircraft values. In the case of g values, the same data are used for the three aerosol input datasets. Multispectral values of g are taken from AERONET columnar values and assumed to be constant with height ( Fig. 4b).
Summing up, for the SW aerosol parametrization in GAME three datasets are tested. In DS1, GRASP-derived spectral profiles at seven wavelengths of the aerosol extinction and SSA are used. In DS2, the Klett retrieved extinction profiles at three wavelengths are used together with the AERONET SSA columnar values at four wavelengths, which are assumed to be constant with height. For DS3, one extinction profile at 550 nm and a column-averaged single-wavelength value of the SSA from the airborne measurements are considered. In the three cases, the columnintegrated AERONET asymmetry parameter at four wavelengths is assumed to be constant with height and used as input.
For the LW calculations, the Mie code is used to obtain α aer (λ, z), SSA(λ, z), and g(λ, z) from the information on the aerosol PSD, complex refractive index (RI), and density, following a similar approach to that used in previous studies (Meloni et al., , 2018Peris-Ferrús et al., 2017). A summary of the aerosol parameters used in the Mie calculations is included in Table 5. Three different datasets are also used for the aerosol parameterization in the LW calculations. In this case, the sensitivity of the model to the PSD used is tested. A similar scheme to that presented for the SW is used, where DS1 relies on GRASP retrievals, DS2 on AERONET products, and DS3 relies on in situ airborne measurements.
Atmos. Chem. Phys., 19, 523-542, 2019 www.atmos-chem-phys.net/19/523/2019/  Table 5. Summary of the data used to obtain α aer (λ, z), SSA(λ, z), and g(λ, z) in the LW spectral range from Mie calculations, i.e., the refractive index, RI; effective radius, r eff ; geometric standard deviation, σ ; and number concentration, N. Three different datasets are used (DS1, DS2, and DS3) based on different particle size distribution (PSD) data used. The indications below the sources of the aerosol parameters indicate whether the parameter is column-integrated (col) or if it is vertically resolved (z) and the number of wavelengths at which it is given (nλ). DB (2017) stands for Di .
The spectral real and imaginary parts of the RI of mineral dust in the LW are obtained from Di , using the Morocco source, and assumed constant with height. The analysis by Di  only covers the spectral range 3-16 µm so an extrapolation assuming the spectral dependence presented in  for shorter and longer wavelengths is performed. This assumption is not exempt from uncertainty, since the refractive index presents a certain variability associated with the different nature of mineral dust properties. For example, the use of the refractive index provided for the Algerian and Mauritanian sources from Di  leads to variations in the ARE of 0.8 and 0.3 W m −2 at the BOA and the TOA, respectively. Additionally, vertical variations in the refractive index are also a source of uncertainty in the obtained radiative fluxes. The mineral dust particle density is assumed to be 2.6 g cm −3 . Regarding the PSD, three parameters (namely the effective radii, r eff ; standard deviation, σ ; and the numeric concentrations, N ) for fine and coarse modes are used. The fine mode comprises particles within the diameter range 0.1-1 µm, whereas for the coarse mode the range 1-30 µm is considered. A third mode at about 30 µm was detected in  for the same dust event using balloon-borne measurements with concentrations up to 10 −4 particles cm −3 . However, this giant mode is not considered in our study due to the lack of data above Granada. Con-sidering the relevance of large particles for the ARE LW (i.e., Sicard et al., 2014a, b;Meloni et al., 2018), neglecting this giant mode may contribute to increase the uncertainties in GAME estimations. However, simulations with GAME assuming the presence of a third mode of similar characteristics to the one observed by  indicate that variations in the ARE are negligible in this case (lower than 0.1 W m −2 ). Even for much higher concentrations (10 −1 particles cm −3 ), variations in the ARE of just 0.3 W m −2 at the BOA and 0.15 W m −2 at the TOA are obtained.
In the case of DS1, N values are obtained from the volume concentration profiles provided by GRASP assuming spherical particles in the range between 0.05 and 15 µm radii (Fig. 5). Values of r eff and σ provided by GRASP (Table 4) are column-integrated and thus assumed to be constant with height. This is also the case for DS2, in which the PSD parameters are column-integrated values provided by the AERONET retrieval in Granada (see Table 4).
For DS3, the volume concentration (or the equivalent N ), r, and σ profiles for the fine and coarse modes ( 42 and retrieved with GRASP, with differences in the total volume concentration profiles for the dust layers lower than 8 µm −3 cm −2 (20 %), which fall within the combined uncertainty. Nonetheless differences are still noticeable, especially in the fine mode. On 17 June, GRASP overestimates the aircraft measurements for the fine mode and underestimates them for the coarse mode, which in turn results in a quite different fine to coarse concentration ratio for DS1 and DS3. Additionally, a slight shift is observed in the vertical structure of the aerosol layers. Differences are mostly technical, i.e., the GRASP retrieval is based on 30 min averaged lidar profiles while the aircraft provide instantaneous measurements, but they can also be partially caused by the discrepancies between the vertical aerosol distribution above Granada (sampled by the lidar) and the concentration measured during the aircraft trajectory as they are not exactly coincident. In addition, for 16 June, there is a 2 h time difference between the sun/sky photometer retrieval used in GRASP calculations and the airborne measurements, which can lead to slight differences in the aerosol properties despite the homogeneity of the dust event during this period. In the following, we quantify the impact these differences may introduce in the calculations of F .

GAME output data
As a result of the simulation, GAME provides vertical profiles of radiative fluxes in the shortwave (F SW ) and longwave (F LW ) spectral ranges. The net flux can be calculated from the obtained profiles for both spectral ranges as where the upward and downward arrows are for upward and downward fluxes, respectively. From the obtained radiative flux profiles, the direct ARE profiles are calculated according to the following equation: where F w and F o are the radiative fluxes with and without aerosols, respectively. The direct ARE can be obtained for the SW (ARE SW ) and the LW (ARE LW ) spectral ranges.
4 Mineral dust effect on shortwave and longwave radiation 4.1 SW radiative fluxes Figure 6 shows the radiative flux profiles for the SW spectral range obtained with GAME using the three different input datasets described in Sect. 3, as well as the Net F SW . The radiative fluxes measured by the pyranometer on board the SAFIRE ATR 42 are also included in the figure. The three GAME simulations show similar values with differences below 8 W m −2 on average, which represents less than 1 % variation. The differences in the obtained fluxes are mostly due to the differences in the aerosol load considered depending on the inputs. Even though the differences in the AOD among the different datasets are small (lower than 0.05), they can lead to differences in F SW and ultimately in the ARE SW . In order to quantify these differences, we performed a sensitivity test by varying the AOD while the other parameters were kept constant. We observed a maximum variation in the F SW of 6.5 W m −2 (0.7 %) at the surface, decreasing with height, for changes in the AOD of up to 0.05, which is the difference we observe between the AOD for DS2 and DS1 on 16 June. This result partly explains the differences among the three datasets. In addition, a sensitivity test performed by exclusively varying the SSA indicates that more absorbing particles are related to less ↓ F SW at the surface, namely a variation of 1 % is observed at the BOA for a decrease in the SSA of 0.03. The influence of the SSA decreases with height and is negligible at the TOA. For the ↑ F SW , a decrease of 0.8 % is observed at the BOA if more absorbing particles are present, but in this case the influence at the TOA is larger (2.2 %). In our case, the larger AOD assumed for DS2 on both days (see Table 4 and Fig. 2), causes the ↓ F SW to be slightly lower compared to DS1. For DS3 the AOD is similar to DS2, but the SSA values used, which are relatively smaller compared to those measured by AERONET (see Fig. 4), lead to lower values of the radiative fluxes than for DS2. The vertical distribution of the SSA also influences the radiative fluxes in the SW component as demonstrated in previous studies , contributing to explain the differences observed among the three datasets analyzed here.
Atmos. Chem. Phys., 19, 523-542, 2019 www.atmos-chem-phys.net/19/523/2019/ The evaluation against the aircraft measurements shows larger differences for altitudes below 2.5 km (∼ 860 mbar) on 16 June, whereas a better agreement is found above. On 17 June, no ↑ F SW aircraft data are available below 2 km. Relative differences between the model and the aircraft measured data (calculated as (F GAME −F aircraft )/F aircraft ) are well below 7 %, being the largest discrepancies observed for the ↓ F SW . Differences among the three GAME outputs and the aircraft pyranometer are lower than 5 % for the Net F SW on both days. Considering the very different approaches followed by the model and the direct measurements by the airborne pyranometer (i.e., vertical resolution, temporal sampling, and data acquisition and processing) -together with the uncertainty in the pyranometer (5 W m −2 ) and the estimated uncertainty in the model outputs that can be as large as 19 W m −2 ) -these differences are quite reasonable. A conclusive result on which input dataset provides a better performance is unlikely because of the comparable results obtained with the three datasets.
The values at the surface (or BOA) and at the TOA for the different radiative fluxes can also be evaluated against different instruments: measurements for the ↓ F SW at the surface are available from the sun photometer; AERONET provides values for the ↓ F SW and ↑ F SW at both the BOA and TOA. The time series for these measurements corresponding to 16-17 June and the results obtained with GAME for the different datasets are shown in Fig. 7. AERONET surface radiative fluxes have been extensively validated at several sites around the world (e.g.,  and, in addition, all AERONET sun photometers are mandatorily calibrated once a year. Thus, in order to compare GAME results with AERONET data, we have performed additional simulations for the time of the closest AERONET measurement on 16 June (at 16:22 UTC), assuming that the aerosol parameterization is constant with time between the flight time and the photometer measurement. The ↓ F SW values at the surface obtained with GAME are 564.8, 551.8, and 547.0 W m −2 for DS1, DS2, and DS3, respectively, close to the 531.4 W m −2 provided by AERONET. On 17 June, GAME simulations at 07:40 UTC (instead of 07:30 UTC, which is the time of the flight) provide ↓ F SW at the surface of 466.3, 468.3, and 456.4 W m −2 , very similar to the AERONET value of 463.7 W m −2 .
At the TOA, the ↑ F SW between GAME and AERONET are in quite good agreement on both days. On 16 June, the ↑ F SW values obtained with GAME simulations are equal to 152.0, 153.0, and 148.5 W m −2 , and with AERONET they are equal to 146.2 W m −2 . On 17 June, the obtained values with GAME are 133.6, 136.6, and 130.9 W m −2 for DS1, DS2, and DS3, and 131.6 W m −2 for AERONET.
The ARE SW profiles, calculated by using Eq.
(2) and GAME simulations for the three input datasets, are shown in Fig. 8, together with the simultaneous values provided by AERONET on 17 June at the BOA and TOA. Comparing the three GAME simulations, we can see that the low discrepancies in the F profiles from Fig. 6 lead to variations in the ARE SW of 10 %-27 % (3-10 W m −2 ) over the averaged profile depending on the input dataset used. The variations in the ARE SW values are tightly connected to differences in the AOD considered as input in the model, as already observed in previous studies Lolli et al., 2018;Meloni et al., 2018). The SSA and the vertical distribution of the aerosol also plays an important role, as observed for DS3, which shows a profile that is quite different from DS2 despite the AOD being quite close for both datasets.
Differences are also observed when comparing ARE SW values obtained from GAME to those retrieved by AERONET. Contrary to GAME simulations, AERONET does not consider the vertical distribution of the aerosols when calculating the ARE SW , and the definition of the ARE SW at the BOA ( BOA ARE SW ) is slightly different. Indeed, AERONET BOA ARE SW is calculated as the difference between the downward fluxes with and without aerosols, the difference between the upward fluxes (reflected by the Earth) being neglected. Considering this, we can correct the BOA ARE SW provided by AERONET by multiplying by a factor 1−alb(λ). The corrected BOA ARE SW value on 17 June is thus −31.9 W m −2 , which is within the range of values provided by GAME at the surface. All discrepancies observed here are mostly intrinsic to the different techniques used for the acquisition of the data and the retrieval algorithms. The effect of the data processing has also been observed in pre- vious studies (Lolli et al., 2018). Moreover, the sensitivity tests performed reveal that an increase in the AOD of 0.05 can lead to a stronger effect of the ARE both at the BOA (up to 6.7 W m −2 ) and the TOA (up to 2.5 W m −2 ), and more absorbing particles (decrease in the SSA of 0.03) lead to more ARE at the BOA and less at the TOA (4 and 2 W m −2 in absolute terms, respectively). Therefore, the differences among the datasets are within the estimated uncertainty.
The ARE SW values obtained at the BOA and TOA for the three datasets and the averaged value, as well as the FE, are included in Table 6. Both at the BOA and TOA, the ARE SW has a cooling effect, as expected for mineral dust in this region according to values obtained in the literature (e.g., Mallet et al., 2016). Differences among the three datasets lead to variations in the ARE SW of up to 30 % (or 20 % for the FE), observing larger variability on 16 June. The values of the ARE SW and the FE SW are highly dependent on the SZA and a straightforward comparison with previous studies is not simple. Nonetheless, the values obtained for this case are within the range of previous values observed in the western Mediterranean region for similar values of SZA, e.g., FE between −263.4 and −157.1 W m −2 at the BOA and −23.8 and −86.2 W m −2 for SZA = 60 • or ARE values ranging between −93.1 and −0.5 W m −2 at the BOA and between −34.5 and +8.5 W m −2 at the TOA for different SZA values (e.g., Sicard et al., 2014a, b;.
Atmos. Chem. Phys., 19, 523-542, 2019 www.atmos-chem-phys.net/19/523/2019/   Figure 9 shows F LW calculated with GAME after obtaining the aerosol properties in the LW spectral range from Mie calculations for the three mentioned datasets (see Sect. 3.2.2). F LW values measured by pyrgeometers located on board the ATR are also shown.

LW radiative fluxes
In general, differences in the F LW are always lower than 6 % (lower than 10 W m −2 on average), with the airborne values being overestimated by the model on 16 June and underestimated on 17 June. On this latter day, larger differences are observed on the Net F LW compared to 16 June, which might be explained by the inaccurate value of LST used due to the lack of precise data. A sensitivity test performed by increasing the air surface temperature measured at the meteorological station by 5 K indicates that ↑ F LW increases its value up to 30 W m −2 at the surface, and around 10 W m −2 from 1 km onwards, which is nonnegligible. This would lead to an overestimation of the aircraft measured values, but still within a 6 % difference. This highlights the need for accurate LST measurements for radiation simulations in the LW spectral range. Additionally, a sensitivity test performed by assuming a 10 % uncertainty in the PSD parameters (r eff , N , and σ ) leads to an estimated uncertainty in the F LW retrieved by GAME of around 1.2 W m −2 . As stated before, the assumption of the refractive index can also introduce variations as large as 0.8 W m −2 . Considering the uncertainty in the pyrgeometer, the fact that the aircraft and the model present dif- Figure 10. Time series of the ↓ F LW at the BOA during the period 16-17 June. Surface measurements of diffuse (red) radiation from the ground-based pyranometer at the Granada station are included. GAME output data for different inputs are represented by the blue circles (DS1), red (DS2), and green (DS3) crosses. ferent vertical resolutions and time samplings, and the uncertainties due to the use of the standard atmosphere or the parameterization of the surface properties, the obtained differences are not significant.
A comparison of GAME results against the observations from the ground-based pyrgeometer at the Granada station is included in Fig. 10. At the BOA, the longwave radiation measured by the pyrgeometer is in quite good agreement with GAME calculations on 16 June, with differences within 1 W m −2 . However, GAME overestimates the pyrgeometer data by 5 W m −2 (1.3 %) on 17 June. This difference on 17 June, even though larger than on 16 June, is still within the uncertainty limits.
As for the ARE LW , Fig. 11 shows the profiles obtained with GAME using the three datasets as inputs. Values at the BOA and TOA for each dataset and the average values are included in Table 7, together with the FE. Opposite to the SW, the ARE LW produces a heating effect both at the BOA and TOA with positive values. The slight differences in the F LW in Fig. 9 due to the use of different aerosol input datasets lead to variations of up to 2 W m −2 in the BOA ARE LW (ranging from 20 % to 26 %), which needs to be considered in the interpretation of the results and reduced for a better estimate of the direct ARE. Despite this, values obtained for this dust event (3.2 W m −2 on average for both days) are in agreement with previous studies performed for mineral dust in the infrared region (Sicard et al., 2014a, b) and the FE obtained are comparable to those reported by Meloni et al. (2018). It is extremely interesting to look at the differences between the two days in terms of AOD ( AOD) and the effective radius for the coarse mode, r eff,c , ( r eff,c ) and their implication for the differences in the ARE LW at the BOA ( BOA ARE LW ). For DS1, AOD ( r eff,c ) is −0.02 (+0.18 µm), which produces a decrease in BOA ARE LW ( BOA ARE LW = −0.5 W m −2 ). For DS2, AOD ( r eff,c ) is −0.04 (+0.18 µm), which produces a decrease in BOA ARE LW ( BOA ARE LW = −1.0 W m −2 ). If we relate these variations to the sensitivity study of , the expected ARE LW increase due to the increase in the coarse mode radii is counterbalanced by the ARE LW decrease when AOD decreases in both cases. Oppositely, for DS3 the AOD ( r c ) is −0.05 (+0.64 µm), producing an increase in BOA ARE LW ( BOA ARE LW = +1.6 W m −2 ). Here, the large increase in the coarse mode radius dominates over the AOD decrease.  indeed show that the largest positive gradient of ARE LW occurs for median radii ranging from 0.1 to 2.0 µm. For DS3 the increase in BOA ARE LW produced by a positive r c is larger than the decrease in BOA ARE LW that would have produced AOD alone. At the TOA, same trends, but much less marked, are observed.

Total mineral dust radiative effect
The total ARE, including both the SW and LW component, is included in Fig. 12 +3.9 (+17.0) +2.9 (+12.6) +2.9 (+15.3) +1.7 (+8.9) −0.04 +0.18 −1.0 DS3 +2.5 (+10.9) +1.3 (+5.7) +4.1 (+22.8) +1.8 (+10.0) −0.05 +0.64 +1.6 Avg. ± SD +3.2 ± 0.7 +2.1 ± 0.8 +3.2 ± 0.8 +1.7 ± 0.1 (+15.0 ± 3.6) (+10.2 ± 3.9) (+18.1 ± 4.1) (+9.6 ± 0.6) Table 8. ARE (and FE indicated between parenthesis) at the BOA and the TOA for the total (SW+LW) spectral range obtained with GAME using DS1, DS2, and DS3 as inputs for 16 and 17 June 2013. The averaged values and standard deviation are also included. imately 20 % of the ARE SW near the surface (except for DS3 on 16 June), and reaches up to 50 % at higher altitudes where the total ARE is quite low (see 16 June in Fig. 12). Overall these ARE LW to ARE SW ratios are in agreement with those found at the BOA in previous studies for the Mediterranean region, which ranged between 9 % and 26 % . As for the TOA, larger ratios are obtained there on 16 June, but it is worthy to note that results are not directly comparable to previous studies because of the differences in SZA and the different vertical resolution in GAME for the SW and LW components above 4 km, which may lead to numeric artifacts in the obtained results.

Conclusions
A moderate Saharan dust event affecting the western Mediterranean region during the ChArMEx/ADRIMED campaign on June 2013 was extensively monitored by ground-based and aircraft instrumentation above the Granada experimental site. Radiative fluxes and mineral dust ARE both in the solar and infrared spectral ranges are calculated for this event with the RTM GAME. Three different aerosol input datasets are used by the GAME RTM in order to evaluate the impact of different input data in GAME calculations.
For the SW, very low variability with the input aerosol data (less than 1 %) is observed for the radiative fluxes. The evaluation of GAME-calculated radiative fluxes against the aircraft data reveals differences between the model fluxes and the measurements below 7 %, with better agreement at altitudes above the planetary boundary layer. The differences between the retrievals with the three aerosol datasets are quite insignificant, especially taking into account the different approaches followed by the model and the pyranometers and the estimated uncertainties for both the measured data (5 W m −2 ) and the model (around 12 W m −2 ). Thus a conclusion on which input dataset provides a better performance is unlikely. The small differences between GAME radiative flux retrievals lead to variations in the ARE SW of up to 33 %, mostly driven by the differences in the aerosol vertical distribution and load, followed by the SSA.
For the LW component, the effect of aerosol on the radiative properties is lower compared to the SW, but certainly nonnegligible and of opposite sign. GAME retrievals using the three aerosol datasets reveal differences in the fluxes lower than 2 W m −2 (less than 1 %). The comparison with the pyrgeometer data measured at the ATR reveals, however, differences around 7 %. The influence of the assumed CO 2 , O 3 , and the water vapor profiles and LST are needed to fully explain this discrepancy between the aircraft and the simulated profiles.
The total ARE, including both the SW and LW components, confirms that mineral dust produces a cooling effect both at the surface and the TOA, as already reported in the literature. On average, the ARE LW represents 20 % of the ARE SW at the surface, therefore clearly indicating that global model estimates need to consider the complete spectrum to avoid an overestimation of the mineral dust cooling effect. Additionally, it is necessary to be aware of the effects of using different measurement techniques and processing methodologies when calculating aerosol radiative properties. Even though the differences observed here when using different aerosol datasets are slight, they still exist and a homogenization of the techniques to feed global models would be beneficial for a better estimate of the ARE and a reduced uncertainty.
Data availability. Part of the data used in this publication were obtained as part of the AERONET and EARLINET networks and are publicly available. For additional data or information please contact the authors.
Author contributions. MJGM and MS designed the study and wrote the manuscript with contributions from all authors. MJGM, RR, JABO, CD, GB, PF, and BT provided data and performed data analysis. RB provided the initial version of the model. All authors have given approval to the final version of the manuscript.
Competing interests. The authors declare that they have no conflict of interest.
Special issue statement. This article is part of the special issue "CHemistry and AeRosols Mediterranean EXperiments (ChArMEx) (ACP/AMT inter-journal SI)". It is not associated with a conference. Abstract. Aerosol radiative properties are investigated in southeastern Spain during a dust event on 16-17 June 2013 in the framework of the ChArMEx/ADRIMED (Chemistry-Aerosol Mediterranean Experiment/Aerosol Direct Radiative Impact on the regional climate in the MEDiterranean region) campaign. Particle optical and microphysical properties from ground-based sun/sky photometer and lidar measurements, as well as in situ measurements on board the SAFIRE ATR 42 French research aircraft, are used to create a set of different levels of input parameterizations, which feed the 1-D radiative transfer model (RTM) GAME (Global Atmospheric ModEl). We consider three datasets: (1) a first parameterization based on the retrievals by an advanced aerosol inversion code (GRASP; Generalized Retrieval of Aerosol and Surface Properties) applied to combined photometer and lidar data, (2) a parameterization based on the photometer columnar optical properties and vertically resolved lidar retrievals with the two-component Klett-Fernald algorithm, and (3) a parameterization based on vertically resolved optical and mi-crophysical aerosol properties measured in situ by the aircraft instrumentation. Once retrieved, the outputs of the RTM in terms of both shortwave and longwave radiative fluxes are compared against ground and in situ airborne measurements. In addition, the outputs of the model in terms of the aerosol direct radiative effect are discussed with respect to the different input parameterizations. Results show that calculated atmospheric radiative fluxes differ no more than 7 % from the measured ones. The three parameterization datasets produce a cooling effect due to mineral dust both at the surface and the top of the atmosphere. Aerosol radiative effects with differences of up to 10 W m −2 in the shortwave spectral range (mostly due to differences in the aerosol optical depth) and 2 W m −2 for the longwave spectral range (mainly due to differences in the aerosol optical depth but also to the coarse mode radius used to calculate the radiative properties) are obtained when comparing the three parameterizations. The study reveals the complexity of parameterizing 1-D RTMs as sizing and characterizing the optical properties of min-eral dust is challenging. The use of advanced remote sensing data and processing, in combination with closure studies on the optical and microphysical properties from in situ aircraft measurements when available, is recommended.

Introduction
The radiative effect by atmospheric aerosol is estimated to produce a net cooling effect of the Earth's climate. However, an accurate quantification of this cooling is extremely difficult. In fact, the aerosol radiative effect (ARE) is affected by large uncertainties. Due to the direct aerosol-radiation interaction, the ARE is estimated to be −0.27 W m −2 on average at the global scale, with an uncertainty range of −0.77 to −0.23 W m −2 ; whereas the radiative effect related to cloud adjustments due to aerosols is −0.55 W m −2 (−1.33 to −0.06 W m −2 ; , which is the largest unknown in the radiative forcing of the atmosphere. The extent to which the ARE uncertainty range reported is due to physical processes or due to the measurement uncertainty itself is still hard to quantify. In previous studies, the AREs in longwave (LW) spectral range were commonly neglected due to the complexity of an accurate quantification of the optical properties in this spectral range Mallet et al., 2008;. However, the contribution of the LW component to the ARE is nonnegligible for large aerosol particles, i.e., marine aerosol or mineral dust (e.g., Markowicz et al., 2003;Otto et al., 2007;Sicard et al., 2014a, b;Meloni et al., 2018).
The contribution of mineral dust to the ARE in the infrared spectral range is especially relevant because of its large size and abundance (Meloni et al., 2018). Mineral dust is estimated to be the most abundant aerosol type in the atmosphere by mass (e.g., , with global emission between 1000 and 3000 Mt yr −1 . The high temporal and spatial variability in dust concentrations and the variability in their microphysical and optical properties present a significant challenge to our understanding of how these particles impact the environment . Many measurements worldwide have been made using different approaches, including satellites, which can provide global coverage of mineral dust properties. However, the retrievals of particle properties are still affected by large uncertainties (Levy et al., 2013) and the information on mineral dust properties is quite scarce .
One of the areas frequently influenced by mineral dust is the Mediterranean Sea region, affected by dust intrusions from the close by Sahara or the Middle East region (Moulin et al., 1998; producing significant perturbations to the shortwave (SW) and the LW radiation balance Papadimas et al., 2012; as well as the regional climate (Nabat et al., 2015). The ARE in the Mediterranean region can be responsible for a strong cooling effect both at the surface (or bottom of the atmosphere, BOA) and at the top of the atmosphere (TOA). The so-called forcing efficiency (FE), which is defined as the ratio between the ARE and the aerosol optical depth (AOD) for the SW spectral component ranges between −150 and −160 W m −2 for solar zenith angles (SZAs) in the range 50-60 • (di , being able to reach values larger than 200 W m −2 at the BOA during strong dust events in the Mediterranean region (Gómez-Amo et al., 2011). The LW component accounts for an effect of up to 53 % of the SW component and with an opposite sign .
The Aerosol Direct Radiative Impact on the regional climate in the MEDiterranean region (ADRIMED) field campaign within the Chemistry-Aerosol Mediterranean Experiment (ChArMEx, http://charmex.lsce.ipsl.fr, last access: 25 January 2018) took place in the Mediterranean region from 11 June to 5 July 2013 . It aimed at characterizing the different aerosol particles and their radiative effects using airborne and ground-based measurements collected in the Mediterranean Basin, with special focus on the western region. In particular, two ChArMEx/ADRIMED flights, F30 and F31 from the French ATR 42 environmental research aircraft of SAFIRE (http://www.SAFIRE.fr, last access: 25 January 2018), took place above southeastern Spain during a Saharan dust episode on 16 and 17 June 2013.
In this paper, we present an analysis of the mineral dust radiative properties during this particular episode and take advantage of the thorough database that is available. Multiple datasets are used as input in a radiative transfer model (RTM) to evaluate the influence of the different measurements and data processing in the retrieved direct ARE. The model used here is the Global Atmospheric ModEl (GAME; Dubuisson et al., 1996, which allows for calculating both the solar and thermal infrared fluxes. An evaluation against aircraft in situ measurements of radiative fluxes is also presented. Two main goals are pursued: (i) the quantification of the direct ARE for two case studies within a dust transport episode and (ii) the evaluation of the model estimate sensitivity to the aerosol input used.
The paper is structured as follows: Sect. 2 includes a description of both the ground-based and in situ aircraft instrumentation and a short description of the retrieval algorithms used for the present study, Sect. 3 is devoted to the description of GAME and the input datasets used here, and results are presented in Sect. 4; finally, a short summary and concluding remarks are included in Sect. 5.

Ground-based measurements
Ground-based measurements used in this work were carried out at the Andalusian Institute for Earth System Research (IISTA-CEAMA) of the University of Granada, Spain (lat 37.16, long −3.61; 680 m a.s.l.) by the Atmospheric Physics Group of the University of Granada (GFAT-UGR). This experiment site is located in the western Mediterranean basin, near the African continent (∼ 200 km). Therefore, long-range transport of mineral dust particles from north Africa is a main source of natural atmospheric aerosol in the region (e.g., Lyamani et al., 2005;. The station is also affected by long-range transported smoke  and fresh smoke from nearby biomass burning . Anthropogenic sources such as pollution from Europe, the Iberian Peninsula, and the Mediterranean Sea (Pérez-Ramírez et al., 2016) also affect the station. Local sources are mainly road traffic and central heating systems . IISTA-CEAMA station is equipped with a CE-318-4 (Cimel Electronique) sun/sky photometer, which belongs to the AERONET network . This instrument performs direct solar irradiance measurements, used to derive AOD, and sky radiance measurements both measured at least at the following nominal wavelengths (λ): 440, 670, 870, and 1020 nm. The AOD product provided by AERONET has uncertainties of ±0.01 for λ > 440 nm and of ±0.02 for λ < 440 nm . AERONET also provides aerosol optical and microphysical properties such as columnar particle size distribution (PSD), real and imaginary parts of the refractive indices (RRI and IRI, respectively), asymmetry factor (g), and single scattering albedo (SSA) using the AOD and sky radiance values in an inversion algorithm . For the present study, AERONET Version 2 Level 1.5 (Level 2.0 when available) data are used. The uncertainty in the retrieval of SSA is ±0.03 for high aerosol load (AOD 440 > 0.4) and SZAs > 50 • ; while for measurements with low aerosol load (AOD 440 < 0.2), the retrieval accuracy of SSA drops down to 0.02-0.07 . For high aerosol load and SZAs > 50 • , errors are about 30 %-50 % for the IRI. For particles in the size range 0.1 < r < 7 µm (r being the aerosol radius), errors in PSD retrievals are around 10 %-35 %, while for sizes lower than 1 µm and higher than 7 µm retrieval errors rise up to 80 %-100 %. The inversion code provides additional variables such as the volume concentration; effective radius, r eff ; and geometric standard deviation of the equivalent lognormal distribution, σ , for fine and coarse modes of the retrieved PSD that will be used in the current study.
The multiwavelength aerosol Raman lidar MULHACEN, based on a customized version of LR331D400 (Raymetrics S.A.) is operated at Granada station as part of EAR-LINET/ACTRIS (European Aerosol Research Lidar Network/Aerosols, Clouds, and Trace Gases Research Infrastructure Network; https://www.actris.eu/default.aspx, last access: 15 June 2018; Pappalardo et al., 2014) since April 2005. The system has a monostatic biaxial configuration, which usually requires an overlap correction to minimize the incomplete overlap effect . The system emits vertically to the zenith by means of a pulsed Nd:YAG laser, with second-and third-harmonic generators, that emits simultaneously at 1064, 532, and 355 nm. The receiving system consists of several detectors, which can split the radiation according to the three elastic channels at 355, 532 (parallel-and perpendicular-polarized;, and at 1064 nm; two nitrogen Raman channels at 387 and 607 nm; and a water vapor Raman channel at 408 nm (Navas-Guzmán et al., 2014). The aerosol backscatter coefficient profiles (β aer (z, λ), z being the vertical height) obtained from the multiwavelength lidar were calculated with the Klett-Fernald method . For the retrieval of the aerosol extinction coefficient profiles (α aer (z, λ)), a height-independent lidar ratio (LR) obtained by forcing the vertical integration of α aer (z, λ) to the AOD from the AERONET photometer  was assumed. The assumption of a constant LRs introduces uncertainty in α aer (z, λ) retrievals, especially when different types of aerosol appear at different layers. In our case, the LR used for the Klett-Fernald retrieval are very similar to those provided by GRASP (see . Considering the different uncertainty sources, total uncertainty in the profiles obtained with the Klett-Fernald method is usually 20 % for β aer (z, λ) and 25 %-30 % for α aer (z, λ) profiles .
Additionally, surface temperature and pressure are continuously monitored at IISTA-CEAMA by a meteorological station located 2 m above the ground. At the same location, the global and diffuse downward radiative fluxes for the SW component are continuously measured with a CM11 pyranometer (Kipp & Zonen) and diffuse downward radiative fluxes for the LW component are measured with a precision infrared radiometer pyrgeometer (Eppley), both being instruments regularly calibrated at the site .

Airborne measurements
The SAFIRE ATR 42 aircraft performed two overpasses above Granada on 16 (flight F30) and 17 June (flight F31) in 2013 during the ChArMEx/ADRIMED campaign. During F30, the SAFIRE ATR 42 descended performing a spiral trajectory from 14:15 to 14:45 UTC; whereas during flight F31, the aircraft ascended in the early morning (from 07:15 to 07:45 UTC) at around 20 km from Granada station (see Fig. 1 from . Additional flight details can be found in previous studies Mallet et al., 2016;. The airborne instrumentation includes a scanning mobility particle sizer (SMPS) and an ultra-high sensitivity aerosol spectrometer (UHSAS) for measuring aerosol number size distribution in the submicron range. The forward-scattering spectrometer probe model 300 (FSSP-300) and the GRIMM optical particle counter (sky-OPC 1.129) were used to measure the optical size distributions in the diameter nominal size range between 0.28 and 20 µm and between 0.3 and 32 µm, respectively. A nephelometer (TSI Inc, model 3563) was used to measure the particle scattering coefficient at 450, 550, and 700 nm, and a cavity attenuated phase shift extinction monitor (CAPS-PMex, Aerodyne Inc.) was employed to obtain the aerosol extinction coefficient (α aer ) at 530 nm. For more details on the aircraft instrumentation see  and references therein. The PLASMA (Photomètre Léger Aéroporté pour la Surveillance des Masses d'Air) system, which is an airborne sun-tracking photometer, was additionally used to obtain AOD with wide spectral coverage (15 channels between 0.34 and 2.25 µm) with an accuracy of approximately 0.01, as well as the vertical profiles of the aerosol extinction coefficient .
Airborne radiative fluxes (F ) were measured with Kipp & Zonen CMP22 pyranometers and CGR4 pyrgeometers. Upward and downward SW fluxes ( ↑ F SW and ↓ F SW ) were measured in the spectral range 297-3100 nm by two instruments located above and below the aircraft fuselage. The same setup was used for the pyrgeometers, which provided the LW upward and downward radiative fluxes ( ↑ F LW and ↓ F LW ) for wavelengths larger than 4 µm. Both pyranometers and pyrgeometers were calibrated in January 2013 and data were corrected for the temperature dependence of the radiometer's sensitivity following .
Radiation measurement data from the aircraft were filtered out for large pitch and roll angles and corrected from the rapid variations in the solar incidence angle around the SZA due to the aircraft attitude (pitch and roll). This correction also depends on aircraft heading angle and solar position. It should be noted that, beforehand, roll and pitch offsets must be determined (the axis sensor is not necessarily vertical on average during a horizontal leg). Cosine errors were taken into account. Finally, data were corrected from variations in the SZA during the flight to ease the comparison with GAME retrievals. After these various corrections, an estimated uncertainty of ±5 W m −2 is considered to affect the data, taking into account the accuracy of the calibration and the acquisition system together with the consistency of airborne measurements (Meloni et al., 2018).

The GRASP code
The GRASP (Generalized Retrieval of Aerosol and Surface Properties) code  provides aerosol optical and microphysical properties in the atmosphere by combining the information from a variety of remote sensors (e.g., . In our case, GRASP was used to invert simultaneously coincident lidar data (range-corrected signal, RCS, at 355, 532, and 1064 nm) and sun/sky photometer measurements (AOD and sky radiances both from AERONET at 440, 675, 870, and 1020 nm) providing a detailed characterization of the aerosol properties, both column-integrated and vertically resolved. It is worthy to note that this GRASP scheme, based on Lopatin et al. (2013), presents the main advantage that it allows for retrieving aerosol optical and microphysical properties for two distinct aerosol modes, namely fine and coarse. The α aer , β aer , SSA (all at 355, 440, 532, 675, 870, 1020, and 1064 nm), and aerosol volume concentration (VC) profiles obtained as output from GRASP will be used as input to GAME in the present study, together with the column-integrated PSD properties (namely r eff and σ for fine and coarse modes). A more in-depth analysis of GRASP output data retrieved using the lidar and sun/sky photometer data at Granada station for the two inversions coinciding with the aircraft overpasses during flights F30 and F31 during ChArMEx/ADRIMED campaign can be found in .
3 GAME radiative transfer model

GAME description
The GAME code is widely described by Dubuisson et al. (2004 and . It is a modular RTM that allows for calculating upward and downward radiative fluxes at different vertical levels from the ground up to 20 km (100 km) in the SW (LW) spectral range. The solar and thermal infrared fluxes are calculated in two adjustable spectral ranges, which in this study were fixed to match those of the aircraft radiation measurements (namely 297-3100 nm for the SW and 4.5-40 µm for the LW) by using the discrete ordinates method . Note that the GAME code has a variable spectral sampling in the SW (depending on the spectral range considered) and a fixed spectral sampling (115 values) in the LW spectral range (Table 1).

GAME input data parameterization
The two considered SAFIRE ATR 42 flights, F30 and F31, took place on 16 and 17 June 2013, respectively, coinciding with ground-based lidar and sun/sky photometer measurements performed at the station. On these days, mineral dust with origin in the Sahara region (southern Morocco near the border with Algeria) reached Granada after ∼ 4 days of traveling, according to back-trajectory analysis (see Supplement Fig. S1) and the results presented in . A homogenous dust layer reaching up to 5 km a.g.l. was ob-Atmos. Chem. Phys., 19, 523-542, 2019 www.atmos-chem-phys.net/19/523/2019/ served on 16 June, whereas on 17 June the dust layer was decoupled from the boundary layer and located between 2 and 4.5 km a.g.l. . A very similar vertical structure was observed for the same dust event above Minorca . Daily maps of Meteosat Second Generation-derived AOD over the Mediterranean from 15 to 18 June during the dust event shown in Fig. 4 of  show the regional extension of the plume over the western Mediterranean region. On 16 June, the F30 flight above Granada site took place between 14:15 and 14:45 UTC coincident with the lidar measurements. The corresponding SZA at 14:30 UTC was 31.49 • . The sun/sky photometer microphysics data were not available until 16:22 UTC, even though the retrieved AOD and its spectral dependence (represented by the Ångström exponent) were very stable between the time of the lidar measurements and the time of the sun/sky photometer inversion. On 17 June, the F31 flight occurred in the early morning (07:15 to 07:45 UTC, with SZA = 61.93 • at 07:30 UTC), and simultaneous lidar and sun/sky photometer data were available. Unfortunately, the airborne vertical profile of extinction by the CAPS measurements was not available during this second flight. Clouds were detected by the lidar on 17 June after 15:00 UTC. Furthermore, a sky camera and the ground-based pyranometer and pyrgeometer data indicate cloud contamination (but not in the zenith) in the radiation data much earlier (around 09:00 UTC), also preventing satellite retrievals in the region. A summary of the experimental data used as input for GAME calculations during these two case studies is presented in Table 2. This input includes surface parameters and atmospheric profiles of meteorological variables, main gas concentrations, and aerosol properties. The aerosol properties used in the present study are parameterized using three different datasets, based on the different instrumentation and retrievals available, i.e., Dataset 1 (DS1), Dataset 2 (DS2) and Dataset 3 (DS3). A more detailed description of the different parameters is provided next.

Surface parameters and profiles of meteorological variables
The surface parameters required for GAME are the surface albedo (alb(λ)) and land-surface temperature (LST). The alb(λ) for the SW range is obtained from the sun/sky photometer data using the AERONET retrieval at 440, 675, 880, and 1020 nm, and for the LW from the integrated emissivity between 4 and 100 µm provided by the Single Scanner Footprint (SSF) Level2 products of the CERES (Clouds and the Earth's Radiant Energy System (http://ceres.larc.nasa. gov/, last access: 10 December 2018) instrument (Table 3). LST values are obtained from MODIS (Moderate Resolution Imaging Spectroradiometer) 1 km daily level-3 data (Wan, 2014) on 16 June. Unfortunately, on 17 June MODIS data were not available due to the presence of clouds and the local surface temperature was obtained from temperature measurements at the Granada site, where the meteorological station is located at 2 m above the ground. LST and alb(λ) values used for the two analyzed cases are included in Table 3. Table 2. Summary of the data sources used to obtain the input data parameterizations for GAME computations both in the SW and LW spectral ranges, including the surface parameters (albedo, alb, and land-surface temperature, LST), profiles of meteorological variables, and main gases and the aerosol parameters. For the aerosol parameters (aerosol extinction, α aer ; single scattering albedo, SSA; and asymmetry parameter, g) three different datasets are used (DS1, DS2, and DS3) based on different instrumentation and retrievals. The indications below the sources of the aerosol parameters indicate whether the parameter is column-integrated (col) or if it is vertically resolved (z) and the number of wavelengths at which it is given (nλ   Figure 1 shows the pressure (P ), temperature (T ), and relative humidity (RH) profiles obtained from the SAFIRE ATR 42 measurements. Data from the meteorological station located at IISTA-CEAMA are used to complete these profiles at the surface level; whereas at altitudes above the aircraft flight, a scaled US standard atmosphere is used for completion. The concentration profiles of the main absorbing gases (O 3 , CH 4 , N 2 O, CO, and CO 2 ) are also taken from the US standard atmosphere, while for the gaseous absorption coefficients the HITRAN database is used (as in Sicard et al., 2014a, b). Variations in the concentration profiles of the main absorbing gases have low impact of the radiative fluxes and the ARE, thus small uncertainty is introduced by this approach. A sensitivity test performed in the present study, varying the O 3 profiles up to double concentrations, indicates maximum differences of 4 W m −2 in the F SW and 3.6 W m −2 in the case of the F LW . For the ARE, differences are negligible (below 0.2 W m −2 ).

Aerosol parameterization
As for the aerosol parameterization, α aer (λ, z), SSA(λ, z), and g(λ, z) are required as GAME input data (Table 2). For the SW wavelengths, these properties can be obtained from the measurements performed with the instrumentation available during the campaign; namely the lidar, the sun/sky pho-tometer, and the in situ instrumentation on board the aircraft. On the other hand, direct measurements of the aerosol properties in the LW are not so straightforward and thus scarce. Hence, the aerosol LW radiative properties are calculated by a Mie code included as a module in GAME. According to , the dust particles nonsphericity effect at the thermal infrared wavelengths is not significant on the LW direct ARE, thus the shape of the mineral dust can be assumed as spherical for the Mie code retrievals introducing negligible uncertainties.
For the SW simulations, we run GAME using three different aerosol input datasets, i.e., DS1, DS2, and DS3 (Table 2), in order to evaluate their influence on the ARE calculations. DS1 relies on a parameterization based on the advanced postprocessing GRASP code, which combines lidar and sun/sky photometer data to retrieve aerosol optical and microphysical property profiles; DS2 relies on Klett-Fernald lidar inversions and AERONET products and corresponds to a reference parameterization (easily reproducible at any station equipped with a single-or multiwavelength lidar and an AERONET sun/sky photometer and without the need of an advanced postprocessing algorithm); and DS3 relies on in situ airborne measurements and corresponds to an alternative parameterization to DS1 and DS2.  For DS1 (Fig. 2a and d), α aer profiles at seven different wavelengths obtained with GRASP are used as input data in GAME. In DS2 (Fig. 2b and e), the α aer profiles are obtained from the lidar data using Klett-Fernald retrievals and adjusting the lidar ratio to the AERONET retrieved AODs, as mentioned in Sect. 2.1. Finally, for DS3 ( Fig. 2c and f) the α aer values are obtained from the aircraft in situ measurements (CAPS and PLASMA data on 16 June and PLASMA on 17 June). A detailed analysis and discussion on the comparison between α aer profiles provided by the aircraft measurements, GRASP, and the lidar system at Granada is already included in . In general, the lidar, GRASP, and the CAPS data are in accordance, observing the same aerosol layers and similar values, with differences within 20 %. GRASP slightly overestimates CAPS data by 3 Mm −1 on average, whereas the differences with PLASMA are larger, reaching 30 % (or 11 Mm −1 ). In the case of the Klett-Fernald retrieval, values are lower than those retrieved with GRASP by up to 19 %. Considering that the uncertainty in α aer is around 30 % for both GRASP and the Klett-Fernald retrieval and 3 % for the CAPS data, this discrepancy is well below the combined uncertainty in the different datasets. Dif-ferences in the α aer profiles translate into differences in the integrated extinction and, hence, into differences in the AOD values used as input in the radiative flux retrievals. The AOD values presented here (included in Table 4) are obtained by integrating the α aer profiles at 550 nm from the surface up to the considered top of the aerosol layer (4.3 km on 16 June and 4.7 km on 17 June). In GRASP retrieved α aer profiles, values above this top of the aerosol layer are slightly larger than zero since GRASP takes into account stratospheric aerosols by an exponential decay (Lopatin et al., 2013), thus the approach used here to calculate the AOD leads to lower values compared to the column-integrated AOD provided by the sun photometer. Differences among the three datasets are more noticeable on 16 June, when the AOD for DS1 is 0.05 lower than for DS2 and DS3; whereas on 17 June the maximum difference is 0.03, obtained between DS1 and DS2. The AOD values at 550 nm reveal that GRASP input data (DS1) and to a lesser extent the aircraft in situ data (DS3) underestimate the aerosol load in the analyzed dust layer compared to AERONET (DS2) due to the differences in the retrieval techniques, e.g., although AERONET provides integrated AOD for the whole column, low α aer values above the aerosol layer are neglected for the AOD calculations in DS1 and DS3. Table 4. Column-integrated number concentration (N), effective radii (r eff ), and standard deviation (σ ) of fine and coarse aerosol modes and AOD at 550 nm for DS1, DS2, and DS3 on 16 and 17 June.
16 June (SZA = 31.49 • ) N f (no. µm −2 ) N c (no. µm −2 ) r eff,f (µm) r eff,c (µm) σ f (µm) σ c (µm) AOD (550 nm Figure 3 presents the SSA values retrieved by the GRASP algorithm, used as input for GAME in DS1, on 16 (F30, Fig. 3a) and 17 June (F31, Fig. 3b). The mean SSA at 440 nm is equal to 0.92 on 15 June, whereas on 17 June it is 0.85. On 17 June the SSA profiles present lower values and more variation with height than on 16 June; the lower SSA values indicate the presence of more absorbing particles on 17 June. The vertical variation on 17 June is associated with the presence of two different layers, whereas a more homogeneous dust layer is observed on 16 June. For DS2, the SSA are taken from AERONET columnar values and assumed to be constant with height (Fig. 4a). The SSA at 440 nm was 0.89 and 0.83 on 16 and 17 June, respectively; as already observed in Fig. 3, SSA values are lower on 17 June due to the intrusion of more absorbing particles. For DS3, SSA values at 530 nm are obtained from the nephelometer and the CAPS or PLASMA on board the ATR. In order to reduce the uncertainty in the measured data, only averaged values for the column will be considered, being 0.88 and 0.83 on 16 and 17 June, respectively (Fig. 4). Therefore, differences of up to 0.04 and 0.02 are observed on 16 and 17 June, respectively, among the SSA values obtained with the three datasets. Despite these differences, the retrieved SSA values obtained here are within the range of typical values for dust aerosols Lopatin et al., 2013) and differences are still within the uncertainty limits, which range between 0.02 and 0.07 depending on the aerosol load for AERONET data  and is 0.04 for the aircraft values. In the case of g values, the same data are used for the three aerosol input datasets. Multispectral values of g are taken from AERONET columnar values and assumed to be constant with height (Fig. 4b).
Summing up, for the SW aerosol parametrization in GAME three datasets are tested. In DS1, GRASP-derived spectral profiles at seven wavelengths of the aerosol extinction and SSA are used. In DS2, the Klett retrieved extinction profiles at three wavelengths are used together with the AERONET SSA columnar values at four wavelengths, which are assumed to be constant with height. For DS3, one extinction profile at 550 nm and a column-averaged single-wavelength value of the SSA from the airborne measurements are considered. In the three cases, the columnintegrated AERONET asymmetry parameter at four wavelengths is assumed to be constant with height and used as input.
For the LW calculations, the Mie code is used to obtain α aer (λ, z), SSA(λ, z), and g(λ, z) from the information on the aerosol PSD, complex refractive index (RI), and density, following a similar approach to that used in previous studies (Meloni et al., , 2018Peris-Ferrús et al., 2017). A summary of the aerosol parameters used in the Mie calculations is included in Table 5. Three different datasets are also used for the aerosol parameterization in the LW calculations. In this case, the sensitivity of the model to the PSD used is tested. A similar scheme to that presented for the SW is used, where DS1 relies on GRASP retrievals, DS2 on AERONET products, and DS3 relies on in situ airborne measurements.   Table 5. Summary of the data used to obtain α aer (λ, z), SSA(λ, z), and g(λ, z) in the LW spectral range from Mie calculations, i.e., the refractive index, RI; effective radius, r eff ; geometric standard deviation, σ ; and number concentration, N. Three different datasets are used (DS1, DS2, and DS3) based on different particle size distribution (PSD) data used. The indications below the sources of the aerosol parameters indicate whether the parameter is column-integrated (col) or if it is vertically resolved (z) and the number of wavelengths at which it is given (nλ). DB (2017) The spectral real and imaginary parts of the RI of mineral dust in the LW are obtained from Di , using the Morocco source, and assumed constant with height. The analysis by Di  only covers the spectral range 3-16 µm so an extrapolation assuming the spectral dependence presented in  for shorter and longer wavelengths is performed. This assumption is not exempt from uncertainty, since the refractive index presents a certain variability associated with the different nature of mineral dust properties. For example, the use of the refractive index provided for the Algerian and Mauritanian sources from Di  leads to variations in the ARE of 0.8 and 0.3 W m −2 at the BOA and the TOA, respectively. Additionally, vertical variations in the refractive index are also a source of uncertainty in the obtained radiative fluxes. The mineral dust particle density is assumed to be 2.6 g cm −3 . Regarding the PSD, three parameters (namely the effective radii, r eff ; standard deviation, σ ; and the numeric concentrations, N ) for fine and coarse modes are used. The fine mode comprises particles within the diameter range 0.1-1 µm, whereas for the coarse mode the range 1-30 µm is considered. A third mode at about 30 µm was detected in  for the same dust event using balloon-borne measurements with concentrations up to 10 −4 particles cm −3 . However, this giant mode is not considered in our study due to the lack of data above Granada. Con-sidering the relevance of large particles for the ARE LW (i.e., Sicard et al., 2014a, b;Meloni et al., 2018), neglecting this giant mode may contribute to increase the uncertainties in GAME estimations. However, simulations with GAME assuming the presence of a third mode of similar characteristics to the one observed by  indicate that variations in the ARE are negligible in this case (lower than 0.1 W m −2 ). Even for much higher concentrations (10 −1 particles cm −3 ), variations in the ARE of just 0.3 W m −2 at the BOA and 0.15 W m −2 at the TOA are obtained.
In the case of DS1, N values are obtained from the volume concentration profiles provided by GRASP assuming spherical particles in the range between 0.05 and 15 µm radii (Fig. 5). Values of r eff and σ provided by GRASP (Table 4) are column-integrated and thus assumed to be constant with height. This is also the case for DS2, in which the PSD parameters are column-integrated values provided by the AERONET retrieval in Granada (see Table 4).
For DS3, the volume concentration (or the equivalent N ), r, and σ profiles for the fine and coarse modes (Fig. 5) are calculated from the data provided by the aircraft in situ measurements in the range between 0.02 and 40 µm diameter.  found a general good accordance between the volume concentration profiles measured by the instrumentation on board the SAFIRE ATR 42 and retrieved with GRASP, with differences in the total volume concentration profiles for the dust layers lower than 8 µm −3 cm −2 (20 %), which fall within the combined uncertainty. Nonetheless differences are still noticeable, especially in the fine mode. On 17 June, GRASP overestimates the aircraft measurements for the fine mode and underestimates them for the coarse mode, which in turn results in a quite different fine to coarse concentration ratio for DS1 and DS3. Additionally, a slight shift is observed in the vertical structure of the aerosol layers. Differences are mostly technical, i.e., the GRASP retrieval is based on 30 min averaged lidar profiles while the aircraft provide instantaneous measurements, but they can also be partially caused by the discrepancies between the vertical aerosol distribution above Granada (sampled by the lidar) and the concentration measured during the aircraft trajectory as they are not exactly coincident. In addition, for 16 June, there is a 2 h time difference between the sun/sky photometer retrieval used in GRASP calculations and the airborne measurements, which can lead to slight differences in the aerosol properties despite the homogeneity of the dust event during this period. In the following, we quantify the impact these differences may introduce in the calculations of F .

GAME output data
As a result of the simulation, GAME provides vertical profiles of radiative fluxes in the shortwave (F SW ) and longwave (F LW ) spectral ranges. The net flux can be calculated from the obtained profiles for both spectral ranges as where the upward and downward arrows are for upward and downward fluxes, respectively. From the obtained radiative flux profiles, the direct ARE profiles are calculated according to the following equation: where F w and F o are the radiative fluxes with and without aerosols, respectively. The direct ARE can be obtained for the SW (ARE SW ) and the LW (ARE LW ) spectral ranges.
4 Mineral dust effect on shortwave and longwave radiation 4.1 SW radiative fluxes Figure 6 shows the radiative flux profiles for the SW spectral range obtained with GAME using the three different input datasets described in Sect. 3, as well as the Net F SW . The radiative fluxes measured by the pyranometer on board the SAFIRE ATR 42 are also included in the figure. The three GAME simulations show similar values with differences below 8 W m −2 on average, which represents less than 1 % variation. The differences in the obtained fluxes are mostly due to the differences in the aerosol load considered depending on the inputs. Even though the differences in the AOD among the different datasets are small (lower than 0.05), they can lead to differences in F SW and ultimately in the ARE SW . In order to quantify these differences, we performed a sensitivity test by varying the AOD while the other parameters were kept constant. We observed a maximum variation in the F SW of 6.5 W m −2 (0.7 %) at the surface, decreasing with height, for changes in the AOD of up to 0.05, which is the difference we observe between the AOD for DS2 and DS1 on 16 June. This result partly explains the differences among the three datasets. In addition, a sensitivity test performed by exclusively varying the SSA indicates that more absorbing particles are related to less ↓ F SW at the surface, namely a variation of 1 % is observed at the BOA for a decrease in the SSA of 0.03. The influence of the SSA decreases with height and is negligible at the TOA. For the ↑ F SW , a decrease of 0.8 % is observed at the BOA if more absorbing particles are present, but in this case the influence at the TOA is larger (2.2 %). In our case, the larger AOD assumed for DS2 on both days (see Table 4 and Fig. 2), causes the ↓ F SW to be slightly lower compared to DS1. For DS3 the AOD is similar to DS2, but the SSA values used, which are relatively smaller compared to those measured by AERONET (see Fig. 4), lead to lower values of the radiative fluxes than for DS2. The vertical distribution of the SSA also influences the radiative fluxes in the SW component as demonstrated in previous studies , contributing to explain the differences observed among the three datasets analyzed here. The evaluation against the aircraft measurements shows larger differences for altitudes below 2.5 km (∼ 860 mbar) on 16 June, whereas a better agreement is found above. On 17 June, no ↑ F SW aircraft data are available below 2 km. Relative differences between the model and the aircraft measured data (calculated as (F GAME −F aircraft )/F aircraft ) are well below 7 %, being the largest discrepancies observed for the ↓ F SW . Differences among the three GAME outputs and the aircraft pyranometer are lower than 5 % for the Net F SW on both days. Considering the very different approaches followed by the model and the direct measurements by the airborne pyranometer (i.e., vertical resolution, temporal sampling, and data acquisition and processing) -together with the uncertainty in the pyranometer (5 W m −2 ) and the estimated uncertainty in the model outputs that can be as large as 19 W m −2 ) -these differences are quite reasonable. A conclusive result on which input dataset provides a better performance is unlikely because of the comparable results obtained with the three datasets.
The values at the surface (or BOA) and at the TOA for the different radiative fluxes can also be evaluated against different instruments: measurements for the ↓ F SW at the surface are available from the sun photometer; AERONET provides values for the ↓ F SW and ↑ F SW at both the BOA and TOA. The time series for these measurements corresponding to 16-17 June and the results obtained with GAME for the different datasets are shown in Fig. 7. AERONET surface radiative fluxes have been extensively validated at several sites around the world (e.g.,  and, in addition, all AERONET sun photometers are mandatorily calibrated once a year. Thus, in order to compare GAME results with AERONET data, we have performed additional simulations for the time of the closest AERONET measurement on 16 June (at 16:22 UTC), assuming that the aerosol parameterization is constant with time between the flight time and the photometer measurement. The ↓ F SW values at the surface obtained with GAME are 564.8, 551.8, and 547.0 W m −2 for DS1, DS2, and DS3, respectively, close to the 531.4 W m −2 provided by AERONET. On 17 June, GAME simulations at 07:40 UTC (instead of 07:30 UTC, which is the time of the flight) provide ↓ F SW at the surface of 466.3, 468.3, and 456.4 W m −2 , very similar to the AERONET value of 463.7 W m −2 .
At the TOA, the ↑ F SW between GAME and AERONET are in quite good agreement on both days. On 16 June, the ↑ F SW values obtained with GAME simulations are equal to 152.0, 153.0, and 148.5 W m −2 , and with AERONET they are equal to 146.2 W m −2 . On 17 June, the obtained values with GAME are 133.6, 136.6, and 130.9 W m −2 for DS1, DS2, and DS3, and 131.6 W m −2 for AERONET.
The ARE SW profiles, calculated by using Eq.
(2) and GAME simulations for the three input datasets, are shown in Fig. 8, together with the simultaneous values provided by AERONET on 17 June at the BOA and TOA. Comparing the three GAME simulations, we can see that the low discrepancies in the F profiles from Fig. 6 lead to variations in the ARE SW of 10 %-27 % (3-10 W m −2 ) over the averaged profile depending on the input dataset used. The variations in the ARE SW values are tightly connected to differences in the AOD considered as input in the model, as already observed in previous studies Lolli et al., 2018;Meloni et al., 2018). The SSA and the vertical distribution of the aerosol also plays an important role, as observed for DS3, which shows a profile that is quite different from DS2 despite the AOD being quite close for both datasets.
Differences are also observed when comparing ARE SW values obtained from GAME to those retrieved by AERONET. Contrary to GAME simulations, AERONET does not consider the vertical distribution of the aerosols when calculating the ARE SW , and the definition of the ARE SW at the BOA ( BOA ARE SW ) is slightly different. Indeed, AERONET BOA ARE SW is calculated as the difference between the downward fluxes with and without aerosols, the difference between the upward fluxes (reflected by the Earth) being neglected. Considering this, we can correct the BOA ARE SW provided by AERONET by multiplying by a factor 1−alb(λ). The corrected BOA ARE SW value on 17 June is thus −31.9 W m −2 , which is within the range of values provided by GAME at the surface. All discrepancies observed here are mostly intrinsic to the different techniques used for the acquisition of the data and the retrieval algorithms. The effect of the data processing has also been observed in pre- vious studies (Lolli et al., 2018). Moreover, the sensitivity tests performed reveal that an increase in the AOD of 0.05 can lead to a stronger effect of the ARE both at the BOA (up to 6.7 W m −2 ) and the TOA (up to 2.5 W m −2 ), and more absorbing particles (decrease in the SSA of 0.03) lead to more ARE at the BOA and less at the TOA (4 and 2 W m −2 in absolute terms, respectively). Therefore, the differences among the datasets are within the estimated uncertainty.
The ARE SW values obtained at the BOA and TOA for the three datasets and the averaged value, as well as the FE, are included in Table 6. Both at the BOA and TOA, the ARE SW has a cooling effect, as expected for mineral dust in this region according to values obtained in the literature (e.g., Mallet et al., 2016). Differences among the three datasets lead to variations in the ARE SW of up to 30 % (or 20 % for the FE), observing larger variability on 16 June. The values of the ARE SW and the FE SW are highly dependent on the SZA and a straightforward comparison with previous studies is not simple. Nonetheless, the values obtained for this case are within the range of previous values observed in the western Mediterranean region for similar values of SZA, e.g., FE between −263.4 and −157.1 W m −2 at the BOA and −23.8 and −86.2 W m −2 for SZA = 60 • or ARE values ranging between −93.1 and −0.5 W m −2 at the BOA and between −34.5 and +8.5 W m −2 at the TOA for different SZA values (e.g., Sicard et al., 2014a, b; Figure 9 shows F LW calculated with GAME after obtaining the aerosol properties in the LW spectral range from Mie calculations for the three mentioned datasets (see Sect. 3.2.2). F LW values measured by pyrgeometers located on board the ATR are also shown.

LW radiative fluxes
In general, differences in the F LW are always lower than 6 % (lower than 10 W m −2 on average), with the airborne values being overestimated by the model on 16 June and underestimated on 17 June. On this latter day, larger differences are observed on the Net F LW compared to 16 June, which might be explained by the inaccurate value of LST used due to the lack of precise data. A sensitivity test performed by increasing the air surface temperature measured at the meteorological station by 5 K indicates that ↑ F LW increases its value up to 30 W m −2 at the surface, and around 10 W m −2 from 1 km onwards, which is nonnegligible. This would lead to an overestimation of the aircraft measured values, but still within a 6 % difference. This highlights the need for accurate LST measurements for radiation simulations in the LW spectral range. Additionally, a sensitivity test performed by assuming a 10 % uncertainty in the PSD parameters (r eff , N , and σ ) leads to an estimated uncertainty in the F LW retrieved by GAME of around 1.2 W m −2 . As stated before, the assumption of the refractive index can also introduce variations as large as 0.8 W m −2 . Considering the uncertainty in the pyrgeometer, the fact that the aircraft and the model present dif- Figure 10. Time series of the ↓ F LW at the BOA during the period 16-17 June. Surface measurements of diffuse (red) radiation from the ground-based pyranometer at the Granada station are included. GAME output data for different inputs are represented by the blue circles (DS1), red (DS2), and green (DS3) crosses. ferent vertical resolutions and time samplings, and the uncertainties due to the use of the standard atmosphere or the parameterization of the surface properties, the obtained differences are not significant.
A comparison of GAME results against the observations from the ground-based pyrgeometer at the Granada station is included in Fig. 10. At the BOA, the longwave radiation measured by the pyrgeometer is in quite good agreement with GAME calculations on 16 June, with differences within 1 W m −2 . However, GAME overestimates the pyrgeometer data by 5 W m −2 (1.3 %) on 17 June. This difference on 17 June, even though larger than on 16 June, is still within the uncertainty limits.
As for the ARE LW , Fig. 11 shows the profiles obtained with GAME using the three datasets as inputs. Values at the BOA and TOA for each dataset and the average values are included in Table 7, together with the FE. Opposite to the SW, the ARE LW produces a heating effect both at the BOA and TOA with positive values. The slight differences in the F LW in Fig. 9 due to the use of different aerosol input datasets lead to variations of up to 2 W m −2 in the BOA ARE LW (ranging from 20 % to 26 %), which needs to be considered in the interpretation of the results and reduced for a better estimate of the direct ARE. Despite this, values obtained for this dust event (3.2 W m −2 on average for both days) are in agreement with previous studies performed for mineral dust in the infrared region (Sicard et al., 2014a, b) and the FE obtained are comparable to those reported by Meloni et al. (2018). It is extremely interesting to look at the differences between the two days in terms of AOD ( AOD) and the effective radius for the coarse mode, r eff,c , ( r eff,c ) and their implication for the differences in the ARE LW at the BOA ( BOA ARE LW ). For DS1, AOD ( r eff,c ) is −0.02 (+0.18 µm), which produces a decrease in BOA ARE LW ( BOA ARE LW = −0.5 W m −2 ). For Figure 11. Direct ARE profiles in the LW spectral range simulated using DS1 (blue line), DS2 (red line), and DS3 (green line) as aerosol input data in GAME for 16 June (a) and 17 June (b).
DS2, AOD ( r eff,c ) is −0.04 (+0.18 µm), which produces a decrease in BOA ARE LW ( BOA ARE LW = −1.0 W m −2 ). If we relate these variations to the sensitivity study of , the expected ARE LW increase due to the increase in the coarse mode radii is counterbalanced by the ARE LW decrease when AOD decreases in both cases. Oppositely, for DS3 the AOD ( r c ) is −0.05 (+0.64 µm), producing an increase in BOA ARE LW ( BOA ARE LW = +1.6 W m −2 ). Here, the large increase in the coarse mode radius dominates over the AOD decrease.  indeed show that the largest positive gradient of ARE LW occurs for median radii ranging from 0.1 to 2.0 µm. For DS3 the increase in BOA ARE LW produced by a positive r c is larger than the decrease in BOA ARE LW that would have produced AOD alone. At the TOA, same trends, but much less marked, are observed.

Conclusions
A moderate Saharan dust event affecting the western Mediterranean region during the ChArMEx/ADRIMED campaign on June 2013 was extensively monitored by ground-based and aircraft instrumentation above the Granada experimental site. Radiative fluxes and mineral dust ARE both in the solar and infrared spectral ranges are calculated for this event with the RTM GAME. Three different aerosol input datasets are used by the GAME RTM in order to evaluate the impact of different input data in GAME calculations.
For the SW, very low variability with the input aerosol data (less than 1 %) is observed for the radiative fluxes. The evaluation of GAME-calculated radiative fluxes against the aircraft data reveals differences between the model fluxes and the measurements below 7 %, with better agreement at altitudes above the planetary boundary layer. The differences between the retrievals with the three aerosol datasets are quite insignificant, especially taking into account the different approaches followed by the model and the pyranometers and the estimated uncertainties for both the measured data (5 W m −2 ) and the model (around 12 W m −2 ). Thus a conclusion on which input dataset provides a better performance is unlikely. The small differences between GAME radiative flux retrievals lead to variations in the ARE SW of up to 33 %, mostly driven by the differences in the aerosol vertical distribution and load, followed by the SSA.
For the LW component, the effect of aerosol on the radiative properties is lower compared to the SW, but certainly nonnegligible and of opposite sign. GAME retrievals using the three aerosol datasets reveal differences in the fluxes lower than 2 W m −2 (less than 1 %). The comparison with the pyrgeometer data measured at the ATR reveals, however, differences around 7 %. The influence of the assumed CO 2 , O 3 , and the water vapor profiles and LST are needed to fully explain this discrepancy between the aircraft and the simulated profiles.
The total ARE, including both the SW and LW components, confirms that mineral dust produces a cooling effect both at the surface and the TOA, as already reported in the literature. On average, the ARE LW represents 20 % of the ARE SW at the surface, therefore clearly indicating that global model estimates need to consider the complete spectrum to avoid an overestimation of the mineral dust cooling effect. Additionally, it is necessary to be aware of the effects of using different measurement techniques and processing methodologies when calculating aerosol radiative properties. Even though the differences observed here when using different aerosol datasets are slight, they still exist and a homogenization of the techniques to feed global models would be beneficial for a better estimate of the ARE and a reduced uncertainty.
Data availability. Part of the data used in this publication were obtained as part of the AERONET and EARLINET networks and are publicly available. For additional data or information please contact the authors.
Author contributions. MJGM and MS designed the study and wrote the manuscript with contributions from all authors. MJGM, RR, JABO, CD, GB, PF, and BT provided data and performed data analysis. RB provided the initial version of the model. All authors have given approval to the final version of the manuscript.
Competing interests. The authors declare that they have no conflict of interest.
Special issue statement. This article is part of the special issue "CHemistry and AeRosols Mediterranean EXperiments (ChArMEx) (ACP/AMT inter-journal SI)". It is not associated with a conference.