the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Impact of formulations of the homogeneous nucleation rate on ice nucleation events in cirrus
Peter Spichtinger
Patrik Marschalik
Manuel Baumgartner
Homogeneous freezing of solution droplets is an important pathway of ice formation in the tropopause region. The nucleation rate can be parameterized as a function of water activity, based on empirical fits and some assumptions on the underlying properties of supercooled water, although a general theory is missing. It is not clear how nucleation events are influenced by the exact formulation of the nucleation rate or even their inherent uncertainty. In this study we investigate the formulation of the nucleation rate of homogeneous freezing of solution droplets (1) to link the formulation to the nucleation rate of pure water droplets, (2) to derive a robust and simple formulation of the nucleation rate, and (3) to determine the impact of variations in the formulation on nucleation events. The nucleation rate can be adjusted, and the formulation can be simplified to a threshold description. We use a stateoftheart bulk ice microphysics model to investigate nucleation events as driven by constant cooling rates; the key variables are the final ice crystal number concentration and the maximum supersaturation during the event. The nucleation events are sensitive to the slope of the nucleation rate but only weakly affected by changes in its absolute value. This leads to the conclusion that details of the nucleation rate are less important for simulating ice nucleation in bulk models as long as the main feature of the nucleation rate (i.e. its slope) is represented sufficiently well. The weak sensitivity of the absolute values to the nucleation rate suggests that the amount of available solution droplets also does not crucially affect nucleation events. The use of only one distinct nucleation threshold function for analysis and model parameterization should be reinvestigated, since it corresponds to a very high nucleation rate value, which is not reached in many nucleation events with low vertical updrafts. In contrast, the maximum supersaturation and thus the nucleation thresholds reached during an ice nucleation event depend on the vertical updraft velocity or cooling rate. This feature might explain some high supersaturation values during nucleation events in cloud chambers and suggests a reformulation of ice nucleation schemes used in coarse models based on a purely temperaturedependent nucleation threshold.
Clouds are one of the most important components in the Earth–atmosphere system. They influence the hydrological cycle and Earth's energy balance via interaction with radiation. Clouds can cool the system by partly scattering and reflecting incoming solar radiation (albedo effect) but also warm the atmosphere by absorbing and reemitting thermal radiation as emitted by the Earth's surface (greenhouse effect). While for liquid clouds a net cooling effect can be derived, the radiative effect for clouds containing ice crystals is still under debate. In particular, for pure ice clouds (socalled cirrus clouds) at high altitudes in the low temperature range (T<235 K), the albedo effect and greenhouse effect are of the same order of magnitude but have different signs, leading to different net effects (see, for example, Fusina et al., 2007; Joos et al., 2014; Gasparini et al., 2017). Thus, details in microphysical properties of ice crystals might decide about a net warming or cooling of cirrus clouds, as can be seen in former model studies (e.g. Zhang et al., 1999). A key aspect of ice crystals is their size which directly affects the scattering and absorption of radiation. Smaller crystals scatter incoming solar light more effectively; thus the optical depth τ is directly dependent on the size, as can be seen in the usual approximation (see, for example, Fu and Liou, 1993)
where D_{e} denotes the effective diameter of the crystal, IWC is the ice water content, Δz represents the vertical extent of the cloud, and a and b are empirically derived constants. Since the available water vapour is mainly determined by thermodynamic conditions, the pathway of ice nucleation often determines the ice crystal number concentration in cirrus clouds and thus their effective size (assuming a certain amount of available water vapour).
Ice crystals can be formed by very different nucleation processes, which can be grouped into two major pathways, namely in situ and liquid origin ice formation (e.g. Krämer et al., 2016; Luebke et al., 2016; Wernli et al., 2016). The overall term “in situ formation” refers to ice formation at humidities below water saturation, whereas “liquid origin formation” subsumes all formation processes where cloud droplets are present and humidity is close to water saturation (e.g. freezing of cloud droplets); see the definition in Wernli et al. (2016). It is well known that the ice crystal number concentration varies crucially as a function of the underlying nucleation process, leading to potentially strong changes in the resulting radiative effect (see, for example, Krämer et al., 2020).
Despite the availability of many observational data and laboratory experiments (e.g. Hoose and Möhler, 2012), and also the development of new theoretical models (e.g. the soccer ball model; see Niedermeier et al., 2011), the details of ice nucleation at the molecular scale are still unknown.
A special situation occurs for the socalled homogeneous freezing of supercooled solution droplets (also short: homogeneous nucleation) at cold temperatures below 235 K. This process describes the spontaneous freezing of supercooled aqueous solution particles containing a small amount of (usually inorganic) substances. Although the details of this freezing process are also not completely understood on a molecular scale, reproducible laboratory experiments allowed for the formulation of an empirical fit for the nucleation rate (Koop et al., 2000). Such a fit bears inherent but maybe also unknown uncertainty, since we have no generally accepted theory for comparison. Other fits or a change in the fit parameters might also lead to different formulations of nucleation rates.
A priori, it is not clear how large the impact of the formulation of nucleation rates might be on simulating nucleation events in models resolving nucleation events in time. This issue is the starting point of our investigation. We want to address three different aspects. First, we want to link the former formulation by Koop et al. (2000) to recent findings on pure water in order to formulate a consistent framework for our models. Second, we want to derive a robust and simple formulation of the homogeneous nucleation rates, which can be used for analytical as well as numerical investigations. Third, we want to investigate the impact of variations of nucleation rates (based on the new formulation) on nucleation events, i.e. on the resulting ice crystal number concentrations.
From theory (e.g. Baumgartner and Spichtinger, 2019) and former idealized box model simulations (e.g. Kärcher and Lohmann, 2002; Ren and Mackenzie, 2005; Spichtinger and Gierens, 2009), we know that ice crystal numbers as produced in homogeneous nucleation events driven by a constant cooling rate (equivalent to a constant vertical velocity) crucially depend on several parameters and, thus, also affect the radiative properties of the formed ice cloud (see, for example, calculations in Krämer et al., 2020; Joos et al., 2009). Therefore, it is of high importance to understand the impact of the formulation of nucleation rates on the resulting ice crystal number concentrations.
We emphasize that all our investigations are meant in a bulk sense; i.e. only integrated quantities such as the ice crystal number and (total) ice crystal mass are considered. Using this approach, we consider the case of a newly forming cirrus cloud and do not focus on the freezing or forming details of single ice crystals.
The study is structured as follows. In the next section, we present the fit by Koop et al. (2000) and its empirical basis, as related to water theories. In Sect. 3 we describe the simple model used for idealized simulations for testing the impact of different formulations of nucleation rates. In Sect. 4 the more compact formulation of the nucleation rate along with several approximations is discussed. The consequences of using the proposed approximations are explored by idealized numerical simulations. In Sect. 5 we investigate the impact of a recently proposed formulation of the saturation vapour pressure over supercooled liquid water on the nucleation events (Nachbar et al., 2019). In Sect. 6 a new formulation of the nucleation rate based on results for freezing of pure supercooled water (Koop and Murray, 2016) is presented, and its impact on the number concentration of nucleated ice crystals is discussed. In Sect. 7 we investigate thresholds of ice nucleation as well as the peak values of supersaturation during nucleation events, Finally, we summarize the results and draw some conclusions in Sect. 8.
Nucleation events are investigated in the phase space spanned by temperature and water activity of the aqueous solution. The latter is defined as the ratio of saturation pressures of water vapour over the solution p_{sol} and pure water p_{liq}, as ${a}_{\mathrm{w}}:=\frac{{p}_{\mathrm{sol}}}{{p}_{\mathrm{liq}}}$. In this representation, the melting curve for different inorganic solutions turns out to be solely temperaturedependent; i.e. ${a}_{\mathrm{w}}^{i}\left(T\right):={a}_{\mathrm{w}}\left({T}_{\mathrm{m}}\right)=\frac{{p}_{\mathrm{ice}}\left(T\right)}{{p}_{\mathrm{liq}}\left(T\right)}$ (see Koop, 2015, his Eq. 5), where p_{ice} denotes the saturation vapour pressure over ice. The important insight here is that the freezing and nucleation events also collapse to a single line in the diagram (see Koop et al., 2000; Koop, 2004, 2015), which can be fitted by shifting the melting curve (deviation Δa_{w}∼0.305). This also means that the nucleation events do not depend on the solute, which is at least true for most inorganic substances (see, for example, Koop, 2004). Thus, the nucleation rate can be solely parameterized as a function of $\mathrm{\Delta}{a}_{\mathrm{w}}={a}_{\mathrm{w}}{a}_{\mathrm{w}}^{i}$. For the fitting procedure in Koop et al. (2000), a polynomial of degree 3 is used and results in the formulation
of the homogeneous nucleation rate coefficient J_{sol}. The nucleation rate coefficient is used to formulate the probability of freezing of aqueous solution droplets. The fit was used in the spirit of the representation of the nucleation rate for pure water as derived by Pruppacher (1995). During this time, three water theories were available, and the nucleation rate (as a cubic polynomial) was chosen according to the stability limit hypothesis (e.g. Mishima and Stanley, 1998), leading to an unlimited increase in the rate (see, for example, Pruppacher, 1995, his Fig. 3). However, meanwhile this water theory can be ruled out by experimental evidence; thus only the two other water theories remain (singularityfree hypothesis vs. liquid–liquid critical point; see Gallo et al., 2019, 2016), which do not imply an unlimited increase in nucleation rates of pure water (see, for example, Koop and Murray, 2016). Thus, the heuristic basis for choosing a cubic polynomial as a fit is not valid anymore.
Note that for atmospheric relevant conditions, both remaining water theories produce essentially the same results. Only at very low temperatures T<150 K, where highly viscous or even glassy states of water occur, is a different behaviour predicted. Such temperatures are not relevant for investigations of ice clouds in the tropopause region, where homogeneous freezing of solution droplets takes place. However, these theories provide the basis for the formulation of the saturation vapour pressure over supercooled water in no man's land (Murphy and Koop, 2005), combining heat capacities of liquid water and amorphous ice.
Finally, using the assumption of solution droplets being in equilibrium with their environment and neglecting size effects, water activity equals the liquid water saturation ratio S_{liq} due to
where p_{v} denotes the partial water vapour pressure. Using this representation of a_{w} together with the ice saturation ratio ${S}_{\mathrm{i}}=\frac{{p}_{\mathrm{v}}}{{p}_{\mathrm{ice}}}$, the computation
shows that Δa_{w} only depends on the ice saturation ratio and temperature.
Note that although recent measurements (Pathak et al., 2021) corroborate the procedure in the study by Murphy and Koop (2005), in a recent study by Nachbar et al. (2019) the combination of liquid water and amorphous ice is called into question, leading to a different formulation of the saturation vapour pressure over supercooled water and thus a different water activity. In the following investigations, we will also use this formulation in order to determine the sensitivity of the nucleation events to the choice of a saturation vapour pressure formulation. Note that for each choice, the water activity ${a}_{\mathrm{w}}^{i}\left(T\right)$ must be recalculated.
We begin with the description of the governing equations for the relevant ice processes in a nucleation event, i.e. homogeneous nucleation and diffusional growth. Both processes are key for determining the properties of the nucleation event, such as the number of nucleated ice crystals and the evolution of the ice saturation ratio (e.g. its peak value). Of course, other processes such as sedimentation and aggregation of ice crystals are important for the evolution of ice clouds but usually act on longer timescales, e.g. when the particles are grown to larger sizes. Thus, we omit these processes and concentrate on nucleation and growth, as in former studies (e.g. Kärcher and Lohmann, 2002; Baumgartner and Spichtinger, 2019).
We formulate the model in terms of averaged quantities for ice crystal mass and number concentration (q_{i}, n_{i}), i.e. as a twomoment scheme. Additionally, the saturation ratio with respect to hexagonal ice, ${S}_{\mathrm{i}}=\frac{{p}_{\mathrm{v}}}{{p}_{\mathrm{ice}}\left(T\right)}$, is used, with the partial water vapour pressure, p_{v}, and the saturation water vapour pressure over hexagonal ice, p_{ice}(T). Thus, the complete set of equations for an adiabatically ascending air parcel can be represented as
including changes of temperature T and pressure p. In these equations, w denotes the vertical velocity of the air parcel, c_{p} is the specific heat capacity of dry air (assumed as a constant; see Baumgartner et al., 2020), L denotes the (constant) latent heat of sublimation, and ρ is the air density. The assumption of an ideal gas is adopted for air and water vapour. The terms Nuc_{n} and Nuc_{q} denote changes due to nucleation, and the terms Dep_{q} and Dep_{s} describe changes due to diffusional growth of ice crystals. The term “Cool” denotes the impact of adiabatic expansion due to upward motion with velocity w; this is also reflected in the change of temperature and pressure, using adiabatic lapse rate and hydrostatic pressure, respectively. For temperature, we would have to consider diabatic changes due to latent heat release in phase changes.
Computing the total derivative of the saturation ratio using the representation ${S}_{\mathrm{i}}=\frac{p\phantom{\rule{0.125em}{0ex}}{q}_{\mathrm{v}}}{{\mathit{\epsilon}}_{\mathrm{0}}{p}_{\mathrm{si}}\left(T\right)}$, where ε_{0} denotes the ratio of molar masses of water and dry air, together with the Clausius–Clapeyron equation yields
and
To a good approximation, for cold temperatures the first term in the brackets in Eq. (11), which describes latent heat release due to phase changes, can be omitted. In the following, we will omit the evolution in Eqs. (8) and (9) for temperature and pressure; i.e. we assume these to be constant during the nucleation event. Thus, we arrive at
As a result of assuming temperature and pressure as being constant, only the vertical velocity w is an external parameter for the supersaturation. For the terms Nuc_{x} and Dep_{x} (x=n, q, s), we have to keep temperature and pressure as fixed parameters: T=T_{env} and p=p_{env}. This approach was also used in former investigations (see, for example, Spreitzer et al., 2017; Baumgartner and Spichtinger, 2019).
The nucleation term can be described as
where V_{d} is the mean volume of a supercooled solution droplet, n_{a} is the number concentration of solution droplets, and m_{0} is the mean mass of a newly frozen solution droplet, which can be set to ${m}_{\mathrm{0}}={\mathrm{10}}^{\mathrm{16}}$ kg. The nucleation rate for the homogeneous freezing of solution droplets is denoted by J_{nuc}. For comparison with former investigations (Kärcher and Lohmann, 2002; Spichtinger and Gierens, 2009), we set the number concentration of the background aerosol to quite a large value of n_{a}ρ=10^{4} cm^{−3} = 10^{10} m^{−3}; since the resulting ice crystal number concentration as produced in nucleation events is usually some orders of magnitude smaller, we do not have to care about a possible consumption of a major part fraction of solution droplets. We will later discuss the impact of this value in terms of nucleation events.
The diffusional growth of ice crystals is determined by the growth rate
with the diffusion constant for water vapour in air ${D}_{\mathrm{v}}^{\ast}={D}_{\mathrm{v}}(p,T){f}_{D}$ as corrected by the factor f_{D} for the kinetic regime, the capacity of ice crystals, C, assuming columnar shape, the Howell factor G_{v}(p,T) describing the impact of latent heat, and the ventilation correction f_{v}, respectively. Note that the capacity also depends on the mean mass of the ice crystal ensemble; i.e. $C=C\left(\stackrel{\mathrm{\u203e}}{m}\right)=C({n}_{\mathrm{i}},{q}_{\mathrm{i}})$. The details of the formulation are given in Appendix A.
Combining the expressions from above, the reduced system of equations reads
Here we make the following remarks:

