The Leipzig Cloud Interaction Simulator ( LACIS ) : operating principle and theoretical studies concerning homogeneous and heterogeneous ice nucleation

Introduction Conclusions References


Introduction
Ice containing clouds, such as cirrus and mixed phase clouds have an impact on the Earth's radiative balance by scattering and absorbing solar and terrestrial radiation (Zuberi et al., 2002;Hung et al., 2003).Ice formation in clouds changes cloud radiative properties (DeMott et al., 2003b), affects cloud dynamics, chemical processes, charge separation in cumulonimbus clouds (Takahashi, 1978), and is the source of effective pathways to form precipitation in mixed phase clouds.Therefore, ice formation processes greatly impact cloud lifetime and Earth's climate (Lohmann, 2006).
Ice formation in the atmosphere takes place via homogeneous and heterogeneous nucleation processes.Homogeneous ice nucleation proceeds from a stochastic event in liquid water or aqueous solution droplets without being catalyzed by a foreign substance.In contrast heterogeneous ice nucleation is induced by foreign substances called ice forming nuclei (IN) (Cantrell and Heymsfield, 2005).In general four different heterogeneous freezing modes are distinguished: deposition, condensation, immersion and contact freezing mode (e.g., Pruppacher and Klett, 1997).In the framework of the present paper, we will mainly concentrate on immersion freezing being defined as: An insoluble aerosol particle becomes immersed in a droplet or acts initially as cloud condensation nucleus (CCN).Due to temperature decrease, ice nucleation takes place directly at the IN surface and induces the freezing of the supercooled droplet (e.g., Megahed, 2007).
In general, the understanding of the physical and chemical processes underlying heterogeneous ice formation is limited.Therefore, more scientific work, both theoretical and experimental, is necessary to elucidate fundamental physical and chemical mechanisms, as well as to develop adequate parameterizations that are suitable for use in cloud and global models (K ärcher and Lohmann, 2003;Cantrell and Heymsfield, 2005).
Various field observations of droplet freezing through heterogeneous ice nucleation show that insoluble substances, especially mineral dust particles, act as effective IN Figures in the atmosphere (DeMott et al., 2003a,b;Sassen et al., 2003;Cziczo et al., 2004;Richardson et al., 2007).Mineral dust particles originate from desert and arid regions and can be lifted into the free troposphere during storm events.Subsequently, the dust particles can be transported over large distances (Prospero, 1999;Sassen et al., 2003;DeMott et al., 2003a).As a result mineral dust particles indirectly influence cloud properties, precipitation, and therefore Earth's climate (Zuberi et al., 2002;DeMott et al., 2003a,b).
Considering laboratory studies, there are numerous methods for investigating homogeneous and heterogeneous ice nucleation.Examples are wind tunnel experiments (Pruppacher and Neiburger, 1968;Diehl and Mitra, 1998), the method of electrodynamic droplet levitation (Davis, 1997;Duft and Leisner, 2004), differential scanning calorimetry (Koop et al., 1999;Chang et al., 1999), optical microscopy in a cold stage cell (Koop et al., 1998;Murray et al., 2010) and cloud chamber methods.Concerning the latter, three different types of cloud chambers are classified according to the mechanism used to achieve supersaturation with respect to water and/or ice.Supersaturation with respect to water and ice can be obtained by adiabatically expanding the gas inside the chamber (expansion cloud chamber, e.g.Aerosol Interaction and Dynamics in the Atmosphere (AIDA, M öhler et al., 2003), using the mixing of warm humidified with cold dry air (e.g.Fast Ice Nucleus CHamber FINCH, Bundke et al., 2008), or by combined heat-vapor diffusion (diffusion cloud chamber, e.g. the Leipzig Aerosol Cloud Interaction Simulator (LACIS), Stratmann et al. (2004), the Continuous Flow thermal gradient Diffusion Chamber (CFDC), Rogers, 1988 and Zurich Ice Nucleation Chamber (ZINC), Stetzer et al., 2008).
During the measurement campaign FROST (FReezing Of duST), which took place in April 2008 at the Leibniz Institute of Tropospheric Research (IfT), the ability of mineral dust particles (Arizona Test Dust, ATD) to function as IN was investigated and quantified at the cloud chamber LACIS (Niedermeier et al., 2010).The aerosol particles used were characterized with respect to shape, chemical composition, hygroscopic growth and droplet activation.The main focus of the experiments performed was the Introduction

Conclusions References
Tables Figures

Back Close
Full immersion freezing of water droplets dispersed in air.It should be noted that each droplet contained just one particle/IN, with the particles all possessing approximately the same size.In order to quantify the immersion freezing behavior, fractions of frozen droplets as a function of temperature were determined over a temperature range from 233.15 K to 240.65 K. Based on these measurements, a parameterization of the ice nucleation rate describing the immersion freezing of ATD particles was derived.
In the framework of the present paper, the operating principle of LACIS for studies of immersion freezing is explained on the basis of model simulations for similar conditions as in Niedermeier et al. (2010).Simulations were carried out with the numerical model FLUENT/FPM (Computational Fluid Dynamics (CFD) code FLUENT, Fluent Inc., 2001) combined with the Fine Particle Model (FPM, Particle Dynamics GmbH, Wilck et al., 2002;Whitby et al., 2003).The numerical model accounts for both, homogeneous and heterogeneous ice nucleation.Concerning homogeneous ice nucleation, Classical Nucleation Theory (CNT) was applied, while for heterogeneous ice nucleation either CNT or the parameterization as given in Niedermeier et al. (2010) were utilized.The latter approaches are evaluated via comparison to the experimental results from the measurement campaign FROST.

Leipzig Aerosol Cloud Interaction Simulator
The laminar flow diffusion cloud chamber LACIS was constructed to investigate cloud microphysical processes like hygroscopic growth and droplet activation of aerosol particles under atmospherically relevant conditions (Stratmann et al., 2004).Basically, LACIS is a laminar flow tube of adjustable length consisting of a variable number of 1 m long tube segments (also called sections).For studying homogeneous and heterogeneous ice nucleation, LACIS was extended to its full length of 7 m by adding sections covering the supercooling temperature range T < T 0 = 273.15K, where ice nucleation can occur.The residence times inside this long version of LACIS range from about 2 to 50 s.The temperature can be varied from 298 down to 223 K under operational Introduction

Conclusions References
Tables Figures

Back Close
Full pressures from 700 hPa to ambient values.Inside LACIS supersaturations with respect to water and/or ice are achieved by a combined heat and vapor diffusion process.For determining suitable experimental conditions and quantitative interpretation of experimental data, the numerical model FLUENT/FPM -FLUENT in combination with the Fine Particle Model developed at IfT.
LACIS as a whole consists of a flow conditioning system, the laminar flow tube itself, and optical particle detectors systems.A schematic of the LACIS instrument setup is given in Fig. 1.

Aerosol particle generation and conditioning
The aerosol flow is generated by dispersing various different aerosol particles (e.g.mineral dust, soot, ammonium sulfate particles) in a particle free air flow.Different coating devices are available for modifying the aerosol particles' surfaces by applying coatings of different atmospheric relevant chemical substances (e.g.sulfuric acid, ammonium sulfate, succinic acid, see Stratmann et al., 2004;Niedermeier et al., 2010).In order to enable size-resolved examination of the aerosol particles, a Differential Mobility Analyzer (DMA, type Vienna Medium, Knutson and Whitby, 1975;Reischl, 1991) is utilized, which selects quasi-monodisperse particles according to their electric mobility.
In addition to the aerosol flow a particle free sheath air flow controlled by a mass flow controller (MFC 1179, MKS, Andover, MA) is provided.Both flows are conditioned with respect to temperature and humidity prior entering the flow tube (Fig. 1).Therefore

Laminar flow tube
The sheath air enters the laminar flow tube isokinetically with the aerosol flow, with the latter forming an approximately 2-mm-diameter aerosol beam at the flow tube center.The volume flow rates of aerosol and sheath air are 0.08 l min −1 and 4.00 l min −1 respectively, (corresponding to a mean flow velocity of u = 0.4 m s −1 ) with the first being adjusted by a mass flow controller (MFC 1179, MKS, Andover, MA) at the LACIS outlet.
The newly developed long version of the laminar diffusion cloud chamber LACIS (Fig. 1) consists of seven linked one-meter laminar flow tube sections with an inner diameter of 15 mm.The wall temperatures of the seven tube sections are adjusted separately by thermostats (TH, Sect.1-5: FP50, Sect.6-7: LH85, Julabo, Seelbach, Germany).The temperature control of the tube walls follows the counter flow principle, i.e., the cooling fluids run in the reverse direction compared to the flow inside the tube.In order to control and monitor these wall temperatures, external resistance thermometers (Pt100, B 1/10 pursuant DIN EN 60751, additional calibration at IfT) are used to control the refrigerant temperatures in the supply (the thermometers are connected to the thermostats' control circuits) and measure them in the return line of the cooling cycle of each tube section.With this configuration a wall temperature accuracy of 0.10 K with a stability of ±0.01 K for section 1 to 5 and for the last two tube sections Introduction

Conclusions References
Tables Figures

Back Close
Full an accuracy of 0.30 K with a stability of ±0.10 K is attained.The flow tube is vertically oriented and operated in a top to bottom flow direction.The flow inside the tube is laminar and axisymmetric with a stable, well-defined aerosol particle beam at the center of the flow tube (Stratmann et al., 2004).The operating parameters of LACIS are summarized in Table 1.
Downstream of the laminar flow tube, a Condensation Particle Counter (CPC 3010, TSI Inc., St. Paul, Minnesota, USA), operating at 1.0 l min −1 , is used to measure the aerosol particle number concentration, and a Dew Point Mirror (DPM 973, MBW Calibration Ltd., Wettingen, Switzerland), operating at 0.7 l min −1 , monitors the outlet dew and frost point temperatures.

Particle detection
To detect seed aerosol particles, water droplets, and ice particles, two different Optical Particle Counters (OPC) are installed at the LACIS setup.The upper white-light optical particle spectrometer (OPC 1 in Fig. 1, described in Kiselev et al., 2005), which is situated between tube Sects. 1 and 2, is used to analyze the size distribution of the aerosol particles and/or water droplets after passing the first section.It was not used in the experiments dealt with here, but is mentioned for completeness.At the outlet of LACIS, the white light aerosol spectrometer (WELAS 1000, Palas GmbH, Karlsruhe, Germany, OPC 2 in Fig. 1) was used during the FROST measurement campaign.When measuring with the WELAS instrument the differentiation between frozen droplets and non-frozen or evaporated droplets, that is needed to determine ice fractions, was realized via evaluation of the measured size distributions (Niedermeier et al., 2010).In order to distinguish directly between ice particles and water droplets/seed particles having identical sizes (via phase state and therefore surface structure), two further optical devices applying different techniques can be employed in future applications: (i) the Thermostabilized Optical Particle Spectrometer (TOPS-ICE, developed at IfT), which allows to distinguish between spherical (droplets) and non-spherical (seeds Introduction

Conclusions References
Tables Figures

Back Close
Full and ice particles) particles by detecting the polarization state of the scattered light and (ii) the LISA instrument (Lacis Ice Scattering Apparatus, University of Hertfordshire, Hatfield AL10 9AB, United Kingdom, Hirst et al., 2001), with which two-dimensional diffraction scattering patterns of the investigated particles are recorded.

Applications and modes of operation
For investigating hygroscopic growth, activation, and heterogeneous ice nucleation behavior of size-segregated well-defined aerosol particles, LACIS can be operated in different modes which depend on the actual boundary conditions.These are the water sub-and supersaturated modes for the temperature range above T 0 = 273.15K, the water sub-, ice supersaturated, and the water super-and ice supersaturated modes below T 0 = 273.15K.If the inlet dew point temperature is lower than the wall temperature, the inner tube walls are in equilibrium with the water vapor of the flow and the dew point temperature remains constant (water sub-saturated mode).Applying this water sub-saturated mode for T > 273.15K with inlet and wall temperature being identical, deliquescence and hygroscopic growth of aerosol particles can be investigated (Wex et al., 2006(Wex et al., , 2007;;Voigtl änder et al., 2008;Ziese et al., 2008;Niedermeier et al., 2008).With higher inlet dew point temperatures compared to the wall temperatures, supersaturation can be achieved as a result of the simultaneous heat and vapor diffusion, which occur at slightly different rates.In this mode, the water vapor condenses on both the aerosol particles and the inner tube walls.As a result a thin water film accumulates on the tube walls and the dew point temperature of the flow converges to the wall-temperature set point.With this method, critical supersaturations for cloud droplet activation and growth, including kinetic effects, can be studied (Wex et al., 2006(Wex et al., , 2007;;Voigtl änder et al., 2008;Niedermeier et al., 2008;Ziese et al., 2008).Homogeneous and heterogeneous ice nucleation processes can be investigated with LACIS.In particular for heterogeneous ice nucleation three different freezing modes can be studied: immersion freezing, deposition nucleation and contact freezing insideout (surface crystallization).In order to analyze deposition nucleation, LACIS can be Introduction

Conclusions References
Tables Figures

Back Close
Full operated in the water subsaturated and ice supersaturated mode at temperatures below 273.15 K.The investigation of contact freezing inside-out (Shaw et al., 2005;Durant and Shaw, 2005), i.e. evaporating the initially formed water droplets followed by freezing due to the contact of the IN with the water-air interface (water sub-and ice supersaturated conditions) is also feasible.In this study, we mainly concentrate on water super-and ice supersaturated conditions to study homogeneous and, especially, immersion freezing processes.A simple and straight forward mode of operation was used for LACIS, i.e., following Niedermeier et al. (2010), Sects.6 and 7 were applied to activate the seed particles to water droplets and subsequently freeze them by further temperature decrease.Other modes are possible but not dealt with here.To achieve reproducible and well-defined measurement conditions, it is ensured that a thin ice layer covers the inner tube walls.This is realized by accumulating liquid water on the tube walls which is then converted to ice by cooling the walls down below 233 K.
3 Numerical model and nucleation rate coefficients

Numerical model
The numerical model is able to simulate the hygroscopic growth of the seed particles, their activation to cloud droplets and subsequent condensational growth, as well as homogeneous and heterogeneous ice nucleation processes under the prevailing thermodynamic conditions inside LACIS.The main fluid dynamical processes to be considered are fluid flow, and heat and mass transfer.Concerning particle dynamics, transport and phase transitions processes need to be accounted for.These processes are mathematically described by the momentum, the vapor mass transport, the energy equation and conservation equations for e.g., particle number and mass.The particle dynamical equations account for transport due to convection, diffusion and external forces (thermophoresis, sedimentation), as well as phase transition processes such as condensation/evaporation and homogeneous and heterogeneous ice nucleation.The Introduction

Conclusions References
Tables Figures

Back Close
Full respective equations are coupled through latent heat release and vapor depletion resulting from the phase transition processes.
To solve the coupled fluid and particle dynamics equations for initial and boundary conditions encountered in LACIS, the Computational Fluid Dynamics code FLUENT (Fluent Inc., 2001) combined with the Fine Particle Model (FPM, particle dynamics GmbH, Wilck et al., 2002;Whitby et al., 2003) are applied.For minimizing computing time, the simulations are realized on a two-dimensional computational grid taking advantage of the system's rotational symmetry.
The fluid momentum equation, assuming steady state conditions, is given by with the density of gas mixture ρ g , velocity vector v , pressure p, dynamic viscosity of the vapor-gas mixture µ, the term V comprising viscosity terms not accounted for in ∇ • (µ∇v ), and the vector of gravitational acceleration g.The vapor mass transport equation has the following form (Bird et al., 1960): where ξ v is the vapor mass fraction, j v represents the mass flux of vapor relative to the mass average velocity, S v specifies the vapor sink due to condensation on particles, droplets or ice particles, D v is the vapor diffusion coefficient in air and α v,g is the thermal diffusion factor of the vapor-gas mixture.The mass flux of vapor j v is governed by two mechanisms, molecular (first term of Eq. 3) and thermal diffusion (second term of Eq. 3).The energy equation for an air-vapor mixture includes heat transport due to convection, conduction and vapor transport.This equation is expressed as (5) Introduction

Conclusions References
Tables Figures

Back Close
Full with the specific enthalpy h, the heat flux q, heat source S h , thermal diffusivity α = k h /ρ g c p with heat conductivity k h and specific heat capacity at constant pressure c p .h v and h g are the specific enthalpies of vapor and gas, M, M v and M g represent the molar weights of the mixture, the vapor and the dry carrier gas, respectively.The Modal Aerosol Dynamics method (MAD) (Whitby and McMurry, 1997) is applied to parameterize the particle size distribution.Therefore, the size distribution is split into modes (represented by j ) (Whitby and McMurry, 1997;Whitby et al., 2002Whitby et al., , 2003)), each representing a distinct particle population.Here the following particle populations are distinguished: seed particle-droplet mode (j = 1), the homogeneously nucleated ice mode (j = 2) and the heterogeneously nucleated ice mode (j = 3).Basically two moments M k i ,j (number and mass) are used for representing each mode, i.e., the total particle number concentration N j (k = 0) and the mass concentrations M i ,j (k = 1), with each chemical component i in the particle being represented by its own mass moment.This corresponds to each mode being internally mixed, monodisperse and moving in size space.From the considered moments, total particle mass and size can be determined for each mode.The moment dynamic equations for the number N j and mass M i ,j concentration for the different modes are given in Table 2 with external particle velocity v , gas velocity v g , particle diffusion coefficient D j and the single particle mass m i ,j of substance i in mode j .The particle diffusion coefficient is computed via , where d p presents the particle diameter assuming spherical shape and C C is the Cunningham correction factor, which can be calculated using where λ specifies the mean free path of the gas molecules.For the description of the dynamic growth of water droplets and ice particles, the single particle growth law according to Barrett and Clement (1988) is used for the droplets and ice particles Introduction

Conclusions References
Tables Figures

Back Close
Full (6) S i and S i ,j describe the saturation ratios in the gas phase and over the particle surface, whereby Kelvin and Raoult effects are accounted for.The equilibrium vapor pressure is given by p i , f mass and f heat are the mass and heat transfer transition functions, k g is the carrier gas heat conductivity and L i represents the latent heat of evaporation and sublimation, respectively.Further quantities in Table 2 are M 2,1 , the mass concentration of liquid water in the seed-particle droplet mode, ρ g , ρ w the gas-mixture and liquid water densities, S IN the total surface area of the seed particles' insoluble core (S IN = N 0 S p ). N 0 is the total particle/droplet number and s p is the surface area of a single particle), and j hom (T ), j het (T ) represent the homogeneous and heterogeneous nucleation rate coefficient with the units m −3 w s and m −2 s, respectively.The newly-developed phase transition model, which transfers particles from the seed particle-droplet mode to either homogeneous or heterogeneous ice mode, is implemented in the moment dynamics equations via the respective sink/source terms S k hom,i and S k het,i .These sink/source terms can also be interpreted as freezing rates.S k hom,i is proportional to the total liquid water volume of the considered droplet population and the temperature dependent homogeneous ice nucleation rate coefficient.S k het,i is a function of total IN surface area (only insoluble core) and the temperature dependent heterogeneous ice nucleation rate coefficient.Both quantities have the same units: number of nucleation events per unit time and gas volume.
The different modes and the particle fluxes between the modes are illustrated in Fig. 2. In the calculations it is assumed, that the seed particles consist of an insoluble core (e.g.ATD) and a soluble coating (e.g.(NH 4 ) 2 SO 4 ).These particles can be either dry or hygroscopically grown or activated droplets.For the latter two cases liquid water is also a component of the seed particle-droplet mode (left solid line in Fig. 2).The homogeneous and heterogeneous ice modes have the same composition, i.e., they Introduction

Conclusions References
Tables Figures

Back Close
Full are made up of the insoluble core, the coating and ice.The material properties for the ice phase used for the computations are specified in the Appendix.Through homogeneous ice nucleation, described by the homogeneous freezing rate Sk hom,i , particles from the seed particle-droplet mode are transferred to the homogeneous ice mode (j = 2).Likewise, the heterogeneous freezing rate Sk het,i determines the transfer to the heterogeneous ice mode (j = 3).The concept outlined above does not depend on any specific homogeneous and/or heterogeneous nucleation rate coefficient, so different coefficients, e.g.those discussed below, can be implemented and tested.

Nucleation rate coefficients
For determining the homogeneous and heterogeneous nucleation rate coefficients to be used in FLUENT/FPM, two different theoretical approaches are adopted: (a) Classical Nucleation Theory is applied for both homogeneous and heterogeneous ice nucleation, and (b) CNT is used for modeling homogeneous nucleation, but immersion freezing is described by implementing a parameterization derived from prior LACIS measurements (Niedermeier et al., 2010).
According to CNT (Pruppacher and Klett, 1997;Seinfeld and Pandis, 1998;Zobrist et al., 2007) the homogeneous nucleation rate coefficient j hom , which is the number of nucleation events per time interval and total liquid water volume of the droplet population, is defined as Gibbs free energy ∆G hom is expressed by where σ w/i represents the interfacial free energy of the water-ice boundary, n i is the number density of molecules in the solid phase and S w/i is the ratio of the saturation vapor pressures over water and ice.The first term of Eq. ( 7) represents the flux of water molecules to the ice phase and the second term describes the equilibrium number of critical embryos in the liquid phase.Values for quantities such as ∆F diff , S w/i and σ w/i are chosen according to Zobrist et al. ( 2007) and references therein.Altogether the homogeneous nucleation rate coefficient is a very steep function of absolute temperature.
It is known that homogeneous ice nucleation rate coefficients from CNT tend to be too small at temperatures lower than 235 K, but as otherwise good agreement with experiments was found (Jeffery and Austin, 1997;Pruppacher and Klett, 1997;Cantrell and Heymsfield, 2005), they can be used in the temperature range considered here.
Heterogeneous ice nucleation is implemented in the CNT by assuming the insoluble foreign substance, i.e. the IN, increases the likelihood to form a critical embryo, but does not disturb the stochastic nature of the freezing process.So heterogeneous ice nucleation can be derived from the homogeneous case by additionally accounting for the energy barrier reduction due to presence of the IN.Usually this is done based on the concept of contact angle (assuming a spherical cap for ice embryo shape) (Mason, 1971;Pruppacher and Klett, 1997;Seinfeld and Pandis, 1998).The reduced critical Gibbs free energy is then given by where f (θ) = 1 4 (2+cosθ)(1−cosθ) 2 is the reduction factor and θ represents the contact angle.θ may vary between 0 and 180 • , where the former case implies that the energy barrier is zero (nucleation occurs as soon as supersaturation is reached) and the latter Introduction

Conclusions References
Tables Figures

Back Close
Full corresponds to homogeneous ice nucleation.The heterogeneous ice nucleation rate coefficient, which is defined as number of nucleation events per time interval and total IN surface, is given by with the number density of liquid molecules in contact with IN surface N s , which has a typical value of 1 × 10 19 m −2 .
As a second approach for heterogeneous ice nucleation, a parameterization derived from LACIS measurements of immersion freezing is implemented into the numerical model.The measurement procedure and the derivation of this parameterization are explained in Niedermeier et al. (2010).The heterogeneous ice nucleation rate coefficient j het,LACIS (T s ) is with a and f het being temperature independent fitting parameters derived from the measured data.The fitting parameter a includes information about the IN surface area of a single particle s p and kinetic effects.

Results and discussion
The numerical model FLUENT/FPM as described above is a suitable tool for exploring LACIS' behavior for a given set of boundary conditions, testing assumptions made during the interpretation of experimental data, and evaluating the feasibility of different theoretical approaches for modeling homogeneous and heterogeneous ice nucleation.
In the framework of this section, (a) the principle behavior of LACIS when operated in immersion freezing mode will be described, and (b) experimental data and their interpretation as given in Niedermeier et al. ( 2010) will be reviewed.In this context, boundary conditions corresponding to those used during the FROST campaign were applied in the numerical simulations.The inlet temperature and the inlet dew point were set to 273.15 K and 265.95 K.The wall temperatures of Sects.6 and 7 were always identical and varied in a range of 233.15 ≤ T w 6,7 ≤ 240.65 K. Furthermore, the wall boundary condition for sections 6 and 7 was always set to ice saturation (S i = 1), which corresponds to ice covered inner tube walls.For both the homogeneous and the heterogeneous studies the seed particle/IN concentrations were set to 300 cm −3 .
As each droplet contained a single seed particle/IN this results in a droplet number concentration of 300 cm −3 as well.When studying homogeneous and heterogeneous ice nucleation, IN were assumed to be spherical with diameters of 187 or 300 nm, internally mixed consisting of an insoluble ATD core and a small amount (mass fraction of 0.019) of ammonium sulfate.The latter was done to reproduce the activation behavior observed in CCN measurements during the FROST campaign.Concerning ATD, the following material properties were assumed: a molecular weight of 65.18 g mol −1 , which is the mass weighted average of the main constituent of ATD and a density of 2600 kg m −3 (M öhler et al., 2006).

Behavior of LACIS operated in immersion mode
In order to study ice nucleation processes in LACIS the thermodynamic conditions such as temperature and saturation with respect to water and ice inside the flow tube have to be known.Due to the coupled heat-vapor diffusion processes taking place, the Introduction

Conclusions References
Tables Figures

Back Close
Full temperature, vapor concentration and saturation profiles are complex and inhomogeneous.As particles/droplets, to good approximation, are confined to the center of the flow tube, the temperature and saturation profiles at the flow tube centerline are of special interest.Typical profiles of temperature and saturation with respect to water and ice are given as a function of time for different wall temperature settings in Fig. 3. Just the last two sections are shown, i.e., the profiles in Sections.6 and 7.The temperature profile (Fig. 3, panel 1) exhibits a steep fall within the first freezing section (t ≤ 1.6 s) and approaches the externally-set wall temperatures in the second freezing section (1.6 ≤ t ≤ 3.2 s).It should be noted that on average after 2.6 s already, the temperature at the centerline T axis reaches its set value within a range of about +0.3 K.This range is defined as temperature error in LACIS.
Within the supercooling temperature range (T < 273.15 K) the vapor pressure over ice is smaller than over supercooled water.Fig. 3, panel 2 shows profiles for the saturation with respect to water (solid lines) and ice (dashed lines).For both water and ice the saturation ratio rises strongly until a maximum is reached, then decreases and approaches a constant value.The actual profiles depend on wall temperature and inlet dew point.For constant inlet dew point, the lower the wall temperature, the higher are the saturation maxima.At the outlet of LACIS the ice saturation approaches 1, while water saturation converges to a subsaturated level.
The behavior of supercooled water droplets and ice particles under these thermodynamic conditions is explained in Fig. 4 the frozen droplets.A mass balance yields that approximately 94% of the water vapor is transported to the ice covered walls.6% of the inlet water vapor flows out of LACIS and only the remaining small amount on the order of 10 −5 participates in microphysical processes described in the following.Additionally, for reference the axial temperature T axis is given by the gray curve and the right ordinate.
At first the seed particles grow hygroscopically.As soon as the water saturation increases above the critical super-saturation, the seed particles become activated to liquid droplets, which is seen in Fig. 4 as an increase in the liquid water mass fraction.Subsequently the droplets grow dynamically by vapor diffusion (continuous rise of liquid water mass fraction in Fig. 4) and reach their maximum droplet sizes, which are approx.
2.1 µm for T w 6,7 = 240.65K, 3.0 µm for T w 6,7 = 239.15K, 3.5 µm for T w 6,7 = 238.15K, 4.7 µm for T w 6,7 = 235.65K and 5.6 µm for T w 6,7 = 233.15K wall temperature settings.At about T axis ≈ 245 K ice nucleation starts taking place and the first droplets freeze.Due to further cooling and passage of time, more droplets freeze.The mass fraction of ice particles increases continuously due to more and more droplets freezing and the depositional growth of the already frozen droplets.Fig. 4 clearly shows that the droplets are first formed at higher temperatures (T axis ≈ 260 K) and require further reduction of temperature to freeze.This indicates that for these conditions, immersion freezing is the only freezing process taking place.Once water saturation falls below 1 (T axis ≈ 240 K, t ≈ 1.6 s) the remaining unfrozen droplets start to evaporate in the ice super-and water subsaturated environment due to the Wegener-Bergeron-Findeisen effect, and evaporate completely before passing the outlet of LACIS.As a result seed particles and different sized ice particles leave the outlet of LACIS.Experimental and theoretical sensitivity studies varying the inlet dew point resulted in dew point temperature ranges to be considered in experimental investigations in which the determined ice fractions are not affected by evaporation of the droplets at too early a stage.
In the immersion freezing parameterization (Eq.11) derived in Niedermeier et al. (2010) is applied.Panel 1 illustrates the homogeneous j hom and heterogeneous j het ice nucleation rate coefficients.First of all j het rises steeply and approaches a nearly constant value (t > 1.5 s corresponds to the beginning of Sect.7 in which the absolute temperature T is nearly constant).j hom increases later and at lower temperatures compared to j het .It is only non-negligible for the wall temperature settings of T w 6,7 = 335.65K and T w 6,7 = 333.15K.But it should be noted that homogeneous and heterogeneous ice nucleation rate coefficients, because of their different dimensions ([j hom ] = m −3 w s and s), can only be compared qualitatively.For quantitative comparison, freezing rates for homogeneous S hom and heterogeneous ice nucleation S het (illustrated in panel 2) are more appropriate, as these quantities feature the same dimensions.S het has a similar shape compared to j het .S hom increases at lower temperatures than S het and is smaller than S het except at T w 6,7 = 333.15K.This implies that first, ice formation takes place via heterogeneous ice nucleation, and only if S hom becomes effective (for T w 6,7 ≤ 335.65 K) ice is also formed due to homogeneous ice nucleation.For T w 6,7 = 333.15K, both the homogeneous and the heterogeneous source terms drop to zero at t = 2.3 s, as all supercooled water droplets have been frozen.This can also be seen in panel 3 depicting the ice fraction f i , which is the ratio of ice particle number N i per total particle number N 0 .At sufficiently low temperatures (T ≈ 243.5 K) immersion freezing takes place and the ice fraction exceeds the experimental detection limit of 10 −4 .The ice fraction f i increases monotonically with decreasing temperature, and for T w 6,7 = 333.15K, f i rises steeply when homogeneous ice nucleation sets in and becomes dominant.In general the lower the wall temperature the higher the ice fraction.Summing up the results of the numerical simulations discussed in this section, it can be stated that over the experimental parameter space (wall temperatures, dew points, residence times, etc.) investigated, ice nucleation is clearly dominated by immersion freezing for all cases apart from T w 6,7 = 333.15K, where homogeneous ice nucleation becomes dominant at about 234.9 K. Homogeneous freezing is negligible Introduction

Conclusions References
Tables Figures

Back Close
Full for T w 6,7 > 235.65 K.The ice nucleation time varies in a range of 1.7 to 2.1 s for the respective wall temperature settings and the vast bulk of ice is formed in Sect.7, where the temperature is almost constant.

Comparison of experiments and model simulations
The numerical model can also be used to interpret experimental data collected with LACIS: e.g.understanding the relative importance of processes involved, evaluating the feasibility of different theoretical approaches, and checking the validity of assumptions made.Specifically, here the relative importance of homogeneous vs. heterogeneous ice nucleation during immersion freezing experiments at LACIS, the applicability of CNT to describe the immersion freezing behavior of ATD particles, and the feasibility of assuming a constant temperature when deriving parameterizations for ice nucleation rate coefficients from LACIS-measured ice fractions, will be explored.
Figure 6 shows a comparison of ice fractions as measured with LACIS and calculated with FLUENT/FPM, as a function of temperature.The experimental data are represented by the orange (ice fractions as measured for 300 nm ATD particles) and the black (ice fractions determined for highly diluted ammonium sulfate solution droplets) squares.Concerning the model simulations, results applying homogeneous (solid line, Eq. 7) and heterogeneous (broken lines, Eq. 10) CNT are given.The temperature on the x-axis corresponds to the wall temperatures of Sects.6 and 7 (T w 6,7 ).
Looking at the experimental data in Fig. 6, it becomes obvious that, when considering ATD particles as IN, around T = 236 K the slope of the f i vs. temperature curve becomes steeper.A similar behavior can be found considering the ice fractions measured for the highly diluted solution droplets.Now comparing theoretical and experimental results, it can be seen that CNT-based homogeneous nucleation theory (solid line) is able to explain the behavior of the ammonium sulfate solution droplets (both slope and absolute values are predicted with reasonable accuracy).The small differences observed, corresponding to a ∆T ≈ 0.5 K, may be a result of the numerical model slightly overpredicting droplet volume and Introduction

Conclusions References
Tables Figures

Back Close
Full consequently ice nucleation rate.However, CNT based heterogeneous nucleation theory, assuming constant contact angles, fails to predict the observed freezing behavior.Neither the ice fractions, nor the slope of the ice fraction vs. temperature curve match (the predicted slope is much steeper), regardless of the actual contact angle.Furthermore, it can be seen that decreasing the total particle surface area by a factor of 2.6 (dotted curve compared to dashed-dotted-dotted curve), the ice fraction decreases slightly, but the shape of the curve remains almost constant.
Interpreting these results the following statements can made for the experimental parameter space considered here: Classical homogeneous ice nucleation theory, utilizing the properties as given in Zobrist et al. (2007), is able to predict the freezing behavior of highly diluted ammonium sulfate solution droplets.Classical heterogeneous ice nucleation theory together with the assumption of constant contact angle fails to predict the experimental observations made of the immersion freezing behavior of ATD particles.
Investigating immersion freezing of water droplets coated with a nonadecanol mono- In Niedermeier et al. (2010) an alternative approach to Zobrist et al. ( 2007) and Marcolli et al. (2007) for parameterizing experimental results based on CNT was suggested (Eq.11).Validity of the assumptions made in this context will be discussed in the following.
Similar to Figs. 6 and 7 depicts ice fractions as function of temperature with the experimental data (orange and black squares) being identical.The orange solid line represents results from the parameterization as given in Niedermeier et al. (2010), i.e.Eq. ( 11) is used for calculating the ice fractions.The other solid lines correspond to results from FLUENT/FPM with different assumptions for calculating the freezing rate.Shown are results determined assuming just homogeneous (blue line), just heterogeneous (green line), and both homogeneous and heterogeneous (red line) ice nucleation being active.The latter curve compares well with the experimental data, indicating that the FLUENT/FPM is a suitable tool for describing the complex fluid/particle dynamical and phase transition processes taking place in LACIS.The results from the model simulations taking either homogeneous (blue line) or heterogeneous (green line, underlying the red curve for T > 235.65 K) ice nucleation into account show clearly that immersion freezing is dominant for T > 234.9 K, while at lower temperatures homogeneous freezing is the main process.It should be noted that in the parameter space investigated here, heterogeneous ice nucleation is not quenching homogeneous ice nucleation.In other words, even for a reasonably effective IN such as pure ATD, there is a temperature range in which either both mechanisms can be active or even homogeneous nucleation can be dominant, although every supercooled droplet has a particle immersed.The main reason is the apparently different temperature dependence of the two freezing mechanisms.The actual explanation of why the two mechanisms feature different temperature dependencies remains the topic of future investigations.Now concentrating on the temperature range in which immersion freezing is dominant and comparing the orange to the green line (underlying the red one in the range of interest) it can be seen that the green line, corresponding to the FLUENT/FPM results, is slightly lower, however still within the level of uncertainty of the experimental data, Introduction

Conclusions References
Tables Figures

Back Close
Full than the orange one representing data generated with the parameterization given in Niedermeier et al. (2010).This parameterization was determined assuming the temperature to be constant during the ice nucleation process and equal to the wall temperature of the last section.The time period where the temperature at the centerline is almost constant (deviation +0.3 K) was taken as ice nucleation time.In contrast, the FLUENT/FPM simulations, even though being based on the same nucleation rate coefficient, account for the detailed temperature variation along the flow tube center line (Fig. 7).The small difference between the two data sets is indicative that the assumptions made in Niedermeier et al. (2010) concerning both nucleation temperature and time are justified.Consequently, the method for determining the fitting coefficients in the parameterization equation, and the parameterization concept itself are justified and valid.

Summary and conclusions
Homogeneous and heterogeneous ice nucleation, in particular the immersion freezing of Arizona Test Dust (ATD) particles, have been studied both theoretically and experimentally with the Leipzig Cloud Interaction Simulator (LACIS).In the framework of the present paper, both the physical LACIS setup as used e.g. during the two FROST measurement campaigns (see also Niedermeier et al., 2010;Sullivan et al., 2010;Reitz et al., 2011;Wex et al., 2011) tal parameter space (temperatures, dew points, seed particle number concentrations, etc.) in which such experiments can be performed.
For the evaluation of different theoretical approaches for the description of homogeneous and heterogeneous ice nucleation, model simulations were carried out applying Classical Nucleation Theory (CNT) for homogeneous and heterogeneous ice nucleation.Concerning the latter, a constant contact angle was assumed.It was found that for the experimental parameter space considered, classical homogeneous nucleation theory, utilizing the properties as given in Zobrist et al. (2007), is able to predict the freezing behavior of highly diluted ammonium sulfate solution droplets.However, classical heterogeneous nucleation theory, together with the assumption of a constant contact angle, fails to predict the experimental observations made concerning the immersion freezing behavior of ATD particles.The main reason for this failure is the apparently too strong temperature dependence of the nucleation rate coefficient.
Assuming CNT for homogeneous ice nucleation and the immersion freezing parameterization according to Niedermeier et al. (2010) for heterogeneous ice nucleation, it was found that the simulated freezing behavior was in good agreement with the measured values.Also, it was found that in the experimental parameter space investigated, heterogeneous ice nucleation is not necessarily quenching homogeneous ice nucleation.Even for a reasonably effective IN such as pure ATD there are temperature regimes where homogeneous ice nucleation is dominant.The main reason are the different temperature dependencies of the two freezing mechanisms.
Finally, reviewing the assumptions concerning constant temperature and ice nucleation time made in Niedermeier et al. (2010) when deriving a parameterization for the nucleation rate coefficient in the immersion freezing mode, the good agreement between parameterization and simulation results shows that both assumptions were highly justified.This underlines the applicability of both the method to determine the fitting coefficients in the parameterization equation, and the validity of the parameterization concept itself.Introduction

Conclusions References
Tables Figures

Back Close
Full To describe the material ice in the numerical model following temperature dependent properties are required: density, heat capacity at constant pressure, interfacial free energy between ice and vapor phase, latent heat of sublimation, thermal conductivity and vapor pressure.For the ice density ρ i the parameterization in Pruppacher and Klett (1997) is applied: with T in K and ρ i in kg m −3 .Heat capacity values of ice are given by Giauque and Stout (1936) in a temperature range of 15 < T < 273 K in the unit J mol −1 K −1 : The experimentally determined interfacial free energy according to Ketcham and Hobbs (1969) is used for the simulations.The latent heat of sublimation L s is derived from experiments of Murphy and Koop (2005): with L s in J mol −1 .The vapor pressure over ice is approximated by the equation of Goff and Gratch 1 .Introduction

Conclusions References
Tables Figures

Conclusions References
Tables Figures
Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | both flows are humidified to defined dew point temperatures by two separate saturators (aerosol: MH-110-12S-4, sheath air: PH-30T-24KS, Perma Pur, Toms River, New Jersey).The saturators consist of semipermeable Nafion tubes surrounded by water jackets, which are temperature controlled by the respective thermostats (TH, aerosol air: F25, sheath air: FP50, Julabo, Seelbach, Germany).The saturator of the aerosol air flow can also be bypassed, so that the aerosol flow remains dry (T d ≈ 233 K).The inlet dew point temperatures of the sheath air flow can be varied in the range between Discussion Paper | Discussion Paper | Discussion Paper | 233 and 293 K by mixing humidified and dry air flows.An inline-connected chilled mirror dew point hygrometer (DPM, DewMaster, EdgeTech, Marlborough, MA) monitors the dew point temperature of the sheath air flow.The aerosol and sheath air flows are combined in the inlet section of LACIS.The inlet serves as heat exchanger (temperature controlled by a thermostat with an accuracy of 0.01 K, TH, F25, Julabo, Seelbach, Germany) for harmonizing the temperatures of both air flows and combining them at the entrance of the laminar flow tube.At this point, initial condition such as inlet temperature T IN , inlet dew point temperature T d,IN and flow velocities u, are well defined and known for consideration in experimental data interpretation and for use in the numerical model.
Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | ) with the Boltzmann constant k B , absolute temperature T , Planck constant h, the diffusion energy across the water-ice interface ∆F diff , and the number density of molecules in the bulk water N v (typical value 3.1 × 10 28 m −3 , Zobrist et al., 2007).The critical Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | f het accounts for IN surface properties and thermodynamic effects.Here T S = T 0 − T axis is the supercooling temperature and C 1 = 5.00×10 5 K 2 and C 2 = 8.24×10 1 K are constants resulting from a simplification of the surface free energy and Gibbs free energy terms.For 300 nm uncoated ATD particles and a ice nucleation time of 1.56 s the fitting parameters have the following values: a = 1.31 and f het = 4.51 × 10 −2 .This formula is valid over a supercooling range of 235 < T ≤ 239 KDiscussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | , whereas model version b in Sect.3.2 is applied for forming ice particles.It depicts the mass fractions of water vapor, liquid water and ice as function of time at the centerline of LACIS for a wall temperature of T w 6,7 = 238.15K. Here, the mass fractions should be interpreted as total mass of species i (water vapor, liquid water and ice) per total mass of gas.The red line in Fig. 4 represents the water vapor mass fraction, which is defined by the inlet dew point temperature.The mass fraction of water vapor decreases along the centerline of LACIS.Sinks are the flux of water vapor to the inner ice covered tube walls (boundary condition S i = 1), condensation of water vapor on the droplets and deposition of water vapor on Discussion Paper | Discussion Paper | Discussion Paper | Fig. 5 the time development of different parameters quantifying homogeneous and heterogeneous ice nucleation are shown.Here homogeneous ice nucleation is calculated according to CNT (Eq.7).For simulating the heterogeneous ice nucleation, Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | layer and deducing the heterogeneous ice nucleation rate coefficients from this measurements,Zobrist et al. (2007) found similar results.The CNT model using constant contact angle cannot reproduce those measurements because the slope the j het curve is too steep.As a resultZobrist et al. (2007) assume a linear temperature dependence of the contact angle in order to get their experimental results and theory into agreement.Marcolli et al. (2007) analyzed the immersion freezing behavior of droplets containing several distinct sized ATD particles with differential scanning calorimeter technique.They also stated, that a regular stochastic model (CNT) cannot explain their experimental results.Simulations assuming a singular model with contact angle distribution, where the contact angles vary between the particles considered, or accounting for a distribution of active sites led also to better agreement.Similar conclusions concerning the insufficiency of a constant contact angle model describing experiments accurately were drawn in (L ü önd et al., 2010), where the immersion freezing behavior of sizeselected kaolinite particles was studied.Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | , and the numerical model developed to design experiments and interpret their results have been presented in detail.The model developed and used for the theoretical investigations is based on the commercially available computational fluid dynamics code FLUENT and the Fine Particle Model (FPM).Both together allow for the coupled solution of the conservation equations for momentum, mass, heat and seed particle, droplet, ice particle number and mass, needed to describe the complex fluid/particle dynamical and phase transition processes taking place inside LACIS.Using this model, the operating principle of LACIS has been presented, (a) outlining its ability to perform both homogeneous and heterogeneous freezing (specifically immersion freezing) experiments, and (b) defining the experimen-Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | ice nucleation active sites in mineral dust particles caused by sulphuric acid condensation, Discussion Paper | Discussion Paper | Discussion Paper |

Fig. 2 .Fig. 3 .Fig. 5 .Fig. 7 .
Fig.2.Schematic of the phase transition in the numerical model with regard of moments (k), modes (j ) and species (i ) given by integral moments M k i ,j .Sink and source term of the seed particle-droplet (solid line), the homogeneous (dashed line) and heterogeneous ice mode (solid line) described by Sk hom,i and Sk het,i , respectively.The modes given in this figure are broader for illustration, whereas monodisperse modes are applied in the numerical model.