Articles | Volume 18, issue 19
Research article
09 Oct 2018
Research article |  | 09 Oct 2018

Ice crystal number concentration estimates from lidar–radar satellite remote sensing – Part 2: Controls on the ice crystal number concentration

Edward Gryspeerdt, Odran Sourdeval, Johannes Quaas, Julien Delanoë, Martina Krämer, and Philipp Kühne

The ice crystal number concentration (Ni) is a key property of ice clouds, both radiatively and microphysically. Due to sparse in situ measurements of ice cloud properties, the controls on the Ni have remained difficult to determine. As more advanced treatments of ice clouds are included in global models, it is becoming increasingly necessary to develop strong observational constraints on the processes involved.

This work uses the DARDAR-Nice Ni retrieval described in Part 1 to investigate the controls on the Ni at a global scale. The retrieved clouds are separated by type. The effects of temperature, proxies for in-cloud updraft and aerosol concentrations are investigated. Variations in the cloud top Ni (Ni(top)) consistent with both homogeneous and heterogeneous nucleation are observed along with differing relationships between aerosol and Ni(top) depending on the prevailing meteorological situation and aerosol type. Away from the cloud top, the Ni displays a different sensitivity to these controlling factors, providing a possible explanation for the low Ni sensitivity to temperature and ice nucleating particles (INP) observed in previous in situ studies.

This satellite dataset provides a new way of investigating the response of cloud properties to meteorological and aerosol controls. The results presented in this work increase our confidence in the retrieved Ni and will form the basis for further study into the processes influencing ice and mixed phase clouds.

1 Introduction

Clouds play a central role in the Earth's energy budget, as they are responsible for large variations in the reflected shortwave and emitted longwave radiation (Stephens et al.2012). The response of clouds to changing greenhouse gases and aerosols remains one of the largest uncertainties in understanding past and future climate changes (Boucher et al.2013). Significant advances have been made into modeling and observing the role of aerosols in liquid clouds (e.g. Wang et al.2011; Wood et al.2011; Seifert et al.2015; Gettelman2015; Ghan et al.2016; Zuidema et al.2016), especially through the use of retrievals of the cloud droplet number concentration (e.g. Quaas et al.2008; Gryspeerdt et al.2016), but the impact of aerosols on high clouds remains uncertain (Gettelman et al.2012; Jensen et al.2016; Zhou et al.2016; Heyn et al.2017). A large part of this uncertainty comes from the difficulty in retrieving cirrus cloud properties at a large enough scale to separate the roles of individual factors controlling the ice crystal number concentration (Ni).