As shown in Spreitzer et al. (2017), it is possible to determine and characterize the steady states of the reduced system, which additionally includes sedimentation. This leads to a nonlinear oscillator with a bifurcation diagram, depending on the updraft velocity w and on the temperature T.

The usefulness of this simple doublemoment scheme depends on the scales of the scenarios. We generally found good agreement with such parcel models and also on an LES (large eddy simulation) scale (and even coarser resolution) with observations, more sophisticated models, and also theory (see, for example, Spichtinger and Gierens, 2009; Spichtinger, 2014; Baumgartner et al., 2022).
Investigations of ice clouds in the cold temperature regime (T<235 K) need to include the nucleation process of homogeneous freezing of aqueous solution droplets. As pointed out in Sect. 1, the formulation by Koop et al. (2000) based on water activity is a meaningful fit to experimental data. However, for theoretical investigations and the use in reduced order models, a simpler but still accurate approximation would be helpful. In the following we present a way to derive such an approximation based on the original fit through measurements by Koop et al. (2000) in addition to recent observations for pure supercooled water.
4.1 Correction of the nucleation rate
In the study by Koop and Murray (2016) a parameterization of the nucleation rate of pure supercooled water J_{pure liq}(T) was derived, based on recent measurements. Thus, in the context of homogeneous freezing of solution droplets, the nucleation rate for pure water particles should coincide with the nucleation rate of solution droplets J_{sol} at water saturation; i.e. the condition
should hold for a value $\mathrm{\Delta}{a}_{\mathrm{w}}^{\ast}$ at water saturation, as claimed and already used in the study by Koop et al. (2000). However, evaluating these two formulations of the nucleation rates at water saturation shows a similar qualitative behaviour down to temperatures T∼235K but a quantitative disagreement; see the blue and black curve in Fig. 1. A reasonable requirement is that the values of both formulations should match in the temperature range 235 K ≤ T ≤ 240 K, since this range is relevant for the freezing of pure water cloud droplets with reasonable sizes. This temperature range at water saturation is equivalent to the range of water activity difference $\mathrm{0.27}\le \mathrm{\Delta}{a}_{\mathrm{w}}\le \mathrm{0.31}$. The offset between the curves can be corrected by shifting the logarithm of the nucleation rate for solution droplets by a constant value. The value of the shift was calculated by minimizing the square distance between the curves in the respective temperature range. Thus, the corrected nucleation rate for aqueous solution droplets reads
with $\mathit{\delta}=\mathrm{1.522}$. The nucleation rates are given in SI units (as used for all quantities throughout this study); i.e. $\left[J\right]={\mathrm{m}}^{\mathrm{3}}\phantom{\rule{0.125em}{0ex}}{\mathrm{s}}^{\mathrm{1}}$ (per cubic metre per second).
We make the following remarks:

The nucleation rate of pure water droplets can be used for a direct parameterization of the nucleation rate of aqueous solution droplets. This will be carried out in Sect. 6.1.

The (new) disagreement (or small shift) of the rates solely stems from the comparison with the new formulation of Koop and Murray (2016), since originally the nucleation rate for solution droplets was chosen in agreement with measurements of nucleation rates for pure water droplets (Koop et al., 2000).

In the following we will refer to the corrected nucleation rate as the “reference” nucleation rate, since, to the best of our knowledge, it provides the best and most recent fit for the homogenous nucleation rate of solution particles.
4.2 Nucleation rate as a function of T and S_{i}
The general strategy of the study is to represent the exponent of the nucleation rate by loworder polynomials in a thermodynamic variable x; i.e.
For instance, the formulation of the nucleation rate for aqueous solution droplets by Koop et al. (2000) is based on a polynomial of degree 3; i.e.
using the thermodynamic quantity $x=\mathrm{\Delta}{a}_{\mathrm{w}}={a}_{\mathrm{w}}{a}_{\mathrm{w}}^{i}$.
Note that the nucleation rate J_{pure liq} for pure water droplets is also based on the same structure; i.e. log _{10}(J_{pure liq}) is a polynomial of order 6 in the thermodynamic variable T (see Koop and Murray, 2016). For analytical investigations of the homogeneous nucleation, it is desirable to represent log _{10}(J) by a polynomial of low degree. As will be shown in the following, the formulation
with a polynomial yields sufficient agreement with the reference. For analytical investigations (e.g. using asymptotic analysis), it is helpful to represent the nucleation rate using a threshold for the humidity to account for the explosive character of nucleation events as used in the analysis by Baumgartner and Spichtinger (2019). Thus, for the nucleation rate for supercooled solution droplets, we make the following ansatz:
where S_{c}=S_{c}(T) is the temperaturedependent threshold value for the saturation ratio. Note that the ansatz is consistent (or even equivalent) with condition Eq. (22). In order to describe J as a function of S_{i} and T, we reformulate Δa_{w} as
using a threshold S_{c}(T) that corresponds to a fixed value J_{0} of the nucleation rate; i.e. $J\left({S}_{\mathrm{c}}\right(T),T)={J}_{\mathrm{0}}$. Taking the logarithm, this equality implies ${p}_{n}\left({x}_{\mathrm{0}}\right)={j}_{\mathrm{0}}={\mathrm{log}}_{\mathrm{10}}\left({J}_{\mathrm{0}}\right)$, with x=Δa_{w}. As in former studies (see, for example, Koop et al., 2000; Kärcher and Lohmann, 2002), we choose ${J}_{\mathrm{0}}={\mathrm{10}}^{\mathrm{16}}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{3}}\phantom{\rule{0.125em}{0ex}}{\mathrm{s}}^{\mathrm{1}}={\mathrm{10}}^{\mathrm{10}}\phantom{\rule{0.125em}{0ex}}{\mathrm{cm}}^{\mathrm{3}}\phantom{\rule{0.125em}{0ex}}{\mathrm{s}}^{\mathrm{1}}$. Note that this choice for the parameterization is quite arbitrary and has no strict physical interpretation, although one can argue with the cooling rates of the underlying experiments and thus with the probability of the freezing of droplets with a given volume within a certain predefined time interval (Koop et al., 2000).
Evaluating Eq. (24) at S_{i}=S_{c}, we arrive at
leading to a description of the threshold
and
if the polynomial p_{n}(x) can be inverted in the relevant range $\mathrm{0.26}\le \mathrm{\Delta}{a}_{\mathrm{w}}\le \mathrm{0.34}$. Combining the equations from above, the nucleation rate can be represented as
which is a threshold description using the thermodynamic variables S_{i} and T. This representation amounts to a reformulation of the original approximation, if the inverse function ${p}_{n}^{\mathrm{1}}\left(x\right)$ exists in the relevant range (i.e. p_{n}(x) is strictly monotonic). In the following, we consider the case of linear and quadratic polynomials, as determined by ansatz (23).

First, we consider the case of a linear polynomial ${p}_{\mathrm{1}}\left(x\right)={a}_{\mathrm{0}}+{a}_{\mathrm{1}}x$.
The inverse function of p_{1}(x)=y is given by ${p}_{\mathrm{1}}^{\mathrm{1}}\left(y\right)=\frac{y{a}_{\mathrm{0}}}{{a}_{\mathrm{1}}}$, implying the threshold$$\begin{array}{}\text{(28)}& {S}_{\mathrm{c}}\left(T\right)={\displaystyle \frac{\mathrm{1}}{{a}_{\mathrm{w}}^{i}\left(T\right)}}{\displaystyle \frac{{j}_{\mathrm{0}}{a}_{\mathrm{0}}}{{a}_{\mathrm{1}}}}+\mathrm{1}.\end{array}$$Substituting Eq. (28) into expression (27) yields
$$\begin{array}{}\text{(29)}& \begin{array}{rl}{\mathrm{log}}_{\mathrm{10}}J({S}_{\mathrm{i}},T)& =j({S}_{\mathrm{i}},T)\\ & ={j}_{\mathrm{0}}+{a}_{\mathrm{1}}{a}_{\mathrm{w}}^{i}\left(T\right)({S}_{\mathrm{i}}{S}_{\mathrm{c}}(T\left)\right)\\ & ={j}_{\mathrm{0}}+A\left(T\right)({S}_{\mathrm{i}}{S}_{\mathrm{c}}(T\left)\right),\end{array}\end{array}$$where $A\left(T\right)={a}_{\mathrm{1}}{a}_{\mathrm{w}}^{i}\left(T\right)$. The coefficients a_{0} and a_{1} can be determined in different ways; see Sect. 4.3. Furthermore, approximations to the functions A(T) and S_{c}(T) can be investigated.

