Determination and analysis of spectral aerosol optical properties by a multi-instrumental approach

Continuous in-situ measurements of aerosol optical properties where conducted from 20 June to 20 July in Granada (Spain) with a 7-wavelength Aethalometer, a Multi Angle Absorption Photometer, and a 3-wavelength integrating Nephelometer. The aim of this work is to describe a methodology to obtain the absorption coe ﬃ cients ( b abs ) for the 5 di ﬀ erent Aethalometer wavelengths. In this way, data have been compensated using algorithms which best estimate the compensation factors needed. Two empirical factors are used to infer the absorption coe ﬃ cients from the Aethalometer measurements: C – the parameter describing the enhancement of absorption by particles in the ﬁlter matrix due to multiple scattering of light in the ﬁlter matrix; and f – the parameter com- 10 pensating for non-linear loading e ﬀ ects in the ﬁlter matrix. Spectral dependence of f found in this study is not very strong. Values for the campaign lie in the range from 1.15 at 370 nm to 1.11 at 950 nm. Wavelength dependence in C proves to be more important, and also more di ﬃ cult to calculate. The values obtained span from 3.40 at 370 nm to 4.35 at 950 nm. Furthermore, the temporal evolution of the Ångström exponent of 15 absorption ( α abs ) and the single scattering albedo ( ω 0 ), is presented. On average α abs is around 1.1 ± 0.3, and ω 0 is 0.78 ± 0.08 and 0.74 ± 0.09 at 370 and 950 nm, respec-tively. These are typical values for sites with a predominance of absorbing particles, and the urban measurement site in this study is such. The b abs average values are of 16 ± 10 Mm − 1 (at 370 nm) and 5 ± 3 Mm − 1 (at 950 nm), respectively. Finally, di ﬀ erences 20 between working days and Sunday have been further analyzed, obtaining higher b abs and lower ω 0 during week than on Sundays as a consequence of the inﬂuence of diesel tra ﬃ c.


Introduction
The radiative forcing in the Earth's atmosphere caused by aerosols is highly uncertain (IPCC, 2013).The direct effect of aerosol influence is exhibited by scattering and/or absorbing solar radiation.By scattering the light they contribute to the cooling of the underlying atmosphere (negative radiative forcing).Furthermore, by absorbing light they contribute to the heating of the atmosphere (positive radiative forcing).The latter is an important component of the radiation budget of the atmosphere which still needs to be better characterized to reduce uncertainties in climate models (Houghton et al., 2001;Ramanathan et al., 2001;Kirchstetter et al., 2004).
The major light-absorbing compounds among aerosols are carbonaceous substances and mineral dust.However, the spectral dependence of aerosol absorption coefficient varies depending on their chemical and physical properties.Light absorption by black carbon (BC) is generally considered to vary weakly with wavelength (Bergstrom et al., 2007).Nevertheless, other aerosol components such as desert aerosols (Dubovik et al., 2002), organic compounds (Jacobson, 1998), or biomass burning (Mukai and Ambe, 1986), exhibit a stronger wavelength dependence.Therefore, the spectral dependence of aerosol absorption coefficient is important to distinguish different aerosol absorbing components (Collaud-Coen et. al, 2004;Fialho et al., 2005;Sandradewi et al., 2008a).
Different methods and instruments are used to measure light absorption (Hansen et al., 1984;Horvath, 1993;Bond and Bergstrom, 2006;Moosmüller et al., 2009;Müller et al., 2011).The most widely used method for in situ measurements is the filter-based technique, which relies on determining the absorption coefficient (b abs ) by measuring the change of light attenuation due by aerosol collected on a filter matrix (Hansen et al., 1982).Instruments currently used for determining the aerosol absorption coefficient (b abs ) using filter-based techniques are the Aethalometer (Hansen et al., 1984), the Particle Soot Absorption Photometer (PSAP) (Bond et al., 1999) and the Multi-Angle Absorption Photometer (MAAP) (Petzold and Schonlinner, 2004).
It is known that all filter-based instruments suffer from several artifacts (Schmid et al., 2006;Müller et al., 2011;Collaud-Coen et al., 2010).Therefore, it is necessary to apply compensation to the data to determine the absorption coefficient.Several studies have proposed different numerical methods to correct these artifacts (Bond et al., 1999 Weingartner et al., 2003;Arnott et al., 2005;Schmid et al., 2006), and have proven to give quite satisfactory results when compared to other non-filter based instruments.
From the instruments mentioned above, the MAAP is the newest filter-based instrument for measuring the aerosol absorption coefficient (Petzold and Schönlinner, 2004;Petzold et al., 2005).Measurements provided by the MAAP are only available for one wavelength, which is a disadvantage since it is important to determine both b abs and its spectral dependence.This information can be obtained by combining the MAAP measurements with those performed by a multi-wavelength Aethalometer AE-31, which performs measurements at 7 different channels covering the range from the UV (370 nm) to NIR (950 nm).
The aim of this work is to obtain the spectral variation of the compensation factors of the Aethalometer data and then apply them to the measurements performed at the different channels of the Aethalometer.In this way, the aerosol absorption coefficients at all the different Aethalometer channels are calculated and then, other optical parameters, such as the spectral single scattering albedo (ω 0 ) or the Ångström exponent of absorption (α abs ), are determined and analyzed.
The data set presented in this study was measured during one month, from 29 June 2012 to 29 July 2012, in Granada (Spain) using a MAAP, an Aethalometer, and a Nephelometer.In this work we will discuss the methodology used for compensating the Aethalometer and MAAP data.Finally, the temporal evolution of the results obtained is presented and compared with results in other works to evaluate the reliability of the data.

