The WeIzmann Supercooled Droplets Observation ( WISDOM ) on a 1 Microarray and application for ambient dust

The WeIzmann Supercooled Droplets Observation on Microarray (WISDOM) is a new setup for studying ice 7 nucleation in an array of monodisperse droplets for atmospheric implications. WISDOM combines microfluidics 8 techniques for droplets production and a cryo-optic stage for observation and characterization of freezing events of individual 9 droplets. This setup is designed to explore heterogeneous ice nucleation in the immersion freezing mode, down to the 10 homogeneous freezing of water (235 K) in various cooling rates (typically 0.1-10 K min). It can also be used for studying 11 homogeneous freezing of aqueous solutions in colder temperatures. Frozen fraction, ice nucleation active surface site (INAS) 12 densities and freezing kinetics can be obtained from WISDOM measurements for hundreds of individual droplets in a single 13 freezing experiment. Calibration experiments using eutectic solutions and previously studied materials are described. 14 WISDOM also allows repeatable cycles of cooling and heating for the same array of droplets. This paper describes the 15 WISDOM setup, its temperature calibration, validation experiments and measurement uncertainties. Finally, application of 16 WISDOM to study the INP properties of size-selected ambient Saharan dust particles is presented. 17


Introduction
In mixed phase clouds, water droplets remain stable in a supercooled state below 273 K and ice nucleates spontaneously as droplets reach the homogeneous freezing temperature, below 236 K (Pruppacher et al., 1998).At warmer temperatures, ice particles may coexist with supercooled droplets, due to heterogeneous nucleation facilitated by the presence of ice nuclei particles (INPs) (Cantrell and Heymsfield, 2005).In cases where INPs are immersed in the droplet before supercooling, referred to as immersion freezing mechanism, the droplets first grow to supercritical size before freezing occurs (de Boer et al., 2011).
Observations and modeling studies suggest that immersion freezing is the prominent mechanism for heterogeneous ice formation in mixed phase clouds (Ansmann et al., 2008;Field et al., 2012;Nagare et al., 2016;Possner et al., 2017;Rosenfeld and Woodley, 2000).
Ice particles affect the radiative and microphysical properties of mixed phase clouds and Earth's hydrological cycle.Therefore, they can influence present and possibly future climate (Hoose and Möhler, 2012;IPCC, 2013).Studying ice formation in clouds is hence important, and yet, due to its complexity, this process is still not fully understood and presents a great challenge to laboratory and field researchers as well as for clouds and climate modelers (DeMott et al., 2010;Schnaiter et al., 2016;Ullrich et al., 2017).
Offline studies of immersion freezing often use cold stage techniques (Budke and Koop, 2015).The basic idea is to place an array of droplets over a cold stage and cool continuously until all are frozen, to obtain a quantitative measurement of their corresponding freezing temperatures (Vali, 1971).The droplets may be microliter-sized and observed with a simple camera.
Smaller droplets, down to the pico-liter range, are usually observed under a microscope.In both cases, freezing events are identified by optical changes in the droplets when they crystalize (Atkinson et al., 2013;Hiranuma et al., 2015;Knopf and Lopez, 2009;Murray et al., 2011).
Cold stage techniques may suffer from technical issues such as droplets evaporation and vapor transfer due to the Wegener-Bergeron-Findeisen process, where ice grows on the expense of supercooled droplets or from seeding of neighboring droplets by formation and surface growth of frost halos (Budke and Koop, 2015).Some cold stages instruments place oil over the droplets or use droplet in oil emulsions to prevent these effects (Murray et al., 2012).Still, results from cold stage experiments may be biased by effects of inhomogeneous temperature of the substrate and the surroundings or by various contaminations caused during droplets' preparation and measurement (Hiranuma et al., 2015).Furthermore, supercooling is limited due to the presence of impurities, which increases with the volume of the droplet.Hence, to allow comprehensive studies down to the homogeneous region, low volumes (<1 µL) are used and generation of these volumes is not trivial and may cause further complications.
Microfluidics is a technology of fluids manipulation in micro-channels array on a small device.Microfluidics is widely used in a range of fields, such as physics, chemistry, biology, life sciences and the food industry (Neethirajan et al., 2011;Sackmann et al., 2014;Whitesides, 2006).Recent studies used microfluidic apparatus to study ice nucleation processes.Riechers et al. (2013) used a microfluidics device to produce and collect monodisperse droplets of water in various sizes, which were subsequently observed under a microscope to study their homogeneous freezing.Stan et al. (2009) recorded nucleation in water droplets and silver iodide seeded droplets, while droplets were flowing during cooling.Schmitz et al. (2009) established 'Dropspots', a static microfluidic array of droplets, later used by Edd et al. (2009) to measure nucleation kinetics.However, in the atmospheric heterogeneous ice nucleation field, microfluidics techniques are not widely adopted, despite many potential advantages.
The WeIzmann Supercooled Droplets Observation on Microarray (WISDOM) is a new instrument combining the cold stage technique with microfluidics technology, and is designed to study immersion freezing of micrometer-sized droplets, while addressing most of the technical issues listed above.The WISDOM setup introduces several advantages of microfluidics to the atmospheric ice nucleation field.WISDOM is based on the 'Dropsopts' static array (Schmitz et al., 2009), which enables the separation and the fixation of the droplets, so that each individual droplet is recorded and studied, and also can be used for repetition of freezing cycles and further exploration of the nucleation process of a specific sample.
In this paper, we present the WISDOM setup, its calibration and validation procedures.For validation experiments, homogeneous and heterogeneous freezing were examined by following homogeneous freezing rates of pure aqueous or solutions or deriving the efficiencies of three types of mineral dust surrogates and collected ambient Saharan dust to validate heterogeneous freezing experiments.Homogeneous nucleation rates are described stochastically using the volume-dependent ice nucleation rate (Jv(T)) in supercooled droplets, given by the frozen fraction (fice) of droplets with volume V at a certain temperature (T) and time intervals (∆t) (Alpert et al., 2011;Murray et al., 2010;Riechers et al., 2013), Heterogeneous freezing is described by a singular approach that assumes that nucleation occurs at a certain temperature due to special nucleation site.Hence, a cumulative number of nucleation sites per unit surface area, ns, is used to describe the heterogeneous nucleation efficiency at a certain temperature, Where fice is the fraction of frozen droplets at temperature T, and A is the specific surface area of the immersed particles in each droplet (Vali, 1971;Vali et al., 2015;Whale et al., 2015).
Validation of WISDOM was further extended below the homogeneous nucleation temperature of pure water, using aqueous solutions as the freezing temperatures of solutions decrease as a function of the solution water activity.

