4.1. Detection of CCN and INP
Estimating concentrations of cloud condensation nuclei (CCN) is crucial for comprehending aerosol-cloud interactions. However, in situ observations of CCN are scarce, and many passive remote sensing methods can only offer proxies like total aerosol optical depth (AOD) for column-effective CCN assessments. [
148].The potential of lidars in the study of aerosol-cloud interactions should be clear, and in this section we will present some findings from observational studies using lidars to characterize the presence of CCN and INP.
A method to extend the CCN spectrum measured at the surface to high altitudes was proposed by Ghan and Collins [
149]. The lidar measured value
at altitude
z where aerosol is exposed to the RH(z), is recalculated to the value it would have in dry conditions using an independent measurement of the vertical profile of RH and surface measurements made with a humidified nephelometer providing the dependence of the extinction on relative humidity. A light scattering hygroscopic growth factor
is defined as the ratio between the extinction coefficient at a various RHs and the extinction coefficient at dry conditions. The same factor is used also to quantify the amount of change in the particle backscattering coefficient due to water uptake:
. Then, surface measurements of the CCN spectrum CCN(
) are scaled by the ratio of the backscatter (or extinction) profile
to the backscatter (or extinction) retrieved at or near the surface,
. This method assumes that the hygroscopic growth factor
measured at the surface is the same as the one
measured at z, and this implies the assumption that the vertical structure of CCN concentration is identical to the vertical structure of dry extinction or backscatter. These assumptions are equivalent to require that both the PSD shape (but not the total aerosol number) and the aerosol composition and particle shape are independent of altitude. In fact, Ghan
et al. [
150] showed that vertical inhomogeneity in the PSD and presumably in particle shape and composition are the mail sources of error in such CCN retrieval, as near-surface CCN properties could be significantly different from CCN properties near the cloud base.
The ability of aerosols to act as CCN depends more on their size rather than chemical composition or mixing state [
151]. This allows to give an assessment on CCN through a determination of the aerosol PSD and an assesment of the number of particles whose dimension is above a certain threshold. It follows that in order to derive CCN concentration one should first use lidar derived optical parameters to retrieve the aerosol PSD. Then an assessment of particle hygroscopicity for different types of aerosol is needed to assess the ability of particles to act as CCN. In should be noted that high RH near clouds can change the aerosol optical properties and complicate the CCN retrieval.
Lv
et al. [
152] developed a method for profiling CCN concentrations that usees backscatter coefficients at 355, 532, and 1064 nm and extinction coefficients at 355 and 532 nm. Three types of aerosols (urban industrial, biomass burning, and dust) are considered, each with a different bimodal PSD. This PSD was inferred by using look-up tables developed based on the aerosol PSD database of the Aerosol Robotic Network (AERONET) database [
153]. The best bimodal size distribution parameters were selected by comparing lidar observations and Mie optical computation [
59] on the aerosol PSD, varying the PSD parameters through typical ranges for the three types of aerosols, until an agreement between measurement and calculations is reached. This "wet" PSD, measured at ambient RH, is then corrected to its "dry" version by using a hygroscopic scattering enhancement factor obtained from Humidified Tandem Differential Mobility Analyzer (HTDMA) or Raman lidar data. Then the aerosol critical radius (
) for CCN activation at selected critical supersaturation is computed from the maximum of the
-Köhler curve [
23]. As, according to
-Köhler theory,
depends on the hygroscopicity parameter
which changes according to aerosol type, such critical radius (
) depends on the aerosol type and for the three types of aerosols in the study, it is taken from the literature [
154,
155,
156]. CCN number concentrations are then determined from the integration of the retrieved PSD, from
upward.
In 2019, Tan
et al. [
157] similarly proposed to retrieve CCN concentrations from multi-wavelength Raman lidar measurements, but instead of using AERONET datasets and estimations of hygroscopicity according to aerosol classification, they relied on in-situ measured PSD, mixing states, and chemical compositions data to define the relationship between CCN concentrations and lidar-derived optical properties. Hygroscopic enhancements of backscatter and extinction with relative humidity have been used to create humidograms, to derive both dry backscatter and extinction and hygroscopicity at different wavelengths. These, together with lidar color ratios, extinctions and backscatters data have been used as input to a Random Forest Regression machine learning algorithm that produces the best estimate for the ratio CCN-to-extinctions. Optical data simulated with Mie computations on in-situ-measured PSD and
-Köhler curve from in-situ measured chemical compositions are used as the training data. A similar in-situ dataset is used as test data.
Lenhardt
et al. [
158] explored the connections between aerosol backscatter and extinction coefficients using the airborne High Spectral Resolution Lidar 2 (HSRL-2) in biomass burning aerosol (BBA)-influenced air masses over the southeast Atlantic. They also examined in situ measurements of cloud condensation nuclei (CCN) concentrations that were spatiotemporally collocated. To ensure accuracy and avoid detecting swollen, highly hygroscopic aerosols that could artificially inflate backscatter and extinction values without corresponding increases in CCN concentration, observations taken at relative humidity (RH) greater than 40% were filtered out. Their findings, illustrated in
Figure 5, demonstrate robust linear relationships between lidar-derived backscatter and extinction coefficients and CCN concentration.
Figure 5.
CCN concentration versus HSRL-2 (a) backscatter and (b) extinction coefficients with blue scatter points representing 355 nm and red scatter points representing 532 nm. This combined data set represents 10 d of observations and 80 total collocated data points (per coefficient), covering all 3 years of ORACLES. Supersaturation for these observations ranges between 0.22 %–0.4 %. The Pearson correlation coefficient is shown, with the Spearman rank correlation coefficient given in parentheses. Error bars are given for a CCN relative uncertainty of 10% and for calculated HSRL-2 uncertainties. Lines of best fit are forced through the origin to represent the practicality of using linear regression equations to quantitatively obtain CCN concentrations using HSRL-2 observables. (Figure and caption from Figure 8 in Lenhardt
et al. [
158] licensed under CC BY 4.0 )
Figure 5.
CCN concentration versus HSRL-2 (a) backscatter and (b) extinction coefficients with blue scatter points representing 355 nm and red scatter points representing 532 nm. This combined data set represents 10 d of observations and 80 total collocated data points (per coefficient), covering all 3 years of ORACLES. Supersaturation for these observations ranges between 0.22 %–0.4 %. The Pearson correlation coefficient is shown, with the Spearman rank correlation coefficient given in parentheses. Error bars are given for a CCN relative uncertainty of 10% and for calculated HSRL-2 uncertainties. Lines of best fit are forced through the origin to represent the practicality of using linear regression equations to quantitatively obtain CCN concentrations using HSRL-2 observables. (Figure and caption from Figure 8 in Lenhardt
et al. [
158] licensed under CC BY 4.0 )
These studies employ multi-wavelength Raman or HRSL lidars, commonly called 3
+ 2
lidars, systems of a certain complexity, coupled with in-situ measurement and RH profiling by means of radiosoundings or remote sensing. A simplified approach was suggested by Mamouri and Ansmann [
159] who proposed an inversion algorithm for data from a more manageable single-wavelength polarization diversity lidar. This lidar is still able to specify aerosol classes (desert, nondesert continental, and marine). The possible presence of external mixing in the aerosol is handled with the procedure outlined in Tesche
et al. [
160]. Following a methodology proposed by Shinozuka
et al. [
161], extensive datasets from AERONET and lidar correlated observations are then used to connect lidar derived extinction to total particle number concentrations (for dry particle dimensions above defined thresholds depending on aerosol class) for the three aerosol classes. An evaluation of the correction for water-uptake [
23,
162] is performed, this time assuming fixed RH values for the typical conditions of observation. Finally, a simple parametrization is used to connect the particle concentrations to the CCN concentrations, obtained from the formers with scaling factors dependent on supersaturation and aerosol class. Thus, height profiles of CCN concentrations can be retrieved from lidar-derived ambient aerosol extinction. This approach has lent itself to derive global climatologies of CCN from the analysis of the satellite borne lidar CALIOP [
163]. A similar approach was also used by Choudhury and Tesche [
164]. In their study, they employed normalized volume size distributions and refractive indices based on the CALIPSO aerosol model [
80] for six aerosol types identified by the CALIPSO lidar. These size distributions were adjusted until agreement was reached between the extinction coefficient inferred from CALIPSO measurements and that calculated through light-scattering computations. These modified size distribution were then used to compute the aerosol number concentration for particles with dimension above defined thresholds depending on aerosol type. Again, the CCN concentration at a certain set of supersaturations was estimated by multiplying the aerosol number concentration with scaling factors which depend on the aerosol type and on the level of supersaturation.
Lidars have been utilized to retrieve profiles of INP concentrations. This is achieved by integrating particle concentration profiles derived from lidar measurements with parameterizations of INP efficiency specific to different aerosol types and freezing mechanisms (such as immersion, condensation, deposition, or contact freezing). Mamouri and Ansmann [
159] apply a regression of lidar derived extinction vs SAD to retrieve the INP concentration from the latter. Indeed, precise knowledge of aerosol type is crucial for using lidar retrievals to estimate INP concentrations. This is because INP concentrations are inferred solely from physical properties such as particle number concentration and/or size distribution (SAD), using parameterizations that have been developed for specific aerosol types. Examples include dust [
165,
166,
167,
168], marine aerosols [
169], soot [
167], and global aerosols [
170]. Therefore, studies of INPs have typically focussed on estimating INP concentrations within these specific aerosol classes.
Using a similar methodology, Haarig
et al. [
171] presented vertical profiles of cloud condensation nuclei (CCN) number concentration, particles with diameter greater than 500 nm, size distribution, mass, and INP concentration. These profiles were derived from measurements of 3
+ 2
Raman lidar with polarization diversity. They compared these measurements with in situ CCN concentrations and INP-relevant aerosol properties collected by aircraft in the Saharan Air Layer (SAL) over the Barbados region.
Extinction coefficient profiles were separately retrieved for mineral dust, marine, and continental aerosols. Empirical conversion factors [
172] were applied to convert these extinction coefficients into particle number concentrations (for particles above a threshold size dependent on aerosol type) and size distributions (SADs). Subsequently, various INP parameterizations [
166,
170] were tested, using particle number concentration and SAD as inputs.
Comparisons with in situ data of mass concentrations and particle number concentrations, which were used as inputs for the INP parameterizations, demonstrated good agreement.
Similarly, Marinou
et al. [
173] have retrieved cloud-relevant particle number concentrations (i.e particles whose linear dimension in dry conditions are above 250 nm) and SAD
using lidar measurements froma a
+ 2
Raman lidar with polarization diversity. INP concentration profiles were estimated using various ice nuclei parameterizations. These lidar-derived results were subsequently compared with direct INP measurements obtained by sampling aerosols along the lidar profile using two UAVs equipped with INP samplers. The collected samples were then analyzed offline using the FRIDGE (FRankfurt Ice nucleation Deposition freezinG Experiment) INP counter [
174]. This approach allowed them to evaluate the effectiveness of different INP parameterizations across different temperature ranges and types of particles.
Wieder
et al. [
175] also provided a direct validation of the INP concentration retrievals. They tested INP retrievals based on data from a
+ 2
Raman lidar with polarization diversity by comparing them with in situ observations of aerosols and INP taken at a nearby mountain site in the Swiss Alps. The sampled air masses predominantly contained Saharan dust and continental aerosols. Various INP parameterizations were also evaluated in this study.
4.2. Impact of Aerosol on Mixed and Cirrus Clouds
Lidars are very effective in studying the impact of particular types of aerosols on the microphysics of clouds. Choi
et al. [
176] demonstrated the ability to distinguish supercooled liquid clouds from ice clouds for the satellite-borne lidar CALIOP from the layer-integrated particle depolarization ratio and backscatter coefficient at 532nm, together with the cloud top and bottom temperatures, and demonstrated an inverse correlation between supercooled cloud presence and dust presence. Utilizing its capability to discriminate cloud phases, Tan
et al. [
177] directly investigated the ice-nucleating potential of dust, polluted dust, and smoke aerosols in mixed-phase clouds. They analyzed vertically-resolved profiles of cloud thermodynamic phase and aerosols from global spaceborne lidar data spanning 5 years. Their findings revealed that the presence of dust aerosols, both in their clean and polluted forms, in various regions globally at temperatures conducive to mixed-phase clouds, reduces the fraction of supercooled clouds compared to total cloud cover in those regions. This reduction was attributed to the ability of dust aerosols to nucleate ice crystals.
The influence of desert dust aerosols on cirrus microphysical properties was investigated using concurrent and co-located datasets from CALIOP, CloudSat, and MODIS over the Taklimakan Desert [
178]. The study highlighted the negative "Twomey effect" under arid conditions [
56], where cloud albedo decreases due to increased heterogeneous nucleation. This process results in fewer but larger ice crystals with higher fall velocities.
The role of wildfire smokes on cirrus formation was investigated by Mamouri
et al. [
179] who showed strong evidence that long range transport of aged smoke (organic aerosol particles) originating from wildfires triggered significant ice nucleation causing the formation of extended cirrus layers, upon gravity wave activity close to the tropopause.
The impact of aerosols on cloud thermodynamic phase was also addressed by Zhang
et al. [
180] over east Asia by combining the 4 years datasets of CloudSat radar and CALIOP lidar measurements. Although temperature differences at cloud top showed to be the most important drivers of the cloud thermodynamic phase, regional differences and seasonal anomalies of glaciated and mixed-phase relative cloud fractions correlated well with variations of dust occurrence frequency. Moreover, the relative frequency of glaciated clouds associated with concomitant presence of mineral dust was higher than the frequency of glaciated clouds associated with presence of polluted dust, smoke, and background aerosols at any given cloud top temperature.
Hofer
et al. [
181] investigated the efficiency of heterogeneous ice formation across stations located in New Zealand (Lauder), Germany (Leipzig), and southern Chile (Punta Arenas), as influenced by cloud-top temperature. They observed that Lauder and southern Chile, which typically experience low concentrations of free-tropospheric aerosols, exhibit lower ice formation efficiency compared to polluted mid-latitude regions like Germany. This suggests that the reduced ice formation efficiency at Lauder and Punta Arenas is linked to low concentrations of INP. Additionally, episodes of continental aerosol outbreaks, more frequent at Lauder than Punta Arenas, moderately enhance the ice formation efficiency at Lauder relative to Punta Arenas.
The diurnal cycle of the supercooled water cloud fraction was investigated by Wang
et al. [
182] using 33 months of lidar observations from the Cloud-Aerosol Transport System aboard the International Space Station. This system provides cloud top phase information at various local times within specific grid and isotherm parameters. The study revealed a strong and statistically significant negative correlation between the dust aerosol extinction coefficient near the location and time of liquid/ice cloud footprints, and the fraction of supercooled water clouds. Specifically, higher dust aerosol extinction coefficients corresponded to lower supercooled water cloud fractions in the observed regions.
4.3. Impact of Aerosol on Warm Clouds
Studies investigating the influence of aerosols on warm clouds predominantly employ dual-FOV lidar techniques. Schmidt
et al. [
183] utilized colocated dual-FOV Raman lidar observations to examine aerosol and cloud optical and microphysical properties beneath and within thin layered liquid water clouds. They complemented these observations with Doppler lidar measurements of updrafts and downdrafts at cloud base to explore the relationship between aerosol concentrations near cloud base and cloud characteristics and dynamics.
The dual-FOV lidar setup enabled the use of two multiple-scattering channels (elastic and nitrogen Raman multiple-backscatter channels) to profile single-scattering extinction coefficient, effective radius, cloud droplet number concentration, and liquid water content [
133]. The study included two case studies that tracked the evolution of altocumulus clouds in clean and moderately polluted conditions.
The impact of updrafts, downdrafts, turbulent mixing, and entrainment of dry air on the microphysical properties of layered clouds was investigated using a combination of Doppler lidar and cloud radar. Significant differences in cloud properties were documented during updraft and downdraft episodes, particularly in droplet extinction, number concentration, and effective radius. Updraft episodes showed a notable increase in extinction coefficient and droplet number concentration, accompanied by lower droplet effective radii attributed to new droplet formation. The liquid water content (LWC) profile retrieved during updraft periods closely resembled the adiabatic LWC profile, whereas during downdraft episodes, higher LWC values indicated an absence of dry intrusions.
Conversely, signs of dry entrainment during downdrafts—such as lower effective radii and reduced LWC values observed around the cloud center—were documented in another case study. This scenario occurred when weaker vertical motions facilitated downward mixing of dry air from above the shallow cloud layer.
Various studies highlight the challenges of studying aerosol-cloud interactions (ACI) in turbulent environments, where turbulent motions can influence processes like evaporation, drop collision, and new droplet formation, leading to opposing effects on cloud microphysical parameters [
184]. Schmidt
et al. [
185] conducted a statistical analysis spanning two years, demonstrating a clear aerosol impact on cloud evolution and cloud droplet number concentration in the lower portions of altocumulus layers during updrafts. Their analysis underscored the importance of considering cloud dynamics in assessing ACI parameters, emphasizing the need to incorporate vertical wind information.
Furthermore, a comprehensive review of contemporary field studies on aerosol-cloud interactions in warm clouds (depicted in
Figure 6) supported findings that Aitken and accumulation-mode particles are activated at cloud base when rapidly uplifted, irrespective of whether clouds form over oceans or continents. This dynamic contrasts with passive satellite remote sensing, which typically yields lower estimates of
compared to ground-based lidar and airborne observations, potentially due to differences in cloud-top dynamics.
Jimenez
et al. [
186] conducted cloud measurements in pristine marine conditions at Punta Arenas in southern Chile using a multiwavelength polarization Raman lidar with dual-field-of-view (FOV), coupled with a Doppler lidar for vertical wind component measurements. They performed a detailed study on aerosol-cloud interactions (ACI), resolving updrafts and downdrafts. The study utilized profiles of aerosol-related parameters near cloud base, cloud microphysical properties just above cloud base, and 1-minute resolution vertical wind data. CCN number concentration was derived from lidar-measured aerosol extinction following Mamouri and Ansmann [
159]. They observed high
values ranging from 0.8 to 1.0. The study highlighted the impact of aerosol water uptake on ACI, noting that the highest
values were obtained when considering aerosol extinction measurements taken approximately 500 meters below cloud base (indicating dry aerosol conditions) in ACI computations.
Wang
et al. [
187] utilized a dual-FOV HSRL to simultaneously profile aerosols and liquid water clouds clouds, focusing on assessing the adiabaticity of these low-level clouds. They found in some case that observed profiles of microphysical properties in these clouds are not perfectly adiabatic, contrary to common assumptions in current retrieval methods [
35,
188,
189]. Additionally, the study confirmed that increased aerosol loading leads to higher droplet number concentrations and reduced droplet effective radius, although no discernible increase in liquid water path (LWP) was detected. This suggests that aerosol-induced water reduction through enhanced evaporation may offset increases caused by suppressed rain formation. The lidar observations were validated through concurrent measurements with cloud radar, Raman lidar, and sun photometer instruments.