Measurement site
Measurements presented in this study were performed in Granada, from 29 June to 29 July 2012.Granada (37.18 • N, 3.58 • W, 680 m a.s.l.) is a non-industrialized medium-size city located in southeastern Spain with a population of around 500 000 inhabitants considering the whole metropolitan area (http://www.juntadeandalucia.es/).Near continental conditions prevail at this site and are responsible for large temperature differences, providing cool winters and hot summers (Lyamani et al., 2010).
The measurement station is located in the southern part of the city, less than 500 m away from a highway that surrounds the city.Local aerosol sources are mainly road traffic (dominated by diesel engines) together with soil re-suspension, especially during the warm-dry season when the reduced rainfall may increase the contribution of local mineral dust.Due to its location in the Mediterranean basin, it is influenced by two major aerosol source regions: Europe, as a source of anthropogenic pollutants, and North Africa as a source of natural mineral dust (Lyamani et al., 2010).

Instrumentation
Sampling for all the different instruments was obtained using a stainless steel tube with 20 cm diameter and 5 m length (Lyamany et al., 2008).The inlet was located about 15 m above the ground.Measurements were carried out without an aerosol size cut-off or heating of the sampled air.From the tube, several stainless steel pipes led the sampled air to each instrument at the appropriate flows.Different diameters were adjusted to maintain the laminar flow in the tubes and minimize particle losses (Baron and Willeke, 2005).
The aerosol light absorption coefficient was measured using two different filter-based instruments.The MAAP (Thermo Scientific) measures the light transmitted through and backscattered from a particle-loaded filter.The b abs, MAAP at 637 nm (Müller et al., 2011) is calculated using radiative transfer model which includes a treatment of the scattering effects of the filter matrix and the light scattered by the aerosol component.A detailed description of the method is provided by Petzold and Schönlinner (2004).The MAAP works at a constant flow rate of 16.7 L min −1 and provides measurements every minute.
The total method uncertainty for the aerosol light absorption coefficient inferred from MAAP measurement is around 12 % (Petzold and Schönlinner 2004;Petzold et al., 2005).
The MAAP was believed to be the instrument which is the least affected by artifacts for obtaining the aerosol absorption coefficient among all the different filter-based methods.Different studies showed that b abs, MAAP measured by MAAP is in good agreement with those measured by photoacoustic spectrometry (e.g.Petzold et al., 2005;Sheridan et al., 2005).However, a recent study by Hyvärinen et al. (2013) showed that MAAP measurements suffer from some artifacts in locations with high concentrations of light absorbing particles.Therefore, to avoid these artifacts, MAAP data have been compensated using the method described in Hyvärinen et al. (2013).A more detailed explanation can be found in Sect.3.1.The other instrument used to measure the b abs is an Aethalometer model AE-31-ER (Magee Scientific).The AE-31 measures light attenuation at 7 different wavelengths covering the UV (370 nm), visible (470, 520, 590, and 660 nm), and NIR (880 nm and 950 nm) ranges.A complete description of the operating principles of Aethalometers can be found in Hansen et al. (2005).
The Aethalometer measures the light attenuation through a quartz filter matrix as aerosols are deposited on the filter.This parameter is defined by: ATN(λ) = − ln where I is the intensity of light that passes though the loaded filter, and I 0 is the intensity of light passing through the unloaded part of the filter.The attenuation coefficient at each wavelength, b ATN (λ) can be obtained by: where A is the filter spot area (1.67 cm 2 ), V the volumetric flow rate, and ∆ATN is the variation in the attenuation measured during the time interval (∆t).The attenuation coefficient is determined from the change of attenuation and is therefore due entirely to the loading of the filter by aerosols.In this work the instrument was operated with a flow rate of 4 L min −1 and a sampling interval of 5 min.The compensation algorithm to obtain the absorption coefficient b abs from b ATN is described in more detail in Sect.3.2.
Aerosol scattering (b scat ) and backscattering (b bsp ) coefficients at ground level were measured using a TSI Model 3563 three-wavelength (450, 550, 700 nm) integrating Nephelometer (Anderson et al., 1996;Heintzenberg et al., 2006).This instrument draws the ambient air through a temperature-monitored inlet at a flow rate of 30 L min −1 , illuminates the sample with a halogen lamp, and measures the scattered light using three photomultiplier tubes.The scattered light is integrated over an angular range which can be adjusted to either 7-170 • or 90-170 • , using a backscatter shutter to give total scatter and backscatter signals.Pressure and temperature are measured in the sample chamber so scattering by air molecules can be calculated and then subtracted from the total scattering to determine scattering by aerosol particles.Relative humidity is also measured inside the Nephelometer and it was 28 ± 10 % during the study period.In this study, non-idealities due to truncation errors were corrected (Anderson and Ogren, 1998).Uncertainties in the nephelometer measurements consist of angular truncation errors, non-lambertian nature of the light source, and wavelength non-idealities.The calibration uncertainty is approximately 7 % (Anderson et al., 1996;Anderson and Ogren, 1998;Esteve et al., 2012).

Data compensation and methods
All filter absorption photometers are affected by loading effects which need to be compensated and the various compensation methods relied mostly on the continuity of data (Weingartner et al., 2003;Virkkula et al., 2007;Hyvärinen et al., 2013) to derive the compensation parameters.Here we propose a more robust method to determine the quality of the compensation for loading effects: the slope of the absorption coefficient vs. the loading of the spot.This method was used to compensate the BC measurements (Park et al., 2010), but here we use it as an independent criterion to evaluate the compensation.This criterion does not use any assumptions or external data, only the measurements from a single instrument -it checks the measured parameters for internal consistency.For a long period of homogeneous measurements in terms of the sources and meteorology, the absorption coefficient should not depend on the loading of the spot.The uncompensated measurements of Black Carbon (BC) in filter photometers are proportional to the time derivative of the attenuation of detected light intensity (MAAP: in reflection; Aethalometer: in transmission), not on the absolute attenuation values.Any dependence of the measured absorption on the loading (attenuation) is evidence of a bias due to the loading of the spot.The parameter, which was shown to be most representative of these effects, is the slope of the attenuation or absorption coefficient or attenuation coefficient as a function of the loading (following the logic in Park et al., 2010).We will demonstrate the usefulness of this method as an independent criterion of the compensation for both absorption photometers used in the study: the MAAP and the Aethalometer AE-31 by plotting the average value of the attenuation or absorption coefficient as a function of the loading of the spot with BC.Hyvärinen et al. (2013) observed that at high BC concentrations the MAAP measurements are also affected by artifacts which are different to those seen in other filterbased instruments, such as the PSAP or the Aethalometer.Two solutions were proposed to avoid or compensate for this problem: either applying a compensation algorithm to the b atn, MAAP given directly by the MAAP's internal algorithm; or using reflectance signal only, as obtained from the photodetector at 165 • , which was found to be relatively free of this artifact, albeit exhibiting more scatter and not completely eliminating the artifact jump at the tape advance (Hyvärinen et al., 2013).

MAAP compensation
In this work, the MAAP was set to give the photodetector raw signals as an output, so the reflected signal at 165 • was chosen to obtain the absorption coefficient at 637 nm.The attenuation coefficient can be determined from the reflected signals by using the following equation (Petzold et al., 2005): where, in our case, (R 0 /R) is the ratio of the photodetector signal at 165 • for a particle free and a particle loaded filter, respectively.The multiplication factor of 0.5 in Eq. (3) has to be applied because the light passes through the layer of sampled aerosol twice before reaching the photodetector.The value obtained from Eq. (3) does not correspond to the absorption coefficient since there is a filter loading effect which influences the measured signals.Petzold et al. (2005) determined an empirical method to compensate this artifact using test aerosols.These test aerosols consisted of pure black aerosol samples from kerosene flame particles, and externally mixed grey and black aerosols of varying single scattering albedo.The obtained relation for the aerosol absorption coefficient for these aerosols using the reflected signal is: Similar to Hyvärinen et al. (2013) values automatically provided by the MAAP (b atn, MAAP ) are lower than those calculated from raw reflected signal for high absorption coefficients (Fig. 1).We do not observe any systematic saturation in the measurements as found previously (Kanaya et al., 2008) or large discontinuities in the data at the tape advance (Brito et al., 2013).This indicates that the use of the reflection signals rather than the default values compensated the data successfully.To check the compensation, we plot (Fig. 2a and b) the average value of the MAAP attenuation and absorption coefficient in a spot loading bin (with width 0.3 µg cm −2 ) as a function of the loading of the spot with BC between a clean filter (no loading) and the BC value with enough data in the campaign to gather enough statistics (just above 5 µg cm −2 ).The slope of the compensated MAAP absorption coefficient is half of that of the uncompensated ones, 1879 proving that the compensation reduces the loading effects.Another feature which is informative is the intersect of the fit: this is the uncompensated MAAP attenuation coefficient (7.83 Mm −1 ; Fig. 2a) or the compensated absorption coefficient (6.94 Mm −1 ; Fig. 2b) characteristic for the campaign.We propose to use the difference between the intersect 6.94 Mm −1 (that is compensated "average" absorption values) and the campaign average 7.25 Mm −1 as an additional internal criterion with an empirical limit of 5 %.The compensation of the MAAP data using the reflection measurements rather than the raw MAAP data satisfies this criterion, and we use the compensated MAAP absorption coefficient for comparison with the Aethalometer attenuation and absorption coefficients.

Data pre-processing
The Aethalometer raw data was first screened with three criteria: data showing excessive noise; and the upper and lower limits of the Aethalometer detection.High noise in the measurements degrades the determination of the absorption spectral dependence.This wavelength dependence can be approximated by a power law expression (Ångström, 1929).A measurement is considered too noisy if the curve cannot fit a spectrum satisfactorily (Rizzo et al., 2011), therefore, only fits with R 2 higher than 0.85 have been chosen to reduce noise.The same data treatment was applied to the Nephelometer data as well.
The second criterion is related to the lower detection limit of the instrument, as the measured attenuation should be above it.According to the Aethalometer manual, this limit corresponds to a mass increase of 5 ng and can be expressed in units of the b abs (Mm −1 ) by the following expression: where m = 5 ng, σ ATN (λ) is the mass-specific attenuation cross-section in m 2 g −1 , ∆t is the time based period (set to 5 min in our case), and Q is the sampling flow in L min −1 (4 L min −1 in our instrument).The attenuation cross section is calculated using the following equation: This equation is based on a calibration at 880 nm using the Malissa-Novakov method, a solvent based thermal desorption method for elemental carbon analysis (Gundel et al., 1984).Also, as a consequence of random voltage fluctuations and from recording discontinuities which take place during the Aethalometer's tape advance a third criterion related to the upper detection limit is applied (deCastro et al., 2008).To identify these bad quality data, the 99th percentile of the attenuation coefficient (i.e., 162.2 Mm −1 at 660 nm) was chosen as the extreme upper limit.
From the whole data set, 4 % of the data did not satisfy any of these three criteria, leaving a total number of 8276 of attenuation measurements.

Data compensation
Aethalometer artifacts, unlike MAAP, are well known and five different algorithms have been proposed to correct them (Weingartner et al., 2003;Arnott et al., 2005;Schmid et al., 2006;Virkkula et al., 2007;Collaud-Coen et al., 2010).The beam crossing the filter suffers: (1) multi-scattering effects due to the filter fibers, (2) single scattering effects, due to the aerosol particles deposited in the filter, and (3) filter loading effects, which are related to the shadowing produced as the particles accumulate on the filter.The result of these optical interactions is the b ATN is generally larger than b abs (Petzold et al., 1997;Ballach et al., 2001).
Due to the aforementioned aerosol-filter interactions, the Aethalometer requires specific site compensation factors.In this sense, we have calculated this compensation 1881 factors for our site and used them to compensate the Aethalometer data by apply-ing the algorithms proposed by Weingartner et al. (2003) and Schmid et al. (2006).
The reason to select these two algorithms out of the five existing ones is that Wein-gartner et al. (2003) and Schmid et al. (2006) found a way to obtain the spectral dependence of the multiple scattering compensation factor from aerosol single scat-tering albedo.Knowing this spectral dependence will allow us to convert b ATN to b abs for all 7 Aethalometer wavelengths.The most important difference between these two algorithms lies in the fact that Weingartner's algorithm does not consider the artifact produced by the single scattering effect, whiles Schmid's does.
In their test with ammonium sulfate, Weingartner et al. (2003) found no significant dependence of b ATN on the scattering component of the aerosol in the filter.There-fore, they proposed an algorithm to obtain b abs without taking into account the single scattering effect of the aerosols: where C compensates for the multiple scattering effects exhibiting values ≥ 1, and R(f , ATN) for the loading-effect with values ≤ 1, respectively.The loading compensation depends on the amount of the sample collected on the filter and is hence dependent on the attenuation measured by the Aethalometer.The parameter can be expressed as proposed in Weingartner et al. (2003): where f is the shadowing factor, which depends on the type of the aerosols.This parameter has been calculated by minimizing the difference between the ratio of b ATN before and after the filter spot change.The calculated median f values were plotted as a function of λ and fitted with a linear equation (Sandradewi et al., 2008b).Values obtained for this parameter are shown in Table 1.
The effect of the compensation of the Aethalometer measurements is shown in Fig. 2c and d.The slope reduces significantly to a value which is non-distinguishable from 0 -the compensation is efficient in eliminating the loading effects.The difference between the intersect 7.14 Mm −1 and the campaign average 7.10 Mm −1 satisfies the empirical 5 % criterion.Additionally, the value lies extremely close to the average and the intersect determined from the analysis of the MAAP compensation (Sect.3.1).
Since the loading effect is small for lightly loaded filters (Weingartner et al., 2003), C can be determined by comparing low loaded Aethalometer measurements (ATN < 10 %) with the ones obtained by a different comparison instrument.In this study we chose the MAAP (b abs, MAAP ) as a comparison and C is determined using the equation proposed in Weingartner et al. (2003) as: C has been calculated at 637 nm, as the MAAP measurements are performed at this wavelength.Since the nearest wavelengths in the Aethalometer are 590 and 660 nm, b ATN (637 nm) has been calculated by approximating the wavelength dependence of the attenuation spectra for each measurement to a power-law expression, such as: where a is a fitting parameter and α atn is the Ångström exponent of attenuation.
After taking logarithm of Eq. ( 10), a "linear fit" was applied to the log-log curve: To calculate C, b ATN at 637 nm for the Aethalometer data has to be calculated since b abs, MAAP in Eq. ( 9) is related to 637 nm.The average C value at the wavelength of 637 nm was determined to be C 637 = 4.22± 0.06 from the arithmetic mean (95 % confidence level of the mean) of the ratios of b ATN (637) and b abs, MAAP (Schmid et al., 1883 Discussion In particular, from the different sites chosen in that paper, our value is close to those obtained (C 660 = 4.12 ± 0.06) at Cabauw (the Netherlands), a background site located near populated and industrialized areas (Collaud-Coen et al., 2010).Some publications (Lack et al., 2009;Nakayama et al., 2010) indicate that the multiple scattering compensation factor might depend on the particle size, as the aerosol penetration depth into the filter varies depending on the size.Since this effect is still difficult to quantify (Rizzo et al., 2011), it is not straightforward to conclude if there will be an overestimation or underestimation in our absorption coefficients due to this effect.A significant part of the aerosol load on the filter tape on the measurement site is always diesel exhaust particles, so at least this part of the loading is extremely homogeneous.Therefore, we did not consider this effect further in this manuscript.Schmid et al. (2006) parameterized (based on Arnott et al., 2005) the spectral dependence of C using the following expression:

Spectral dependence of C and aethalometer absorption coefficients
where C * is the multiple scattering compensation factor which includes the effects of aerosol scattering, m s is the fraction of the aerosol scattering coefficient erroneously interpreted as absorption, ω 0 is the aerosol single scattering albedo, and λ is the wavelength.Both C * (λ) and m s were calculated for ammonium sulfate particles by Arnott et al. (2005).They however showed that these values gave unsatisfactory results when used to correct ambient measurements.
Equation ( 12) has been used to obtain C * 637 using C 637 calculated in the previous section with Eq. ( 9).Once it is calculated, C * values for the other wavelengths will be obtained assuming the same wavelength dependence as in Arnott et al. (2005).Values of m s used here were taken from Arnott et al. (2005) since with our data it is not possible to determine them.ω 0 (637) was calculated using the following expression (Schmid et al., 2006): where b abs, MAAP is the MAAP absorption coefficient measured at 637 nm and b scat is the scattering coefficient from the Nephelometer, interpolated to 637 nm by using Eqs.( 10) and ( 11) for b scat (λ).
To obtain the different C values for the remaining Aethalometer wavelengths, ω 0 (λ), C * (λ), and m s (λ), are required (Eq.12).Equation ( 13) is not useful in this case since MAAP absorption is only obtained at one wavelength and cannot be applied for the other Aethalometer channels.Therefore, for calculating ω 0 (λ) we used the following expression (Schmid et al., 2006): where ω 0,ref is the reference single scattering albedo calculated at 637 nm (λ ref ) using Eq. ( 13), α scat is the Ångström exponent of scattering and α abs is the Ångström exponent of absorption.α scat depends on the b scat measured at different wavelengths and can be obtained from the nephelometer measurements (450-700 nm).However, α abs is calculated from the b abs measured by the Aethalometer at different wavelengths.from our values in Granada during July 2012, α scat,450−700 = 1.3 ± 0.3 and ω 0 (637) = 0.76 ± 0.08, consistent with the prevalence of diesel traffic on the highway near the measurement site.These differences between both sites are significant and make the parameterization function of C(λ) given by Schmid et al. (2006) not applicable to our measurements.Therefore, based on their work, the iterative procedure has been used to calculate α abs and, then, C(λ) and b abs (λ), following these steps: 1.The Ångström exponent of attenuation (α atn ) is calculated by fitting the logarithm of the spectral attenuation coefficients with a linear fit as presented in Eq. ( 11).
2. The α atn is used in Eq. ( 14) as a first estimation of α abs , to obtain ω 0 (λ) and then, the compensation factor C for each Aethalometer wavelength is calculated using Eq. ( 12).
3. The compensation factor C(λ) at different wavelengths is used in Eq. ( 7) to obtain the new compensated absorption coefficients (b abs ).
4. This new b abs (λ) is used again in step 1 to get a better estimation for α abs .

Non-compensated/compensated Aethalometer vs. MAAP data
In Fig. 3a, non-compensated Aethalometer data (b ATN ) at 637 nm obtained using Eq. ( 11) are compared to b abs, MAAP values.The Aethalometer and the MAAP are well correlated, with R = 0.917, while the slope 3.66 shows the relationship between the Aethalometer attenuation coefficient and the MAAP absorption.After Weingartner's compensation (Weingartner et al., 2003)  1.005 and there appears a slight increase of the correlation coefficient: from 0.917 to 0.926 (Fig. 3b).The compensation factors used (C 637 = 4.22 ± 0.06 and f 637 ≈ f 660 = 1.131) were already obtained in Sect.3.2.2.Dispersion of some points is due to some noise in the data which has not been removed by the filtration algorithm described above in Sect.3.2.1.We can further reduce this dispersion by averaging the data over a more extended period of time.Furthermore, the compensated b abs, W values are on average 73 % lower than corresponding non-corrected b ATN values and differ from b abs, MAAP on average by around 2 %.
The statistics of the comparison between the two data sets are shown in Table 2. Average values of b abs, MAAP and b abs, W with their standard deviations are 9.2±5.8Mm −1 and 9.4 ± 6.1 Mm −1 , respectively.The 25th and 95th percentile for both data sets are very similar indicating a similar distribution of the data.The statistical analysis shows a good agreement between b abs, MAAP and b abs, W . Therefore it is concluded that Weingartner's compensation, with our site specific compensation parameters, compensates satisfactorily the Aethalometer b ATN to obtain absorption coefficients which agree well with those measured by the MAAP.

Compensation parameters at different wavelengths and final compensation
The spectral dependence of C can be obtained using Eq. ( 12).To start, C * and m s as proposed in Arnott et al. (2005) for 660 nm were used (see of the two absorption photometers, C * 637 can be calculated by using Eq. ( 12) and assuming m s (637) ≈ m s (660) = 0.0713 (Table 1).The value obtained is 3.93 which differs from the value obtained by Arnott et al. (2005) in the laboratory (C * Arnott of 2.182 at 660 nm), as can be seen in Table 1.On the other hand, this value is consistent with their proposal of C * (521) = 3.69 for ambient measurements.Following this, C * (λ) values at other wavelengths have been calculated assuming the same spectral dependence of C * Arnott (λ) given by Arnott et al. (2005) and normalized to C * Arnott (660).As can be seen from Table 1, the obtained C * (λ) are much higher than those estimated by Arnott et al. (2005) or by Weingartner et al. (2003) in the laboratory.
Once we have retrieved the C * (λ) values, the iterative method explained in Sect.3.2.3 has been applied for each measurement to calculate the C(λ), which will be used to compensate b ATN (λ).Almost 93 % of the data converges at the fourth iteration while the remaining 7 % do at the fifth iteration.The averaged values for C(λ) are shown in Table 1.It is evident that C increases with λ, from 3.42 in the UV (370 nm) channel to 4.59 in the infrared (950 nm).Comparing these values with our reference C 637 , calculated using Eq. ( 9), a higher wavelength dependence is observed at lower wavelengths (∼ 23 % for C 370 ) that at near infra-red wavelengths (∼ 4 % for C 950 ).Weingartner et al. (2003) assumed that there was no wavelength dependence on C since the difference in their measurements between 450-660 nm was smaller than 10 %.This is consistent with our data since the difference in C between 470-660 nm is ≈ 8 %, and the only wavelength differing more is 370 nm (see Table 1).Furthermore, the differences in our calculated C values at each wavelength are similar to the ones reported in Schmid et al. (2006) for a α abs ≈ 1.In addition, the average α abs obtained with the iterative method was found to be 1.09 with a standard deviation of around 0.25.We performed a sensitivity analysis on the effect of the compensation on the C values: a 15 % change of the parameter f causes a change in the C values smaller than 1 % for all wavelengths.
The obtained absorption coefficients are, on average, around a 72 % and 79 % lower than the corresponding b ATN values measured at 370 and 950 nm, respectively.Rizzo et al. ( 2011) reported a decrease of a 75 % at 450 nm, which is similar to the findings in this study.Figure 4 shows a boxplot of the compensated absorption coefficients b abs at each wavelength obtained from the Aethalometer and the MAAP.On one hand, it can be seen that the MAAP absorption coefficient is in reasonable agreement with the boxes from the compensated Aethalometer data at 590 and 660 nm.This is expected when the data are properly compensated, since the spectral behavior of the absorption coefficient follows a power law as shown in Eq. ( 10).Rizzo et al. (2011) performed sensitivity tests and showed that the main source of error due to b abs and α abs is driven by the choice of α scat .Since they used averaged values over the whole measurement period, their sensitivity tests consisted in varying these values and see how this affected the final result for b abs and α abs .They obtained a maximum deviation in the results of ±10 % and ±40 % on b abs and α abs , respectively.Following their tests, in our measurements, only ±1 % for b abs and ±2 % for α abs deviations are observed.This is due to the fact that we have used concurrent absorption and scattering measurements at multiple wavelengths and have thus reduced deviations in b abs and α abs significantly.

Spectral absorption coefficient, spectral single scattering albedo, and
Ångström exponent of absorption

Temporal evolution of aerosol spectral properties
The statistics obtained for b abs (λ), and ω 0 (λ) from the whole measurement campaign are shown in Table 3.For α abs, 370-950 an average value of 1.09 is calculated with a standard deviation of 0.25.Values of α abs around 1 are related to the presence of traffic aerosols (Sun et al., 2007) and are typical of urban areas (Bergstrom et al., 2007).This shows a good agreement with the characteristics of the site.Maxima are greater than 1.55 (P95) and reach values up to 2.5.These high values are related to special , respectively.Lowest values of b abs were measured at near infrared channels, 880 nm and 950 nm.In Table 3 it is obvious that for the 7 channels the difference between P25 and P95 is considerable.This is due to the fact that our site is strongly affected by traffic which generates high peaks in the absorption coefficients causing strong differences between b abs at low traffic hours and rush hours (see Fig. 5d).
Finally, the statistics for ω 0 obtained for the seven Aethalometer channels were calculated.Averaged ω 0 values from 370 nm to 950 nm lie between 0.78 and 0.74, with a standard deviation of 0.08-0.09.Lyamani et al. (2010) reported an average value of ω 0 (670) of 0.73 ± 0.06 in summer, which is close to the value obtained in this work, ω 0 (670) of 0.76 ± 0.08.P25 values are 0.68 and 0.74 at 950 nm and 370 nm, and P95 values are 0.88 and 0.90 at 950 nm and 370 nm, respectively.In general, the ω 0 (λ) average values show that during the measurement period the atmosphere in Granada at surface level contained a large fraction of absorbing particles.
Figure 5 shows the temporal evolution of hourly and daily averaged data for the different optical aerosol parameters calculated in the previous section during the period from 29 June to 29 July.Figure 5a and b display the hourly and daily average values of ω 0 (λ) obtained for the 370 nm and 950 nm channels of the Aethalometer.In Fig. 5a it can be observed that the ω 0 presents minima of 0.5-0.6 correspond to maximum peaks of b abs (Fig. 5d).These maxima in b abs, are consequence of traffic emissions, this suggests that lower values of ω 0 are caused by traffic emissions which increase the absorbing component at surface level (Lyamani et al., 2010(Lyamani et al., , 2011)).During the whole period ω 0,950 is higher than ω 0,370 , except for ordinal days 181 (29 June) and 202-203 (20 and 21 July) (Fig. 5).A more detailed analysis of these days will be discussed below.In the case of the α abs (Fig. 5c), except for the mentioned special episodes, values of this parameter show point to point variations no greater than 0.1.Hourly average values of α abs are mostly between 0.9 and 1.2.The spectral dependence of b abs does not vary a lot, as well as the daily average of α abs , which exhibits values around 1.1 for most of the month.This corroborates the fact that the site in Granada is dominated by urban aerosols whose major source is traffic (Lyamani et al., 2008(Lyamani et al., , 2010;;Titos et al., 2012).
Figure 5d shows the hourly average values of the b abs (370) and b abs (950).Both channels exhibit very similar temporal behavior, although absorption coefficients in the UV are considerably higher than those in the near infrared.The same evolution is seen for b abs at the other five channels, with values located between these two channels (not shown).The values of b abs on Sundays are in general considerably lower than those obtained during the rest of the week.This can be related to the fact that the traffic emissions on Sunday are considerably lower than during the rest of the week.
For the two events already mentioned -ordinal days 181 (29 June) and 202-203 (20 and 21 July), an extended analysis of surface level characteristics was made using back trajectories and results from the Navy Aerosol Analysis and Prediction System (NAAPS, http://www.nrlmry.navy.mil/aerosol/).The model used for calculating back trajectories is the Hybrid Single-Particle Lagrangian Integrated Trajectory Model (HYS-PLIT) developed by the National Oceanic and Atmosphere Administration (Draxler and Rolph, 2003).The model combines a Lagrangian approximation for resolving the air mass transport with an Eulerian approximation for the diffusion of pollutants.On the other hand, the NAAPS model is a global model that predicts concentrations of sulfates, dust and smoke aerosols in ground level and in the atmospheric column.The NAAPS global model is an offline model that uses meteorological analysis and forecasts fields from the Navy Operational Global Analysis and Prediction System (Hogan and Brody, 1993) developed by Christensen (1997).LIDAR images obtained from Granada's website (http://atmosfera.ugr.es)have been included as well to complement the analysis of surface level characteristics.Lidar measurements were taken with a Raman lidar model LR331D400 (Raymetrics S.A., Greece) included in EARLINET (Boesenberg et al., 2001).The Raman lidar system is configured in a monostatic biaxial alignment pointing vertically to the zenith.The transmitter of the lidar system is a pulsed Nd:YAG laser with fundamental emission at 1064 nm and additional emissions at 532 and 355 nm by using second and third harmonic generators.Output energies are 110, 65 and 60 mJ at 1064, 532 and 355 nm respectively, and pulses of 7 to 9 ns can be fired with a pulse repetition frequency of 1, 2, 5 and 10 Hz.The receiving system consists of a 0.4 m diameter Cassegrain telescope and a wavelength separation unit with dichroic mirrors, interference filters and a polarization cube, that discriminates seven channels corresponding to elastic wavelengths (1064, 532 parallel-polarized, 532 perpendicular-polarized and 355 nm) and to nitrogen and water vapor Raman-shifted wavelengths (387, 408 and 607 nm).Raman signals were not used in this study.The optics set-up is such that the maximum overlap is reached at about 400 m above the instrument.More information on the instrument can be found in Guerrero-Rascado et al. (2008) The days selected for a further analysis are 29-30 June and 20-21 July, since α abs values are considerable higher in comparison with the rest of the month (Fig. 5c). Figure 6 shows the 5 days back trajectories arriving at Granada on 29 June and 20 July 2012.Figure 7, shows the concentration of different types of particles (dust, biomass smoke) at ground level predicted by the NAAPS model for the 29 June and 20 July, respectively, at 12.00 and 18.00 UTC.On the 29 June, the trajectories arriving at 1500 and 3000 m a.g.l.came from the West Sahara, while the trajectory at 500 m a.g.l.came from the Mediterranean Sea between Spain and North Africa (Fig. 6a).In addition, from Fig. 7a it can be seen that on 29 June, NAAPS predicted a significant surface concentration of dust particles (80-160 µg m −3 ) at Granada.Furthermore, no significant contributions from biomass burning and sulfate aerosols were obtained from the NAAPS model.Figure 8a shows the temporal evolution of the lidar range corrected signal (RCS) at 532 nm, on 29 July.Early in the morning it is evident the presence of an aerosol layer decupled of the Planetary Boundary Layer.The analysis of the backward trajectories suggests that the decoupled layer extending from 2000 to 3000 m a.s.l. was originated in North Africa.This is a clear evidence of the Saharan dust outbreak over our station.
To study whether there is any noticeable effect in the optical properties at ground level or not, the 24 h average values of both ω 0 (λ) and b abs (λ) obtained on 29-30 June and 20-21 July are shown in Fig. 9a and b, respectively.Figure 9a shows a strong spectral dependence of ω 0 (λ) on 29 June with considerably lower values (∼ 0.80) for the UV channel than in the IR, which is a well-known characteristic of dust.On 30 June, the spectral dependence of ω 0 (λ) is not strong, and can be related to fact that the dust intrusion is over and urban aerosols are dominant again.Collaud-Coen et al. (2004) proposed a method for detecting Saharan dust events based on the change of the Ångström exponent of single scattering albedo.Negative exponent values of the ω 0 (α SSA ) are due to the large size of mineral aerosols.In our case, on 29 June a α SSA = −0.09 is obtained, which is related to the presence of dust in the atmosphere.On 30 June this wavelength dependence decreased to α SSA = −0.007,indicating that the mineral dust intrusion was over.In the case of the b abs (λ) a stronger spectral dependence for the values obtained on 29 than for the 30 of June is found (Fig. 9b).This spectral dependence leads to the values of α abs of 1.34 and 1.09 on 29 and 30 June, respectively.The return of the exponent to a value close to 1 is consistent with diesel aerosols becoming predominant again.Figure 5c shows that the α abs reached values up to nearly 2.0 during the dust intrusion on 29 June.To study the consistency of the α abs values obtained during this dust event, they were compared with those obtained in other studies.Collaud-Coen et al. (2004) reported α abs values between 1.5 and 1.8 during Saharan dust events at Jungfraujoch, as well as a non-negligible spectral dependence of ω 0 .Values of α abs of 2.2 were reported by Bergstrom et al. (2002) for mixed urban pollution and desert dust aerosols in ACE Asia program.Therefore, it can be concluded that our determination of the Saharan dust episodes and the associated values of the Ångström exponent of single scattering albedo are in good agreement with the published values.Another event was detected on 20 (day 202) and 21 (day 203) July.HYSPLIT back trajectories arriving on the 20 July at 1500 and 300 m a.g.l.reached Granada proceeding from the Atlantic Ocean, entering the Iberian Peninsula from the Northwest, while at 500 m a.g.l.they came again from the Mediterranean Sea between Spain and North Africa (Fig. 6b).NAAPS predictions, presented in Fig. 7b and c, show influence of both dust and smoke aerosols at surface level on 20 July.The temporal evolution of RCS at 532 nm obtained from lidar measurements (Fig. 8b) shows the presence of aerosol particles up to 4000 m a.s.l.According to the back trajectories analysis, air masses with different origin contribute to this aerosol layer.This includes North African air masses, loaded with mineral dust, and air masses from regions affected by fires.
Values on the 20 July were daily averaged starting at 12.00, when the event is first detected.The analysis of the wavelength dependence of ω 0 leads to α SSA = −0.08 on the 20 July, showing again presence of dust, and to α SSA = −0.011 on 21 July, when the episode is over.In Fig. 9a a stronger spectral dependence of ω 0 during the event is observed.For 20 July, differences between ω 0 in the UV and in the IR are around 0.06, while on 21 July are lower than 0.01.Similar to the event on 29 July, variations in the spectral dependence indicate a change in the dominant source and origin of the aerosols which is corroborated with NAAPS predictions of dust and smoke aerosols at surface level.Figure 9b shows that the spectral dependence of b abs does not show significant differences from one day to another, the values of α abs are 1.28 and 1.31 for the 20 and 21 July, respectively.
Although there are similarities in the spectral dependencies of the optical parameters, both events differ: on 29 June the event is dominated by dust with larger ω 0 and lower b abs ; on 20 and 21 July the aerosol is a mixture of smoke and dust, with lower ω 0 and higher b abs .

Comparison between working days and Sundays
In Fig. 5d, based on the daily averages (blue dots), it was already observed that b abs was generally higher during working days than on Sundays.Due to these differences, In Fig. 10a it can be seen that during daytime ω 0 is in general slightly higher on Sundays than working days.At night-time, ω 0 values are very similar, the difference being only 0.05.The fact that a less absorptive atmosphere is found on weekends than on working days results directly from the b abs behavior.During the day time, b abs for working days is considerably higher than during Sundays (Fig. 10b), which also causes lower ω 0 from Monday to Friday.It is also seen that on working days and Sundays b abs has two maxima, one in the morning (between 04-09 h) and another in the afternoon (between 18-21 h).Those maxima are strongly related to minima in ω 0 .Usually, the morning b abs maximum is considerably higher on working days than on Sundays due to its relation with traffic volume, which is more intense from Monday to Friday.This causes a more pronounced decrease in ω 0 during the working days.On the other hand, the evening maximum appears to be very similar for working days and Sundays.
The diurnal variation of the exponent α abs shows no important differences between working days and Sundays.The α abs values obtained range from 1.00 to 1.25, which are typical values for urban aerosols (Sun et al., 2007).Unlike ω 0 , this parameter appears to be not affected by maxima in b abs .This corroborates the fact that α abs is a parameter related to the type of absorbing particles and their source rather than to their concentration.

Conclusions
The aim of this study has been to obtain the compensated absorption coefficient for all seven Aethalometer AE-31 channels generalizing the method proposed by Schmid et al. (2006).Once the Aethalometer coefficients have been compensated, the temporal evolution of ω 0 (λ), b abs (λ), and α abs (λ) was analyzed.In addition, during the measurement period two different sources of aerosolized particulate matter affecting Granada were investigated and the differences between working days and Sundays were presented.
The attenuation coefficient at 637 nm was compensated applying the method pro-posed by Weingartner et al. (2003) and then compared to b abs obtained from MAAP measurements.Furthermore, site specific parameters for the loading (R(f , ATN)) and multiple scattering (C) were obtained.The values for f range from 1.15 at 370 nm to 1.11 at 950 nm, which is consistent with previous values reported by Weingartner et al. (2003), Sandradewi et al. (2008b), andCollaud-Coen et al. (2010), among oth-ers.Compensated data from both instruments were checked for any remaining loading effects and none were found.
For the multiple scattering effect, the C value at 637 nm was found to be 4.22 ± 0.06 and it is in agreement with values published in the literature for sites with similar characteristics.Results of the comparison with the MAAP give a slope of 1.005 and a correlation coefficient of 0.926.Therefore, the compensated factors obtained in this study provide a satisfactory compensation of the data and are used to calculate the C factor at other wavelengths.
To obtain the wavelength dependence of C(λ), the additional parameters m s and C * are needed, which were obtained by Arnott et al. (2005) in laboratory experiments.These values proved to be quite low when applied to ambient measurements, so new C * (λ) values were calculated in this study.The C * value for ambient data calculated at 637 nm for this campaign exhibits a value of 3.93, while m s used is the value pro-vided by Arnott et al. (2005).Differences in C * agree with those expected from Arnott et al. (2005) for ambient measurements.The same spectral dependence as in Arnott et al. (2005) has been used to obtain C * at other wavelengths.
Applying the iterative procedure proposed by Schmid et al. (2006) Once the different parameters are corrected and calculated at each wavelength, a study of their temporal evolution was performed.Averaged values of ω 0 range between 0.784 (at 370 nm) and 0.737 (at 950 nm) indicating that during July 2012 Granada's atmosphere was significantly burdened with absorbing particles.This is in agreement with results obtained by Lyamani et al. (2010) for the summer period.An average value of 1 for α abs is explained by the fact that Granada's site is dominated by urban aerosols (Sun et al., 2007) whose major source is traffic.
Finally, the difference between working days and Sundays shows that the absorption coefficients are higher during working days than during Sundays as a consequence of the traffic intensity on the highway close to the station.These values cause decreases in ω 0 leading to an atmosphere with more absorbing particles during the working days.However, α abs is not affected and exhibits no differences between work days and Sundays as this value is more related to the aerosol type and their source, which does not strongly change at this site.
Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | ; 1873 Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper |

1875
Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper |

1877
Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper |

Discussion
Paper | Discussion Paper | Discussion Paper | Discussion Paper |

Discussion
Paper | Discussion Paper | Discussion Paper | Discussion Paper | Paper | Discussion Paper | Discussion Paper | Discussion Paper | 2006).Collaud-Coen et al. (2010) obtained different C values ranging between 2.8 and 7.8 at several sites.
Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Schmid et al. (2006) used fixed values for α scat,450−700 and ω 0 (545) measurements given byChand et al. (2006) for the Amazonia, and a range of different α abs based on Kirchstetter et al. (2004) to obtain a parameterization function of C(λ) and α abs .Then they applied an iterative procedure over this function to obtain C(λ) and α abs .The fixed values used were α scat = 2.0 ± 0.4 and ω 0 (532) = 0.92 ± 0.02, which differ 1885 Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | is applied to the Aethalometer attenuation coefficient, we obtain the compensated Aethalometer absorption coefficient (b abs, W ). The slope between b abs, W and b abs, MAAP changes to Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper |

1889(
Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper |aerosol episodes (dust and biomass burning particles) affecting Granada during the measurement period, which will be described in more detail below.Average values obtained for the compensated b abs ranged from 16 Mm −1 nm) with standard deviations of 10 and 3 Mm −1 Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper |

1891
Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper |

1893
Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | the hourly average values in UTC of ω 0 (370), ω 0 (950), b abs (370), b abs (950), and α abs, 370-950 have been calculated for working days (Monday to Friday) and Sundays.Average diurnals variations of ω 0 , b abs , and α abs for working days and Sundays are shown in Fig. 10a-c, respectively.

1903
Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper |

Figure 4 .
Figure 4. Boxplot of the Aethalometer (black) and MAAP (red) corrected absorption coefficients (babs).The square inside the boxes represents the average; the central line corresponds to the median; the edges of the box are the 25 th and 75 th percentiles; the whiskers correspond to 5 th and 95 th percentiles; and horizontal lines outside the boxes, maxima and minima.

Fig. 4 .
Fig. 4. Boxplot of the Aethalometer (black) and MAAP (red) corrected absorption coefficients (b abs ).The square inside the boxes represents the average; the central line corresponds to the median; the edges of the box are the 25th and 75th percentiles; the whiskers correspond to 5th and 95th percentiles; and horizontal lines outside the boxes, maxima and minima. 37
Table 2).The C 660 value of 2.21 was obtained, which is only ∼ 50 % of the value obtained at 637 nm (see previous Sect.3.2.2).Applying this compensation factor leads to unsatisfactory comparison with b abs, MAAP, since the new compensated data is overestimating the MAAP data (slope of 1.92).The statistics of the comparison of this two data sets are shown in Table 2 as b abs, Arnott .These results corroborate that the C Arnott et al. (2005)values fromArnott et al. (2005)are not applicable for this campaign.To get a better estimation of the spectral dependence on C, specific ambient measurements were calculated to obtain C * (λ).As we only have available C 637 for the direct comparison 1887 Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper |

Table 1 .
Schmid et al. (2006) obtained for the 7 Aethalometer wavelengths: f corresponds to the "shadowing effect" in the loading compensation; 100α is the value given inArnott et al. (2005)to the scattering fraction of particles (m s inSchmid et al., 2006); M is the value given inArnott et al. (2005)to the multiple scattering compensation factor (C * inSchmid et al. (2006)and in this work); C * and C are the values specifically site-calculated in this work.

Table 3 .
Statistics of data for α