Droplets production and trapping
The WISDOM setup, shown in Figure 1, is made of a microfluidic setup which include a pressure controlled pump with 4 independent flow channels (OB1 MK3 by Elveflow), a stereoscope (SMZ-171 by Motic) that permits a full view of all channels and inlets, and a CCD camera (GS3 by Point Grey) that enable real time monitoring of the droplets production.The flows in the channels are continuous and controlled by the pressure pump.One channel is connected to the continuous (oil) phase, and a second channel contains the sample (aqueous solution that can contain INPs).The two phases meet in a narrow junction where monodisperse droplets are generated due to the pressure exerted by one phase over the other.The ratio between the flows determines the size of the emerging droplets; the volume increases with increasing flow rate of the sample.In this setup, droplets are suspended in an oil mixture, consisting of mineral oil (Sigma Aldrich) and a 2 weight percent (wt%) nonionic surfactant (span80®, Sigma Aldrich), added for droplets stabilization (Riechers et al., 2013).Hence, an array of picoliter (micrometer-size) droplets is generated directly on a device.
The principle of the Schmitz et al. (2009) design is that the droplets flow into round chambers that are connected by constriction channel.At a certain flow, the droplets are squeezed through the constriction channel and the array fills up with droplets.
When the flow is too weak or stopped, the constriction channel stops the droplets' movement and they are trapped in the chambers.The droplets are isolated and stable in the chambers and it is safe to move the device from the generation stage to the cold stage for the freezing experiments.
Devices are fabricated following the Schmitz et al. (2009) protocol.Briefly, the device pattern is imprinted on a polydimethylsiloxane (PDMS) polymer, later glued to a 1mm thick microscope glass slide using air plasma treatment.After the plasma treatment the PDMS surfaces are hydrophilic (Eddings, 2008).Therefore, the devices were used only in the following day, after their surfaces became hydrophobic, following their exposure to the atmosphere, or after their annealing at 60C for half an hour.