Second, we consider the case of a quadratic polynomial
${p}_{\mathrm{2}}\left(x\right)={b}_{\mathrm{0}}+{b}_{\mathrm{1}}x+{b}_{\mathrm{2}}{x}^{\mathrm{2}}=a(xb{)}^{\mathrm{2}}+c$.
Since a quadratic function is not strictly monotonic in general, inverting the quadratic polynomial leads to two functions; i.e.$$\begin{array}{}\text{(30)}& {p}_{\mathrm{2}}^{\mathrm{1}}\left(y\right)=b\pm \sqrt{{\displaystyle \frac{yc}{a}}}.\end{array}$$If one solution can be ruled out (e.g. due to physical constraints), we can formulate
$$\begin{array}{}\text{(31)}& {\mathrm{log}}_{\mathrm{10}}J={p}_{\mathrm{2}}\left(({S}_{\mathrm{i}}{S}_{\mathrm{c}}(T\left)\right){a}_{\mathrm{w}}^{i}\left(T\right)+{p}_{\mathrm{2}}^{\mathrm{1}}\left({j}_{\mathrm{0}}\right)\right)\end{array}$$using the threshold description
$$\begin{array}{}\text{(32)}& {S}_{\mathrm{c}}\left(T\right)=\left(b\pm \sqrt{{\displaystyle \frac{{j}_{\mathrm{0}}c}{a}}}\right){\displaystyle \frac{\mathrm{1}}{{a}_{\mathrm{w}}^{i}\left(T\right)}}+\mathrm{1}.\end{array}$$Equivalently, we can derive a formulation
$$\begin{array}{}\text{(33)}& \begin{array}{rl}{\mathrm{log}}_{\mathrm{10}}J={c}_{\mathrm{0}}& +{q}_{\mathrm{1}}\left(T\right)({S}_{\mathrm{i}}{S}_{\mathrm{c}}(T\left)\right)\\ & +{q}_{\mathrm{2}}\left(T\right)({S}_{\mathrm{i}}{S}_{\mathrm{c}}(T){)}^{\mathrm{2}}\end{array}\end{array}$$with appropriate functions q_{1} and q_{2}, which might be useful for analytic investigations.
We will use this quadratic ansatz again for a direct approximation of the nucleation rate of pure water droplets (see Sect. 6.1).
4.3 Linear polynomial fit for the nucleation rate
In this section we investigate approximations of the exponent of the nucleation rate of aqueous solution droplets J_{sol} and their impact on nucleation events in an idealized scenario. We concentrate on the reference formulation (Koop et al., 2000). Since the polynomial p_{3}(x) in the original formulation
nearly behaves as a linear polynomial in the relevant range $\mathrm{0.26}\le \mathrm{\Delta}{a}_{\mathrm{w}}\le \mathrm{0.34}$, it can be easily approximated by a linear relation; i.e. ${p}_{\mathrm{3}}\left(x\right)\approx {b}_{\mathrm{0}}+{b}_{\mathrm{1}}x$. For this we can use two different approaches: (i) using a leastsquares fit and (ii) a Taylor expansion at a prescribed value y_{0}. While the first approach is just a fitting procedure in the relevant range $\mathrm{0.26}\le \mathrm{\Delta}{a}_{\mathrm{w}}\le \mathrm{0.34}$, the second approach relies on an a priori choice for the evaluation point ${y}_{\mathrm{0}}\in [\mathrm{0.26},\mathrm{0.34}]$, and it is not evident from the outset which value should be used to provide an accurate approximation. For this, we investigate the sensitivity of p_{3} to a small perturbation $\mathit{\epsilon}=y{y}_{\mathrm{0}}$; i.e. we consider
with the coefficients
and
The Taylor approximation provides a range for the slope of the linear approximation; these values motivate the sensitivity analysis in Sect. 4.5.2. In the relevant range $\mathrm{0.26}\le y\le \mathrm{0.34}$, for y=Δa_{w} we obtain slopes in the range $\mathrm{221}\le {b}_{t\mathrm{1}}\le \mathrm{453}$. This investigation gives us a hint about possible variations in the slope, which will be used later for the sensitivity analysis in Sect. 4.5.2.
In contrast, using a leastsquares fitting routine for $\mathrm{0.26}\le \mathrm{\Delta}{a}_{\mathrm{w}}\le \mathrm{0.34}$, we obtain a linear function
with ${b}_{\mathrm{ls},\mathrm{0}}=\mathrm{62.19267}$ and ${b}_{\mathrm{ls},\mathrm{1}}=\mathrm{254.7749}$. For further investigations, we only use the linear fit from Eq. (38). We observe that the linear fit p_{ls}(x) best approximates p_{3} close to the inflection point x_{infl}≈0.30756 (see Fig. 2, left panel).
For each linear approximation $p\left(x\right)={b}_{\mathrm{0}}+{b}_{\mathrm{1}}\cdot x$ of p_{3}(x), the exponent of the nucleation rate and the saturation ratio threshold become, as demonstrated in Sect. 4.2,
Since ${a}_{\mathrm{w}}^{i}$ is a rather complicated function of temperature, it is particularly useful in the context of analytical investigations to have simpler approximations of this quantity. This motivates us to approximate ${a}_{\mathrm{w}}^{i}$ and its inverse $\frac{\mathrm{1}}{{a}_{\mathrm{w}}^{i}}$ in the relevant temperature range $\mathrm{190}\le T\le \mathrm{230}$ K by polynomials q(T) of degree deg(q)≤2. Similarly, we can approximate the nucleation threshold S_{c}(T) by polynomials s(T) of degree deg(s)≤2. For the approximations we use a leastsquares procedure within the temperature range $\mathrm{190}\le T\le \mathrm{230}$ K. The results are presented in Fig. 2 (middle and right panels).
Combining the approximations q(T) and s(T) yields the formulation
of log _{10}(J). As can be seen in Fig. 2, the nucleation threshold is accurately approximated by a linear relation (deviation is smaller than 0.3 %). In former studies (e.g. Kärcher and Lohmann, 2002; Ren and Mackenzie, 2005) linear fits were derived for the nucleation thresholds; however, these fits deviate significantly more from the reference in comparison to ours (see Fig. 2). The deviation depends on the respective formulation (or approximation) of ${a}_{\mathrm{w}}^{i}$.
4.4 Thresholds for prescribed nucleation rate values
The threshold description in Sect. 4.3 was based on the choice j_{0}=16, corresponding to a nucleation rate $J={\mathrm{10}}^{\mathrm{16}}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{3}}\phantom{\rule{0.125em}{0ex}}{\mathrm{s}}^{\mathrm{1}}$. As already mentioned, the choice of j_{0} is quite arbitrary, and these high values of J are very often not reached in the numerical simulations (see Sect. 4.5). For a better diagnostics of the nucleation events and the relative strength of nucleation events, we introduce a similar concept for nucleation thresholds, based on a prescribed nucleation rate value $J\sim {\mathrm{10}}^{{x}_{\mathrm{0}}}$. For this purpose we use Eq. (40) of the nucleation threshold based on the linear approximation of the nucleation rate with a fixed but arbitrary value x_{0}>0 for the nucleation rate value; hence, we can write
where the function $A\left(T\right)={b}_{\mathrm{1}}{a}_{\mathrm{w}}^{i}\left(T\right)$ only depends on the linear approximation of J, as stated in Sect. 4.2. Note that obviously S_{cx0}(T)=S_{c}(T) for x_{0}=j_{0}. This leads to the formulation of the nucleation rate
with a general nucleation value x_{0} and its associated threshold function S_{cx0}(T). The threshold function is just shifted by the value $\frac{{x}_{\mathrm{0}}{j}_{\mathrm{0}}}{A\left(T\right)}$; i.e. the type of the threshold function remains the same. This formulation will be used for the theoretical investigations using small perturbations (see Sect. 4.6)
4.5 Numerical simulations of nucleation events for different approximations
In the following we investigate the impact of our approximations of log _{10}(J) on nucleation events. The setup is as follows: we use the simple bulk ice physics model as described by the set of ordinary differential equations (15), (16), and (17) in Sect. 3. A nucleation event is ensured by assuming a constant vertical velocity, which directly translates into a constant adiabatic cooling of the air parcel and, thus, an initially increasing saturation ratio. Instead of changing the temperature adiabatically, we directly control the supersaturation as described in Sect. 3; this allows us to control the nucleation event without the need to disentangle the different contributions of temperature and supersaturation.
The nucleation event always shows the same qualitative behaviour: due to the supersaturation source ∼wS_{i} with a constant updraft w, the variable S_{i} increases, and the nucleation term produces ice crystals, which can grow by water vapour diffusion, constituting a sink for supersaturation. The peak value of S_{i} is reached once the source and sink of supersaturation balance. Afterwards the variable S_{i} decreases due to diffusional growth and thus shuts off the nucleation term. The peak value depends crucially on the number of nucleated ice crystals that are needed, to balance the source for S_{i} by the diffusional growth (depending on the product of number concentration and mean radius of ice crystals). The number concentration of ice crystals produced in the nucleation event clearly depends on the vertical velocity w (source term) and the environmental conditions (diffusion depends on temperature and pressure). For details of the time evolution of nucleation events, see Appendix B.
4.5.1 Standard approximation
We compare the following four different representations of the nucleation rate using numerical simulations:

We represent the nucleation rate in the water activity formulation by Koop et al. (2000) with the correction as described in Sect. 4.1 (reference nucleation rate).

We represent the water activity approximated by the linear fit as described in Sect. 4.3 (see Eq. 38, linear regression).

We represent the nucleation rate as a function of S_{i} an T as described in Sect. 4.2, based on the formulation
$$\begin{array}{}\text{(43)}& {\mathrm{log}}_{\mathrm{10}}J={j}_{\mathrm{0}}+A\left(T\right)({S}_{\mathrm{i}}{S}_{\mathrm{c}}(T\left)\right)\end{array}$$of the exponent of the nucleation rate. We compare the following two sets of approximations for A(T) and S_{c}(T):
 a.
A linear approximation is used for A(T) and a quadratic approximation for S_{c}(T), respectively.
 b.
A constant approximation is used for A(T) and a linear approximation for S_{c}(T), respectively.
 a.
