A First Case Study of CCN Concentrations from Spaceborne Lidar Observations

Abstract

We present here the first cloud condensation nuclei (CCN) concentration profiles derived from measurements with the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) aboard the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO), for different aerosol types at a supersaturation of 0.15%. CCN concentrations, along with the corresponding uncertainties, were inferred for a nighttime CALIPSO overpass on 9 September 2011, with coincident observations with the Facility for Airborne Atmospheric Measurements (FAAM) BAe-146 research aircraft, within the framework of the Evaluation of CALIPSO’s Aerosol Classification scheme over Eastern Mediterranean (ACEMED) research campaign over Thessaloniki, Greece. The CALIPSO aerosol typing is evaluated, based on data from the Copernicus Atmosphere Monitoring Service (CAMS) reanalysis. Backward trajectories and satellite-based fire counts are used to examine the origin of air masses on that day. Our CCN retrievals are evaluated against particle number concentration retrievals at different height levels, based on the ACEMED airborne measurements and compared against CCN-related retrievals from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensors aboard Terra and Aqua product over Thessaloniki showing that it is feasible to obtain CCN concentrations from CALIPSO, with an uncertainty of a factor of two to three.

1. Introduction

Aerosols impact the formation, the optical properties, and the life cycle of clouds. These aerosol-cloud interactions (ACI) mark the largest source of uncertainty in climate change studies [1]. Cloud droplets emerge from the condensation of supersaturated water vapor on aerosol particles that can nucleate liquid cloud droplets: so-called cloud condensation nuclei (CCN). Particles can act as CCN depending on, e.g., their size, solute content, the presence of surface-active or slightly soluble substances, the wettability and shape of insoluble particles, as well as the ambient supersaturation. The CCN concentration hence has a strong link to the cloud droplet number concentration (CDNC) [2].

High aerosol concentrations increase CDNC, which, at constant cloud water contents, leads to smaller cloud droplet effective radii, to an increase in cloud albedo, and, hence, to a negative radiative forcing (first indirect effect) [3]. This increase in CDNC is expected to increase cloud lifetime and cloud coverage reducing precipitation formation efficiency (second indirect effect) [4]. For polluted convective clouds, the suppression of precipitation may allow cloud droplets to reach the freezing level, where they can release latent heat which further invigorates cloud updrafts. This is expected to suppress precipitation and the development of downdrafts, and so prolong the growth of convective clouds [5,6]. The presence of absorbing aerosols leads to local heating that can evaporate existing clouds, cools the surface, stabilizes the temperature profile, and, hence, suppresses convection and cloud formation. This is the so-called semi-direct aerosol effect [7,8]. However, absorbing aerosols affect clouds also thermodynamically and may increase or decrease cloud coverage over an area, depending on the relative position of the aerosol and cloud layers [9].

All these effects are competing and the difficulty of discriminating between aerosol effects and effects of the meteorological background undermine the precise quantification of the influence that aerosols exert on clouds, radiation, and precipitation, and by extension on climate [10]. On top of that, the efforts to constrain model ACI estimates by observational data are based on the assumption that the aerosol optical depth (AOD) or Aerosol Index (AI: AOD multiplied by the Ångström exponent) measured by satellite sensors constitute a good proxy for the number of CCN that reside in the atmosphere [11,12,13,14,15,16,17,18,19,20].

Shinozuka et al. [21] suggest that AOD or AI are not always a good proxy of CCN concentrations in the atmosphere (see also Stier et al., [22]), and that parameterizations that connect AOD with CCN concentrations involve an uncertainty of a factor of two to three. Based on these findings and previous studies [23,24,25,26], more recent studies show the potential of using ground-based polarization lidar measurements to obtain estimates of CCN concentration [27,28,29,30,31]. In addition, satellite-based polarimeters, such as the Polarization and Directionality of the Earth’s Reflectances (POLDER) instrument, may allow for obtaining a better proxy for the columnar CCN concentrations compared to AOD and AI [32]. The use of polarimeter-based products instead of AOD or AI points towards a more negative radiative forcing, due to aerosol-cloud interactions (RFaci) than the one reported in IPCC [1].