Freezing experiments and detection
The droplets array is placed in a commercial cryostage (Linkam, THMS600) coupled to an optical microscope (Olympus, BX-51 with 10X magnification, transmission mode).Experiments are monitored by a microscope mounted CCD camera (Allied Vision Technologies, Oscar F-510C) for automatic identification of droplets and their freezing events.Both the device and the cooling stage are cleaned with 2-propanol.Then the device is placed over the stage together with a thin layer of oil on its bottom to provide good thermal conductivity.Each freezing experiment starts with dry N2 purging to replace the moist atmosphere inside the cryostage to prevent condensation.During the experiment, N2 flow prevents water condensation on the cryostage window.Freezing experiments are conducted with a cooling rate of 1 K min -1 which is relevant for atmospheric conditions and also allows good thermalization of the droplets, as will be shown in the calibration section (section 3.1).Each cooling cycle is followed by a heating cycle, where melting is observed.Analysis of the melting onset is then used to verify that the thermal conductivity is good and thus validate the measurement.
In-house LabVIEW software is used to record a freezing experiment movie file and analyze it offline.The temperature readings by the Linkam cryostage temperature sensor (<±0.25 K for the operated temperature range) and the movie frames are synchronized and integrated.In most cases, 1 second (or 0.017 K at 1 K min -1 ) per frame is used.Currently, the WISDOM setup operates with two types of devices that differ in their droplets' trap diameter: 40 µm and 100 µm.Approximately 550 and 120 droplets can be monitored per experiment in the smaller and larger diameter devices, respectively.Statistically, for the same sample, larger droplets encompass more INP surface area within each droplet, which can be more sensitive for detecting rare active sites.The device can be reused for the same sample, if it is not clogged or destroyed during the experiment.
However, because the channels of the 40 µm device are smaller they tend to clog faster (for instance by large particles).

Automatic detection of phase transitions
The optical brightness of a droplet changes during a phase transition (freezing or melting) due to the different interaction of light with the liquid and the solids.For phase transition detection, an in-house image processing LabVIEW program monitors automatically the optical brightness change.The program detects the droplets using a spherical shape criterion and sets a square surrounding the droplet that defines an array of pixels that are attributed to that specific droplet.A change in the optical brightness is represented by the gray level value of the image's pixels, ranging from 0 to 255.Freezing is calculated per movie frame and is defined as the subtraction of the brightness mean value for each droplet in two consecutive frames (ΔGL), thus allowing derivation of freezing rates.At the beginning of the analysis, the first 15 frames are used to identify the noise level of the signal by calculating its standard deviation std(ΔGL).The program then searches for the maximal freezing signal that is also greater than 5 times the noise level.The temperature associated with this freezing signal is assigned as the freezing temperature for that droplet.
In this algorithm, the program can distinguish successfully between a phase transition event and noise that arises from the camera signal, droplet movement or any other interruption.Figure 2 presents a spectral analysis for different types of phase transitions observed in WISDOM.Since WISDOM operates in transmission microscopy mode, the light is scattered more efficiently by ice crystals in comparison with a liquid droplet and a freezing event involves droplet darkening and a negative signal.The negative signal of a freezing event of a single droplet is shown in Figure 2a.In comparison, during melting, the droplet becomes brighter until all the crystals melt, and the signal is positive.In Figure 2b+c the analysis of a melting signal and a eutectic melting signal are presented for the entire frame.instruments.Therefore, we performed a thorough temperature calibration using the known eutectic melting points and the melting points of several aqueous solutions as calibration reference points.Although ice nucleation experiments are performed while cooling, the calibration experiments were done while heating to improve the calibration precision and to avoid biases associated with supercooling of the liquids (Budke and Koop, 2015).