These are specific cases; however arbitrary combinations of approximations for A(T) and S_{c}(T) might be used.
Figure 3 shows the approximated exponents of the nucleation rate together with the (corrected) reference formulation by Koop et al. (2000) for the three temperatures T=196, 216, and 236 K as functions of Δa_{w}. These temperatures are chosen for direct comparison with former studies (Kärcher and Lohmann, 2002; Spichtinger and Gierens, 2009). Evidently, the linear fit with respect to water activity is very close to the reference, and the same is true for the case of a linear function A(T) and a quadratic approximation S_{c}(T). For the simplest approximation (constant function A(T) and linear approximation S_{c}(T)), larger deviations from the reference nucleation rate can be seen. At T=196 K, there is a strong underestimation in the lower range of Δa_{w}, whereas for T=236 K the underestimation is most pronounced for higher values of Δa_{w} (green vs. black curves). In both cases, we expect deviations in the number concentrations of nucleated ice crystals during the nucleation event and the maximum saturation ratio attained.
We investigate standard nucleation events in terms of (i) the resulting ice crystal number concentration at the end of the simulation as in former studies (e.g. Kärcher and Lohmann, 2002; Spichtinger and Gierens, 2009) and (ii) the maximum (peak) supersaturation, which was reached during the nucleation event. Although the latter is usually not considered, it is of interest for comparisons with real measurements, for example, in cloud chambers.
Figure 4 shows the results of the numerical simulations, i.e. the number of nucleated ice crystals (top panel) and the maximum saturation ratio (bottom panel) at environmental pressure p=200 hPa (the results are similar for other environmental conditions).
Comparing the number of nucleated ice crystals and the maximum saturation ratio, it is evident that the difference between the reference calculation, based on the corrected nucleation rate by Koop et al. (2000), and the simulations using the approximated nucleation rates is rather small.
For almost all nucleation events, the deviation from the reference simulations is not larger than ±15 %. To assess these deviations, one should keep in mind that measurements of ice crystal number concentrations are quite difficult, and the uncertainties are usually larger than 15 %. For instance, for the forwardscattering spectrometer probe instrument, which was used in many flight campaigns (e.g. Voigt et al., 2017), the uncertainty is estimated by about ∼ 10 % (de Reus et al., 2009). Thus, the deviations in our simulations and the uncertainties of realistic measurements are roughly of the same order. This fact renders it presumably impossible to decide on the correctness of any of the different formulations and approximations of the nucleation rate based on the available observations.
Finally, we conclude that a linear approximation of the reference nucleation rate by Koop et al. (2000) is accurate enough to represent nucleation events in a physically meaningful way. Thus, we can use this description as well as the derived formulations of the nucleation rate as a function of temperature T and saturation ratio S_{i} in order to investigate which parameters of the nucleation rates significantly affect the outcome of nucleation events. This will be carried out in the next section.
4.5.2 Impact of the parameters of the linear approximation
Generally, we are interested in the impact of the formulation of the nucleation rate on nucleation events. The original parameterization by Koop et al. (2000) is based on a cubic polynomial, which has slopes in the range $\mathrm{221}\le b\le \mathrm{453}$; see Sect. 4.3. The linear approximation is sufficiently good for representing the reference rate; thus, we now use this simple linear representation ${\mathrm{log}}_{\mathrm{10}}J={b}_{\mathrm{0}}+{b}_{\mathrm{1}}\cdot \mathrm{\Delta}{a}_{\mathrm{w}}$ in order to test the sensitivity of nucleation events on the two parameters b_{0}, b_{1}. Parameter b_{0} controls the absolute value of the nucleation rate, while parameter b_{1} accounts for its steepness, i.e. the slope.
In a first step, we investigate the impact of the slope of the nucleation rate given by coefficient b_{1}. One should keep in mind that during the nucleation event, the value of $\mathrm{\Delta}{a}_{\mathrm{w}}=({S}_{\mathrm{i}}\mathrm{1}){a}_{\mathrm{w}}^{i}\left(T\right)$ is increasing as S_{i} increases; thus the exponent of the nucleation rate basically grows linearly. Consequently, an increase in the saturation ratio immediately translates into an increase in Δa_{w}; hence the abscissa in Fig. 5 may be thought of as representing the saturation ratio. If high values of the nucleation rate are already reached at lower supersaturation values, the nucleation is triggered earlier in comparison to the reference scenario.
However, an earlier onset of ice nucleation implies that the newly nucleated ice crystals already start to grow by diffusion. Consequently, the growing ice crystals tend to decrease the saturation ratio and, if they are sufficiently numerous, prematurely stop the nucleation event. In this case, fewer ice crystals will nucleate, and a smaller maximum saturation ratio will be reached compared to the reference. The opposite mechanism is expected for smaller values of the nucleation ratio in comparison to the reference; i.e. higher ice crystal concentrations will occur together with a larger maximal saturation ratio.
In order to illustrate this mechanism more quantitatively, we artificially changed the slope of the linear function. The reference slope b_{1}≈255 is either reduced to a value of b_{1}=100 or enhanced to a value b_{1}=500, which is motivated by the values of the Taylor approximation, derived in Sect. 4.3. In both cases, the parameter b_{0} of the linear function is adapted such that the inflection point of the polynomial p_{3}(Δa_{w}) at Δa_{w}∼0.311 is met for better comparison with the reference simulations. The resulting nucleation rates are displayed in Fig. 5, while the number of nucleated ice crystals and the maximum ice saturation ratio during the nucleation event are summarized in Fig. 6: the top panel shows the concentrations of nucleated ice crystals, and the bottom panel shows the maximum saturation ratio during the nucleation events.
In the case of the enhanced or reduced slope, as indicated in Fig. 5, we see the theoretically proposed behaviour in the ice crystal number concentration exactly: the values are reduced for reduced slopes and enhanced for enhanced slopes. The reductions are by up to a factor of 0.4, the enhancements are by up to a factor of 2.4, and the largest changes can be seen at the highest temperature T=236 K.
In the bottom panel of Fig. 6, a dependency on temperature and vertical velocity is seen. For very low vertical velocities, the maximum supersaturation behaves as expected, i.e. reduced values for the reduced slope and enhanced values for the enhanced slope. For very high vertical velocities, this behaviour is reversed; i.e. we see reduced values of S_{i,max} for enhanced slopes and enhanced values of S_{i,max} for reduced slopes. The transition slightly depends on the temperature. This can be explained as follows: for low vertical velocities, Δa_{w} (and thus the supersaturation) is always below the inflection point Δa_{w}∼0.311. Thus the nucleation rate is always smaller for the enhanced slope in comparison to the reference, while it is always larger in comparison to the reference for the reduced slope. Therefore, in the case of an enhanced slope, the nucleation starts later compared to the reference. This leads to the behaviour as described above. However, beyond the inflection point, the behaviour is reversed, and thus the resulting maximum supersaturation is now enhanced for a reduced slope, and it is reduced for an enhanced slope. The inflection point is reached at different vertical velocities for different temperatures, i.e. for lower temperatures at lower values of w and for higher temperatures at higher values of w. Note that only the maximum supersaturation is affected upon Δa_{w} crossing the inflection point, while no influence on the number concentration of ice crystals is seen.
After having varied the slope of the nucleation rate, we now turn to its absolute values and modify coefficient b_{0}, which translates into a change of values of J by ${\mathrm{10}}^{{b}_{\mathrm{0}}}$. In order to investigate the sensitivity, we add a constant value $\mathrm{\Delta}b\in \mathit{\{}\mathrm{6},\phantom{\rule{0.125em}{0ex}}\mathrm{3},\phantom{\rule{0.125em}{0ex}}\mathrm{3},\phantom{\rule{0.125em}{0ex}}\mathrm{6}\mathit{\}}$ to the coefficient b_{0}, resulting in an increase or decrease in the absolute value of the nucleation rate by a factor of 10^{Δb}. In Fig. 7 the results in terms of ice crystal number concentration and maximum supersaturation are displayed.
Maybe surprisingly, the absolute values of the number concentrations of ice crystals in comparison to the reference formulation are not crucially affected (see Fig. 7, top panel), although some deviations occur (up to a factor of 2). The strongest deviations can be seen for warm temperatures (T=236 K) at very low vertical velocities. Overall, the relative deviations from the reference events in variables n_{i} and peak values of S_{i} are within the interval [0.4, 2], but for vertical velocities in the range $w\ge \mathrm{0.05}\phantom{\rule{0.125em}{0ex}}\mathrm{m}\phantom{\rule{0.125em}{0ex}}{\mathrm{s}}^{\mathrm{1}}$, the relative deviation is within the interval [0.8, 1.4].
Comparing the influence of a scaling of the absolute values of the nucleation rate and the steepness of the rate, we conclude that the correct steepness of the nucleation rate is much more important than the absolute value of J. Even changes by orders of magnitude in the values of the nucleation rate have a minor impact on the number of nucleated ice crystals. A similar conclusion was also drawn in the theoretical study by Baumgartner and Spichtinger (2019). In that study, the authors investigated a slightly simplified system of equations by means of asymptotic analysis. The simplified system describes the temporal evolution of the number concentration of ice crystals and the saturation ratio, and an approximate asymptotic solution was constructed. To leading order, the approximate solution for the number concentration of ice crystals was completely independent of the precise values of the nucleation rate, but the steepness contributed directly. The only necessary condition on the values of the nucleation rate was that it attains large values, i.e. significantly larger than the other coefficients within the equations.
For the maximum supersaturation values, the impact of the absolute value of J is much more pronounced. As expected, upon reduction of the nucleation rate by a factor of 10^{Δb} with $\mathrm{\Delta}b\in \mathit{\{}\mathrm{6},\phantom{\rule{0.125em}{0ex}}\mathrm{3}\mathit{\}}$ the supersaturation reaches much higher values of S_{i} until the values of the rescaled nucleation rate become large enough to initiate the nucleation of ice crystals. For the enhancement of the absolute values of the nucleation rate, the results are reversed: the maximum supersaturation is reduced, since the enhanced nucleation rate attains values that allow for the production of ice crystals for smaller saturation ratios. This behaviour is represented in the bottom panels of Fig. 7.
We remark that this idealized enhancement of the nucleation rate can also be seen in the connection with the aerosol number concentration n_{a}. A change of n_{a} by some orders of magnitudes while no changes in J are applied has the same effect as changing the absolute value of the nucleation rate (or the parameter b_{0} in the argument of the exponential function). Thus, a strong reduction or enhancement of the available solution droplets will only slightly change the amount of ice crystals in a nucleation event. Therefore, we can conclude that for a meaningful approximation of the nucleation rate, the exact number concentration of available aerosols is also not crucial for the strength of the homogeneous nucleation event but perhaps for the starting time of the event. Including size effects of the solution droplets might additionally change the picture quantitatively (see, for example, Baumgartner et al., 2020).
4.6 Impact of perturbations in S_{i} and T on the nucleation rate
In this section we investigate the impact of changes in S_{i} and/or T on the nucleation rate by employing a perturbation analysis. A short explanation of this technique is given in Appendix D. In the real atmosphere, variations of the temperature due to dynamical processes will introduce such changes, for example, from a passing or even breaking gravity wave. In numerical simulations, these variations (also often called fluctuations) are often artificially introduced (e.g. Jensen and Pfister, 2004). In any case, the impact of such changes is investigated using perturbation analysis (also called asymptotic analysis).
We start with the linear approximation of the nucleation rate as formulated in Eq. (40) with $A\left(T\right)={b}_{\mathrm{1}}{a}_{\mathrm{w}}^{i}\left(T\right)$. We can estimate the usual values of the function A(T) in the temperature range $\mathrm{190}\phantom{\rule{0.125em}{0ex}}\mathrm{K}\le T\le \mathrm{230}$ K using $\mathrm{0.51}\le {a}_{\mathrm{w}}^{i}\left(T\right)\le \mathrm{0.66}$ such that $\mathrm{129}\le A\left(T\right)\le \mathrm{169}$. For a very simple but still sufficiently accurate constant approximation of ${a}_{\mathrm{w}}^{i}\left(T\right)$, we can set ${a}_{\mathrm{w}\mathrm{0}}^{i}=\mathrm{0.574312}$ (see Fig. 2, pink line), such that $A\left(T\right)\approx {A}_{\mathrm{0}}={b}_{\mathrm{1}}{a}_{\mathrm{w}\mathrm{0}}^{i}=\mathrm{146.32}$. Finally we can state $A\left(T\right)=\mathcal{O}\left({\mathit{\epsilon}}^{\mathrm{2}}\right)$ with the usual perturbation approach ε∼0.1, such that we set $A\left(T\right)={A}^{\ast}{\mathit{\epsilon}}^{\mathrm{2}}$ with ${A}^{\ast}=\mathcal{O}\left(\mathrm{1}\right)$ as ε→0. For the nondimensionalization of the threshold function in the linear approximation ${S}_{\mathrm{c}}\left(T\right)\approx {s}_{\mathrm{0}}+{s}_{\mathrm{1}}T$, we have to estimate the order of the coefficients for the relevant temperature range. Using $\mathrm{190}\phantom{\rule{0.125em}{0ex}}\mathrm{K}\le {T}_{\mathrm{ref}}\le \mathrm{230}$ K and the definition T=T_{ref}ϑ with the nondimensional temperature ϑ, we find
with ${\mathit{\sigma}}_{\mathrm{1}}={s}_{\mathrm{1}}{T}_{\mathrm{ref}}$. Obviously, ${s}_{\mathrm{0}}={\mathit{\sigma}}_{\mathrm{0}}=\mathrm{2.27697}=\mathcal{O}\left(\mathrm{1}\right)$ and $\mathrm{0.66}\le {\mathit{\sigma}}_{\mathrm{1}}\le \mathrm{0.8}$, such that σ_{1}=𝒪(1). Using the simplest approximation A(T)=A_{0} and ${S}_{\mathrm{c}}\left(T\right)={s}_{\mathrm{0}}+{s}_{\mathrm{1}}T$ for the general formulation of the threshold function S_{cx0} (see Eq. 41), we can simplify the expression as
Using nondimensionalization, we end up with the following representation:
where σ_{x0}=s_{x0} and ${\mathit{\sigma}}_{\mathrm{1}}={s}_{\mathrm{1}}{T}_{\mathrm{ref}}$. Finally, we use the estimation ${A}_{\mathrm{0}}={A}^{\ast}{\mathit{\epsilon}}^{\mathrm{2}}$ to obtain
Since ${j}_{\mathrm{0}}=\mathcal{O}\left({\mathit{\epsilon}}^{\mathrm{1}}\right)$ and x_{0}=𝒪(ε^{β}) with $\mathit{\beta}\ge \mathrm{1}$, we find ${\mathit{\sigma}}_{x\mathrm{0}}=\mathit{\delta}+{\mathit{\sigma}}_{\mathrm{0}}=\mathcal{O}\left(\mathit{\epsilon}\right)+\mathcal{O}\left(\mathrm{1}\right)=\mathcal{O}\left(\mathrm{1}\right)$. After nondimensionalizing the argument in the nucleation rate, we can now investigate the response of the nucleation rate upon a perturbation (i) in the saturation ratio (i.e. in the same way as the numerical simulations are set up), (ii) in temperature, and (iii) in adiabatic changes of temperature driving changes in the saturation ratio simultaneously. In reality, case (iii) is almost exclusively relevant.
First, we estimate the increase of J due to variations of S_{i} at a constant temperature T=T_{ref}. For this purpose we start at a given value of the saturation ratio S_{i}, which corresponds to a certain threshold x_{0} via the relation (41). We choose this value as a reference value ${S}_{\mathrm{ref}}={S}_{\mathrm{c}x\mathrm{0}}\left({T}_{\mathrm{ref}}\right)={s}_{\mathrm{c}x\mathrm{0}}\left(\mathrm{1}\right)={\mathit{\sigma}}_{\mathrm{0}}{\mathit{\sigma}}_{\mathrm{1}}$; this corresponds to a reference value of the nucleation rate $J={J}_{\mathrm{ref}}={J}_{\mathrm{unit}}\cdot {\mathrm{10}}^{{x}_{\mathrm{0}}}$ (with ${J}_{\mathrm{unit}}=\mathrm{1}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{3}}\phantom{\rule{0.125em}{0ex}}{\mathrm{s}}^{\mathrm{1}}$). Assuming the expansion
for the saturation ratio, where ${S}_{\mathrm{0}}={S}_{\mathrm{ref}}{\mathit{\sigma}}_{x\mathrm{0}}{\mathit{\sigma}}_{\mathrm{1}}$, we investigate the impact of such a perturbation on the exponent j. Keeping the temperature fixed as in the numerical simulations, we arrive at
We are interested in the relative change of the nucleation rates $\frac{J({S}_{\mathrm{i}},T)}{{J}_{\mathrm{ref}}}$, which translates into $j({S}_{\mathrm{i}},\mathrm{1})j({S}_{\mathrm{ref}},\mathrm{1})$. By definition, we have ${x}_{\mathrm{0}}=j({S}_{\mathrm{ref}},\mathrm{1})$; thus we obtain
Inspecting Eq. (50), it is evident that a nonzero perturbation term S_{α} in Eq. (48) is connected with the factor ε^{α−2}; hence a change of order 𝒪(ε^{α}) in supersaturation translates into a change of order 𝒪(ε^{α−2}) in the exponent of J. For instance, a change by S_{i}∼0.01 translates into a change of 𝒪(1) in j and thus in a change by a factor of 10 in the nucleation rate J.
Second, we consider perturbations of temperature without changing the saturation ratio, although this might not happen in the atmosphere. Using the approach above with a constant reference value of saturation, i.e. ${S}_{\mathrm{ref}}={s}_{\mathrm{c}x\mathrm{0}}\left(\mathrm{1}\right)={\mathit{\sigma}}_{\mathrm{0}}{\mathit{\sigma}}_{\mathrm{1}}$ and temperature perturbations $\mathit{\vartheta}=\mathrm{1}+\mathit{\epsilon}{\mathit{\vartheta}}_{\mathrm{1}}+{\mathit{\epsilon}}^{\mathrm{2}}{\mathit{\vartheta}}_{\mathrm{2}}+{\mathit{\epsilon}}^{\mathrm{3}}{\mathit{\vartheta}}_{\mathrm{3}}+\mathcal{O}\left({\mathit{\epsilon}}^{\mathrm{4}}\right)$, we find the following expression:
The relative change of the nucleation rate is then given by
Thus, a temperature perturbation ϑ_{α} of order 𝒪(ε^{α}) leads to a relative change in j of order 𝒪(ε^{α−2}). Note the sign of the perturbations, which turns into the opposite sign in the change of j. Because of the strictly monotonic decrease of the threshold function S_{cx0}(T), a negative temperature change leads to a higher threshold and in turn to a lower nucleation rate at a given saturation ratio.
Instead of perturbing the saturation ratio and the temperature individually, these quantities are connected in the real world. To a good approximation, their joint variation is through a purely adiabatic change. Therefore, we finally investigate the impact of adiabatic temperature changes on the saturation ratio and in turn on the nucleation rate. For this purpose, we have to consider the dependence of S_{i} on adiabatic temperature changes. We start with the cooling source term of the saturation ratio
The term $\mathit{\gamma}\left(T\right)=\frac{\mathrm{1}}{\mathit{\kappa}}\frac{L}{{R}_{\mathrm{v}}T}$ within the brackets has the values $\mathrm{28.8}\le \mathit{\gamma}\left(T\right)\le \mathrm{23.2}$ for $\mathrm{190}\phantom{\rule{0.125em}{0ex}}\mathrm{K}\le T\le \mathrm{230}$ K, such that we find $\mathit{\gamma}\left(T\right)=\mathcal{O}\left({\mathit{\epsilon}}^{\mathrm{1}}\right)={\mathit{\gamma}}^{\ast}{\mathit{\epsilon}}^{\mathrm{1}}$ and ${\mathit{\gamma}}^{\ast}\phantom{\rule{0.25em}{0ex}}\sim \mathrm{2.5}<\mathrm{0}$. Approximating the total differential in Eq. (53) with finite differences ΔS_{i} and ΔT, we arrive at
We set as an approximation S_{i}=S_{ref} and T=T_{ref}, such that we can set
with S_{l}=𝒪(1). We assume l≥1 since we do not consider changes of the saturation ratio of order 𝒪(1). The analogous expansion for the temperature reads
Combining these expansions, Eq. (53) becomes
or equivalently
The only nontrivial balance is achieved for $l=k\mathrm{1}$; i.e.
Note that $l=k\mathrm{1}\ge \mathrm{1}$; i.e. we have to consider k≥2 for the perturbation of temperature. This is a meaningful restriction since we are interested in small changes of temperature in the cold temperature regime, i.e. a change in the temperature in the order of ∼1 K in physical units. Hence, we would not expect adiabatic temperature changes of order 𝒪(ε), corresponding to changes of order ∼10 K. Thus, we assume an asymptotic expansion
for the temperature perturbation. We are generally interested in adiabatic expansions due to vertical upward motion, which in turn leads to decreasing temperatures; hence we conclude ϑ_{k}<0 for k≥2. Since ${\mathit{\gamma}}^{\ast}<\mathrm{0}$, Eq. (59) leads to positive changes in the saturation ratio s_{k}>0 for ϑ_{k}<0. Generally, warming due to adiabatic compression can be studied in the same way by setting ϑ_{k}>0.
Now we consider the nucleation rate in the formulation for arbitrary thresholds x_{0} in the nucleation rate using ${S}_{\mathrm{ref}}={\mathit{\sigma}}_{x\mathrm{0}}{\mathit{\sigma}}_{\mathrm{1}}={s}_{\mathrm{c}x\mathrm{0}}\left(\mathrm{1}\right)$:
Thus, for k≥2 we find terms of the form ${A}^{\ast}\left({s}_{\mathrm{ref}}{\mathit{\gamma}}^{\ast}{\mathit{\vartheta}}_{k+\mathrm{1}}+{\mathit{\sigma}}_{\mathrm{1}}{\mathit{\vartheta}}_{k}\right){\mathit{\epsilon}}^{k\mathrm{2}}$ of order 𝒪(ε^{k−2}). Comparing the nucleation rates, we find for the relative change
For the relative impact of these terms, we use the estimations ${\mathit{\gamma}}^{\ast}<\mathrm{2.3}$ and σ_{1}≤0.8. We have to distinguish two scenarios for perturbations ϑ_{k}<0:

ϑ_{k}<0 for all k≥2. In this case we can assume
$$\begin{array}{}\text{(63)}& {S}_{\mathrm{ref}}{\mathit{\gamma}}^{\ast}{\mathit{\vartheta}}_{k+\mathrm{1}}+{\mathit{\sigma}}_{\mathrm{1}}{\mathit{\vartheta}}_{k}>\mathrm{0}.\end{array}$$Therefore, an adiabatic temperature perturbation ϑ_{k} of order 𝒪(ε^{k}) (k≥2) leads to relative changes in j of order 𝒪(ε^{k−3}). Note that the changes in saturation ratio are always dominant and larger than the changes in the threshold, which changes j by order 𝒪(ε^{k−2}) in the opposite direction.

ϑ_{k}<0 and ${\mathit{\vartheta}}_{k+\mathrm{1}}=\mathrm{0}$ for a distinct k≥2. In this case, the previously discussed temperature effect can be seen; i.e. the nucleation threshold is changed, leading to a reduction of the nucleation rate exponent. This effect is merely academic, since we have to switch off higher perturbations in temperature, which is quite unlikely.
One should keep in mind that we investigated the relative increase in the exponent of the nucleation rate. A relative change of order 𝒪(ε^{k}) in the exponent translates into a relative change of order $\mathcal{O}\left({\mathrm{10}}^{{\mathit{\u03f5}}^{k}}\right)$ in the nucleation rate J, thus ranging over several orders of magnitudes. For instance, in the first scenario, changes of temperature of order ∼1K lead to changes in j of about ∼10, which in turn translate into a change of the nucleation rate J by a factor of 10^{10}.
Overall, we can state that changes in S_{i} are most important for changing j, either stemming from adiabatic temperature changes or driven directly as in our numerical studies.
Since the formulation of the nucleation rate by Koop et al. (2000) relies on the water activity, and thus on the function ${a}_{\mathrm{w}}^{i}\left(T\right)=\frac{{p}_{\mathrm{ice}}\left(T\right)}{{p}_{\mathrm{liq}}\left(T\right)}$, the saturation vapour pressure over liquid water (i.e. in no man's land) plays an important role. In this section we investigate the impact of choosing another formulation for p_{liq}(T) on the nucleation rate and thus on the nucleation events.
5.1 New representation of saturation water vapour
In the formulation by Murphy and Koop (2005), the extrapolation of the saturation vapour pressure into the no man's land of the phase diagram of water is based on the assumption that the state of amorphous ice is thermodynamically equivalent to supercooled liquid water. Therefore, the specific heat of liquid water can be extended in the supercooled regime using measurements of amorphous ice. This leads to the established formulation in Murphy and Koop (2005).
Recently, a new representation of the saturation vapour pressure over supercooled liquid water was proposed by Nachbar et al. (2019). In this study, the authors consider different states of water in the low temperature range. They conclude that amorphous ice is thermodynamically different from supercooled water; thus they provide a different extrapolation for the saturation vapour pressure (Nachbar et al., 2019).
Although the deviation between the two curves is very small – even in the low temperature range less than 10 % – its impact on saturation ratios as well as on the nucleation thresholds is quite large, as can be seen in Fig. 8.
The curves of water saturation as well as the nucleation thresholds are systematically shifted to higher values. In addition, the new curves have a more linear shape than the curves resulting from Murphy and Koop (2005). The ratio of the saturation pressures over ice and liquid (i.e. the functions ${a}_{\mathrm{w}}^{i}\left(T\right)$ and $\frac{\mathrm{1}}{{a}_{\mathrm{w}}^{i}\left(T\right)}$ behave differently), ${a}_{\mathrm{w}}^{i}\left(T\right)$, is much closer to a quadratic curve, as can be seen in the top panel of Fig. 9. These new fits were used for the formulation of the approximated nucleation rate. Thus, we do not change the general approach for approximating the nucleation rate, etc.; we only use a different representation of the function ${a}_{\mathrm{w}}^{i}\left(T\right)$.
5.2 Numerical simulations of nucleation events
In Fig. 10 the results of the nucleation events using the new representation of the saturation vapour pressure due to Nachbar et al. (2019) are shown. As for former experiments, the ice crystal number concentration (top panel) and the maximum supersaturation values (bottom panel) are shown.
For the ice crystal number concentration, the impact of the new formulation of p_{liq} is small; the relative deviation from the reference simulations using the original vapour pressure formulation by Murphy and Koop (2005) is always smaller than 15 %. The deviation increases with decreasing temperature and is most prominent for lower vertical updrafts ($w<\mathrm{1}\phantom{\rule{0.125em}{0ex}}\mathrm{m}\phantom{\rule{0.125em}{0ex}}{\mathrm{s}}^{\mathrm{1}}$).
For the maximum saturation ratio, the change as compared to the reference simulations is much more prominent. As can be seen in Fig. 8 the nucleation thresholds for a value of $J={\mathrm{10}}^{\mathrm{16}}{\mathrm{m}}^{\mathrm{3}}\phantom{\rule{0.125em}{0ex}}{\mathrm{s}}^{\mathrm{1}}$ are increasing with decreasing temperature with a larger slope compared to the reference case. This behaviour can clearly be seen in the maximum supersaturation; for decreasing temperature the maximum supersaturation is increasing to higher values in comparison to the reference simulations (Fig. 10, bottom panel). The increase does not depend on the vertical velocities.
We remark that at the moment it is not clear which thermodynamic hypothesis and thus which resulting approximation for the saturation vapour pressure over liquid water is physically correct; however the formulation by Murphy and Koop (2005) seems to agree with recent measurements (Pathak et al., 2021). In particular, it is not clear if the formulation of Nachbar et al. (2019) can be extrapolated to values T<200 K. Thus, we cannot recommend using a certain formulation.
Up to now we have always employed the reference nucleation rate in our computations, i.e. the formulation as in Koop et al. (2000) but corrected by a constant offset (see Sect. 4.1), in order to match the nucleation rate for pure water droplets by Koop and Murray (2016) in a certain temperature range. In this section we take a different point of view, assuming that we can just directly adopt the formulation by Koop and Murray (2016) for the nucleation rate of aqueous solution droplets, providing an exact match of both curves by definition. In the following we discuss the consequences of using such a direct approach in terms of nucleation events.
6.1 Direct fit to nucleation rate of pure water
In order to arrive at a direct fit, we assume that at water saturation, the freezing of pure water droplets should behave as the freezing of solution droplets at supercooled states. To avoid a complicated reformulation of the formula from Koop and Murray (2016), in terms of the water activity Δa_{w}, we use a quadratic polynomial fit to the original formulation J_{hom}(T) at water saturation; see Appendix C for details. Figure 11 presents the original data (black curve) together with the quadratic fit (red curve) and the corrected formulation of Koop et al. (2000) (blue curve) from Sect. 4.1. In contrast to the (corrected) formulation by Koop et al. (2000), the nucleation rate reaches a maximum at Δa_{w}∼0.345 and decreases afterwards. As a result, there is a significant deviation between the two nucleation rates (J_{K2000}>J_{KM2016}) for the range Δa_{w}>0.32. Thus, we can expect that for cold temperatures and/or high upward motions, there will be large deviations in the ice crystal number concentrations produced within nucleation events.
It should be kept in mind that the range of the parameterization of the nucleation rate as given in Koop et al. (2000) is restricted to the interval $\mathrm{0.26}\le \mathrm{\Delta}{a}_{\mathrm{w}}\le \mathrm{0.34}$. As a result, it is not clear if the parameterization works well for values Δa_{w}>0.34. However, there are measurements (see Laksmono et al., 2015) for the freezing of pure water droplets that also show a kind of plateau at cold temperatures (corresponding to high values of Δa_{w}). Thus, for higher values Δa_{w}>0.34, we use the value J_{KM2016}(Δa_{w})=J_{KM2016}(0.34) to (a) mimic the plateau in the measurements and (b) avoid numerical issues in the simulations. Assuming that the nucleation rate does not depend on quantities other than water activity, it may now be used in numerical simulations of homogeneous nucleation events.
6.2 Numerical simulations of nucleation events
After having obtained the direct formulation of the nucleation rate of Koop and Murray (2016), we now investigate its impact on nucleation events using numerical simulations as before. For completeness, two different types of simulations are done: (1) simulations using the standard formulation of p_{liq} by Murphy and Koop (2005) and (2) simulations using the new formulation of p_{liq} by Nachbar et al. (2019). The results of the simulations are shown in Fig. 12.
First we consider the ice crystal number concentrations (top panel). For low vertical updrafts, the values of n_{i} are only slightly affected in case of using the adapted nucleation rate. For higher vertical velocities, there is a reduction in the produced ice crystal number concentrations; this reduction increases with increasing vertical updrafts. This effect can be explained as follows. The nucleation rates differ significantly for higher values Δa_{w}≥0.31; i.e. the slope of the adapted rate is (much) smaller than the original nucleation rate by Koop et al. (2000). For higher updrafts, the supersaturation reaches higher values, which is equivalent to higher values of Δa_{w}. Thus, the nucleation rates differ for these high updraft events, and fewer ice crystals are produced when using the adapted nucleation rate. Apart from the influence at high vertical velocities, there is almost no difference in the ice crystal number concentrations between the nucleation events using the different formulations.
Considering the values of maximum supersaturation (bottom panel), there is a similar behaviour as for n_{i}. At low vertical velocities there is almost no difference between the reference nucleation rate and the newly adapted rate. In the case of using the saturation vapour pressure according to Nachbar et al. (2019), the observed shift in the maximum supersaturation values stems from the increased difference between the values of the saturation vapour pressures at low temperatures; see Sect. 5. At higher updrafts ($w>\mathrm{0.5}\phantom{\rule{0.125em}{0ex}}\mathrm{m}\phantom{\rule{0.125em}{0ex}}{\mathrm{s}}^{\mathrm{1}}$), the maximum supersaturation values increase nonlinearly. For the coldest temperature (T=196 K), we note a dramatic increase up to very high values (${S}_{\mathrm{i},\mathrm{max}}\sim \mathrm{1.8}$). However, note that in all cases the values of the maximum supersaturation stays below water saturation; hence no liquid origin ice formation would occur.
For the evaluation of measurements of ice clouds, the possible range of supersaturation is often estimated using the socalled Koop line, i.e. the supersaturation threshold S_{c}(T) which corresponds to a nucleation rate value $J={\mathrm{10}}^{\mathrm{16}}{\mathrm{m}}^{\mathrm{3}}\phantom{\rule{0.125em}{0ex}}{\mathrm{s}}^{\mathrm{1}}={\mathrm{10}}^{\mathrm{10}}{\mathrm{cm}}^{\mathrm{3}}\phantom{\rule{0.125em}{0ex}}{\mathrm{s}}^{\mathrm{1}}$. In many investigations (see, for example, Krämer et al., 2009), this function is used as an upper bound for possible values of S_{i} inside and also outside of ice clouds. However, from our investigations in this study so far, we have to carefully consider two different aspects from a purely theoretical point of view:

The nucleation threshold assigned to the frequently used value j_{0}=16 is arbitrary chosen; there is no convincing physical justification for using this particular value. In Koop et al. (2000) different values $J={\mathrm{10}}^{{j}_{\mathrm{0}}}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{3}}\phantom{\rule{0.125em}{0ex}}{\mathrm{s}}^{\mathrm{1}}$, with ${j}_{\mathrm{0}}\in [\mathrm{1},\phantom{\rule{0.125em}{0ex}}\mathrm{17}]$, are used, but for testing the impact of droplet sizes, they used the value j_{0}=16. Nucleation of ice crystals is not a switching process; it occurs gradually and smoothly, although the nucleation rates are very steep functions of the supersaturation. The size or strength of the nucleation event cannot be determined just by the maximum of the supersaturation; the amount of ice crystals as formed in the nucleation event is determined by the integral over the supersaturation curve (see, for example, the discussion in Dinh et al., 2016). Thus, it is possible to form many crystals in lower updrafts, even if the high nucleation threshold is not reached. From our simulations, we observe that the peak supersaturation for nucleation events depends crucially on the vertical velocity, i.e. on the temperature rate, which is prescribed during the event. This is quite obvious from the differential equation determining the change of S_{i}: the peak value is given by $\frac{\mathrm{d}{S}_{\mathrm{i}}}{\mathrm{d}t}=\mathrm{0}$, i.e. when source and sink terms balance each other. Since the source includes the vertical velocity linearly, the dependence of the peak supersaturation on w is obvious, although it is not linear.