Following these developments, we present a first case study in which CCN concentrations estimated from measurements with the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) aboard the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO; Winker et al., [33]) satellite are compared to airborne in-situ measurements. Specifically, CCN concentrations for a nighttime CALIPSO track that crossed Eastern Europe, the Eastern Mediterranean, and Africa on 9 September 2011 are evaluated with airborne retrievals from the Evaluation of CALIPSO’s Aerosol Classification scheme over Eastern Mediterranean (ACEMED) research campaign. ACEMED offers a unique opportunity for a spaceborne lidar CCN evaluation study in the Eastern Mediterranean, as it made available CCN-related airborne data, which are spatially and temporally coincident with CALIPSO. In addition, our results are compared against data from the unvalidated CCN-related Moderate Resolution Imaging Spectroradiometer (MODIS)/Terra and Aqua product. The used data and methodology are presented in Section 2. Results are presented and discussed in Section 3. Section 4 summarizes the main findings and conclusions of this research.

2. Data and Methods

2.1. Data

2.1.1. CALIOP/CALIPSO

CALIPSO is part of the A-Train formation of satellites and has flown at a 705-km altitude sun-synchronous polar orbit since 2006. It crosses the equator at around 13:30 LT and has a 16-day repeat cycle [33]. CALIOP, the core instrument on CALIPSO, is an Nd:YAG elastic backscatter and polarization lidar, which transmits linear polarized laser light at 532 nm and 1064 nm [34]. CALIOP provides profiles of aerosol- and cloud-related properties [33] from the ground up to 30 km, by measuring the backscattered signals and the state of polarization of the backscattered light. The CALIPSO algorithm breaks down the measured signal into distinct atmospheric layers, e.g., clean air, aerosols, clouds, surface, [35,36], and attributes subtypes to cloud and aerosol layers [37]. Each aerosol subtype is associated with a pre-set lidar ratio (extinction-to-backscatter ratio) that is needed to calculate the aerosol extinction coefficient [36]. For this work, CALIPSO version 4.10 (v4.10) level-2 products [36] were acquired through the ICARE Data and Services Center (http://www.icare.univ-lille1.fr/). More specifically, we make use of aerosol and cloud backscatter coefficient profiles at 532 nm, the vertical feature mask (clean air, cloud, tropospheric aerosol, stratospheric aerosol, surface, subsurface, totally attenuated, low/no confidence), the aerosol subtype mask (marine, dust, polluted continental/smoke, clean continental, polluted dust, elevated smoke, dusty marine, polar stratospheric cloud (PSC) aerosols, volcanic ash, sulfate/other), and other auxiliary quality-related (e.g., cloud–aerosol discrimination score) and meteorological (relative humidity) parameters. These data are available at a horizontal resolution of 5 km, and a vertical resolution of 60 m for altitudes below 20 km.

2.1.2. ACEMED EUFAR Campaign

Two CALIPSO underflights were performed with the Facility for Airborne Atmospheric Measurements (FAAM, http://www.faam.ac.uk/) BAe-146 research aircraft within the framework of the ACEMED campaign over Thessaloniki, Greece (40.6° N, 22.9° E, 60 m a.s.l.): one on 2 September 2011 and one on 9 September 2011. Airborne in-situ and lidar remote-sensing data between 00:00 and 02:00 UTC on 9 September 2011 were used to retrieve ambient aerosol microphysical, optical, and hygroscopic properties at different altitudes [38], by applying the In situ/Remote sensing aerosol Retrieval Algorithm (IRRA). In order to derive the ambient aerosol properties in high-humidity conditions, IRRA utilizes the ISORROPIA II thermodynamic equilibrium model [39]. The BAe-146 aircraft performed several transect flights towards the SSW-NNE direction between 40 and 41.2° N, taking measurements at different altitudes over land and ocean (coastline located at 40.6° N) and a box-shape flight around the area of Thessaloniki. A detailed description of the aircraft measurements and the methodology is given in Tsekeri et al. [38]. For the scope of our study, we used the available dry particle number concentrations (in cm⁻³) at different height levels (1.8, 2.7 and 3.2 km a.s.l. over land and 1.3, 2.1, 2.7, 3.2 km a.s.l. over the sea) from Tsekeri et al. [38].

2.1.3. MODIS/Terra and Aqua

The two Moderate Resolution Imaging Spectroradiometer (MODIS) sensors [40] aboard NASA’s Terra and Aqua satellites are retrieving aerosol optical properties from space with daytime Equator crossing times at 10:30 LT (morning) and 13:30 LT (noon), respectively. For the area of Thessaloniki examined here, the two instruments take measurements at around 09:30 and 11:30 UTC [41]. MODIS has a viewing swath of 2330 km, and measures backscattered radiation in 36 spectral bands between 0.415 and 14.235 μm at a spatial resolution of 250, 500, and 1000 m. The instrument achieves almost global coverage on a daily basis. The Dark Target (DT) algorithms, one for vegetated/dark-soil land areas [42,43,44,45] and one for oceanic areas [43,44,45,46], provide a CCN product, which has rarely been validated or used so far. This product relies on the DT algorithm for oceanic areas, and is based on the retrieval of the size distribution of aerosols [43,47]: the MODIS CCN product is the integral of the small-mode aerosol size distribution from a fixed radius of 30 nm. Starting from Collection 6, the MODIS CCN (Cloud_Condensation_Nuclei) scientific data set (SDS) was renamed to PSML003_Ocean (Particles of the Small Mode Larger than 0.03 μm) [45]. In this work, we utilized MODIS/Terra and Aqua level-2 Collection 6.1 (MYD04_L2 and MYD04_L2) PSML003_ocean best solution data (in particles/cm² or else cm⁻²) to evaluate our CCN estimates.

2.1.4. CAMS Reanalysis

Aerosol mass mixing ratio data from the Copernicus Atmosphere Monitoring Service (CAMS) reanalysis product [48] are used to support our analysis. The CAMS reanalysis includes 3-dimensional fields of both chemical species and aerosols for the time period from 2003 to 2016. It is based on European Centre for Medium-Range Weather Forecasts (ECMWF) Integrated Forecast System (IFS) (CY42R1 release) and a 4D-Var data assimilation system [48,49,50]. Satellite retrievals of total column CO, tropospheric column NO₂, aerosol optical depth, and total column, partial column, and profile ozone retrievals are assimilated. The emissions include the MACCity (MACC and CityZEN EU projects) anthropogenic emissions, the GFAS (Global Fire Assimilation System) fire emissions, the MEGAN (Model of Emissions of Gases and Aerosols from Nature) biogenic emissions, and several greenhouse gas emission datasets (Inness et al., [48] and references therein). Online parameterizations are used to calculate sea salt and dust surface fluxes using as input surface wind speed and other factors [51,52]. The CAMS reanalysis aerosol dataset includes five basic aerosol types: sea salt, dust, organic matter, black carbon, and sulfates. These are further subdivided into eleven aerosol types: three different size bins are used for sea salt particles (0.03–0.50, 0.50–5.00, and 5.0–20.0 μm) and dust (0.03–0.55, 0.55–0.90, and 0.9–20.0 μm), while organic matter and black carbon are differentiated into a hydrophilic and a hydrophobic mode. The data used here cover the time period from 8 to 10 September 2011, with a temporal resolution of three hours. The spatial resolution is approximately 80 km with 60 hybrid pressure levels (top level at 0.1 hPa). The data were obtained from ECMWF’s web repository (https://apps.ecmwf.int/data-catalogues/cams-reanalysis/?class=mc&expver=eac4) through the corresponding WebAPI.

2.1.5. HYSPLIT Backward Trajectories

The Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model [53], developed by NOAA’s Air Resources Laboratory, is used to calculate 96-h backward trajectories that end over Thessaloniki, Greece, at six arrival heights of 0.5, 1.5, 2.5, 3.5, 5.5, and 7.5 km above ground level at 00:00 UTC on 9 September 2011. The calculations are driven by Global Data Assimilation System (GDAS) reanalysis data with a resolution of 1° × 1°.

2.1.6. FIRMS

Data from NASA’s Fire Information for Resource Management System (FIRMS) are utilized to investigate if the air masses arriving over Thessaloniki likely contained traces of biomass-burning smoke from active fires in the greater region. The fire location data as retrieved from MODIS spectral measurements [54] are available at a 1 km spatial resolution. The data used here are part of the Collection 6 MCD14ML dataset and cover the period from 5 to 9 September 2011.

2.2. Methodology—CCN Estimates

Prior to estimating CCN concentrations, the CALIPSO level-2 data are separated according to the contributions of dust and non-dust aerosols [55,56,57,58,59,60]. Specifically, the pure-dust-separation methodology used here was developed within the framework of the EARLINET-ESA LIVAS (LIdar climatology of Vertical Aerosol Structure for space-based lidar simulation studies) research project [61]. This method obtains an improved typing compared to the standard CALIPSO product [57,58]. An overview of the workflow for obtaining CCN concentration profiles from CALIPSO measurements is shown in detail in Figure 1.

First of all, data quality control filters are applied, following Winker et al. [62] and Marinou et al. [58]. The pure dust separation is applied only when the CALIPSO aerosol classification indicates the presence of dust, polluted dust, or dusty marine aerosols. The pure dust backscatter component (β₁) is calculated from the measured particle linear depolarization ratio (δ_p), the theoretical depolarization ratio of the pure dust (δ₁), and non-dust components (δ₂), and the total backscatter (β_t) of the aerosol layer [55] as

$${\mathsf{\beta }}_1={\mathsf{\beta }}_{\mathrm{t}}\frac{({\mathsf{\delta }}_{\mathrm{p}}-{\mathsf{\delta }}_2)(1+{\mathsf{\delta }}_1)}{({\mathsf{\delta }}_1-{\mathsf{\delta }}_2)(1+{\mathsf{\delta }}_{\mathrm{p}})}$$

(1)

where values of δ_p < 0.31 are assumed to represent a mixture of pure dust and non-dust aerosol components. Values of 0.31 ± 0.04 and 0.05 ± 0.03 are used for δ₁ and δ₂, respectively [55,63,64]. The pure dust backscatter coefficient profile is multiplied by an appropriate Saharan dust lidar ratio (LR) of 55 ± 11 sr to calculate the corresponding pure dust extinction coefficients at 532 nm [58].

The non-dust mixture backscatter coefficients (β_t–β₁) are dedicated to marine aerosols and polluted continental aerosol, depending on the CALIPSO classification, as dusty marine and polluted dust, respectively. Non-dust aerosol subtypes are considered as provided in the CALIPSO aerosol typing. As in this work we are interested in tropospheric aerosols which act as CCN, we focus on the tropospheric aerosol subtypes (marine, polluted continental/smoke, and clean continental), and the elevated smoke, which is also frequently observed in the lower troposphere. The LRs used to convert the respective backscatter coefficients to extinction coefficients are given in Kim et al. [36].

Here, we consider three major aerosol types: mineral dust, continental aerosol (polluted continental/smoke, clean continental, and elevated smoke), and marine aerosol. The corresponding CCN concentrations are obtained following recently suggested parameterizations (Equations (2)–(5)). The method followed for the CCN estimates was originally suggested by Mamouri and Ansmann [27] for ground-based lidar measurements in Cyprus, Eastern Mediterranean, and was recently applied over areas with diverse aerosol typologies, e.g., Barbados, Caribbean [29], Germany, Central Europe [30,65] and Tajikistan, Central Asia [31]. First, using the extinction coefficients (σ_d, σ_c and σ_m) for the three major aerosol types (dust, continental and marine), we calculated profiles with particle number concentrations (n_100,d,dry, n_50,c,dry and n_50,m,dry), following Equations (2)–(4). n_i,j,dry expresses the concentration of dry particles with a radius larger than i nm, which is the reservoir of favourable CCN for aerosol type j. For dust (index d), i equals to 100 nm, while for continental (index c) and marine aerosols (index m), i equals to 50 nm.

$${\mathrm{n}}_{100,\mathrm{d},\mathrm{dry}}={\mathrm{c}}_{100,\mathrm{d}}\hspace{0.17em}\cdot \hspace{0.17em}{{\mathsf{\sigma }}_{\mathrm{d}}}^{{\mathrm{x}}_{\mathrm{d}}}(\mathrm{z})$$

(2)

$${\mathrm{n}}_{50,\mathrm{c},\mathrm{dry}}={\mathrm{c}}_{60,\mathrm{c}}\hspace{0.17em}\cdot \hspace{0.17em}{{\mathsf{\sigma }}_{\mathrm{c}}}^{{\mathrm{x}}_{\mathrm{c}}}(\mathrm{z})$$

(3)

$${\mathrm{n}}_{50,\mathrm{m},\mathrm{dry}}={\mathrm{c}}_{100,\mathrm{m}}\hspace{0.17em}\cdot \hspace{0.17em}{{\mathsf{\sigma }}_{\mathrm{m}}}^{{\mathrm{x}}_{\mathrm{m}}}(\mathrm{z})$$

(4)

Following Table 3 in Mamouri and Ansmann [27], c_100,d was set to 6.5 ± 1.8 Mm cm⁻³ (based on measurements at Cape Verde and Barbados), c_60,c to 25.3 ± 3.3 Mm cm⁻³ (Germany), c_100,m to 7.2 ± 3.7 Mm cm⁻³ (Barbados), x_d to 0.70 ± 0.05 (Cape Verde and Barbados), x_c to 0.94 ± 0.03 (Germany), and x_m to 0.85 ± 0.03 (Barbados). Despite that the CCN calculation method used here is the same with Mamouri and Ansmann [27] their aerosol types (dust, continental and marine) were identified in a different way. They also used the particle linear depolarization ratio to separate between dust and non-dust particles, but they separated non-dust particles into a continental and marine component, using a combination of backward trajectory analysis, Ångström exponent information and marine backscatter estimates.

In fact, the particle number concentrations in Equations (2)–(4) express the CCN concentration values (n_CCN) for the different aerosol types for a supersaturation (ss) of 0.15%. For higher supersaturation, smaller particles may also act as CCN (i.e., particles with radius larger than 40 and 30 nm at supersaturations of 0.25% and 0.40%, respectively). The n_30,dry/n_50,dry and n_40,dry/n_50,dry ratios are estimated at 1.7 ± 0.8 and 1.35 ± 0.7 by Mamouri and Ansmann [27], while similar values are found, e.g., in Ji and Shaw [66] and Shinozuka et al. [21]. These ratios can be used as an enhancement factor f_ss in a generalized CCN equation

$${\mathrm{n}}_{\mathrm{C}\mathrm{C}\mathrm{N}\mathrm{j}}=\hspace{0.17em}{\mathrm{f}}_{\mathrm{s}\mathrm{s}}\cdot \hspace{0.17em}{\mathrm{n}}_{\mathrm{i},\mathrm{j},\mathrm{dry}}$$

(5)

where f_ss is equal to 1.0, 1.35, and 1.7 for supersaturations of 0.15%, 0.25%, and 0.40%, respectively. The application of these enhancement factors may be more complex for natural aerosols [67]_. For dust particles with soluble coatings at a supersaturation of 0.15–0.20%, even particles with a dry radius up to 50 nm may act as CCN. This would lead to four times larger CCN concentrations than the ones estimated using n_100,d,dry. The uncertainties introduced by the chemical composition of aerosols (unknown supersaturation and minimum radius of particles that act as CCN) are expected to lead to an overall uncertainty of a factor of two to three in the CCN estimates [21,27]. Here, we report CCN estimates for a supersaturation of 0.15%. The corresponding measurement uncertainties are calculated using error propagation, and are much lower than the overall uncertainty discussed above.

3. Results and Discussion

3.1. CALIPSO Aerosol Typing

For the scope of this work, a CALIPSO track for 9 September 2011 starting at 00:31 UTC was utilized. This track spans from ~57° N to 0° N, and passes over Thessaloniki, Greece (around 00:45 UTC) where the ACEMED campaign took place. The CALIPSO track (dotted red line) can be seen in Figure 2, along with the different areas of interest that crosses. To have a closer look into the greater Thessaloniki area, a part of the track was examined separately (solid red line in Figure 2).

The backscatter coefficients at 532 nm, along with the corresponding vertical feature mask and aerosol subtypes, can be seen in Figure 3. Unfortunately, over Eastern Europe and Central Africa there are only a few valid measurements (Figure 3a), mostly due to cloudiness, as shown in Figure 3c (the filtering applied on the data prior to the application of our pure dust separation algorithm allows only cloud free scenes). On the contrary, a strong aerosol signal is observed over the Eastern Mediterranean, the Sahara, and the Sahel region (Figure 3b). Mixed aerosols are found over the Eastern Mediterranean, while dusty marine and marine aerosols dominate over water. Over the North African coast (Libya), polluted continental/smoke and polluted dust aerosols prevail, while over the Balkan Peninsula (Greece-Bulgaria), the aerosol situation includes a mixture of polluted continental/smoke, polluted dust, and elevated smoke. Over the Sahara, dust prevails and shifts to polluted dust going south towards the Sahel. Over the Thessaloniki area, elevated smoke and dust dominate over land and elevated smoke, dusty marine, and marine over the sea (coastline at 40.6° N).

Independent data from the CAMS reanalysis are utilized to qualitatively evaluate the CALIPSO aerosol typing. The CAMS provides mass mixing ratio vertical cross sections for different aerosol types such as sea salt, dust, carbonaceous (organic matter and black carbon), and sulfates, along the CALIPSO track. The mass mixing ratio for these aerosol types together and the dominant aerosol type by means of mass mixing ratio are shown in Figure 4. In line with CALIPSO, a mixture of sea salt, dust, and carbonaceous aerosols is found over the Eastern Mediterranean; dust dominates over the Sahara and the Sahel, and carbonaceous aerosols from biomass-burning activities prevail over Central Africa. Over Eastern Europe, where there are only a few CALIPSO observations, we see a mixture of carbonaceous and sulfates, while over ~45° N, sea salt dominates, probably due to transport from the North Sea. Over Thessaloniki, the CAMS shows a plume of carbonaceous aerosols over a layer of dust and sulfates. The presence of carbonaceous aerosols over Thessaloniki and the broader Balkan region is probably due to the existence of several fires in the area the days before. In Figure 5, we see the fire spots which were detected from MODIS/Terra and Aqua satellite sensors (NASA FIRMS) over the area, along with the HYSPLIT calculated 4-day back-trajectories ending over Thessaloniki at 00:00 UTC on 9 September 2011. Indeed, the air masses that arrived over the site had previously passed over locations with fires.

3.2. CALIPSO Extinction Coefficients and CCN Concentrations

Vertical cross sections of CALIPSO extinction coefficients at 532 nm and the corresponding uncertainties for continental, dust, and marine aerosols are shown in Figure 6 for the examined CALIPSO track.

In line with the discussion in the previous section, continental aerosols dominate (values larger than 100 Mm⁻¹) over Eastern Europe, land areas of the Eastern Mediterranean, areas of the Sahara region across the North African coast and close to the Sahel, and over the Sahel. Dust dominates over the Sahara Desert (most values between 30–70 Mm⁻¹) and at altitudes higher than ~1 km a.s.l. over the Mediterranean Sea, while marine aerosols dominate over the Mediterranean Sea (most values between 30–70 Mm⁻¹) at lower altitudes within the boundary layer.

The corresponding vertical cross sections with our calculated CCN concentrations can be seen in Figure 7. Continental CCN concentrations larger than 2000 particles/cm³ (usually reported in cm⁻³) at a supersaturation of 0.15% appear over the areas where continental aerosols dominate. Dust CCN concentrations of 80 to 140 cm⁻³ appear over the areas where dust dominates, while marine CCN concentrations larger than 200 cm⁻³ appear at low altitudes over the Mediterranean Sea, where marine aerosols dominate.

Profiles of CCN concentrations at a supersaturation of 0.15% for the different regions of interest are presented in Figure 8. The uncertainty range of a factor of two is adapted from, e.g., Mamouri and Ansmann [27], Ansmann et al., [28], and Haarig et al., [29]. Only profiles for the Eastern Mediterranean, the Sahara and the Sahel regions are shown as, for the other two regions (Eastern Europe and Central Africa), there were large gaps and, hence, the spatial representativeness of their profiles is very limited. Continental CCNs clearly dominate for the three regions. Over the Eastern Mediterranean, the total CCN concentration is about 1900 cm⁻³ at altitudes around 0.5 km a.s.l., and about 1150 cm⁻³ at altitudes around 2.8 km a.s.l. For altitudes between 1.0 and 2.5 km a.s.l., CCN concentrations are around 1000 cm⁻³. Above 3.0 km a.s.l., CCN concentrations decrease and reach zero at about 4.0 km a.s.l. Continental aerosols contribute 94%, dust 1.5% (values up to 36 cm⁻³), and marine aerosols 4.5% (values up to 158 cm⁻³) to the total CCN concentration. Over the Sahara, the total aerosol CCN concentrations are much larger, reaching values up to 3900 cm⁻³ close to the surface. The CCN concentrations decrease sharply reaching zero at 5.3 km a.s.l. Continental aerosols account for 94% of the total CCN concentrations and dust for the remaining 6% (values up to 79 cm⁻³). Over the Sahel region, the CCN concentration of all aerosol types together reaches values of up to 1080 cm⁻³ at an altitude of 0.6 km a.s.l. and decrease to 540 cm⁻³ at 3.0 km a.s.l., and zero above 4.2 km a.s.l. Continental aerosols account for 91% of the total CCN concentrations and dust for the remaining 9% (values up to 75 cm⁻³).

Our findings are in line with recent ground-based lidar studies. Mamouri and Ansmann [27] derived CCN concentrations up to 2000 cm⁻³ at a supersaturation of 0.15% at low altitudes close to the surface at Cyprus, Eastern Mediterranean. Haarig et al. The authors of [29] found CCN concentrations up to 500 cm⁻³ at a supersaturation of 0.2% at 2.0 to 3.0 km a.s.l. over Barbados, Caribbean. Düsing et al. [65] inferred CCN concentrations up to 2000 cm⁻³ at 1.0 km a.s.l. at a supersaturation of 0.15% at Melpitz, Germany, Central Europe, while Genz et al. [30] retrieved CCN concentrations up to 1000 cm⁻³ close to the surface at the same location. Hofer et al. [31] found concentrations up to 1000 cm⁻³ at a supersaturation of 0.2% close to the surface in Tajikistan, Central Asia. This is a first confirmation that our satellite-based CCN estimates are within a reasonable range of values.

3.3. Evaluation of CALIPSO-Derived CCN Concentrations

The CALIPSO CCN concentration profiles for all the aerosol types together, for dust, continental and marine aerosols at a supersaturation of 0.15% (solid lines), and the corresponding factor of two uncertainty range (transparent lines) for the land (spatial average for 40.6°–41.2° N) and sea areas (spatial average for 40°–40.6° N) around Thessaloniki where the ACEMED took place are shown in Figure 9a,b, respectively. Dry particle number concentrations from the IRRA algorithm retrieved using in situ and lidar remote sensing data from the ACEMED flight [38], are overplotted. The flight tracks used for the IRRA retrievals (spatial representativeness) relative to the CALIPSO track are show in Figure 6 of Tsekeri et al. [38]. To compare our CCN profiles with the IRRA retrievals, the three CALIPSO height levels which are closer to each IRRA height level are defined, and an average CCN value is calculated. It has to be noted that the CALIPSO overpass time (around 00:45 UTC) falls within the IRRA retrievals time (between 00:00 and 02:00 UTC). Over land CALIPSO overestimates CCN concentrations up to four times, relative to the IRRA retrievals, while over the sea it underestimates by ~43%. In the first case, all the IRRA retrievals fall outside the CALIPSO factor of two uncertainty range while in the second case two values out a total of four fall outside. The calculated CCN values along with the IRRA retrievals for the ACEMED flight altitudes appear in

Table 1. Original and hydration-corrected CALIPSO-derived CCN concentrations (in cm⁻³) and dry particle number concentrations (in cm⁻³), retrieved with the IRRA algorithm using ACEMED flight measurements at different height levels over Thessaloniki, Greece. The factor of two uncertainty range is given in parentheses. Columnar CCN concentrations (including particles with radius larger than 30 nm) based on the original and the corrected CALIPSO data and collocated MODIS PSML003_Ocean data (in cm⁻²) are also given.

Area	Altitude	CALIOP (Original)	CALIOP (Hydration-Corrected)	In Situ	MODIS
Land	2.1 km	1816 (908,3632)	1504 (752,3008)	727	-
Land	2.7 km	4505 (2253,9010)	2851 (1426,5702)	1318	-
Land	3.2 km	6370 (3185,12740)	2086 (1043,4172)	779	-
Sea	1.3 km	609 (305,1218)	508 (254,1016)	1427	-
Sea	2.1 km	1683 (842,3366)	1405 (703,2810)	1834	-
Sea	2.7 km	1264 (632,2528)	912 (456,1824)	1501	-
Sea	3.2 km	794 (397,1588)	459 (230,918)	2814	-
Sea	column	5.51 × 108 (2.76 × 108,11.02 × 108)	4.37 × 108 (2.19 × 108,8.74 × 108)	-	7.27 × 108

.

As the parameterizations for the calculation of CCN concentrations assume dry particles, it is necessary to correct for hydration effects, i.e., the uptake of water by hydrophilic aerosols. Aerosols observed during the ACEMED campaign experienced hygroscopic growth that affected their optical properties [38]. Over land areas, the relative humidity (RH) reached values larger than 90%, whereas, over the sea, RH reached values of 80% (not shown here), as seen from the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) meteorological profiles, which are included in the CALIPSO v4.1 product [68]. Significant deviations can be expected between the extinction coefficients of ambient and dry particles. In order to correct for the hydration effect, we use dry-to-ambient extinction coefficient ratios (DAR) for different RH values from Tsekeri et al. [38]. A third-degree polynomial is used to fit the DAR and RH data and applied thereafter on the MERRA-2 RH data. The original (ambient) CALIPSO extinction coefficients at 532 nm are multiplied by the emerging DAR values to calculate “dry” CALIPSO extinction coefficients at 532 nm, which are subsequently used in Equations (2)–(5) to obtain the CCN concentrations that are now corrected for the hydration effect.

The original and the corrected total CCN concentrations are shown in Figure 9 for the land and sea areas around Thessaloniki, respectively. Over land, the CALIPSO estimates overestimate the in-situ retrievals by a factor of two. This changes to an underestimation of 57% over sea. In the first case, one IRRA retrieval falls inside the CALIPSO factor of two uncertainty range and a second one marginally falls outside this range. In the second case, two values out of a total of four fall outside (one of them marginally). The corrected CCN values for the ACEMED flight altitudes are also given in Table 1. Overall, the correction leads to a significant decrease of 110% of the retrieved CCN concentrations over land and a smaller decrease of 22% over the sea, due to the overall lower RH. Six out of seven points from the in-situ measurements agree with the CALIPSO estimates within a factor of three.

We also compare our results with MODIS-derived columnar concentration (in cm⁻²) of fine mode particles with radius larger than 30 nm. More specifically, the spatial average of MODIS/Terra and Aqua pixels, whose center falls within a half-degree window centered at the ACEMED sea area is used for the comparison. To calculate columnar concentrations, CALIPSO-derived CCN concentrations are integrated over the atmospheric column. However, as our results are reported at a supersaturation of 0.15% assuming that particles with a radius larger than 50 nm may act as CCN, the columnar concentrations were multiplied by a factor of 1.7 that expresses the ratio of the n_30,dry and n_50,dry concentrations [27] in the comparison to the merged MODIS/Terra and Aqua data. The CALIPSO-derived column CCN concentrations are smaller than the ones inferred from MODIS. However, the MODIS values lie within the factor of two uncertainty range of the CALIPSO values as listed in Table 1. This comparison cannot be considered a validation as the sample is very limited, the MODIS dataset itself is not yet validated and the two datasets are not temporally collocated (MODIS/Terra and Aqua sensors pass over Thessaloniki around 09:30 and 11:30 UTC, respectively). However, the fact that two completely different satellite-based retrieval schemes applied on measurements from a passive and an active sensor give comparable results is promising.

4. Summary and Conclusions

This work investigates the potential of using spaceborne CALIPSO lidar measurements to obtain profiles of CCN concentrations as suggested by Mamouri and Ansmann [27] for ground-based lidar measurements. A CALIPSO measurement on 9 September 2011 over the Thessaloniki region (in northern Greece) is used for an assessment of the CCN retrieval based on independent in-situ observations with the FAAM Bae-146 research aircraft during the ACEMED campaign. The CALIPSO extinction coefficient is separated into contributions of dust and non-dust aerosols. The presence of biomass burning aerosols over Thessaloniki is corroborated by data from the CAMS aerosol reanalysis, and the combined inspection of MODIS fire counts and HYSPLIT backward trajectories.

The extinction coefficients related to mineral dust, continental aerosols, and marine aerosols are used to estimate CCN concentrations at a supersaturation of 0.15%. Continental CCN concentrations larger than 2000 cm⁻³ appear over areas where continental aerosols dominate, dust CCN concentrations of 80 to 140 cm⁻³ appear over dust dominating areas, while marine CCN concentrations larger than 200 cm⁻³ appear at low altitudes over the Mediterranean Sea. Continental aerosols account for more than 90% of the total CCN concentrations over all the examined regions. The values appearing here are of the same magnitude with values appearing in recent ground-based lidar studies giving confidence to our calculations.

Our CCN concentrations are evaluated against dry particle number concentrations retrieved using in situ and lidar remote sensing measurements from a FAAM BAe-146 flight over Thessaloniki during the ACEMED campaign. We suggest here that it is necessary to apply a correction for the hydration effect. To do so we used dry-to-ambient extinction coefficient ratios for different RHs from Tsekeri et al. [38]. Without such a correction, CALIPSO will tend to overestimate CCN concentrations, especially at RHs larger than 70–80%. Most of the CALIPSO-derived CCN concentrations are found to match the aircraft measurements within a factor of two. A comparison to the unvalidated MODIS CCN product over Thessaloniki shows smaller values for the CALIPSO estimate, again with an agreement within a factor of two.

Overall, this first effort shows that it is feasible to obtain CCN concentrations using CALIPSO observations with an uncertainty of a factor of two to three. The need for more detailed validation studies with measurements at various altitudes from other campaigns and improved parameterizations is acknowledged in order to produce in the future a robust climatological product that will be suitable for the evaluation of climate models, for intercomparison with other satellite-based products, and to constrain global ACI estimates.