Droplets thermalization
The temperature of the Linkam stage was measured at the upper center part of the cooling stage and hence may differ from the actual temperature of the droplets in the device due to thermal effects such as temperature gradients and temperature lag.
During cooling or heating, a vertical temperature gradient may develop between the top of the device, in contact with the inner ambient of the cryostage, and the bottom of the device, which is in contact with the cooling silver block.This gradient is expected to increase in magnitude, as the temperature of the stage decreases or increases below or above ambient temperature.
Edd et al. ( 2009) used a similar setup and found a difference between the top temperature and the bottom temperature of about 2 K around 237 K and 3 K around 227 K. Stan et al. (2009) also reported a vertical gradient of 1-2 K, that was reduced to 0.5 K with a flow of cooled N2 over their device.In addition, a thermal lag may arise during cooling or heating as the rate of temperature change is high and precludes proper temperature equilibration.Hence, a more accurate measurement of the droplet temperature is taken as a sum of the stage temperature with the contributions of both thermal gradient and lag.
Figure 3 demonstrates the combined effects of temperature change rate and device properties on the thermalization of pure water droplets (double distilled, 18.2 MΩ cm).Specifically, freezing and melting experiments at different rates were performed.The temperature difference (ΔT) is the difference between the measured values and the extrapolated temperature at equilibrium conditions (0 K min -1 ).As expected, at slower temperature cooling (heating) rates, the droplets are more equilibrated with the stage temperature and ΔT is negligible.However, ΔT increases at higher temperature cooling (heating) rates (e.g.; 10 K min -1 ).We observed that during cooling (heating) the droplet is warmer (colder) than the stage and will freeze (melt) at colder (warmer) temperature at higher cooling (heating) rates.We also found that because ΔT is higher, in absolute value, for devices of thicker PDMS and/or in devices which hold larger droplets, it should be considered in the final temperature calibration for these scenarios.Furthermore, ΔT was found to be almost symmetric for higher temperature cooling (heating) rates.However, for 1 K min -1 , ΔT during cooling is higher than that for heating.Our conjecture is that this can be an effect of the higher thermal gradient that develops as the temperature decreases well below ambient (236 K).

Melting of aqueous solutions
Figure 4 presents the measured melting points of NaCl solutions with different water activities.Reported melting points represent the temperature in which all ice crystals in the droplets completely melted, in contrast with melting temperatures reported for pure liquids such as water, where the onset of melting is defined as the melting point.Melting temperature results were consistent with theoretical melting temperatures reported in Koop and Zobrist (2009).This provides support to our conclusion that droplets thermalize with the cooling stage when using a heating rate of 0.1-1 K min -1 .For faster heating rates (i.e. 10 K min -1 ), the thermal lag was more pronounced, leading to a melting point shift of about 2-3 K.For more concentrated solutions, faster heating rates shifted the melting points more.

Melting of eutectic solutions
Some aqueous solutions, such as NaCl and MgCl2, arrange in a super-lattice at a certain wt% to form a solid with a welldefined melting point (eutectic) (252.05K for NaCl and at 239.95 K for MgCl2) (Borgognoni et al., 2009;Farnam et al., 2016).
Interestingly, this type of melting has a smaller optical signature compared to that of melting points of pure substances, as can be seen in Figure 2b.We have set a specific water activity for a solution by determining its quantitative composition using the extended aerosol inorganic model (E-AIM) (Clegg et al., 1998) at room temperature (298 K).For calibration purposes, because eutectic melting had a negligible variation for different water activities used in the range of 0.99 to 0.95, we decided to take their average to achieve a single melting value.These eutectic melting temperatures are colder than the melting point of pure water and, therefore, are used for expanding WISDOM calibration range.
The final calibration is obtained for a device with a specific PDMS thickness and at a specific cooling (heating) rate.For example, devices with 100 μm diameter sized droplet and of 4mm PDMS thickness have a linear calibration curve of Tdrop=0.97Tstage-0.46at 0.1 K min -1 .

Measurement reproducibility and device variability
Device's inter-variability was determined from 20 devices by comparing their corresponding homogeneous freezing temperatures of pure water.Specifically, each device was recycled three times with freshly prepared droplets.Our results showed high reproducibility in the median freezing temperature, where 50% of the probed droplets froze (T50), and high reproducibility in the melting point temperature.Variation within devices was always smaller than ±0.2K at 1 K min -1 and 0.1 K min -1 (variation within the devices over the whole freezing range is presented in Appendix A).