As described above in Sect. 5, it is still not clear which formulation of the saturation vapour pressure is physically correct. However, the use of the formulation by Nachbar et al. (2019) leads to a higher saturation vapour pressure and thus to a higher nucleation threshold, even for arbitrary values j_{0} and its associated nucleation threshold S_{cx0}(T).
Taking these two aspects into account, we can observe the following behaviour. In Fig. 13 (top panel), we compare the nucleation thresholds for the saturation vapour pressure according to Murphy and Koop (2005) for j_{0}=10 (red curve) and j_{0}=16 (dark blue curve), with the range of peak supersaturations for vertical velocities $\mathrm{0.01}\phantom{\rule{0.125em}{0ex}}\mathrm{m}\phantom{\rule{0.125em}{0ex}}{\mathrm{s}}^{\mathrm{1}}\le w\le \mathrm{2}\phantom{\rule{0.125em}{0ex}}\mathrm{m}\phantom{\rule{0.125em}{0ex}}{\mathrm{s}}^{\mathrm{1}}$ (black vertical bar) and the maximum value for a very unrealistic value $w=\mathrm{10}\phantom{\rule{0.125em}{0ex}}\mathrm{m}\phantom{\rule{0.125em}{0ex}}{\mathrm{s}}^{\mathrm{1}}$ (black crosses). For comparison, the wellknown Koop line as fit and proposed by Kärcher and Lohmann (2002) is plotted (light blue curve). It is quite obvious that for typical vertical velocity values the “classical” Koop line is not reached; i.e. the peak supersaturation is below the threshold. Nevertheless, for strong cooling rates (very high vertical velocities), as are used in experiments in cloud chambers, high supersaturations are reached, which still partly remain below the Koop line. If we change the saturation vapour pressure to the formulation by Nachbar et al. (2019), the qualitative picture remains the same (bottom panel in Fig. 13): even for high vertical updrafts the high nucleation rates are reached; for moderate and small updrafts, the peak supersaturation stays well below the classical nucleation threshold. However, the nucleation thresholds are generally shifted to higher values of supersaturation due to the different saturation vapour pressure formulation. It seems that these values fit better to the experiments in the AIDA cloud chamber as reported in Baumgartner et al. (2022) and Schneider et al. (2021). This might be interpreted as a hint that the formulation by Nachbar et al. (2019) might be the more appropriate formulation for the saturation vapour pressure, although the formulation by Murphy and Koop (2005) agrees well with recent measurements (Pathak et al., 2021). In any case, one has to consider the impact of the cooling rate on the peak supersaturation in a nucleation event. Therefore, the use of the Koop line in the currently applied way is misleading and does not correspond to the actual physics of nucleation events. Note that the temperaturedependent threshold is used in some parameterizations of ice clouds in climate and numerical weather prediction models (see, for example, Kärcher et al., 2006; Köhler and Seifert, 2015). A simple but albeit more realistic extension of such schemes would be a threshold depending on both vertical velocity w and temperature T; a 2D fit to the maximum supersaturation data from our simulations might be a first attempt in this direction.
Finally, we can also investigate the peak supersaturation values for the new empirical nucleation rate formulation, as derived in Sect. 6.1. Generally, we see the same behaviour as for the reference simulations, with a monotonic increase of peak supersaturation values with increasing vertical velocity (see Fig. 14). The use of the saturation vapour pressure formulation by Nachbar et al. (2019) additionally enhances the peak values as seen before. However, the peak values for cold temperatures and very high vertical velocities are strongly enhanced in comparison with the reference simulations. Also these high values are still in line with the measurements in the AIDA chamber, as reported by Baumgartner et al. (2022) and Schneider et al. (2021).
We have investigated the impact of the representation of nucleation rates and diffusional growth on idealized nucleation events, as driven by a constant vertical updraft (i.e. a constant cooling rate). In a first step, we have investigated the original formulation of the nucleation rate for homogeneous freezing of aqueous solution droplets in the formulation by Koop et al. (2000); for a better agreement with the nucleation rate of pure water droplets, a simple shift could be applied. For analytical purposes and simple model calculations, a less complicated formulation is desired. We showed that a linear fit to the original formulation depending on the difference in water activity $\mathrm{\Delta}{a}_{\mathrm{w}}={a}_{\mathrm{w}}{a}_{\mathrm{w}}^{i}$ is accurate enough to reproduce the ice crystal number concentrations quantitatively. Based on this linearization approach, we derived a threshold formulation of the nucleation rate, which can be used for analytical investigations, as already presented in Baumgartner and Spichtinger (2019). Again, the new formulations are good enough to represent nucleation events quantitatively as compared to the reference nucleation formulation.
Using the linear approximation as a starting point, we investigated the impact of different formulations on idealized nucleation events, changing the two relevant parameters (slope and constant offset). These investigations led to the first major results:

The absolute values of the nucleation rate only have a marginal impact on the resulting ice crystal number concentrations in a nucleation event. Even a scaling by up to 6 orders of magnitudes did not severely affect the resulting number concentrations. However, the maximum supersaturations changed, and the resulting deviations range up to a few percent of relative humidity. In addition, the time of nucleation onset is slightly shifted.