A key microphysical property of ice clouds, the Ni links the aerosol environment to dynamic effects driving cloud updrafts and the generation of supersaturation (Pruppacher and Klett1997). Through changes in the ice crystal size, changes in the Ni can have far-reaching implications for a cloud, impacting the radiative (Liou1986; Fusina et al.2007), precipitation and cloud lifetime properties (Lindsey and Fromm2008). The Ni is often used as a prognostic variable in two-moment cloud microphysics schemes (e.g. Lohmann et al.2007; Salzmann et al.2010). This highlights a requirement to understand the controls on the Ni in order to improve our understanding and parameterization of cloud processes. While aircraft measurements of the Ni exist, they are restricted in space and time. They can also be affected by shattering of ice crystals at the instrument inlet (McFarquhar et al.2007; Jensen et al.2009; Korolev et al.2013) and difficulties in measuring the smallest crystals (O'Shea et al.2016). In this paper, the new DARDAR-Nice satellite dataset described in Part 1 (Sourdeval et al.2018a) allows the processes that control the Ni to be investigated globally.

It is known that the temperature plays a strong role in determining the ice crystal nucleation rate. The homogeneous nucleation rate is a strong function of temperature and supersaturation (Koop et al.2000), with atmospherically relevant nucleation only taking place at temperatures colder than 235 K. This strong temperature dependence in the nucleation rate does not necessarily correspond to a strong temperature dependence in the Ni (Heymsfield and Miloshevich1993). A weak Ni temperature dependence was found by Gayet et al. (2004). Krämer et al. (2009) found similar results, with a slight reduction in the Ni for the coldest measurements. Higher Ni values have been observed at cold temperatures during ATTREX (Jensen et al.2013a, 2016) than in Krämer et al. (2009), leading to a weak combined temperature dependence. However, using different datasets targeting different cloud types, Muhlbauer et al. (2014) and Jensen et al. (2013b) both showed an increase in Ni with decreasing temperature, demonstrating that there is still considerable uncertainty regarding the Ni temperature dependence.

The in situ homogeneous nucleation of ice crystals is also dependent on the supersaturation (Koop et al.2000; Lohmann and Kärcher2002), which is often generated through cooling due to vertical air motion. Large-scale updrafts cannot reproduce observed cirrus properties on their own, the smaller scale variation in updraft provided by gravity waves is necessary (Kärcher and Ström2003) and is occasionally able to produce cirrus in regions of large-scale subsidence (Muhlbauer et al.2014). These small-scale updrafts can produce Ni values as high as 50 000 L−1 (Hoyle et al.2005), highlighting the important role that vertical motion can play in determining the Ni.

Although ice can form by in situ nucleation, many ice crystals also form through the freezing of liquid condensate. This liquid-origin cirrus often originates from high updraft regions in mixed-phase clouds, forming thicker cirrus than those composed of in situ ice (Krämer et al.2016; Luebke et al.2016). Synoptic-scale updrafts can also produce liquid-origin cirrus in the midlatitudes (Wernli et al.2016). The Ni formed through these liquid-origin processes is also strongly dependent on the cloud-scale updrafts, with higher updrafts maintaining higher ice supersaturations and producing larger Ni values (Kärcher and Seifert2016; Kärcher2017).

Aerosol also plays a role in determining the Ni, although its impact is complicated by variations in ice crystal nucleation pathways and aerosol properties. While any particle can theoretically act as a homogeneous nucleation center given a high enough supersaturation, in practice these aerosols are often hydrophilic liquid aerosols (Kojima et al.2004). Increases in the aerosol number can result in an increase in Ni through increased homogeneous nucleation. However, in many situations, the Ni is limited by dynamical concerns, which limits the impact of aerosols (DeMott et al.1997; Lohmann and Kärcher2002; Kay and Wood2008; Barahona and Nenes2008; Jensen et al.2013a, 2016).

A second class of aerosols, known as ice nucleating particles (INP) are able to nucleate ice heterogeneously and can freeze liquid water droplets at temperatures warmer than 235 K. At these warmer temperatures, the presence of INP will often control the Ni near nucleation locations (Kärcher and Lohmann2003), as they form the sites for creating an ice crystal, although the nucleating ability of these INP is also a strong function of temperature (Hoose and Möhler2012). As heterogeneous nucleation can take place at a lower supersaturation than homogeneous nucleation, the introduction of INP has the ability to prevent homogeneous nucleation by depressing the supersaturation. As the Ni produced through homogeneous nucleation events is typically higher than the INP concentration (and so the Ni from heterogeneous nucleation), this implies that the introduction of INP into a clean atmosphere can reduce the Ni (DeMott et al.1997; Kärcher and Lohmann2003). In situ (Gayet et al.2004) and satellite studies (Chylek et al.2006) have provided some evidence for a possible decrease in Ni with increasing aerosol based on regional and hemispheric differences in ice crystal properties, although it has proved difficult to conclusively link these Ni changes to a change in INP.

The relative role of heterogeneous and homogeneous nucleation in the atmosphere is unclear, making it difficult to develop observational constraints on the impact of aerosols on the Ni (e.g. Cziczo et al.2013; Gasparini and Lohmann2016; Kärcher and Seifert2016; Jensen et al.2016). In addition, changing conditions over the life cycle of a cloud can result in a switch between nucleation mechanisms (Krämer et al.2016); nucleation is also not the only control on the Ni. The rarity of INP suggests that other processes, such as ice multiplication, are required to explain the Ni observed in the lower atmosphere (Heymsfield et al.2017).

These four factors (temperature, supersaturation/updraft, ice origin and aerosol environment) are all thought to influence the Ni in high clouds, but significant uncertainties remain in assessing the role of these factors on the Ni. First, although they have been investigated using aircraft and balloon measurements (Podglajen et al.2016, 2017), the ice origin and in-cloud updraft are difficult to determine from observations at a global scale and over a significant period of time. A recent classification (Gryspeerdt et al.2018a) has shown some skill at determining these quantities when compared to a convection permitting model and is used in this work to account for this issue.

Second, the Ni is a difficult property to measure at a global scale. Aircraft measurements are limited in space and time and have been afflicted by shattering of crystals on the tips of measurement probes, casting doubt on some earlier measurements of the Ni (McFarquhar et al.2007; Jensen et al.2009; Korolev et al.2013). Additionally, due to the highly variable nature of cirrus clouds and their strong sensitivity to environmental conditions, it can be difficult to separate the relative roles of aerosol and dynamics (Gayet et al.2004).

The retrieval presented in Part 1 of this work (Sourdeval et al.2018a) has demonstrated that the Ni can be retrieved using a combined radar–lidar retrieval, and that this compares well to new in situ aircraft measurements where shattering is accounted for. Combined with simultaneous retrievals of the ice water content, this allows the global distribution of the Ni and the factors that control it to be investigated. Using reanalysis aerosol concentrations and a proxy for the INP concentration, the impact of aerosols on high clouds is also investigated, highlighting avenues for future research into cirrus cloud processes.

2 Methods

The Ni dataset used in this work (DARDAR-Nice) has been described in detail in Part 1 of this work (Sourdeval et al.2018a), so only the main features are outlined here. The DARDAR-Nice product is based on the DARDAR retrieval (Delanoë and Hogan2010), a combined raDAR–liDAR retrieval of ice cloud water content (IWC) and ice crystal effective radius using data from the CloudSat and CALIPSO satellites at approximately 13:30 LST (local solar time). Only daytime data from the period from 2006 to 2013 is used in this work due to constraints in the reanalysis data availability. The properties are retrieved at a horizontal resolution of 1.7 km and a vertical resolution of 60 m. Both the DARDAR IWC and the DARDAR-Nice Ni retrieval compare favorably to in situ aircraft data (Deng et al.2013; Sourdeval et al.2018a). The best agreement is at temperatures below −30C, where the retrievals are more accurate due to the dominance of unimodal ice crystal size distributions and reduced ambiguity in the cloud phase.

To investigate the controls on ice crystal nucleation, a more in-depth study is performed of the Ni near the cloud top (Ni(top)). As the cloud top is the location of the coldest temperature in the cloud, it has the highest theoretical nucleation rates. Although the cloud top is close to the nucleation region in wave clouds (e.g. Spichtinger and Gierens2009a), this is not always the case and the Ni(top) can be rapidly reduced by differential sedimentation and entrainment (e.g. Jensen et al.2013a). However, as the coldest temperature, it provides a limitation on the maximum nucleated Ni within the cloud, restricting the impact of temperature variability due to the vertical extent of the cloud. The cloud top is taken to be the top 120 m of the cloud and only the uppermost cloud layer (in multilayer situations) is used to avoid issues with ice being seeded by ice crystals falling from overlying layers. The data are also restricted to locations where the retrieval has gone through at least two iterations, limiting the impact of prior assumptions about the cloud structure.

Four main controls on the Ni(top) are considered in this work: temperature, cloud-scale updraft, ice origin and aerosol. Temperature data in this study are taken from the ECMWF ERA-Interim reanalysis (Dee et al.2011). Information about the cloud-scale updraughts and the ice-origin (liquid/ice) cannot yet be directly obtained at a global scale using satellites. To provide an indication of these cloud properties, the IC-CIR classification from Gryspeerdt et al. (2018a) is used. This classification is based on the assumed cirrus source (orographic, frontal, convective or synoptic) and determined at 13:30 LST using cloud retrievals from the MODIS instrument (Platnick et al.2017) and reanalysis data. This classification selects orographic clouds by assuming the product of the topographic variation and the wind speed is related to the in-cloud updraft, similar to global climate model parameterizations (e.g. Joos et al.2008). The Oro 2 and Oro 1 regimes are the regimes with the respective highest and second highest sextiles of the parameterized in-cloud updraft. Frontal and convective regimes are selected as connected regions of high level cloud that intersect with reanalysis fronts and regions of large-scale updraft, respectively. Finally, the synoptic regime is taken as a residual, with clouds being considered synoptic if they do not fit any of the other classes. Through comparisons with convection permitting model data and classifications based on reanalysis data, this classification has been shown to provide useful information on the cloud scale updrafts and the ice origin. While not a direct retrieval of these properties, it does allow some inferences to be made regarding the response of the Ni to these factors. The results in this paper are restricted to daytime data, which in turn restricts it to 13:30 LST due to the orbit of the satellites used to construct the DARDAR-Nice dataset.

To investigate a possible aerosol link to Ni, we use the “monitoring atmospheric composition and climate” (MACC) aerosol reanalysis product (Morcrette et al.2009; Benedetti et al.2009), which assimilates MODIS aerosol optical depth (AOD) retrievals into the ECMWF integrated forecast system. The MACC product provides altitude information for aerosols along with speciation information. Although the MACC speciation has not yet been validated, the MODIS cloud droplet number concentration shows a stronger sensitivity to hydrophilic aerosol types than hydrophobic aerosols, suggesting that the MACC speciation conveys useful information about the aerosol type (Gryspeerdt et al.2016). Further from sources, where ageing and other assumptions come into play, the speciation may be less accurate. In the upper troposphere, liquid aerosol is thought to be the dominant aerosol component (Kojima et al.2004), although recent studies have noted that glassy organic aerosols are abundant in the upper atmosphere and can act as INP at temperatures below −55C (Murray2008; Wilson et al.2012). This work takes the MACC total mass concentration as a measure of the liquid aerosol concentration – high (low) aerosol is more (less) than 6 µg m−3, with the caveat that this measure of aerosol may shift towards a measure of INP at the coldest temperatures considered in this work.

The response of the Ni to aerosol is likely to vary by aerosol type (Pruppacher and Klett1997). Although MACC provides a dust speciation, it is not clear whether this can be used to determine the presence of INP. Previous studies have suggested that high altitude aerosol may be responsible for glaciating clouds between 0 and −35C (Choi et al.2010; Kanitz et al.2011; Zhang et al.2012; Tan et al.2014). Based on this previous work, the glaciation of clouds between 0 and −35C is used as a proxy for the occurrence of INP.

Cloud phase is determined using the DARDAR phase detection algorithm (v1.1.4; Delanoë and Hogan2010). This algorithm uses the different response of lidar backscatter and radar reflectivity to liquid and ice hydrometeors to identify glaciated clouds. Clouds with a peak in lidar backscatter at the cloud top are treated as liquid or mixed phase and those with only a strong radar return are treated as glaciated. Experience suggests that the retrieved phase can be unreliable for clouds less than 600 m thick, so these are excluded from this part of the analysis.

Using this cloud top phase product, a “glaciation index” is developed using the phase retrievals between 0 and −35C. As approximately half of cloud tops are glaciated at −20C, glaciated clouds with a top temperature warmer than −20C are identified as “warm ice”, and liquid topped clouds colder than −20C are identified as “cold liquid”. The warm ice pixels are taken to indicate the presence of INP within 100 km (the approximate spatial scale of aerosol variability from Weigum et al.2012), whilst the cold liquid pixels are taken to indicate a relatively INP-free environment. If both (or neither) are detected within 100 km, that pixel is excluded from the analysis. The cloud phase is only used for the uppermost cloud layer when determining this INP proxy to reduce the impact of overlying ice clouds seeding ice in lower layers. In addition, regions with nearby higher cloud layers (those within a 10:1 glide slope) are also excluded from the glaciation index. To reduce the impact of random errors in the phase retrieval, two or more neighboring pixels are required for a detection of warm ice or cold liquid. This glaciation index is produced for each DARDAR vertical profile and is then used as a proxy for the occurrence of INP at temperatures colder then −35C in that profile. The use of cloud glaciation as an INP proxy for colder temperatures in the atmosphere assumes that cloud glaciation is correlated to INP between 0 and −35C and that there is a sufficient vertical correlation in INP occurrence. These assumptions are discussed in Sect. 4.2.1.

This combination of reanalysis temperature and aerosol data, along with previously determined clouds regimes and a proxy for INP are used globally for daytime data over the period from 2006 to 2013 to investigate the role of different processes on the Ni.

3 Results

3.1 Global Ni(top) distribution

The global Ni(top) distribution for crystals larger than 5 µm (Ni(top)5µm) displays several features that highlight the role of different cloud processes in controlling the Ni(top)5µm. The zonal mean Ni(top)5µm (Fig. 1a) shows a strong temperature dependence, with significant increases in the Ni(top)5µm as the temperature decreases from a mean of around 40 L−1 at −35C to almost 140 L−1 at −70C. This temperature dependence is particularly strong at temperatures colder than −40C in both the northern and the southern midlatitudes. There is also a strong Ni(top)5µm increase in the tropics, although the initial increase in Ni(top)5µm at −40C is weaker. Although the Ni produced by heterogeneous nucleation should also increase as temperatures decrease due to increasing INP concentrations (DeMott et al.2010), this strong increase in Ni(top)5µm at −40C along with a continuing Ni increase at colder temperatures is indicative of homogeneous nucleation, which is only significant at temperatures below around −35C.

At temperatures warmer than −35C, the mean Ni(top)5µm is relatively small, especially in the Northern Hemisphere where it averages less than 50 L−1. This is expected from heterogeneous nucleation, where the Ni(top)5µm is limited by available INP. The mean Ni(top)5µm is much larger in the Southern Hemisphere and the tropics, although this is skewed by the long tail of the Ni(top)5µm distribution (Fig. 2). A phase misclassification, with liquid topped cloud being mistaken for ice cloud might explain this hemispheric contrast, due to the large amounts of supercooled water in the Southern Hemisphere (Choi et al.2010). A strong lidar backscatter at the cloud top would lead to a large retrieved Ni(top)5µm (if it was misclassified as an ice cloud). As liquid topped clouds at sub-zero temperatures are more common in the Southern Hemisphere, this would result in an erroneously large mean Ni in the Southern Hemisphere.

Figure 1(a) The zonal mean DARDAR-Nice cloud top Ni (Ni(top)5µm) for crystals larger than 5 µm as a function of temperature from DARDAR-Nice data for the period from 2006 to 2013. Temperatures warmer than −35C are in greyscale. (b) The mean Ni(top)5µm at −50C. Grey indicates missing data. (c) The number of cloud top retrievals at −50C. Zonal means and maps of Ni are available in Part 1.


Figure 2The conditional Ni(top) (L−1) for 5 C temperature bins for each of the main cloud regimes (orographic, frontal, convective and synoptic) from Gryspeerdt et al. (2018a, O2, F, C and S). The top row shows the Ni(top) for particles greater than 5 µm (Ni(top)5µm) and the bottom row shows the Ni(top) for particles greater than 100 µm (Ni(top)100µm). The columns are normalized so that they sum to 100 %. The vertical line is drawn at −35C – approximately the homogeneous nucleation threshold temperature. At temperatures warmer than −35C, the grid lines show the INP numbers predicted by the DeMott et al. (2010) parameterization for 0.1, 1, 10 and 100 L−1 aerosols >0.5µm. The grid lines at temperatures below −35C are the Ni values following Koop et al. (2000) for 1, 10 and 100 cm s−1 updraft speeds for a mean pressure and an aerosol particle number of 300 cm−3, following Krämer et al. (2009). The regime names and definitions are given in Table 1 of Gryspeerdt et al. (2017). Note the log scale for Ni(top).


There are large geographical variations in Ni(top)5µm. At −50C, the Ni(top)5µm is strongly affected by the topography (Fig. 1b). High Ni(top)5µm values are retrieved in mountainous regions over land and around the edge of the Antarctic ice sheet, similar to results from orographic cirrus parameterizations in global climate models (e.g. Joos et al.2008; Barahona et al.2017) and other satellite retrievals (Mitchell et al.2016, 2018). This is consistent with a high Ni(top)5µm being generated through orographic uplift, which can generate the strong updrafts and high supersaturations required for homogeneous nucleation. While it is possible that the increased Ni(top)5µm is due to an increase in the vertical transport of INP, the lack of a similar pattern in the cloud supercooled fraction at −20C (Choi et al.2010) makes this explanation unsatisfactory. The Ni(top)5µm in the tropics is comparatively low, even in regions of significant topography such as the Ethiopian highlands. This is due to low wind speeds in the tropics reducing the in-cloud orographic updraft, similar to the GCM results of Joos et al. (2008). The high orographic Ni(top)5µm also partially explains the hemispheric asymmetry in Ni(top) in the midlatitude and polar regions, due to the high Ni(top)5µm generated by orographic clouds over the Andes and around the edge of the East Antarctic ice sheet. The Ni(top)5µm in the tropics is significantly lower than the Ni5µm at −50C (Part 1). This is partly due to the low number of cloud tops at this temperature in the tropics, meaning that there is a clear sampling bias. Additionally, the cloud top temperature plays an important role in determining the Ni5µm, giving it a much weaker temperature dependence. This temperature dependence is investigated further in the following two subsections.

3.2 Cloud regime dependence

The location map and temperature dependence of the Ni(top)5µm (Fig. 1) and the results from Part 1 hint that there may be a significant regime dependence in the Ni(top), in particular a strong role for orographic clouds and a possible role for convective clouds, given the low Ni(top) in the tropics. Separating the Ni(top) data by regime using the classification of Gryspeerdt et al. (2018a) allows this dependence to be independently investigated. Due to the strong temperature dependence and the large variability of the Ni(top), joint probability histograms, showing the probability of a Ni(top) retrieval at a given temperature are shown in Fig. 2. Following the results of Part 1, the Ni(top) is investigated for crystals bigger than 5 µm (Ni(top)5µm) and 100 µm (Ni(top)100µm). With a minimum size of 5 µm, Ni(top)5µm typically lines up with the smallest sizes measured by in situ instruments, while with a larger minimum size, Ni(top)100µm covers a size range where less shattering is expected and where the normalized size distribution performs well (Delanoë et al.2005; Sourdeval et al.2018a). As noted in the previous section, the skewed Ni(top) distribution makes a simple linear average complicated to interpret. For the remainder of this work, normalized joint histograms are used, showing the probability of finding a particular Ni, given that a certain temperature has been observed.

There are a number of broad similarities between the regimes. Each regime shows a very similar increase in Ni(top)5µm with decreasing temperature, with the decrease becoming weaker at very colder temperatures, rising to around 100 L−1 at −75C. While this is larger than the Ni values reported from many measurements of tropical tropopause layer cirrus (e.g. Jensen et al.2013a), this may be due to sampling differences between the satellite and in situ measurements, with some of the thinnest clouds being missed by the DARDAR retrieval. It is also possible that uncertainties in the shape of the particle size distribution (PSD) can lead to overestimations of Ni(top)5µm by as much as a factor of 2 (Sourdeval et al.2018a). The temperature dependence is similar to that observed during the SPARTICUS and MACPEX campaigns (Jensen et al.2013b; Muhlbauer et al.2014), although the temperature dependence is stronger than that observed in Krämer et al. (2009) where the Ni was sampled in cloud, rather than at the cloud top. However, if the satellite data are sampled in a manner similar to previous work, it reproduces the in situ results (Sourdeval et al.2018a), giving confidence in the magnitude and temperature dependence of the results presented here. There is evidence of possible retrieval errors, as both the orographic and convective regimes have a small number of retrievals of over 30 L−1, with some as high as 50 L−1, around −15C. This suggests that the possible phase misclassification and the high Ni(top) values observed in the zonal mean are more common in certain regimes.

All of the regimes also show a peak in the highest Ni(top)5µm percentiles at temperatures just colder than −35C. The strength varies by regime, with the orographic regime showing a stronger peak and only a weak peak being observed in the frontal and convective regimes. The peak is barely present in the synoptic regime. An increase in the largest Ni(top)5µm values at this temperature is consistent with homogeneous nucleation, either through an increase in the freezing of liquid droplets or by increased homogeneous nucleation through the freezing of unactivated aqueous haze particles. At these colder temperatures, the Ni(top) is roughly parallel with the contours of Ni expected through homogeneous nucleation at a constant updraft (Krämer et al.2009), especially in cases where the peak in Ni(top)5µm at −40C is small (Fig. 2). Warmer than −35C, the Ni(top)5µm is broadly consistent with the number of INP predicted by the DeMott et al. (2010) parameterization, but the Ni(top)5µm becomes increasingly large compared to the number of INP as the temperature increases. As the DARDAR-Nice retrieval has not been evaluated at this temperature, it is unclear if this is a real effect, or if it is due to the possible phase classification issue mentioned previously.

Figure 3(a) The difference in the Ni(top)5µm as a function of temperature between the Oro 2 and Oro 1 regimes (highest and second highest sextiles of estimated in-cloud updraft, respectively). Red above blue at a given temperature indicates an increased Ni(top)5µm in the Oro 2 regimes compared to Oro 1. (b) The difference between the frontal and synoptic regimes. Note the different color scale from (a).


The variation in the strength of this peak is clearly seen when comparing the Oro 2 and Oro 1 regimes (the highest and second highest sextiles of the estimated in-cloud updraft) in Fig. 3a. While there is little difference between the regimes at warmer temperatures, below −35C there is a strong increase in the Ni(top)5µm in the Oro 2 regime. This increase peaks at about −50C, reducing and then almost disappearing at the coldest temperatures studied. The high Ni(top)5µm retrieved in these clouds and the strong dependence on the in-cloud updraft explain the geographical pattern shown in Fig. 1b, where high Ni(top)5µm are observed in mountainous regions. A high Ni(top)5µm in these regimes is supported by results from previous in situ studies, where large Ni values were recorded in orographic clouds (Jensen et al.1998; Field et al.2001; Baker and Lawson2006).

It is possible that this increased Ni(top)5µm is the result of increased aerosol concentrations carried to lower temperatures in the stronger updrafts of the Oro 2 regime. However, the lack of a difference in Ni(top)5µm between the regimes at temperatures warmer than −35C indicates that an increase in INP is not driving this change in Ni(top)5µm, which in turn makes it less likely that this change is due to a change in liquid aerosol. As the Oro regimes are defined by the estimated in-cloud updraft (Gryspeerdt et al.2018a), the difference between the regimes shown in Fig. 3 is likely due to a change in the updraft environment impacting freezing processes.

A change in the updraft environment could modify the Ni(top)5µm by changing the likelihood of homogeneous nucleation, either through allowing more liquid droplets to reach temperatures where they can freeze homogeneously or by increasing the nucleation of haze droplets (Koop et al.2000). These processes cannot be easily distinguished in the current study, although the lack of a significant occurrence of liquid-topped cloud in orographic regions (Tan et al.2014) suggests that an increased cloud droplet number concentration is not the leading contributor to the increase in Ni(top)5µm. This strong response to updraft changes would support previous studies that highlighted the updraft limited nature of many cirrus clouds (Kay and Wood2008; Barahona and Nenes2008). A larger difference exists between the frontal and synoptic regimes (Fig. 3b), indicating that the magnitude of this updraft effect could be stronger than is shown here. However, the difference between the frontal and synoptic regimes cannot be easily attributed to updraft variations.

The temperature dependence of crystals larger than 100 µm (Ni(top)100µm) at the cloud top displays a different pattern (Fig. 2, bottom row). While Ni(top)100µm and temperature are negatively correlated at warmer temperatures, the Ni(top)100µm reaches a peak at around −50C and there is a decrease in the Ni(top)100µm as temperatures reduce further, with the strongest decrease observed in the orographic regime. This is consistent with a shift towards smaller ice crystals at the cloud top with colder temperatures. The synoptic regime shows a weaker decrease in Ni(top)100µm, indicating a slightly larger role for larger ice crystals in this regime. This shift towards smaller crystals at the cloud top is expected due to slower depositional growth and aggregation of ice crystals at colder temperatures resulting in crystals precipitating from the cloud top region before they grow larger than 100 µm.

3.3 The Ni within clouds

The behavior of the Ni within clouds as a function of temperature displays some significant contrasts to the Ni(top)5µm (Fig. 4). While all of the regimes show an increase in the Ni5µm with decreasing temperature, this increase is much weaker than the increase in the Ni(top)5µm. Similarly, although the peak that is visible in the Ni(top)5µm at about −40C is still visible in Ni5µm in the orographic and frontal regimes, it is much weaker than the peak observed at the cloud top. The synoptic regime has the strongest temperature dependence of all of the regimes. One explanation is the lower average cloud depth, such that the Ni5µm retrieved is often closer to the cloud top than in the other regimes. In all of the regimes, the Ni5µm is much larger in-cloud than at the cloud top for temperatures warmer than −30C, with values up to 100 L−1 being commonly observed. The smaller Ni5µm values that are more typically observed in the synoptic regime than the other regimes suggest that seasonal variations of the cloud classes (Gryspeerdt et al.2018a) are likely responsible for the seasonal variations in Ni5µm observed in Part 1.

Figure 4As in Fig. 2, but using the Ni from throughout the cloud, rather than just the cloud top. The temperature scale is the temperature of the Ni retrieval, rather than that of the cloud top.


Similar to the Ni5µm, the temperature dependence of the Ni100µm is very different internally within clouds compared to at cloud tops. The temperature dependence is much weaker, with almost no temperature dependence at temperatures warmer than −50C. There is a decrease in the Ni100µm at the lowest temperatures, similar to the decrease in the Ni(top)100µm seen in Fig. 2 and is explained by the retrievals at these temperatures being closer to the cloud top than at warmer temperatures. The synoptic regime has the lowest Ni100µm at warmer temperatures, which may again be due to the lower geometrical thickness of clouds in this regime, such that the Ni100µm is typically located closer to the cloud top, resulting in a lower Ni(top)100µm for any given temperature inside a cloud.

The larger Ni100µm values at warmer temperatures mean that larger crystals comprise a larger proportion of the Ni5µm, with a reduced contribution of small crystals to the Ni5µm. A weaker temperature dependence of the Ni5µm, especially at temperatures colder than −35C, is in better agreement with the results from Krämer et al. (2009), although a temperature dependence remains. It is possible that the weak temperature dependence in previous results could be due to a lack of measurements near the cloud top, where the temperature dependence is strongest. This may also explain the apparent mismatch between the INP and Ni concentration in aircraft data (e.g. Kanji et al.2017), as the retrieved Ni(top) values are a much closer match to the INP concentrations predicted by the DeMott et al. (2010) parameterization than the internal Ni at temperatures warmer than −35C (Fig. 2). Further sampling differences between the satellite and in situ studies due to the detection limits of satellite instruments and the structuring of flight campaigns may explain the remaining differences between Ni determined using different methods.

3.4 Vertical structure of Ni

The Ni5µm, Ni100µm and ice water content (IWC) all change significantly as a function of depth through the cloud (Fig. 5). For clouds with a top temperature between −40 and −50C (Fig. 2), the Ni5µm continues to increase until about 500 m from the cloud top, at which point it starts to decrease again (Fig. 5, top row panels). The Ni5µm distribution width stays approximately constant from about 1 km into the cloud until around 2–3 km from the cloud top, when it reaches a temperature of around −30C where liquid water can form more easily. At this point the Ni5µm distribution broadens significantly. Similar to the Ni5µm, the Ni100µm also grows quickly when moving down through the cloud, moving to a slower growth regime after the first 500 m from the cloud top. This shift in the Ni100µm growth regime is roughly coincident with the location of the Ni5µm peak. All of the regimes also show an increase in the IWC (Fig. 5, bottom row) with increasing depth in the cloud, which displays a sharp increase over 2.5 km from the cloud top. This sharp increase is consistent with a possible increase in ice through liquid water processes in warmer parts of the cloud.

Figure 5Retrieved properties as a function of the distance from the cloud top. This is for clouds with tops between −40 and −50C. Note the nonlinear scale on the horizontal axis.


There are some differences between the regimes. The synoptic regime has a much weaker peak in Ni5µm below the cloud top and a consequently lower Ni5µm throughout the depth of the regime. Despite having similar values at the cloud top to the other regimes, the Ni100µm and the IWC in the synoptic regime remain lower than the other regimes through the cloud, possibly due to lower in-cloud updrafts. At about 2.5 km from the cloud top, both Ni100µm and IWC increase until they are roughly comparable to the other regimes, suggesting that the signal from liquid water swamps any signal based on ice nucleation.

The peak is also temperature dependent, almost disappearing in clouds with colder tops (see supplementary information) and varying in size and location between the regimes. When the peak occurs in the synoptic regime, it is within 300 m of the cloud top in 67 % of cases, compared to only 48 % of cases for the frontal regime. These distances are comparable to the thickness of nucleation regions of between 20 and 500 m noted in Jensen et al. (2016). The enhancement of the Ni5µm within this peak in the synoptic regime is also smaller, with an average peak of 130 L−1, compared to 270 L−1 in the frontal regime and 325 L−1 for the orographic regime. The increased strength of this peak in regimes expected to have a stronger updraft along with its location close to the cloud top may indicate homogeneous nucleation. Model studies of cirrus clouds suggest that homogeneous nucleation can produce peaks in Ni cloud to the cloud top (Spichtinger and Gierens2009a, b), with an increased Ni at higher updraft velocities. The disappearance of the peak at colder temperatures gives it a similar temperature dependence to the peak in the Ni(top)5µm (Fig. 2) providing further supporting evidence of the impact of homogeneous nucleation on Ni in this temperature range.

Figure 6The difference in the conditional histograms between cases with high MACC total aerosol mass concentration (>6µg m−3) and low total mass concentration (<6µg m−3) for the four main regimes. The grid lines are the same as Fig. 2. The upper set of plots show the difference in Ni(top)5µm and the lower in Ni(top)100µm. The changes sum to zero vertically; red over blue indicates an increase in the Ni(top)5µm/Ni(top)100µm for a given temperature and regime.


It is possible that the varying sensitivities of the CloudSat radar and the CALIOP lidar to crystal size and the attenuation of the CALIOP lidar in the upper levels of the cloud could be generating this vertical structure. The lower vertical resolution and sensitivity to small crystals of the radar could result in it missing the cloud top, which would generate a peak in the Ni5µm at the level where the retrieval starts to include radar information. However, the results in Part 1 show no evidence of a bias in the Ni retrieval as a function of the instruments contributing to the retrieval (Sourdeval et al.2018a). This is primarily due to the sensitivity of the instruments to different ice crystal size distributions. Although the lidar-only retrievals have a higher expected error, they usually only occur in cases where there is a monomodal size distribution dominated by small crystals that can be accurately constrained by the lidar alone. Additionally, the disappearance of the peak at colder temperatures indicates that it is a physical property of the clouds, rather than a property of the retrieval, as the instrument sensitivities would not be expected to strongly vary with temperature.

4 The relationship to aerosol

4.1 The relationship to liquid aerosol

Figure 6 shows how the Ni(top)5µm distribution changes as a function of temperature and MACC reanalysis aerosol (used to indicate high concentrations of liquid aerosol). In most of the regimes, there is a positive relationship between MACC aerosol and Ni(top)5µm at temperatures below −35C (shown by red above blue in Fig. 6). In the synoptic regime, this positive aerosol–Ni(top)5µm relationship only exists for temperatures warmer than −60C – at temperatures colder than this, the relationship becomes weak and noisy. In the other regimes, the positive relationship is maintained to very cold temperatures. At temperatures warmer than −35C, the relationship becomes a lot weaker, with almost no aerosol–Ni(top)5µm relationship existing in the orographic and convective regimes. In the frontal regime, there is a slight negative relationship, with a stronger negative relationship in the synoptic regime. It is possible that this negative relationship is related to a misclassification of ice and liquid at these warmer temperatures being a function of the MACC aerosol, particularly in regions where INP rich aerosol constitute a majority of the aerosol population.

The aerosol–Ni(top)100µm relationship shows a weaker pattern than the aerosol–Ni(top)5µm relationship, with the smaller enhancement of the Ni(top)100µm at colder temperatures in most regimes indicating a shift to smaller crystal sizes. The change in the synoptic regimes is the strongest, likely related to the strong relationship for the Ni(top)5µm. A negative relationship between the aerosol environment and the crystal size has been noted in previous work (Jiang et al.2011; Zhao et al.2018) and often corresponds to an increase in Ni(top), although positive relationships have been observed over the Indian Ocean (Chylek et al.2006).

It is difficult to demonstrate causality with observed aerosol–cloud relationships, to the extent that it is not clear that this relationship is a change in Ni(top) due to a change in aerosol. However, this strong relationship between MACC aerosol and Ni(top) is consistent with an increased ice crystal nucleation through homogeneous nucleation, which can be sensitive to the concentration of liquid aerosol (e.g. Kärcher and Lohmann2002). In situations where the Ni(top) is primarily determined by the freezing of liquid droplets, an increase in cloud droplets in high aerosol regions could also lead to an increased Ni(top), although the number of droplets frozen is relatively insensitive to the total number of liquid droplets (Kärcher and Seifert2016). As with the impact of in-cloud updraft on Ni(top), further investigation is required to determine if one of these mechanisms is dominant. As liquid water has been found in clouds at temperatures as cold as −40C, increased droplet freezing cannot be ruled out, even though many clouds are frozen before reaching this temperature (Choi et al.2010). At colder temperatures, it seems likely that homogeneous nucleation plays a role, as liquid droplets cannot form at these temperatures. In this case, the stronger updrafts in the frontal and convective regimes are important for generating the high supersaturations in which homogeneous nucleation can occur. Changing aerosol types may also play a role at temperatures colder than −60C, where the increasing impact of glassy aerosols may lead the aerosol to nucleate ice heterogeneously. A combination of the weak expected updrafts and the increasing ability of glassy aerosol to act as an INP at low temperatures may explain the weak aerosol–Ni(top) relationship in the synoptic regime below −60C. While there is a clear relationship in Fig. 6, the change in the mean Ni(top) is small, even for this large aerosol perturbation. At −50C, the mean Ni(top)5µm increases from around 140 to 175 L−1, an increase of 25 %. Much of this change is driven by changes in the high updraft tail of the distribution, and producing a 25 % change in Ni(top) at −50C would require an updraft in excess of 1 m s−1 (Kärcher and Lohmann2002). While plausible for the convective and orographic regimes (Gryspeerdt et al.2017), the large updrafts required to generate such a sensitivity may indicate that this relationship is affected by an updraft mediated covariation.

4.2 The relationship to an INP proxy

The sparse nature of INP measurements (e.g. Mamouri and Ansmann2016) and the high sensitivity of the Ni to low INP concentrations means that it is difficult to use retrieved aerosol properties to investigate the effect of INP on the Ni(top). To avoid this issue, the glaciated fraction of clouds lower in the atmosphere (−20C) is used as a proxy for the presence of INP at other levels in the atmosphere.

Figure 7(a) The DARDAR supercooled fraction at −20C, defined as the fraction of the DARDAR cloud top phase retrievals between −17.5 and −22.5C from 2006 and 2013 that are classed as liquid. (b) The conditional probability of observing a daily mean supercooled fraction, given a specified MACC dust mixing ratio for the regions specified in (a). The black line shows the mean supercooled fraction for each aerosol bin.


4.2.1 Cloud glaciation and INP

The addition of the CloudSat data in the DARDAR phase retrieval allows smaller quantities of ice to be detected than in the lidar-only studies, but it produces a very similar pattern of cloud glaciation (Fig. 7) to the previous CALIOP studies (Choi et al.2010; Tan et al.2014). The supercooled fraction is calculated as the number of DARDAR liquid phase retrievals divided by the total number of liquid and ice retrievals between −17.5 and −22.5C from 2006 to 2013. Only the phase of the top cloud layer is considered and only where this layer is more than 600 m thick.

There is a strong hemispheric contrast with a higher glaciated fraction over the Northern Hemisphere and a high supercooled fraction over the Southern Ocean and Antarctica, as observed in previous aircraft (Huang et al.2012) and satellite studies (Choi et al.2010; Tan et al.2014). High glaciated fractions are observed over desert locations in central Asia and Iran, stretching across the North Pacific to the Americas. This is consistent with previous studies suggesting that dust is a good INP. Previous studies have found Asian dust over California, suggesting that transport across the Pacific is not unexpected (Creamean et al.2013). There is also a significant proportion of glaciated cloud downwind of the Andes, which appears to originate near the Altiplano and Patagonia. These are sources of high altitude dust (Ginoux et al.2012) and would support the hypothesis that high altitude dust is able to glaciate clouds. While glaciated cloud in this region has previously been noted (Choi et al.2010), the lower resolution of the previous study made it difficult to determine the source of possible INP. The longer dataset and increased spatial resolution of Fig. 7a make the source in the upper Andes much clearer. Although southern Africa and Australia are also sources of dust (Ginoux et al.2012), this dust is emitted at lower altitudes, which would explain the lower glaciated fractions downwind of these regions.

The origin of the glaciated region over the north Atlantic is less clear, as there are not many local sources of high level dust in the region. It is possible that the dust here has been transported across the Sahara and lofted by cyclone systems crossing the Atlantic. Furthermore, it is possible that the black carbon or ash (Grawe et al.2016) from North America may act as an INP. This might explain the lower supercooled fraction over Siberia, where black carbon from fires typically occurs without the other aerosols that are found in industrial pollution, allowing it to act as an INP (Rosenfeld et al.2011) despite the low amounts of high level dust in this region.

The role of the ice nucleation impact of dust for driving the cloud glaciated fraction is supported by comparing the cloud glaciated fraction to reanalysis aerosol fields (Fig. 7b). Strong negative correlations between the occurrence of supercooled liquid cloud at −20C and the mass concentration of reanalysis dust (Fig. 7b) are observed in some regions, with glaciated cloud dominating at high mass concentrations of reanalysis dust. However, this correlation varies by region. A stronger relationship is found in regions close to dust sources, such as over the N. Pacific (B) and central Asia (D); the relationship is much weaker in the N. Atlantic (A) and the Southern Ocean (C) where the dust is further from source (Fig. 7b).

The stronger dust–glaciation relationship close to the dust source, where the MACC aerosol speciation is best, suggests that the supercooled fraction of clouds at −20C is strongly related to the occurrence of INP. The weaker relationship further from source suggests that although the MACC speciation has been shown to provide useful information on aerosol type, this speciation is less reliable further from source. This is supported by results in liquid clouds, where the dust optical depth–cloud droplet number concentration relationship becomes stronger further from dust sources (Gryspeerdt et al.2016).

Due to the reduced speciation skill from MACC far from dust sources, the occurrence of glaciated cloud at −20C is used as a proxy for the occurrence of INP instead of the reanalysis aerosol. This relies on the following two assumptions: (1) cloud glaciation at −20C is related to INP at −20C; and (2) INP at −20C are correlated to INP at other temperatures.

Based on the relationship to MACC dust aerosol, the first assumption holds in many cases. Although the second assumption is tenuous, previous studies have found similar relationships between cirrus cloud properties and both column and layer aerosol optical depth (AOD) (Zhao et al.2018), similar to model results showing a significant correlation between high altitude cloud condensation nuclei concentration and column AOD (Stier2016). Significant vertical aerosol autocorrelation has also been observed in global climate models (Weigum2014). Additionally, there is very unlikely to be a negative correlation between the INP at the two temperature levels, with the worst case being no correlation. As such, the relationship between the Ni and the INP proxy is unable to give a quantitative result for the impact of INP on the Ni(top)5µm, but it is able to provide a qualitative indication of the sign of the INP impact.

4.2.2 The INP relationship to Ni(top)

The relationship of the Ni(top)5µm to the proxy for INP is shown in Fig. 8. There are a number of features that are similar between the regimes, in particular the strong negative relationship between INP occurrence and Ni(top)5µm at temperatures warmer than −35C. As with the large mean Ni(top)5µm values shown in Fig. 1a, this may be due to liquid clouds being misclassified as ice, resulting in large Ni(top)5µm values being retrieved. The requirement for “warm-ice” means that supercooled liquid occurs less frequently in the high INP cases, and as such it is less likely to be misclassified as ice. The lower frequency of this misclassification then reduces the Ni(top)5µm and Ni(top)100µm in cases of high INP. The weaker Ni(top)5µm response in the synoptic and orographic regimes supports this, as the misclassification in these regimes is weaker (Fig. 2). The warmer temperatures are shaded out in Fig. 8 due to the impact of this potential misclassification.

Figure 8As in Fig. 6, but showing the difference in the conditional histograms as a function of the INP proxy. Red indicates an increase in the occurrence of a particular bin at a higher inferred INP and blue indicates a decrease, meaning that red above blue indicates an increase in Ni(top) with increased INP for a given temperature. The shaded regions are likely affected by a phase misclassification at warmer temperatures.


At colder temperatures, the INP–Ni(top)5µm relationship starts to vary between the regimes. All of the regimes show a decrease in the Ni(top)5µm between around −35 and −50C, the temperatures where the peak in Ni(top)5µm is observed connected with in-cloud updraft (Fig. 3). This decrease is strongest in the orographic regime and weakest in the synoptic regime, similar to the Ni(top)5µm peak observed in the different regimes (Fig. 2). At temperatures colder than −50C, the relationship becomes different again. In all of the regimes, there is an increase in Ni(top)100µm with increasing INP. This is consistent with an increasing number of INP shifting the size distribution towards a smaller number of larger ice crystals. In the orographic and synoptic regimes, this increase also appears in the Ni(top)5µm, generating a positive relationship between the INP proxy and the occurrence of small ice crystals.

As with the previous section, the impact of meteorological covariations cannot be ruled out when interpreting these plots. However, they are consistent with a reduction in Ni(top)5µm due to a suppression of homogeneous nucleation by INP at around −50C. This relationship has previously been found in satellite relationships between the aerosol environment and the ice crystal size, with an increase in the crystal radius in situations where heterogeneous nucleation controls the Ni (Chylek et al.2006; Zhao et al.2018). This would fit with the results in previous sections, suggesting that the Ni(top) at this temperature range just slightly colder than −35C is influenced by homogeneous nucleation. This effect would only be expected in a narrow range of updrafts (Kärcher and Lohmann2002), so further work is necessary to understand the cause of this relationship.

The increase in large crystals at the coldest temperatures (below −60C) is consistent with an INP effect on Ni(top) if heterogeneous nucleation was dominant at these temperatures. This would fit with the results from Fig. 6, where at the coldest temperatures, there was a relatively small response of the Ni(top)5µm to MACC total (liquid) aerosol, suggesting that homogeneous nucleation was not controlling the Ni(top)5µm in synoptic cirrus. At these coldest temperatures, dust can act as an INP at very low supersaturations (as low as 105 %; Möhler et al.2006) and organic aerosol can occur in a glassy state allowing it to act as an INP. This may explain relationships consistent with heterogeneous nucleation and a classical Twomey effect at these temperatures. It is important to note that this proxy for INP relies upon the correlation between cloud glaciation at −20 and INP at −50C, but the absence of this correlation would produce no relationship in Fig. 8, giving some confidence to the qualitative nature of these results.

If the peak in Ni(top)5µm at temperatures colder than −35C is primarily due to droplet freezing, an increase glaciated fraction at warmer temperatures could also result in this reduction of Ni(top)5µm with increasing INP. As the number of INP and Ni warmer than −35C is much lower than the cloud droplet number concentration, the increase in cloud glaciation could result in a reduction in the number of cloud droplets available to form ice crystals at −35C. This would result in a negative relationship between cloud glaciation at −20C and the Ni(top)5µm at colder temperatures as observed in Fig. 8. As with the relationship of Ni(top) to updraft (Fig. 3) and aerosol (Fig. 6), the difference between an aerosol impact on homogeneous nucleation, a change in droplet freezing or an updraft-mediated covariation (no causal effect of aerosol) cannot be distinguished by this analysis.

Figure 9As in Fig. 5, but showing the difference in the retrieved properties depending on the cloud top properties. Red over blue indicates that clouds with above median properties at the cloud top (Ni(top)5µm, IWC(top)) have higher values of the retrieved properties at a specified depth from the cloud top. Note the nonlinear scale on the horizontal axis.


5 Vertical information propagation

The changes in Ni(top)5µm observed in the previous section have impacts throughout the depth of the cloud. Figure 9 shows how Ni and IWC information propagates vertically within a cloud. The cloud profiles are split into two categories, based on whether they have above or below median values of the cloud top properties (Ni(top), IWC(top)). The difference in the vertical structure of the clouds (in a similar manner to Fig. 5) is shown, with red over blue indicating an increase in the retrieved quantity at a given distance from the cloud top for profiles that were above median in that property at the cloud top.

The top row of Fig. 9 shows that Ni5µm information propagates a significant distance through the cloud. Clouds with an increased Ni(top)5µm maintain a higher Ni5µm at distances at least 3km from the cloud top in all regimes. However, as shown in the second row, vertical information about the Ni100µm does not propagate nearly as far through the cloud. The vertical propagation is the highest in the synoptic regime. The vertical propagation of IWC information is very similar to the Ni100µm, with the relationship to the cloud top IWC being significantly reduced more than 500 m from the cloud top.

The large vertical propagation of the Ni5µm indicates that the changes in Ni5µm at the cloud top found in the previous section can have considerable impacts at lower levels in the cloud. However, the lower vertical propagation of the information about the larger crystals (Ni100µm, IWC) would support the suggestion that the growth of the ice crystals after nucleation is primarily controlled by meteorological factors that do not play a large role in the nucleation processes that control Ni100µm. Note that the temperature of the cloud top and the distance from the cloud top can still play a large role in determining the Ni (Figs. 4 and 5).

6 Discussion

These results show that the Ni(top)5µm is strongly affected by several factors including temperature (Fig. 1), cloud type (Fig. 2) and updraft (Fig. 3), and that changes in the Ni(top)5µm can be maintained at large distances from the cloud top (Fig. 9). The dependence of the Ni(top)5µm on the in-cloud updraft and the relationship to reanalysis liquid aerosol (Fig. 6) at temperatures between −35 and −60C is consistent with the impact of homogeneous nucleation processes on the Ni(top)5µm. This is supported by the relationship of the Ni(top)5µm to the INP proxy (Fig. 8), where a reduction in Ni(top)5µm with increasing INP could be indicative of an INP suppression of homogeneous nucleation (Kärcher and Lohmann2003). The relationship with INP is also consistent with heterogeneous nucleation having a strong role to play in determining the Ni(top)5µm in synoptic cirrus clouds at temperatures colder than 60 C.

Uncertainties in the retrieval have been covered in Part 1 of this work (Sourdeval et al.2018a). However, there are a few points to note with regards to the relationship of the Ni to other cloud and meteorological properties. Although there is significant uncertainty in the Ni retrieval, many of these uncertainties are random errors and not systematic functions of the meteorological properties investigated here. Even ice crystal shape, which can be a major issue in ice cloud retrievals, is a function of temperature (to first order) and so does not impact the majority of the results which are presented in this work stratified by temperature. The geographical variations in Fig. 1b show a similar pattern to those from Mitchell et al. (2016, 2018), with high Ni(top) observed in mountainous regions and a reduced Ni(top) in the tropics. The similarity of the results from these two different retrieval products, each with a different physical basis supports the conclusions drawn from these datasets regarding the global distribution of Ni. There is also little evidence to suggest that there are large biases caused by the retrieval only being able to use one instrument (radar or lidar). Cases where only the lidar detects a cloud are often characterized by monomodal ice distributions, which are well represented by the Delanoë et al. (2005) parameterization. As such, these cases are retrieved with similar accuracy to the full radar–lidar retrieval (Sourdeval et al.2018a).

The cloud phase classification is of critical importance to the warmer clouds included in this study and there is evidence of the misclassification of a small number of cases at temperatures warmer than −35C (Fig. 8). This can make it difficult to interpret results at these temperatures, so they are not a focus of this work. The change in phase of these clouds as a function of aerosol is likely to dominate the radiative response of clouds to aerosols at these temperatures.

There are a number of limitations of this study that could be addressed in future work. The lack of information about the location of INP is a serious issue when investigating the impact of aerosol on Ni. While the INP proxy in this work is able to provide a qualitative estimate of the role of INP in determining the Ni, for a quantitative estimate a better proxy or measure of the global INP concentration is required.

Additionally, the impact of meteorological covariations makes it difficult to assign causality to the aerosol–Ni(top) relationships observed in Fig. 6. The lack of a complete picture of the atmosphere makes it difficult to directly control for meteorological variability. The causal link between aerosol and Ni(top) is thought to be strong (e.g. Kärcher and Lohmann2002), but the lack of observations of in-cloud updrafts also limits how accurately the impact of aerosol on the Ni can be determined. Although the cloud regimes used have some ability to constrain the cloud-scale updraft (Gryspeerdt et al.2018a), the updraft is a critical component in determining the Ni through its influence on the supersaturation. The in-cloud updraft is assumed to be largely independent of the aerosol properties in this work, but it is possible that the reanalysis aerosol is related to the in-cloud updraft, such that more aerosol is vertically transported in conditions with high in-cloud updrafts. In this case, a positive correlation between the Ni and MACC reanalysis aerosol could be generated. However, as MACC does not explicitly simulate in-cloud updrafts, the impact of this confounding issue is likely to be small.

It is also possible that using cloud glaciation as a proxy allows other meteorological covariations, which could generate apparent relationships between the INP and Ni. However, the in-cloud updraft is of a second-order importance in determining the cloud top phase compared to the INP concentration (Bühl et al.2013). The inclusion of a glide-slope test when determining the INP proxy means that it is also unlikely that clouds are being glaciated by undetected ice falling from higher cloud layers. The separation into cloud regimes also limits the impact of these kind of meteorological covariations, which might be expected between different regimes, but would be weaker within them.

The behavior of the Ni retrieval in this work follows the expected behavior of the Ni determined in several previous studies based on satellite remote sensing, in situ, theoretical and modeling results. This provides further evidence that the DARDAR-Nice Ni retrieval described in Sourdeval et al. (2018a) is able to retrieve the Ni in a variety of situations.

7 Conclusions

Few global studies exist of the controls on the ice crystal number concentration (Ni), especially regarding the role of aerosols. In this study, the DARDAR-Nice Ni retrieval from Part 1 (Sourdeval et al.2018a) is used to investigate possible controls on the Ni at a global scale for the period from 2006 to 2013. A special emphasis is placed on the Ni at the cloud top (Ni(top)), due to the close proximity to ice crystal nucleation locations within many high clouds (Spichtinger and Gierens2009a; Diao et al.2017).

Strong relationships between the Ni(top) and updraft, cloud type and particularly temperature are observed (Figs. 1 and 2), with a higher Ni(top) for crystals larger then 5 µm (Ni(top)5µm) being found at colder temperatures in all regimes, which is consistent with an increased nucleation rate at lower temperatures. Fewer crystals larger than 100 µm, (Ni(top)100µm) are found at the coldest temperatures, possibly due to the reduced depositional growth rate meaning that they sediment from the cloud top region before they can grow to a sufficient size.

Many of the regimes show an increase in the Ni100µm and a decrease in the Ni5µm with increasing distance from the cloud top (Fig. 5) due to the size sorting impact of sedimentation. The rate of change of the Ni moving away from the cloud top depends on the regime, with much slower changes in the synoptic regime indicating a role of meteorological factors in determining ice crystal growth rates. This is supported by the weaker temperature dependence of the Ni within clouds compared to the Ni(top) (Fig. 4), which may also explain the apparent weak dependence of Ni on temperature (Krämer et al.2009) and INP (Kanji et al.2017) found in previous studies. Given the large difference between the Ni(top) and the Ni deeper in the cloud, this may suggest that the cloud top would make a good target for future in situ campaigns examining the controls on ice nucleation.

There are indications of homogeneous nucleation or possibly the freezing of liquid droplets determining the Ni(top). At temperatures just colder than −35C, there is a peak in the upper quantiles of the Ni(top)5µm (Fig. 2). This peak is related to the updraft strength in the cloud, with the reliably high updrafts in the orographic regime giving it the strongest peak (Fig. 3). This is further supported by the relationship between the Ni(top)5µm and the MACC reanalysis aerosol (Fig. 6), with an increased Ni(top)5µm being observed in high aerosol environments. This indicates a possible dependence on the liquid aerosol concentration, particularly for smaller crystals, although this analysis cannot make a conclusive statement about the causality of this relationship. An investigation into the covariances between the MACC reanalysis aerosol, the DARDAR-Nice Ni and meteorological factors is an important target for future work.

As previous work has suggested that INP occurrence is related to cloud glaciation, the glaciated fraction at −20C is used as a qualitative proxy for INP occurrence (Fig. 7). At temperatures between −50 and −35C, there is a reduction in Ni(top)5µm with increasing INP (Fig. 8), which may indicate an INP suppression of the homogeneous nucleation (Kärcher and Lohmann2003), providing further supporting evidence for the role of homogeneous nucleation in determining the Ni(top). At colder temperatures, some regimes show an increasing Ni(top), and the Ni(top)100µm in particular, which may be evidence of heterogeneous nucleation controlling the Ni(top) and shifting the size distribution towards larger crystals. However, as with the relationship to liquid aerosol, meteorological covariations could be generating these relationships. Further studies are required to separate the role of these different mechanisms in controlling the Ni(top) and to isolate the role of aerosols in these relationships.

While changes to the Ni(top) are important for radiative considerations, changes in the Ni(top) can have implications for the cloud many kilometers below the cloud top (Fig. 9). This far reaching impact into the life cycle of ice and mixed-phase clouds demonstrates the importance of developing strong observational constraints on the controlling factors of the Ni. The results presented in this work provide a global context for existing theory and in situ measurement based hypotheses about cloud properties, highlighting areas for future research to further constrain ice and mixed-phase cloud processes.

Data availability

The DARDAR data product was retrieved from the ICARE data center (, Sourdeval et al.2018b). The IC-CIR cloud classification is available at (, Gryspeerdt et al.2018b).


The supplement related to this article is available online at:

Author contributions

EG, OS and JQ conceived the study, EG, OS and PK performed the analysis, and EG and OS wrote the paper. All of the authors assisted with the interpretation of the results and commented on the paper.

Competing interests

The authors declare that they have no conflict of interest.


This work was supported by funding from the European Research Council under the European Union's Seventh Framework Programme (FP7/2007-2013)/ERC grant agreement no. FP7-306284 (“QUAERERE”); the Bundesministerium für Bildung und Forschung, grant numbers 01LK1210D, 01LK1503A and 01LK1505E; and the Deutsche Forschungsgemeinschaft, grant number QU 311/14-1. Edward Gryspeerdt is supported by an Imperial College Junior Research Fellowship.

Edited by: Matthias Tesche
Reviewed by: two anonymous referees


Baker, B. A. and Lawson, R. P.: In Situ Observations of the Microphysical Properties of Wave, Cirrus, and Anvil Clouds. Part I: Wave Clouds, J. Atmos. Sci., 63, 3160–3185,, 2006. a

Barahona, D. and Nenes, A.: Parameterization of cirrus cloud formation in large-scale models: Homogeneous nucleation, J. Geophys. Res., 113, D11211,, 2008. a, b

Barahona, D., Molod, A., and Kalesse, H.: Direct estimation of the global distribution of vertical velocity within cirrus clouds, Sci. Rep., 7, 6840,, 2017. a

Benedetti, A., Morcrette, J.-J., Boucher, O., Dethof, A., Engelen, R. J., Fisher, M., Flentje, H., Huneeus, N., Jones, L., Kaiser, J. W., Kinne, S., Mangold, A., Razinger, M., Simmons, A. J., and Suttie, M.: Aerosol analysis and forecast in the European Centre for Medium-Range Weather Forecasts Integrated Forecast System: 2. Data assimilation, J. Geophys. Res., 114, D13205,, 2009. a

Boucher, O., Randall, D., Artaxo, P., Bretherton, C., Feingold, G., Forster, P., Kerminen, V.-M., Kondo, Y., Liao, H., Lohmann, U., Rasch, P., Satheesh, S., Sherwood, S., Stevens, B., and Zhang, X.: Clouds and Aerosols, in: book section 7, Cambridge University Press, Cambridge, UK and New York, NY, USA, 571–658,, 2013. a

Bühl, J., Ansmann, A., Seifert, P., Baars, H., and Engelmann, R.: Toward a quantitative characterization of heterogeneous ice formation with lidar/radar: Comparison of CALIPSO/CloudSat with ground-based observations, Geophys. Res. Lett., 40, 4404–4408,, 2013. a

Choi, Y.-S., Lindzen, R. S., Ho, C.-H., and Kim, J.: Space observations of cold-cloud phase change, P. Natl. Acad. Sci. USA, 107, 11211–11216,, 2010. a, b, c, d, e, f, g

Chylek, P., Dubey, M. K., Lohmann, U., Ramanathan, V., Kaufman, Y. J., Lesins, G., Hudson, J., Altmann, G., and Olsen, S.: Aerosol indirect effect over the Indian Ocean, Geophys. Res. Lett., 33, L06806,, 2006. a, b, c

Creamean, J. M., Suski, K. J., Rosenfeld, D., Cazorla, A., DeMott, P. J., Sullivan, R. C., White, A. B., Ralph, F. M., Minnis, P., Comstock, J. M., Tomlinson, J. M., and Prather, K. A.: Dust and biological aerosols from the Sahara and Asia influence precipitation in the western U.S., Science, 339, 1572–1578,, 2013. a

Cziczo, D. J., Froyd, K. D., Hoose, C., Jensen, E. J., Diao, M., Zondlo, M. A., Smith, J. B., Twohy, C. H., and Murphy, D. M.: Clarifying the dominant sources and mechanisms of cirrus cloud formation, Science, 340, 1320–1324,, 2013. a

Dee, D. P., Uppala, S. M., Simmons, A. J., Berrisford, P., Poli, P., Kobayashi, S., Andrae, U., Balmaseda, M. A., Balsamo, G., Bauer, P., Bechtold, P., Beljaars, A. C. M., van de Berg, L., Bidlot, J., Bormann, N., Delsol, C., Dragani, R., Fuentes, M., Geer, A. J., Haimberger, L., Healy, S. B., Hersbach, H., Hólm, E. V., Isaksen, L., Kållberg, P., Köhler, M., Matricardi, M., McNally, A. P., Monge-Sanz, B. M., Morcrette, J.-J., Park, B.-K., Peubey, C., de Rosnay, P., Tavolato, C., Thépaut, J.-N., and Vitart, F.: The ERA-Interim reanalysis: configuration and performance of the data assimilation system, Q. J. Roy. Meteorol. Soc., 137, 553–597,, 2011. a

Delanoë, J. and Hogan, R. J.: Combined CloudSat-CALIPSO-MODIS retrievals of the properties of ice clouds, J. Geophys. Res., 115, D00H29,, 2010. a, b

Delanoë, J., Protat, A., Testud, J., Bouniol, D., Heymsfield, A. J., Bansemer, A., Brown, P. R. A., Forbes, R. M., and Delanoë, J.: Statistical properties of the normalized ice particle size distribution, J. Geophys. Res., 110, D10201,, 2005. a, b

DeMott, P. J., Rogers, D. C., and Kreidenweis, S. M.: The susceptibility of ice formation in upper tropospheric clouds to insoluble aerosol components, J. Geophys. Res., 102, 19575,, 1997. a, b

DeMott, P. J., Prenni, A. J., Liu, X., Kreidenweis, S. M., Petters, M. D., Twohy, C. H., Richardson, M. S., Eidhammer, T., and Rogers, D. C.: Predicting global atmospheric ice nuclei distributions and their impacts on climate, P. Natl. Acad. Sci. USA, 107, 11217–11222,, 2010. a, b, c, d

Deng, M., Mace, G. G., Wang, Z., and Lawson, R. P.: Evaluation of Several A-Train Ice Cloud Retrieval Products with In Situ Measurements Collected during the SPARTICUS Campaign, J. Appl. Meteorol. Clim., 52, 1014–1030,, 2013. a

Diao, M., Bryan, G. H., Morrison, H., and Jensen, J. B.: Ice Nucleation Parameterization and Relative Humidity Distribution in Idealized Squall-Line Simulations, J. Atmos. Sci., 74, 2761–2787,, 2017. a

Field, P. R., Cotton, R. J., Johnson, D., Noone, K., Glantz, P., Kaye, P. H., Hirst, E., Greenaway, R. S., Jost, C., Gabriel, R., Reiner, T., Andreae, M., Saunders, C. P. R., Archer, A., Choularton, T., Smith, M., Brooks, B., Hoell, C., Bandy, B., and Heymsfield, A.: Ice nucleation in orographic wave clouds: Measurements made during INTACC, Q. J. Roy. Meteorol. Soc., 127, 1493–1512,, 2001. a

Fusina, F., Spichtinger, P., and Lohmann, U.: Impact of ice supersaturated regions and thin cirrus on radiation in the midlatitudes, J. Geophys. Res., 112, D24S14,, 2007. a

Gasparini, B. and Lohmann, U.: Why cirrus cloud seeding cannot substantially cool the planet, J. Geophys. Res., 121, 4877–4893,, 2016. a

Gayet, J.-F., Ovarlez, J., Shcherbakov, V., Ström, J., Schumann, U., Minikin, A., Auriol, F., Petzold, A., and Monier, M.: Cirrus cloud microphysical and optical properties at southern and northern midlatitudes during the INCA experiment, J. Geophys. Res., 109, D20206,, 2004. a, b, c

Gettelman, A.: Putting the clouds back in aerosol–cloud interactions, Atmos. Chem. Phys., 15, 12397–12411,, 2015. a

Gettelman, A., Liu, X., Barahona, D., Lohmann, U., and Chen, C.: Climate impacts of ice nucleation, J. Geophys. Res., 117, D20201,, 2012. a

Ghan, S., Wang, M., Zhang, S., Ferrachat, S., Gettelman, A., Griesfeller, J., Kipling, Z., Lohmann, U., Morrison, H., Neubauer, D., Partridge, D. G., Stier, P., Takemura, T., Wang, H., and Zhang, K.: Challenges in constraining anthropogenic aerosol effects on cloud radiative forcing using present-day spatiotemporal variability, P. Natl. Acad. Sci. USA, 113, 5804–5811,, 2016. a

Ginoux, P., Prospero, J. M., Gill, T. E., Hsu, N. C., and Zhao, M.: Global-scale attribution of anthropogenic and natural dust sources and their emission rates based on MODIS Deep Blue aerosol products, Rev. Geophys., 50, 3005,, 2012. a, b

Grawe, S., Augustin-Bauditz, S., Hartmann, S., Hellner, L., Pettersson, J. B. C., Prager, A., Stratmann, F., and Wex, H.: The immersion freezing behavior of ash particles from wood and brown coal burning, Atmos. Chem. Phys., 16, 13911–13928,, 2016. a

Gryspeerdt, E., Quaas, J., and Bellouin, N.: Constraining the aerosol influence on cloud fraction, J. Geophys. Res., 121, 3566–3583,, 2016. a, b, c

Gryspeerdt, E., Quaas, J., Ferrachat, S., Gettelman, A., Ghan, S., Lohmann, U., Morrison, H., Neubauer, D., Partridge, D. G., Stier, P., Takemura, T., Wang, H., Wang, M., and Zhang, K.: Constraining the instantaneous aerosol influence on cloud albedo, P. Natl. Acad. Sci. USA, 114, 4899–4904,, 2017. a, b

Gryspeerdt, E., Quaas, J., Goren, T., Klocke, D., and Brueck, M.: An automated cirrus type classification, Atmos. Chem. Phys., Atmos. Chem. Phys., 18, 6157–6169,, 2018a. a, b, c, d, e, f, g

Gryspeerdt, E., Quaas, J., Goren, T., Klocke, D., and Brueck, M.: Identification and classification of Cirrus (IC-CIR): A cirrus classification based on satellite and reanalysis data (2003–2013), Centre for Environmental Data Analysis, date, 2018b. a

Heymsfield, A. J. and Miloshevich, L. M.: Homogeneous Ice Nucleation and Supercooled Liquid Water in Orographic Wave Clouds, J. Atmos. Sci., 50, 2335–2353,<2335:HINASL>2.0.CO;2, 1993. a

Heymsfield, A. J., Krämer, M., Luebke, A., Brown, P., Cziczo, D. J., Franklin, C., Lawson, P., Lohmann, U., McFarquhar, G., Ulanowski, Z., and Van Tricht, K.: Cirrus Clouds, Meteorol. Monogr., 58, 2.1–2.26,, 2017. a

Heyn, I., Block, K., Mülmenstädt, J., Gryspeerdt, E., Kühne, P., and Quaas, J.: Assessment of simulated aerosol effective radiative forcings in the terrestrial spectrum, Geophys. Res. Lett, 44, 1001–1007,, 2017. a

Hoose, C. and Möhler, O.: Heterogeneous ice nucleation on atmospheric aerosols: a review of results from laboratory experiments, Atmos. Chem. Phys., 12, 9817–9854,, 2012. a

Hoyle, C. R., Luo, B. P., and Peter, T.: The Origin of High Ice Crystal Number Densities in Cirrus Clouds, J. Atmos. Sci., 62, 2568–2579,, 2005. a

Huang, Y., Siems, S. T., Manton, M. J., Protat, A., and Delanë, J.: A study on the low-altitude clouds over the Southern Ocean using the DARDAR-MASK, J. Geophys. Res., 117, D18204,, 2012. a

Jensen, E. J., Toon, O. B., Tabazadeh, A., Sachse, G. W., Anderson, B. E., Chan, K. R., Twohy, C. W., Gandrud, B., Aulenbach, S. M., Heymsfield, A., Hallett, J., and Gary, B.: Ice nucleation processes in upper tropospheric wave-clouds observed during SUCCESS, Geophys. Res. Lett., 25, 1363–1366,, 1998. a

Jensen, E. J., Lawson, P., Baker, B., Pilson, B., Mo, Q., Heymsfield, A. J., Bansemer, A., Bui, T. P., McGill, M., Hlavka, D., Heymsfield, G., Platnick, S., Arnold, G. T., and Tanelli, S.: On the importance of small ice crystals in tropical anvil cirrus, Atmos. Chem. Phys., 9, 5519–5537,, 2009. a, b

Jensen, E. J., Diskin, G., Lawson, R. P., Lance, S., Bui, T. P., Hlavka, D., McGill, M., Pfister, L., Toon, O. B., and Gao, R.: Ice nucleation and dehydration in the Tropical Tropopause Layer, P. Natl. Acad. Sci. USA, 110, 2041–2046,, 2013a. a, b, c, d

Jensen, E. J., Lawson, R. P., Bergman, J. W., Pfister, L., Bui, T. P., and Schmitt, C. G.: Physical processes controlling ice concentrations in synoptically forced, midlatitude cirrus, J. Geophys. Res., 118, 5348–5360,, 2013b. a, b

Jensen, E. J., Ueyama, R., Pfister, L., Bui, T. V., Lawson, R. P., Woods, S., Thornberry, T., Rollins, A. W., Diskin, G. S., DiGangi, J. P., and Avery, M. A.: On the Susceptibility of Cold Tropical Cirrus to Ice Nuclei Abundance, J. Atmos. Sci., 73, 2445–2464,, 2016. a, b, c, d, e

Jiang, J. H., Su, H., Zhai, C., Massie, S. T., Schoeberl, M. R., Colarco, P. R., Platnick, S., Gu, Y., and Liou, K.-N.: Influence of convection and aerosol pollution on ice cloud particle effective radius, Atmos. Chem. Phys., 11, 457–463,, 2011. a

Joos, H., Spichtinger, P., Lohmann, U., Gayet, J.-F., and Minikin, A.: Orographic cirrus in the global climate model ECHAM5, J. Geophys. Res., 113, D18205,, 2008. a, b, c

Kanitz, T., Seifert, P., Ansmann, A., Engelmann, R., Althausen, D., Casiccia, C., and Rohwer, E. G.: Contrasting the impact of aerosols at northern and southern midlatitudes on heterogeneous ice formation, Geophys. Res. Lett., 38, L17802,, 2011. a

Kanji, Z. A., Ladino, L. A., Wex, H., Boose, Y., Burkert-Kohn, M., Cziczo, D. J., and Krämer, M.: Overview of Ice Nucleating Particles, Meteorol. Monogr., 58, 1.1–1.33,, 2017. a, b

Kärcher, B.: Cirrus Clouds and Their Response to Anthropogenic Activities, Curr. Clim. Change Rep., 3, 45,, 2017. a

Kärcher, B. and Lohmann, U.: A parameterization of cirrus cloud formation: Homogenous freezing of supercooled aerosols, J. Geophys. Res., 107, 4010,, 2002. a, b, c, d

Kärcher, B. and Lohmann, U.: A parameterization of cirrus cloud formation: Heterogeneous freezing, J. Geophys. Res., 108, 4402,, 2003. a, b, c, d

Kärcher, B. and Seifert, A.: On homogeneous ice formation in liquid clouds, Q. J. Roy. Meteorol. Soc., 142, 1320–1334,, 2016. a, b, c

Kärcher, B. and Ström, J.: The roles of dynamical variability and aerosols in cirrus cloud formation, Atmos. Chem. Phys., 3, 823–838,, 2003. a

Kay, J. E. and Wood, R.: Timescale analysis of aerosol sensitivity during homogeneous freezing and implications for upper tropospheric water vapor budgets, Geophys. Res. Lett., 35, L10809,, 2008. a, b

Kojima, T., Buseck, P. R., Wilson, J. C., Reeves, J. M., and Mahoney, M. J.: Aerosol particles from tropical convective systems: Cloud tops and cirrus anvils, J. Geophys. Res., 109, D12201,, 2004. a, b

Koop, T., Luo, B., Tsias, A., and Peter, T.: Water activity as the determinant for homogeneous ice nucleation in aqueous solutions, Nature, 406, 611–614,, 2000. a, b, c, d

Korolev, A. V., Emery, E. F., Strapp, J. W., Cober, S. G., and Isaac, G. A.: Quantification of the Effects of Shattering on Airborne Ice Particle Measurements, J. Atmos. Ocean. Tech., 30, 2527–2553,, 2013. a, b

Krämer, M., Schiller, C., Afchine, A., Bauer, R., Gensch, I., Mangold, A., Schlicht, S., Spelten, N., Sitnikov, N., Borrmann, S., de Reus, M., and Spichtinger, P.: Ice supersaturations and cirrus cloud crystal numbers, Atmos. Chem. Phys., 9, 3505–3522,, 2009. a, b, c, d, e, f, g

Krämer, M., Rolf, C., Luebke, A., Afchine, A., Spelten, N., Costa, A., Meyer, J., Zöger, M., Smith, J., Herman, R. L., Buchholz, B., Ebert, V., Baumgardner, D., Borrmann, S., Klingebiel, M., and Avallone, L.: A microphysics guide to cirrus clouds – Part 1: Cirrus types, Atmos. Chem. Phys., 16, 3463–3483,, 2016. a, b

Lindsey, D. T. and Fromm, M.: Evidence of the cloud lifetime effect from wildfire-induced thunderstorms, Geophys. Res. Lett., 35, L22809,, 2008. a

Liou, K.-N.: Influence of Cirrus Clouds on Weather and Climate Processes: A Global Perspective, Mon. Weather Rev., 114, 1167–1199,<1167:IOCCOW>2.0.CO;2, 1986. a

Lohmann, U. and Kärcher, B.: First interactive simulations of cirrus clouds formed by homogeneous freezing in the ECHAM general circulation model, J. Geophys. Res., 107, 4105,, 2002. a, b

Lohmann, U., Stier, P., Hoose, C., Ferrachat, S., Kloster, S., Roeckner, E., and Zhang, J.: Cloud microphysics and aerosol indirect effects in the global climate model ECHAM5-HAM, Atmos. Chem. Phys., 7, 3425–3446,, 2007. a

Luebke, A. E., Afchine, A., Costa, A., Grooß, J.-U., Meyer, J., Rolf, C., Spelten, N., Avallone, L. M., Baumgardner, D., and Krämer, M.: The origin of midlatitude ice clouds and the resulting influence on their microphysical properties, Atmos. Chem. Phys., 16, 5793–5809,, 2016. a

Mamouri, R.-E. and Ansmann, A.: Potential of polarization lidar to provide profiles of CCN- and INP-relevant aerosol parameters, Atmos. Chem. Phys., 16, 5905–5931,, 2016. a

McFarquhar, G. M., Um, J., Freer, M., Baumgardner, D., Kok, G. L., and Mace, G.: Importance of small ice crystals to cirrus properties: Observations from the Tropical Warm Pool International Cloud Experiment (TWP-ICE), Geophys. Res. Lett., 34, L13803,, 2007. a, b

Mitchell, D. L., Garnier, A., Avery, M., and Erfani, E.: CALIPSO observations of the dependence of homo- and heterogeneous ice nucleation in cirrus clouds on latitude, season and surface condition, Atmos. Chem. Phys. Discuss.,, in review, 2016. a, b

Mitchell, D. L., Garnier, A., Pelon, J., and Erfani, E.: CALIPSO (IIR-CALIOP) Retrievals of Cirrus Cloud Ice Particle Concentrations, Atmos. Chem. Phys. Discuss.,, in review, 2018. a, b

Möhler, O., Field, P. R., Connolly, P., Benz, S., Saathoff, H., Schnaiter, M., Wagner, R., Cotton, R., Krämer, M., Mangold, A., and Heymsfield, A. J.: Efficiency of the deposition mode ice nucleation on mineral dust particles, Atmos. Chem. Phys., 6, 3007–3021,, 2006. a

Morcrette, J.-J., Boucher, O., Jones, L., Salmond, D., Bechtold, P., Beljaars, A., Benedetti, A., Bonet, A., Kaiser, J. W., Razinger, M., Schulz, M., Serrar, S., Simmons, A. J., Sofiev, M., Suttie, M., Tompkins, A. M., and Untch, A.: Aerosol analysis and forecast in the European Centre for Medium-Range Weather Forecasts Integrated Forecast System: Forward modeling, J. Geophys. Res.-Atmos., 114, D06206,, 2009. a

Muhlbauer, A., Ackerman, T. P., Comstock, J. M., Diskin, G. S., Evans, S. M., Lawson, R. P., and Marchand, R. T.: Impact of large-scale dynamics on the microphysical properties of midlatitude cirrus, J. Geophys. Res., 119, 3976–3996,, 2014. a, b, c

Murray, B. J.: Inhibition of ice crystallisation in highly viscous aqueous organic acid droplets, Atmos. Chem. Phys., 8, 5423–5433,, 2008. a

O'Shea, S. J., Choularton, T. W., Lloyd, G., Crosier, J., Bower, K. N., Gallagher, M., Abel, S. J., Cotton, R. J., Brown, P. R. A., Fugal, J. P., Schlenczek, O., Borrmann, S., and Pickering, J. C.: Airborne observations of the microphysical structure of two contrasting cirrus clouds, J. Geophys. Res., 121, 13510–13536,, 2016. a

Platnick, S., Meyer, K. G., King, M. D., Wind, G., Amarasinghe, N., Marchant, B., Arnold, G. T., Zhang, Z., Hubanks, P. A., Holz, R. E., Yang, P., Ridgway, W. L., and Riedi, J.: The MODIS Cloud Optical and Microphysical Products: Collection 6 Updates and Examples From Terra and Aqua, IEEE T. Geosci. Remote, 55, 502–525,, 2017. a

Podglajen, A., Hertzog, A., Plougonven, R., and Legras, B.: Lagrangian temperature and vertical velocity fluctuations due to gravity waves in the lower stratosphere, Geophys. Res. Lett., 43, 3543–3553,, 2016. a

Podglajen, A., Bui, T. P., Dean-Day, J. M., Pfister, L., Jensen, E. J., Alexander, M. J., Hertzog, A., Kärcher, B., Plougonven, R., and Randel, W. J.: Small-Scale Wind Fluctuations in the Tropical Tropopause Layer from Aircraft Measurements: Occurrence, Nature, and Impact on Vertical Mixing, J. Atmos. Sci., 74, 3847–3869,, 2017. a

Pruppacher, H. R. and Klett, J. D.: Microphysics of Clouds and Precipitation, Kluwer Acad., Norwell, Mass., 1997. a, b

Quaas, J., Boucher, O., Bellouin, N., and Kinne, S.: Satellite-based estimate of the direct and indirect aerosol climate forcing, J. Geophys. Res., 113, D05204,, 2008. a

Rosenfeld, D., Yu, X., Liu, G., Xu, X., Zhu, Y., Yue, Z., Dai, J., Dong, Z., Dong, Y., and Peng, Y.: Glaciation temperatures of convective clouds ingesting desert dust, air pollution and smoke from forest fires, Geophys. Res. Lett., 38, 21804,, 2011. a

Salzmann, M., Ming, Y., Golaz, J.-C., Ginoux, P. A., Morrison, H., Gettelman, A., Krämer, M., and Donner, L. J.: Two-moment bulk stratiform cloud microphysics in the GFDL AM3 GCM: description, evaluation, and sensitivity tests, Atmos. Chem. Phys., 10, 8037–8064,, 2010. a

Seifert, A., Heus, T., Pincus, R., and Stevens, B.: Large-eddy simulation of the transient and near-equilibrium behavior of precipitating shallow convection, J. Adv. Mod. Earth Syst., 7, 1918–1937,, 2015. a

Sourdeval, O., Gryspeerdt, E., Krämer, M., Goren, T., Delanoë, J., Afchine, A., Hemmer, F., and Quaas, J.: Ice crystal number concentration estimates from lidar–radar satellite remote sensing – Part 1: Method and evaluation, Atmos. Chem. Phys., 18, 14327–14350,, 2018a. a, b, c, d, e, f, g, h, i, j, k, l

Sourdeval, O., Gryspeerdt, E., Krämer, M., Goren, T., Delanoë, J., Afchine, A., Hemmer, F., and Quaas, J.: Ice crystal number concentration from satellite lidar-radar observations (DARDAR-Nice), AERIS,, 2018b. a

Spichtinger, P. and Gierens, K. M.: Modelling of cirrus clouds – Part 1a: Model description and validation, Atmos Chem Phys, 9, 685–706,, 2009a. a, b, c

Spichtinger, P. and Gierens, K. M.: Modelling of cirrus clouds – Part 1b: Structuring cirrus clouds by dynamics, Atmos. Chem. Phys., 9, 707–719,, 2009b. a

Stephens, G. L., Li, J., Wild, M., Clayson, C. A., Loeb, N., Kato, S., L'Ecuyer, T., Stackhouse, P. W., Lebsock, M., and Andrews, T.: An update on Earth's energy balance in light of the latest global observations, Nat. Geosci., 5, 691–696,, 2012. a

Stier, P.: Limitations of passive remote sensing to constrain global cloud condensation nuclei, Atmos. Chem. Phys., 16, 6595–6607,, 2016. a

Tan, I., Storelvmo, T., and Choi, Y.-S.: Spaceborne lidar observations of the ice-nucleating potential of dust, polluted dust, and smoke aerosols in mixed-phase clouds, J. Geophys. Res., 119, 6653–6665,, 2014. a, b, c, d

Wang, M., Ghan, S., Ovchinnikov, M., Liu, X., Easter, R., Kassianov, E., Qian, Y., and Morrison, H.: Aerosol indirect effects in a multi-scale aerosol-climate model PNNL-MMF, Atmos Chem Phys, 11, 5431–5455,, 2011. a

Weigum, N. M.: Scales of variability of atmospheric aerosols, PhD thesis, University of Oxford, Oxford, 2014.  a

Weigum, N. M., Stier, P., Schwarz, J. P., Fahey, D. W., and Spackman, J. R.: Scales of variability of black carbon plumes over the Pacific Ocean, Geophys. Res. Lett., 39, L15804,, 2012. a

Wernli, H., Boettcher, M., Joos, H., Miltenberger, A. K., and Spichtinger, P.: A trajectory-based classification of ERA-Interim ice clouds in the region of the North Atlantic storm track, Geophys. Res. Lett., 43, 6657–6664,, 2016. a

Wilson, T. W., Murray, B. J., Wagner, R., Möhler, O., Saathoff, H., Schnaiter, M., Skrotzki, J., Price, H. C., Malkin, T. L., Dobbie, S., and Al-Jumur, S. M. R. K.: Glassy aerosols with a range of compositions nucleate ice heterogeneously at cirrus temperatures, Atmos. Chem. Phys., 12, 8611–8632,, 2012. a

Wood, R., Mechoso, C. R., Bretherton, C. S., Weller, R. A., Huebert, B., Straneo, F., Albrecht, B. A., Coe, H., Allen, G., Vaughan, G., Daum, P., Fairall, C., Chand, D., Gallardo Klenner, L., Garreaud, R., Grados, C., Covert, D. S., Bates, T. S., Krejci, R., Russell, L. M., de Szoeke, S., Brewer, A., Yuter, S. E., Springston, S. R., Chaigneau, A., Toniazzo, T., Minnis, P., Palikonda, R., Abel, S. J., Brown, W. O. J., Williams, S., Fochesatto, J., Brioude, J., and Bower, K. N.: The VAMOS Ocean-Cloud-Atmosphere-Land Study Regional Experiment (VOCALS-REx): goals, platforms, and field operations, Atmos. Chem. Phys., 11, 627–654,, 2011. a

Zhang, D., Wang, Z., Heymsfield, A., Fan, J., Liu, D., and Zhao, M.: Quantifying the impact of dust on heterogeneous ice generation in midlevel supercooled stratiform clouds, Geophys. Res. Lett., 39, L18805,, 2012. a

Zhao, B., Gu, Y., Liou, K.-N., Wang, Y., Liu, X., Huang, L., Jiang, J. H., and Su, H.: Type-Dependent Responses of Ice Cloud Properties to Aerosols From Satellite Retrievals, Geophys. Res. Lett., 45, 3297–3306,, 2018. a, b, c

Zhou, C., Penner, J. E., Lin, G., Liu, X., and Wang, M.: What controls the low ice number concentration in the upper troposphere?, Atmos. Chem. Phys., 16, 12411–12424,, 2016. a

Zuidema, P., Redemann, J., Haywood, J., Wood, R., Piketh, S., Hipondoka, M., and Formenti, P.: Smoke and Clouds above the Southeast Atlantic: Upcoming Field Campaigns Probe Absorbing Aerosol's Impact on Climate, B. Am. Meteorol. Soc., 97, 1131–1135,, 2016. a

Short summary
The concentration of ice crystals in a cloud affects both the properties and the life cycle of the cloud. This work uses a new satellite retrieval to investigate controls on the ice crystal concentration at a global scale. Both temperature and vertical wind speed in a cloud have a strong impact on the concentration of ice crystals. The ice crystal number is also related to the aerosol environment; defining this relation opens up new ways to investigate human impacts on clouds and the climate.
Final-revised paper