Homogeneous freezing rates of pure water
Homogeneous nucleation in supercooled water occurs in WISDOM between 238 and 237 K for a cooling rate of 1 K min -1 and droplets diameter of 100 μm. Figure 5 shows WISDOM nucleation rates in comparison with other similar instruments.It is seen that the slope of the rate and temperatures are similar to the slopes reported for other instruments.The temperature where 50% of droplets froze (T50) is also in the expected range according to model results of Hoffer (1961).WISDOM rates are slightly slower, but within the uncertainty of the instruments used by Riechers et al. (2013) and Stan et al. (2009).Stöckel et al. (2005) show a higher nucleation rate.This discrepancy can be explained by a decrease in the number of surface nucleation events due to the oil phase surrounding our droplets, whereas in Stöckel et al. (2005), droplets are suspended in air which allows surface nucleation may occur.

Homogeneous and heterogeneous freezing of aqueous solutions
The water-activity-based ice nucleation theory by Koop et al. (2000) describes the dependence of the freezing temperature depression on the water activity (aw) of the solution, regardless of the solute nature.Figure 6 presents the theoretical freezing and melting temperature curves from Koop et al. (2000) with homogeneous ice nucleation results measured in WISDOM, for four solutions with atmospheric relevance.Water activities for NaCl, ammonium sulfate (AS), glucose and levoglucosan mixtures were derived from the AIM model and were corrected for glucose and levoglucosan, for which water activity is temperature dependent (Knopf and Lopez, 2009;Zobrist et al., 2008).The experiments were conducted at 1 K min -1 for 40 and 100 μm droplet diameters.The results follow the theoretical curves of the water-activity-based ice nucleation, and the dependence of the homogeneous freezing on the droplet volume is as expected (Hoffer, 1961;Kuan-Ting and Wood, 2016) as the curve of the smaller diameter droplets (green curve) is slightly colder compared with the larger volume droplets (dark green curve).
Similar experiments were conducted for 0.1 wt% of Arizona Test Dust particles (ATD, Powder Technology Inc.) immersed in glucose solution droplets.The ATD particles facilitate the ice nucleation at warmer temperatures, in agreement with similar studies (Hartmann et al., 2011;Niedermeier et al., 2010), and the freezing depression follow the water-activity-based ice nucleation curves.Here, the dependence of the freezing point on the droplet volume is more pronounced, as the surface area of the immersed particles is higher, hence they contain higher number of nucleation sites (Marcolli et al., 2007) as will be shown in the next section for two more types of dust.
Below 223 K, ice nucleation occurs at slightly lower temperatures than expected by the theoretical freezing curve.As the WISDOM temperature calibration is not valid in this temperature range, we cannot conclude if this is due to a change of the thermal conductivity of the device or an effect of the high concentration of the solute in the water.

Standard dust powder
Heterogeneous freezing efficiencies of suspended mineral dusts K-Feldspar and Illite-NX in supercooled water droplets are presented in Figure 7 and summarized in Table 1, and are compared to recent published data.The particles are suspended at different wt% and the frozen fraction of each suspension is derived as a function of temperature as represented by the color bar.To examine the freezing efficiency and compare the different mineral dust types, the results are normalized to the surface area within each droplet.Experiments were performed at 1 K min -1 for 40 and 100 μm droplets diameters.Suspension preparation and evaluation of the surface area are described in the appendix B.
The results demonstrate the effect of dust surface area immersed in the droplets on the freezing parameters.The freezing temperatures increase with increasing surface area and are also reflected in the warming of the median frozen fraction (T50) colored in yellow.The spectra of the number of nucleation sites per unit surface area (ns) also support surface area dependence because all spectra converge to a single line.The ns results show the increase of nucleation sites at colder temperatures.Results from WISDOM are in good agreement with similar analyses from other instruments.In particular, ns is in best agreement with the Leeds-NIPI (Broadley et al., 2012;Murray et al., 2011) results both for K-Feldspar and for Illite-NX particles.Results of Illite-NX particles are also in good agreement with the Binary instrument (Budke and Koop, 2015) and reside within the uncertainty of both instruments.The linear trend of few wt% support the assumption that particles in suspension are uniformly distributed and the droplets contain approximately the same surface area.