The slope of the nucleation rate (or more precisely in the argument of the exponential function) has a much larger impact on the resulting nucleation event and the ice crystal number concentration. Variations in the slope change the number concentrations in the nucleation events by up to a factor of 2.5 (in both directions). Also, the maximum supersaturation is affected by a deviation of a few percent of relative humidity.
As a final conclusion of this part of our work, we can state that the shape of the nucleation rate is of high importance for the representation of the nucleation process, whereas the absolute strength of the rate is almost negligible, if the values are high enough. This shows that the nucleation process (homogeneous freezing of solution droplets) itself is a quite robust process; thus the accurate formulation is maybe less critical than we thought. Also the amount of available solution droplets as controlled by the background aerosol does not affect the nucleation events itself; it can be seen as a scaling factor of the nucleation rate, in the same sense as in the sensitivity analysis of changing the absolute values of nucleation rates. As long as the amount of aerosol particles is some orders of magnitude larger than the ice crystal number concentration as predicted for a nucleation event, this does not play a role for the nucleation events, and we do not have to care about exhausting the reservoir of solution droplets.
We also investigated the impact of a recently published formulation of the water saturation pressure based on a thermodynamic assumption of different phases of water in the very low temperature range (Nachbar et al., 2019). This new formulation leads to changes in the function ${a}_{\mathrm{w}}^{i}$, which directly affected the nucleation rate based on Δa_{w}. Following the derivations of the threshold description, approximations were constructed. The new resulting functions ${a}_{\mathrm{w}}^{i}\left(T\right)$ and S_{c}(T) can be accurately approximated with polynomials of smaller degrees, as compared to the standard formulation. The new formulation of p_{liq} only marginally changed the resulting ice crystal number concentrations. However, the impact on the maximum supersaturations increased with decreasing temperature up to few percent of relative humidity. Overall, the two different representations of the saturation vapour pressure over liquid water produced very similar, even almost identical, results. Thus, a decision about the validity of a certain formulation must be left to extensive experimental measurements.
In a more speculative part of the study, we adapted the nucleation rate of homogeneous freezing of pure water droplets (Koop and Murray, 2016) as a new parameterization for homogeneous freezing of aqueous solution droplets. This representation is quite similar for low values of Δa_{w} to the original formulation by Koop et al. (2000) and its approximations. However, for very high water activities (i.e. high supersaturations as driven by large vertical updrafts), there is a significant deviation from the reference nucleation rate. Thus, for some cases in the parameter space (high updrafts and low temperatures), there is a significant deviation in the number concentrations and, more obviously, in the maximum supersaturations, which almost reach water saturation in some cases. This approach showed that the shape of the nucleation rate is important for the resulting nucleation events; strong deviations of the shape from its reference affect the results of the nucleation event significantly. Whether this representation of the nucleation rate is a more accurate approximation to the actual physics of ice nucleation remains an open question and might be an objective for experimental investigations.
Finally, we investigated the commonly used threshold for homogeneous nucleation (Koop line) in the light of peak supersaturation values during nucleation events. This threshold corresponds to a nucleation rate of $J={\mathrm{10}}^{\mathrm{16}}{\mathrm{m}}^{\mathrm{3}}\phantom{\rule{0.125em}{0ex}}{\mathrm{s}}^{\mathrm{1}}$ but is only rarely reached during nucleation events. Nucleation itself starts usually at much lower values of S_{i}, corresponding to lower values of the nucleation rate. The peak supersaturation during a nucleation event, characterized as an equilibrium between sources and sinks of supersaturation, depends on temperature and vertical velocity. The peak supersaturation is a much more physical quantity to investigate the strength of a nucleation event. The peak supersaturation as diagnosed from the numerical simulations might be a more physical representation of ice nucleation in coarseresolution models in comparison to the frequently used nucleation threshold.
It should be emphasized that all the results and conclusions are meant in a bulk sense, i.e. for a large collection of ice crystals such as a newly forming cirrus cloud. If one is interested in the details of ice formation for a single or only a small number of particles, then all details of the nucleation rate might be equally important. In that respect, our study shows that homogeneous cirrus formation is a robust physical process.
In this appendix, we present the details of the model as used for the numerical simulations of the nucleation events. Note that we use the mathematical (and also programming) notation of logarithms; i.e. log denotes the natural logarithm (to base e).
A1 Background aerosol
For the aqueous solution droplets in the tropopause region, we assume a size distribution of lognormal type:
with a modal radius ${r}_{\mathrm{sol}}=\mathrm{75}\cdot {\mathrm{10}}^{\mathrm{9}}\phantom{\rule{0.125em}{0ex}}\mathrm{m}$ and a geometric standard deviation σ_{r}=1.5. These values are adapted from the more complex model by Spichtinger and Gierens (2009), using the fact that the dry aerosol population, as used in Spichtinger and Gierens (2009), has grown to larger sizes by water vapour uptake (i.e. assuming Köhler theory; see, for example, Köhler, 1936). The mean volume of the solution droplets
is calculated from the third moment of the lognormal distribution.
A2 Mass distribution for ice crystals
For the ice crystals, we assume a mass distribution of lognormal type
with a parameter
representing the width of the distribution, as described in Spichtinger and Gierens (2009). This distribution is used for the derivation of the rates in the system of ordinary differential equations for the mean quantities of ice mass and number concentration. The integration of weighting functions of the type ${m}^{k},k\in {\mathbb{R}}_{+}$ leads to general moments, which can be computed analytically:
Note that for the averaged quantities, we obtain n_{i}=μ[m]_{0}, q_{i}=μ[m]_{1}, respectively. Thus, we use a doublemoment scheme in our model.
A3 Diffusion constant
For the diffusion of water vapour in dry air, we use the following expression
which is an empirical fit to measurement data (Hall and Pruppacher, 1976). Note that the valid temperature range is different in the book (Pruppacher and Klett, 2010) and in the original article (Hall and Pruppacher, 1976). For analytical investigations, a representation using a quadratic temperature dependence constitutes a good approximation for a restricted temperature range. For the kinetic correction, we use the function
where r denotes the radius of the ice crystal (using a bulk density of ice ${\mathit{\rho}}_{\mathrm{b}}=\mathrm{0.81}\phantom{\rule{0.125em}{0ex}}\mathrm{kg}\phantom{\rule{0.125em}{0ex}}{\mathrm{m}}^{\mathrm{3}}$), and the parameters are given by
using the mean free path of water molecules in air λ (acc. to Pruppacher and Klett, 2010), the Cunningham correction factor C_{cunn}=0.7, and the mean velocity of water molecules ${\stackrel{\mathrm{\u203e}}{c}}_{\mathrm{v}}$. We set the accommodation coefficient α_{m}=0.5 for comparison with former investigations (Kärcher and Lohmann, 2002); this value is also within the range as recommended in recent work by Skrotzki et al. (2013).
For representing the growth rates for the ensemble of ice crystals, by comparison with numerical integration we find that using a shifted mean mass ${m}_{\mathrm{1}}={c}_{\mathrm{1}}\cdot \stackrel{\mathrm{\u203e}}{m}$, c_{1}≈0.819 in the kinetic correction function $f({r}_{\mathrm{1}},a,b)$ is a good approximation.
A4 Howell factor
Latent heat release due to phase changes during diffusional growth changes the surface temperature of the ice crystal. For taking this into account, we use the Howell factor
In the approximation, we neglect the kinetic corrections for diffusion coefficient D_{v} and heat conductivity of air K_{T}.
A5 Capacity of ice crystals
For ice crystals, we assume spherical shape for small crystals and columnar shape for large crystals as in Spichtinger and Gierens (2009); thus the shape factor, or capacity, can be determined exactly using the electrostatic analogy (McDonald, 1963), using a prolate spheroid with semi axes a and b; the capacity can be analytically expressed by
using the eccentricity ${\mathit{\u03f5}}^{\prime}=\sqrt{\mathrm{1}{\left(\frac{b}{a}\right)}^{\mathrm{2}}}$ and the length L of the crystal, which in turn is a function of the crystal mass. Note that the eccentricity changes with crystal growth since the aspect ratio is changing (see Spichtinger and Gierens, 2009, their Eq. 17). We find a very good approximation to the piecewise definition of the capacity by Spichtinger and Gierens (2009) depending on the ice crystal mass
with constants
The representation of the capacity in the ice crystal ensemble is given by the integration, leading to general moments $\mathit{\mu}[m{]}_{{b}_{i}}$.
A6 Ventilation correction
The empirical ventilation corrections usually depend on the use of two dimensionless numbers, i.e. the Schmidt number N_{Sc} and the Reynolds number N_{Re}
using the dynamic viscosity of air μ (e.g. Dixon, 2007). Thus, the size of the ice crystal L is influencing the Reynolds number via the product v_{t}(m)L, using the terminal velocity v_{t} for an ice crystal of mass m. The effect of ventilation, i.e. the additional uptake of water vapour by the airflow around the particle, crucially depends on the shape of the particles. For columnarshaped ice crystals, we adapt the empirical quadratic fit by Liu et al. (2003) to the simulation data (Ji and Wang, 1999) as follows:
For the formulation of the terminal velocity of columnarshaped ice crystals, v_{t}(m), we use the formulation by Spichtinger and Gierens (2009), including also the correction for temperature and pressure, respectively. For representing the ensemble of ice crystals, by comparison with the numerical integration, we find that using a shifted mean mass ${m}_{\mathrm{2}}={c}_{\mathrm{2}}\cdot \stackrel{\mathrm{\u203e}}{m}$ with c_{2}=1.5 in the formulation of the Reynolds number leads to a very good agreement.
In this section we report on the results of the reference simulations, using the corrected formulation of the nucleation rate for supercooled aqueous solution droplets by Koop et al. (2000). For evaluating the quality of the simplified model, we compare the number concentration of ice crystals as obtained from standard nucleation events with results from literature, i.e. with a model using sophisticated particle physics (Kärcher and Lohmann, 2002) and a complex bulk physics scheme (Spichtinger and Gierens, 2009). In Fig. B1 the results are represented for the temperatures T=196, 216, and 236 K at pressure p=200 hPa, as prescribed in Kärcher and Lohmann (2002).
In comparison, we see an overall good agreement of our simple model with the more sophisticated models (Kärcher and Lohmann, 2002; Spichtinger and Gierens, 2009). However, we have to remark here that the deviation in the results for temperature T=236 K at low vertical velocities is the result of the neglect of the ventilation correction in the model by Kärcher and Lohmann (2002). In summary, our simplified approach compares very well with the results of the other studies.
In Fig. B2 a typical nucleation event is shown. Here, two different nucleation parameterizations are used, the reference by Koop et al. (2000) (black line) and the linear fit (red line). There are small differences in the time evolution of the variables saturation ratio S_{i} (left panel), number concentration n_{i} (middle panel), and mean mass m (right panel), but in general there is the same behaviour in both cases.
The source of supersaturation (i.e. cooling by vertical updraft and adiabatic expansion) leads to an increase in S_{i} until nucleation starts at about t_{start}∼40 s, i.e. at very low values of the nucleation rate. S_{i} is still increasing since the sink of depositional growth is not strong enough to reduce water vapour efficiently; thus, the ice crystal number concentration is further increasing due to permanent ice nucleation. At the peak supersaturation, the source and sink of supersaturation are balanced (t_{peak}∼110 s); after this time, S_{i} is decreasing due to the dominant growth term. The number concentration does not change much from this time on, but as long as the values of S_{i} are large enough, ice nucleation still takes place. At about t∼125 s the nucleation event is complete, and no further nucleation takes place, since the nucleation rate is too small. Note that during the time interval [t_{start},t_{peak}] the mean mass m is almost constant (this feature is more prominent in the linear fit case), whereas for t>t_{peak} the mass increases. For t<t_{peak} the nucleation is dominant; thus diffusional growth just compensates for the number increase by mass, whereas afterwards crystal growth is dominant over nucleation. This feature was already seen in former investigations, which leads to a model reduction for analytical investigations (Baumgartner and Spichtinger, 2019). The different nucleation parameterizations agree qualitatively for a nucleation event; however, the nonlinear reference rate leads to some variations. While for the linear fit case, the increase in n_{i} is approximately an exponential growth n_{i}(t)∼exp (αt), and in turn the mean mass is almost constant in the relevant time interval, for the reference case the change deviates slightly from exponential growth.
Note that the thresholds of constant nucleation rates in Fig. B2 (left panel) can be calculated from Eq. (39) using the respective values for j_{0} (i.e. ${j}_{\mathrm{0}}\in \mathit{\{}\mathrm{10},\mathrm{12},\mathrm{14},\mathrm{16}\mathit{\}}$) in the formulation of the supersaturation threshold.
In Koop and Murray (2016) a polynomial of degree 6 is used for fitting the experimental values of the nucleation rate for pure supercooled water. Since polynomials of high degree are difficult to evaluate numerically, we present fits with polynomials of lower degrees, which are still accurate in the relevant temperature range. The original formulation of the nucleation rate is
with a polynomial p_{n}(x) of degree n=6 using the melting temperature of pure water T_{m}=273.15 K. The coefficients c_{i} are reported in Koop and Murray (2016, Table VII), where the nucleation rate is given in units of cubic centimetres per second (${\mathrm{cm}}^{\mathrm{3}}\phantom{\rule{0.125em}{0ex}}{\mathrm{s}}^{\mathrm{1}}$). We reformulate the nucleation rate in SI units (i.e. $\left[J\right]={\mathrm{m}}^{\mathrm{3}}\phantom{\rule{0.125em}{0ex}}{\mathrm{s}}^{\mathrm{1}}$) by a factor of 10^{6} and approximate the logarithmic values log _{10}(J) by polynomials of degree 2 and 4, respectively; i.e.
The coefficients are given in Table C1. For this purpose, we use a leastsquares fit for the temperature range $\mathrm{225}\le T\le \mathrm{245}$ K, for which supercooled water droplets can still exist (see, for example, Fig. 4 in Koop and Murray, 2016). In Fig. C1 (top panel) the approximations are shown in comparison with the original fit, while the ratio $r=\frac{{p}_{i}\left(T\right)}{p\left(x\right)}$ is shown in the bottom panel.
As can be seen, the relative error for the polynomial fit p_{4}(T) is less than 0.25 %, while even for the quadratic fit p_{2}(T), the error is smaller than 2 %. For practical applications in the relevant temperature range $\mathrm{225}\le T\le \mathrm{240}$ K, the quadratic fit might be sufficient. If the original polynomial is used, a sophisticated evaluation of the polynomial is recommended (e.g. Horner scheme).
The perturbation analysis or asymptotic analysis as applied in Sect. 4.6 is a wellknown technique to investigate the impact of perturbations on a mathematical object such as a mathematical expression or the solution of an equation. A good general introduction is given in Holmes (2013), and an application in meteorology is explained in Klein (2010).
The basic idea is to introduce a small parameter ε, to expand the quantity of interest in powers of ε, i.e. ε, ε^{2}, …, and to substitute this expansion into the mathematical object of interest (see Eq. 48 for such an expansion). Since the resulting expression should hold for any value of ε and even in the limit ε→0, all the contributions from the various powers of ε may be considered individually. Given that the parameter ε is assumed as being small, effects that stem from terms with higher powers of ε will only have a small impact, whereas effects with a lower power of ε will be dominant.
For practical applications, it is common to also scale coefficients and parameters of the mathematical expression in powers of ε. This step ensures that the mutual magnitudes of the parameters stay consistent, even in the limit ε→0. This task involves usually some free choices, and is known as distinguished limit.
In the spirit of the works of Klein and Majda, the scaling of the parameters is often done by assuming ε∼0.1; i.e. substituting this value of ε into the scaling of the parameters yields a realistic value of the parameters (see, for example, Hittmeir and Klein, 2018; Baumgartner and Spichtinger, 2019; Klein and Majda, 2006). As an example, the parameter $A\left(T\right)\approx {A}_{\mathrm{0}}={b}_{\mathrm{1}}{a}_{w\mathrm{0}}^{i}\approx \mathrm{149.32}$ is written as ${A}_{\mathrm{0}}={A}^{\ast}{\mathit{\epsilon}}^{\mathrm{2}}$, with ${A}^{\ast}=\mathcal{O}\left(\mathrm{1}\right)$, where the latter may be understood as A^{∗} being independent of ε. With ${A}^{\ast}\approx \mathrm{1.4932}$, the value ε∼0.1 restores the original value of A_{0}.
In essence, the goal is to determine the nonzero parts of the expansion of the mathematical expression, and the respective power of ε indicates how strong this contribution is. As an example, if the final expansion is found to be ${\mathit{\epsilon}}^{\mathrm{1}}{\mathit{\omega}}_{\mathrm{1}}+{\mathit{\omega}}_{\mathrm{0}}+\mathit{\epsilon}{\mathit{\omega}}_{\mathrm{1}}+\mathcal{O}\left({\mathit{\epsilon}}^{\mathrm{2}}\right)$, then the term ω_{−1} will be dominant since this term is associated with the lowest power of ε. If the result would be the expansion ${\mathit{\omega}}_{\mathrm{0}}+\mathit{\epsilon}{\mathit{\omega}}_{\mathrm{1}}+\mathcal{O}\left({\mathit{\epsilon}}^{\mathrm{2}}\right)$, then we may conclude that ω_{0} is the dominant part, and all ω_{k} values for k≥1 only contribute small corrections (since ε is small).
Another fruitful use of perturbation analysis is to allow an equation to determine the possible matching powers of ε, i.e. to answer the following question: which powers of ε are needed to achieve a balance in the equation at hand? As an example, from physical considerations we know that equation Eq. (59) holds. As a consequence, after having expanded both sides of the equation in powers of ε, the expansions on both sides must agree (otherwise the equation would not hold). This is only possible if the powers of ε match; hence the appropriate powers for ε may be inferred.
The model output data and the model configurations are available from the corresponding author upon request.
PM carried out calculations and approximations of the nucleation rates. PS ran the numerical simulations. MB and PS carried out the investigations using asymptotics. All authors were involved in the preparation of the manuscript and have read and approved the final paper.
The contact author has declared that none of the authors has any competing interests.
Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Patrik Marschalik and Peter Spichtinger acknowledge support by the German Bundesministerium für Bildung und Forschung (BMBF) within the HD(CP)^{2} initiative, project S4 (01LK1216A). Manuel Baumgartner acknowledges support by the Deutsche Forschungsgemeinschaft (DFG) within the Transregional Collaborative Research Centre TRR165 Waves to Weather, project Z2. Peter Spichtinger acknowledges support by the DFG within the Transregional Collaborative Research Centre TRR301 TPChange, project B7, and within the research unit Multiscale Dynamics of Gravity Waves (MSGWaves) through grant SP 1163/52. We thank Martina Krämer for fruitful discussions during her stay as a Gutenberg Research College fellow at JGU Mainz and Katrin AmannWinkel for providing us a study about recent measurements of supercooled water. Finally, we thank two anonymous reviewers for their helpful comments and suggestions, which led to a significantly improved manuscript.
This research has been supported by the Bundesministerium für Bildung und Forschung (grant no. 01LK1216A) and the Deutsche Forschungsgemeinschaft (grant no. TRR165, project Z2; grant no. TRR301, project B7; grant no. SP 1163/52).
This paper was edited by Daniel Knopf and reviewed by two anonymous referees.
Baumgartner, M. and Spichtinger, P.: Homogeneous nucleation from an asymptotic point of view, Theor. Comp. Fluid Dyn., 33, 83–106, https://doi.org/10.1007/s00162019004840, 2019. a, b, c, d, e, f, g, h
Baumgartner, M., Weigel, R., Harvey, A. H., Plöger, F., Achatz, U., and Spichtinger, P.: Reappraising the appropriate calculation of a common meteorological quantity: potential temperature, Atmos. Chem. Phys., 20, 15585–15616, https://doi.org/10.5194/acp20155852020, 2020. a, b
Baumgartner, M., Rolf, C., Grooß, J.U., Schneider, J., Schorr, T., Möhler, O., Spichtinger, P., and Krämer, M.: New investigations on homogeneous ice nucleation: the effects of water activity and water saturation formulations, Atmos. Chem. Phys., 22, 65–91, https://doi.org/10.5194/acp22652022, 2022. a, b, c
de Reus, M., Borrmann, S., Bansemer, A., Heymsfield, A. J., Weigel, R., Schiller, C., Mitev, V., Frey, W., Kunkel, D., Kürten, A., Curtius, J., Sitnikov, N. M., Ulanovsky, A., and Ravegnani, F.: Evidence for ice particles in the tropical stratosphere from insitu measurements, Atmos. Chem. Phys., 9, 6775–6792, https://doi.org/10.5194/acp967752009, 2009. a
Dinh, T., Podglajen, A., Hertzog, A., Legras, B., and Plougonven, R.: Effect of gravity wave temperature fluctuations on homogeneous ice nucleation in the tropical tropopause layer, Atmos. Chem. Phys., 16, 35–46, https://doi.org/10.5194/acp16352016, 2016. a
Dixon, J.: The Shock Absorber Handbook, Wiley, ISBN 9780470510209, https://doi.org/10.1002/9780470516430, 2007. a
Fu, Q. and Liou, K.: Parameterization of the Radiative Properties of Cirrus Clouds, J. Atmos. Sci., 50, 2008–2025, https://doi.org/10.1175/15200469(1993)050<2008:POTRPO>2.0.CO;2, 1993. 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, D24514, https://doi.org/10.1029/2007jd008449, 2007. a
Gallo, P., AmannWinkel, K., Angell, C. A., Anisimov, M. A., Caupin, F., Chakravarty, C., Lascaris, E., Loerting, T., Panagiotopoulos, A. Z., Russo, J., Sellberg, J. A., Stanley, H. E., Tanaka, H., Vega, C., Xu, L., and Pettersson, L. G. M.: Water: A Tale of Two Liquids, Chem. Rev., 116, 7463–7500, https://doi.org/10.1021/acs.chemrev.5b00750, 2016. a
Gallo, P., Loerting, T., and Sciortino, F.: Supercooled water: A polymorphic liquid with a cornucopia of behaviors, J. Chem. Phys., 151, 210401, https://doi.org/10.1063/1.5135706, 2019. a
Gasparini, B., Münch, S., Poncet, L., Feldmann, M., and Lohmann, U.: Is increasing ice crystal sedimentation velocity in geoengineering simulations a good proxy for cirrus cloud seeding?, Atmos. Chem. Phys., 17, 4871–4885, https://doi.org/10.5194/acp1748712017, 2017. a
Hall, W. and Pruppacher, H.: Survival of Ice Particles Falling From Cirrus Clouds in Subsaturated Air, J. Atmos. Sci., 33, 1995–2006, https://doi.org/10.1175/15200469(1976)033<1995:TSOIPF>2.0.CO;2, 1976. a, b
Hittmeir, S. and Klein, R.: Asymptotics for moist deep convection I: refined scalings and selfsustaining updrafts, Theor. Comput. Fluid Dyn., 32, 137–164, https://doi.org/10.1007/s001620170443z, 2018. a
Holmes, M. H.: Introduction to Perturbation Methods, vol. 20 of Texts in Applied Mathematics, SpringerVerlag, New York, 2nd edn., https://doi.org/10.1007/9781461454779, 2013. 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, https://doi.org/10.5194/acp1298172012, 2012. a
Jensen, E. and Pfister, L.: Transport and freezedrying in the tropical tropopause layer, J. Geophys. Res., 109, D02207, https://doi.org/10.1029/2003JD004022, 2004. a
Ji, W. and Wang, P.: Ventilation coefficients for falling ice crystals in the atmosphere at lowintermediate Reynolds numbers, J. Atmos. Sci., 56, 829–836, https://doi.org/10.1175/15200469(1999)056<0829:VCFFIC>2.0.CO;2, 1999. a
Joos, H., Spichtinger, P., and Lohmann, U.: Orographic cirrus in a future climate, Atmos. Chem. Phys., 9, 7825–7845, https://doi.org/10.5194/acp978252009, 2009. a
Joos, H., Spichtinger, P., Reutter, P., and Fusina, F.: Influence of heterogeneous freezing on the microphysical and radiative properties of orographic cirrus clouds, Atmos. Chem. Phys., 14, 6835–6852, https://doi.org/10.5194/acp1468352014, 2014. a
Kärcher, B. and Lohmann, U.: A parameterization of cirrus cloud formation: Homogeneous freezing of supercooled aerosols, J. Geophys. Res.Atmos., 107, 4698, https://doi.org/10.1029/2001JD000470, 2002. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o, p
Kärcher, B., Hendricks, J., and Lohmann, U.: Physically based parameterization of cirrus cloud formation for use in global atmospheric models, J. Geophys. Res., 111, D01205, https://doi.org/10.1029/2005JD006219, 2006. a
Klein, R.: ScaleDependent Models for Atmospheric Flows, Annu. Rev. Fluid Mech., 42, 249–274, https://doi.org/10.1146/annurevfluid121108145537, 2010. a
Klein, R. and Majda, A. J.: Systematic multiscale models for deep convection on mesoscales, Theor. Comput. Fluid Dyn., 20, 525–551, https://doi.org/10.1007/s0016200600279, 2006. a
Köhler, C. G. and Seifert, A.: Identifying sensitivities for cirrus modelling using a twomoment twomode bulk microphysics scheme, Tellus B, 67, 24494, https://doi.org/10.3402/tellusb.v67.24494, 2015. a
Köhler, H.: The nucleus in and the growth of hygroscopic droplets, T. Faraday Soc., 32, 1152–1161, 1936. a
Koop, T.: Homogeneous ice nucleation in water and aqueous solutions, Z. Phys. Chem., 218, 1231–1258, 2004. a, b
Koop, T.: Atmospheric water, in: Water: Fundamentals as the Basis for Understanding the Environment and Promoting Technology, edited by: Debenedetti, P. G., Ricci, M. A., and Bruni, F., Proceedings of the International School of Physics “Enrico Fermi”, 187, 45–75, https://doi.org/10.3254/978161499507445, 2015. a, b
Koop, T. and Murray, B. J.: A physically constrained classical description of the homogeneous nucleation of ice in water, J. Chem. Phys., 145, 211915, https://doi.org/10.1063/1.4962355, 2016. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o, p, q, r, s, t
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, e, f, g, h, i, j, k, l, m, n, o, p, q, r, s, t, u, v, w, x, y, z, aa, ab, ac, ad, ae, af, ag, ah
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, https://doi.org/10.5194/acp935052009, 2009. a
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, https://doi.org/10.5194/acp1634632016, 2016. a
Krämer, M., Rolf, C., Spelten, N., Afchine, A., Fahey, D., Jensen, E., Khaykin, S., Kuhn, T., Lawson, P., Lykov, A., Pan, L. L., Riese, M., Rollins, A., Stroh, F., Thornberry, T., Wolf, V., Woods, S., Spichtinger, P., Quaas, J., and Sourdeval, O.: A microphysics guide to cirrus – Part 2: Climatologies of clouds and humidity from observations, Atmos. Chem. Phys., 20, 12569–12608, https://doi.org/10.5194/acp20125692020, 2020. a, b
Laksmono, H., McQueen, T. A., Sellberg, J. A., Loh, N. D., Huang, C., Schlesinger, D., Sierra, R. G., Hampton, C. Y., Nordlund, D., Beye, M., Martin, A. V., Barty, A., Seibert, M. M., Messerschmidt, M., Williams, G. J., Boutet, S., AmannWinkel, K., Loerting, T., Pettersson, L. G. M., Bogan, M. J., and Nilsson, A.: Anomalous Behavior of the Homogeneous Ice Nucleation Rate in “NoMan’s Land”, J. Phys. Chem. Lett., 6, 2826–2832, https://doi.org/10.1021/acs.jpclett.5b01164, 2015. a
Liu, H., Wang, P., and Schlesinger, R.: A numerical study of cirrus clouds. Part I: Model description, J. Atmos. Sci., 60, 1075–1084, https://doi.org/10.1175/15200469(2003)60<1075:ANSOCC>2.0.CO;2, 2003. 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, https://doi.org/10.5194/acp1657932016, 2016. a
McDonald, J.: Use of Electrostatic Analogy in Studies of Ice Crystal Growth, Z. Angew. Math. Phys., 14, 610–620, https://doi.org/10.1007/BF01601268, 1963. a
Mishima, O. and Stanley, H.: The relationship between liquid, supercooled and glassy water, Nature, 396, 329–335, 1998. a
Murphy, D. and Koop, T.: Review of the vapour pressure of ice and supercooled water for atmospheric applications, Q. J. Roy. Meteor. Soc., 131, 1539–1565, 2005. a, b, c, d, e, f, g, h, i, j, k, l, m, n
Nachbar, M., Duft, D., and Leisner, T.: The vapor pressure of liquid and solid water phases at conditions relevant to the atmosphere, J. Chem. Phys., 151, 064504, https://doi.org/10.1063/1.5100364, 2019. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o, p, q, r
Niedermeier, D., Shaw, R. A., Hartmann, S., Wex, H., Clauss, T., Voigtländer, J., and Stratmann, F.: Heterogeneous ice nucleation: exploring the transition from stochastic to singular freezing behavior, Atmos. Chem. Phys., 11, 8767–8775, https://doi.org/10.5194/acp1187672011, 2011. a
Pathak, H., Späh, A., Esmaeildoost, N., Sellberg, J. A., Kim, K. H., Perakis, F., AmannWinkel, K., LaddParada, M., Koliyadu, J., Lane, T. J., Yang, C., Lemke, H. T., Oggenfuss, A. R., Johnson, P. J. M., Deng, Y., Zerdane, S., Mankowsky, R., Beaud, P., and Nilsson, A.: Enhancement and maximum in the isobaric specificheat capacity measurements of deeply supercooled water using ultrafast calorimetry, P. Natl. Acad. Sci. USA, 118, e2018379118, https://doi.org/10.1073/pnas.2018379118, 2021. a, b, c
Pruppacher, H. R.: A new look at homogeneous ice nucleation in supercooled water drops, J. Atmos. Sci., 52, 1924–1933, 1995. a, b
Pruppacher, H. R. and Klett, J. D.: Microphysics of Clouds and Precipitation, vol. 18 of Atmospheric and Oceanographic Sciences Library, Kluwer Academic Publishers, Dordrecht, ISBN 9780792342113, 2010. a, b
Ren, C. and Mackenzie, A. R.: Cirrus parametrization and the role of ice nuclei, Q. J. Roy. Meteor. Soc., 131, 1585–1605, https://doi.org/10.1256/qj.04.126, 2005. a, b, c, d
Schneider, J., Höhler, K., Wagner, R., Saathoff, H., Schnaiter, M., Schorr, T., Steinke, I., Benz, S., Baumgartner, M., Rolf, C., Krämer, M., Leisner, T., and Möhler, O.: High homogeneous freezing onsets of sulfuric acid aerosol at cirrus temperatures, Atmos. Chem. Phys., 21, 14403–14425, https://doi.org/10.5194/acp21144032021, 2021. a, b
Skrotzki, J., Connolly, P., Schnaiter, M., Saathoff, H., Möhler, O., Wagner, R., Niemand, M., Ebert, V., and Leisner, T.: The accommodation coefficient of water molecules on ice – cirrus cloud studies at the AIDA simulation chamber, Atmos. Chem. Phys., 13, 4451–4466, https://doi.org/10.5194/acp1344512013, 2013. a
Spichtinger, P.: Shallow cirrus convection – a source for ice supersaturation, Tellus A, 66, 19937, https://doi.org/10.3402/tellusa.v66.19937, 2014. a
Spichtinger, P. and Gierens, K. M.: Modelling of cirrus clouds – Part 1a: Model description and validation, Atmos. Chem. Phys., 9, 685–706, https://doi.org/10.5194/acp96852009, 2009. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o
Spreitzer, E. J., Marschalik, M. P., and Spichtinger, P.: Subvisible cirrus clouds – a dynamical system approach, Nonlin. Processes Geophys., 24, 307–328, https://doi.org/10.5194/npg243072017, 2017. a, b
Voigt, C., Schumann, U., Minikin, A., Abdelmonem, A., Afchine, A., Borrmann, S., Boettcher, M., Bucuchholz, B., Bugliaro, L., Costa, A., Curtius, J., Dollner, M., Doernbrack, A., Dreiling, V., Ebert, V., Ehrlich, A., Fix, A., Forster, L., Frank, F., Fuetterer, D., Giez, A., Graf, K., Grooss, J.U., Gross, S., Heimerl, K., Heinold, B., Hueneke, T., Jaervinen, E., Jurkat, T., Kaufmann, S., Kenntner, M., Klingebiel, M., Klimach, T., Kohl, R., Kraemer, M., Krisna, T. C., Luebke, A., Mayer, B., Mertes, S., Molleker, S., Petzold, A., Pfeilsticker, K., Port, M., Rapp, M., Reutter, P., Rolf, C., Rose, D., Sauer, D., Schaefer, A., Schlage, R., Schnaiter, M., Schneider, J., Spelten, N., Spichtinger, P., Stock, P., Walser, A., Weigel, R., Weinzierl, B., Wendisch, M., Werner, F., Wernli, H., Wirth, M., Zahn, A., Ziereis, H., and Zoger, M.: MLCIRRUS – The Airborne Experiment on Natural Cirrus and Contrail Cirrus with the HighAltitude LongRange Research Aircraft HALO, B. Am. Meteorol. Soc., 98, 271–288, https://doi.org/10.1175/BAMSD1500213.1, 2017. a
Wernli, H., Boettcher, M., Joos, H., Miltenberger, A. K., and Spichtinger, P.: A trajectorybased classification of ERAInterim ice clouds in the region of the North Atlantic storm track, Geophys. Res. Lett., 43, 6657–6664, https://doi.org/10.1002/2016GL068922, 2016. a, b
Zhang, Y., Macke, A., and Albers, F.: Effect of crystal size spectrum and crystal shape on stratiform cirrus radiative forcing, Atmos. Res., 52, 59–75, https://doi.org/10.1016/S01698095(99)000265, 1999. a
 Abstract
 Introduction
 Empirical fit of the nucleation rate
 Model description
 Investigations of the nucleation rates
 Impact of saturation vapour pressure formulation
 Another approach to formulate the nucleation rate
 Thresholds of ice nucleation
 Summary and outlook
 Appendix A: Model description – details
 Appendix B: Reference simulation results
 Appendix C: Simple fit for nucleation rate of pure water droplets
 Appendix D: Perturbation analysis
 Code and data availability
 Author contributions
 Competing interests
 Disclaimer
 Acknowledgements
 Financial support
 Review statement
 References
 Abstract
 Introduction
 Empirical fit of the nucleation rate
 Model description
 Investigations of the nucleation rates
 Impact of saturation vapour pressure formulation
 Another approach to formulate the nucleation rate
 Thresholds of ice nucleation
 Summary and outlook
 Appendix A: Model description – details
 Appendix B: Reference simulation results
 Appendix C: Simple fit for nucleation rate of pure water droplets
 Appendix D: Perturbation analysis
 Code and data availability
 Author contributions
 Competing interests
 Disclaimer
 Acknowledgements
 Financial support
 Review statement
 References