Ambient mineral dust
WISDOM can also be used for analyzing collected ambient particles.Mineral dust particles were collected in Rehovot, Israel (31.9N, 34.8E about 80m AMSL), during dust storm event on 12-13 March 2017).The dust was transported from the Sahara Desert and North Africa.Size-segregated ambient dust particles were collected on cyclopore polycarbonate filters using a Micro-orifice Uniform deposit Impactor (MOUDI; MSP Corporation model 110-R, (Marple et al., 1991)), that operated at 30 L min -1 and for 24 hrs, similarly to Huffman et al. (2013) and Mason et al. (2015).MOUDI has eleven stages with cut points (D50) of 0.056, 0.10, 0.18, 0.32, 0.56, 1.0, 1.8, 3.2, 5.6, 10, and 18 µm.The size distribution of the particles was obtained by Optical Particle Counter (OPC; GRIMM Technologies model 1.109) in the range of 0.25-32 µm, and used for estimations of surface area immersed in the droplets (further details in Appendix C).
For heterogeneous freezing experiments, a quarter of each filter is placed with 300 µL DDW in 1.5 ml Eppendorf vial and particles were extracted by intensive dry sonication (Hielcher; model UP200St VialTweeter).In Figure 8, the spectra of the nucleation sites per unit surface area (ns) of three super-micron stages (D50 of 1.0, 1.8, 3.2 µm) are presented and summarized in T1.It is also seen that there are slightly more active sites for the larger particles (3.2 µm), as their surface area is higher and there is a higher probability to contain an active site.In Figure 9, ns curves of the collected dust is compared to references of K-feldspar standard particles, analyzed in different instruments (the Leeds-NIPI (Atkinson et al., 2013), LACIS (Niedermeier et al., 2015)) and to measurements of ambient dust samples, from different locations around the world, including Israeli settled dust, that was analyzed in the AIDA chamber.Moreover, the freezing of the size resolved mineral dust analyzed in this study by WISDOM (slope in the temperature range) is consistent with the (grey) polygon that represents the estimated freezing efficiency for natural concentrations of K-feldspar in internally mixed mineral types (Atkinson et al., 2013).The results are also in agreement with Niemand et al. (2012), especially between 243 and 249 K.At warmer temperatures, ns of ambient dust in this study showed lower efficiency than in Niemand et al. (2012).This difference can extend to one order of magnitude in ns, and is more pronounced at smaller particles that were analyzed (around 1-1.8µm diameter).For the larger particles, more nucleating sites are observed.Both the current study and Niemand et al. ( 2012), suggest that K-feldspar is involved with the warmer part of their data (>248 K) as the results are consistent with the Atkinson et al. (2013) scale for ambient samples.The slope of the ns derived in this study is similar to the slope from Atkinson et al. (2013) at warmer temperatures.For the colder regime, the slope is similar to the slope presented by standard K-feldspar particles in Niedermeier et al., (2015).

WISDOM in comparison to other cold stage instruments
The microfluidics technology used in WISDOM solves some substantial issues inherent in other currently used instruments: (1) good control of the size and number of monodisperse droplets, (2) fast production of hundreds of nearly monodisperse droplets minimizes sample sedimentation or agglomeration that may occur in a suspension, leading to a good estimation of the surface area of the suspended material.Moreover, several droplet diameters can be employed in the same device without its modification, (3) good statistics achieved by individual analysis of hundreds of droplets, (4) monodisperse droplets individually analyzed, in contrast to some emulsion techniques (such as Differential Scanning Calorimeter (DSC)), allow to obtain the frozen fraction at each temperature, and to achieve detailed information about active sites and freezing rates, (5) the use of oil minimizes possible artefacts from droplets' evaporation, neighbor seeding or vapor transfer due to the Wegener-Bergeron-Findeisen processes, (6) the small volumes decrease freezing artefacts by impurities, thus allowing to reach the homogeneous freezing threshold (-37ºC), ( 7) possible investigation of several freezing cycles for the same droplets, (8) the microfluidics method and the small droplet volumes enable working with small sample volumes which can be an advantage when working with atmospheric samples.
WISDOM has a very accurate temperature calibration that spans a wide temperature range, using the eutectic freezing method.
WISDOM most resembles the instrument used by in Edd et al. (2009).However, it seems that issues with temperature calibration in Edd et al. (2009) led to a temperature offset, and hence different freezing rates.Stan et al. (2009) achieved better temperature accuracy and high statistics.However, the freezing experiment was conducted in a flow mode, which is more complicated than in the WISDOM setup and requires complicated modeling.In addition, the cooling rates that were used were very fast, which induces additional errors.Riechers et al. (2013) had high temperature accuracy as they also used a DSC.
However, they had to collect the droplets from the device as there was no static array option and this may add further complication and contamination.
The microfluidics technology has also disadvantages.These may include: (1) oil may interact with some of the analyzed particles, possibly leading to biased data, (2) the microchannels are susceptible to clogging, (3) it is not possible to perform any post analysis to the droplets content after the experiment, (4) the small droplets' volumes reduce the sensitivity to rare active sites.This may be solved by performing many experiments or by using larger droplets with more surface area within the droplets.

Summary and conclusions
The new setup WISDOM is based on microfluidics technology and its detailed validation is presented.Based on a set of validation measurements and a good agreement with other instruments, we conclude that WISDOM is a suitable tool for studying atmospheric ice nucleation, both in homogeneous and heterogeneous immersion freezing modes.Results of homogeneous freezing correspond to water-activity-based nucleation theory in supercooled droplets and represent well volume nucleation rates.Heterogeneous freezing in supercooled droplets also agrees well with literature data.Furthermore, freezing efficiency dependence on the particles surface area within the droplets is clearly observed.Using microfluidics allows a mass production of picoliter monodisperse droplets using low volumes of suspensions, which can be beneficial for immersion freezing studies over a wide range of supercooling down to homogenous temperature region.The good reproducibility of the devices, proved using pure water freezing cycles, enables the recycling of the same device for few freezing cycles.It is also shown that the temperature uncertainty can be reduced if the temperature calibration includes the microfluidic devices properties in the working temperature change rates, especially for melting experiments.In this work we have also demonstrated how WISDOM can be applied for studying the ice nucleation properties of ambient samples that contain very small quantity of sample.The particles were collected using the MOUDI during Saharan dust storm event.Results are in correspondence with literature data of ambient dust and further support Atkinson et al. (2013) and the possible importance of K-feldspar for ice nucleation in clouds, but further analysis of the mineralogy is still needed in order to verify that.

Appendix C: Collection of ambient particles during dust storm events in Rehovot
The GRIMM measurement was synchronized to the MOUDI stages for the estimation of the total surface area that was collected on the filter for droplets surface area estimation.For that, two base assumptions were made: (1) all the particles that were collected are extracted to the water later used for the freezing experiments, (2) sphericity of the particles.The GRIMM bins are synchronized to the MOUDI stages based on collection efficiency of the MOUDI, obtained from Marple et al. (1991).
For example, on certain MOUDI stage, all the particles that own diameter that is larger than the D50 have high chance to be impacted on that stage.All the rest of the sizes, that are smaller in their diameter, will continue to the next stage and will have high chance to deposit there.Hence, the GRIMM's bins were synchronized to the MOUDI D50 stages.For the ns calculations, the surface area was based on the number of particles that were measured in a certain bin and their total surface area.To calculate the surface area of a particle (assuming sphericity) in a certain bin, the midpoint of that bin was used as a radius.To calculate the total mass of the particles in each filter, dust density of Quartz was used (2.65 g cm -3 ), as this is usually the dominant mineral (Mahowald et al., 2014).The error of the ns data is propagated from the error in the frozen fraction, the error of the droplet's volume and the error of the MOUDI's collection efficiency in the different stages, the later was the dominant one.
For control, analysis of blank filter was done.The blanks were sonicated before analysing them and freezing was mostly colder than the freezing temperatures that are presented here and hence no special reduction of the final active sites was done.         is also presented.References of K-feldspar standard particles activated in WISDOM, Leeds-NIPI (Atkinson et al., 2013) and LACIS (Niedermeier et al., 2015) instruments are presented, as well as ambient dust particles that were analyzed in AIDA and included Israeli dust (Niemand et al., 2012).
Temperature accuracy is a most important parameter in ice nucleation experiments.An error propagation analysis byRiechers et al. (2013) demonstrated how the temperature uncertainty may lead to a distribution of temperatures between different

Figure B1 .
Figure B1.Accumulated active site density spectra (ns) of K-feldspar (a) and Illite-NX (b) particles with different surface areas suspended in water at cooling rate of 1°C per minute.The error bars are located at three representative frozen fractions, 0.1, 0.5 and 0.9.

Figure 1 :
Figure 1: The WISDOM setup.a) The design of the microfluidic device is based on Schmitz et al. (2009).Aqueous solutions (including the sample) and oil are connected through the inlets and merge in a junction to generate monodisperse droplets.Subsequently, droplets flow into a trap array and settle in them as the flow is stopped.The device is transferred into a cooling stage for subsequent freezing experiments.b) upper and c) side views of the device, which is made of PDMS, plasma glued to a microscope glass slide, placed over the cooling silver block.

Figure 2 :
Figure 2: Spectra of different phase transition events as observed in WISDOM.a) freezing, b) eutectic melting, and c) melting onset and clear point (liquefaction) are the mean of all sampled droplets in a single experiment.The phase transition is defined optically by the brightness information obtained by the gray level of the image pixels.std(ΔGL) describes the standard error of the difference in mean GL for two consecutive frames.At the beginning of the experiment the noise level is studied and freezing, or melting are detected only if std(ΔGL) is as least 5 times greater than the noise std level.Freezing and melting examples are for pure water droplets and the eutectic melting example is for aqueous solution droplets of NaCl.Eutectic melting point of NaCl and pure water melting point are marked by the yellow and red lines in b and c, correspondingly.In all cases the droplets diameter was 100 μm.

Figure 3 .
Figure 3.The temperature difference (ΔT), defined as the temperature difference between the stage temperature and the droplet extrapolated temperature at equilibrium conditions at different cooling (heating) rates.Freezing and melting points of pure water are represented by circles and squares (40 and 100 μm droplet diameter, in correspondence) for different PDMS thicknesses and are represented by different colors.C denotes cooling and H denotes heating.Droplets are close to equilibrium with the stage temperature at rates <0.1 K min -1 and ΔT increases with increasing temperature change rate and with the PDMS height.

Figure 4 .
Figure 4. Temperature calibration by melting points of eutectic solutions and pure water droplets for different heating rates.Calibration is presented for 100 μm droplets with 4 mm PDMS thickness.The onsets of pure water droplets are also considered.Eutectic melting is used for the colder temperatures range (<253 K) while clear point (liquefaction) at various water activities is taken for the warmer temperature range.The upper panel presents the temperature difference between the reference value and the cooling stage temperature after calibration.Most of the differences are within the range ±0.2 K.

Figure 5 .
Figure 5.The volume-dependent homogeneous freezing of pure water, derived for 100 μm droplets with 4 mm PDMS height.WISDOM rates are compared to relevant literature data.The obtained fit from WISDOM is Jv(T)=exp(-3.4T+817.6).Temperature uncertainty for WISDOM is ±0.3 K.

Figure 6 .
Figure 6.Homogeneous and heterogeneous ice nucleation temperatures for 40 and 100 μm aqueous solution droplets as a function of solution water activity.Freezing and melting curves are derived fromKoop et al. (2000).Heterogeneous ice nucleation is performed with 0.1 wt % ATD particles immersed in the droplets.

Figure 7 .
Figure 7. Accumulated active site density spectra (ns) of K-feldspar and Illite-NX particles as a function of temperature from validation experiments of immersion freezing in WISDOM.Frozen fraction values are represented by a color bar, for few surface area values that are exposed in 40 and 100 μm droplets.The dependence of the nucleation site density on the surface area is illustrated here.WISDOM uncertainties, propagated from surface area estimation and measured frozen fraction errors, are included within the size of the markers.For validation, previous immersion freezing measurements are also presented (Hiranuma et al. (2015) and Atkinson et al. (2013)).T-binned data (1°C) normalized by the BET surface area from Hiranuma et al. (2015) is presented in the color squares, only for wet suspension analysis (Binary (red), CSU-IS (orange), Leeds-NIPI (purple), M-AL (green), M-WT (black), NC-State-CS (brown) and CU-RMCS (blue)).Hiranuma et.al. log fit and ns (BET) parameterization are also presented.

Figure 8 .
Figure 8. Accumulated active site density spectra (ns) of ambient super-micron mineral dust particles collected in Israel during Saharan dust events in 2017, for three different sampling stages of the MOUDI; D50 of 1, 1.8 and 3.2 µm.