^{1}Institute for Terrestrial and Planetary Atmospheres/School of Marine and Atmospheric Sciences, Stony Brook University, Stony Brook, New York, USA ^{*}These authors contributed equally to this work. Abstract. Immersion freezing of water and aqueous (NH_{4})_{2}SO_{4} droplets containing leonardite (LEO) and Pahokee peat (PP) serving as surrogates for humiclike substances (HULIS) has been investigated. Organic aerosol containing HULIS are ubiquitous in the atmosphere; however, their potential for ice cloud formation is uncertain. Immersion freezing has been studied for temperatures as low as 215 K and solution water activity, a_{w}, from 0.85 to 1.0. The freezing temperatures of water and aqueous solution droplets containing LEO and PP are 5–15 K warmer than homogeneous ice nucleation temperatures. Heterogeneous freezing temperatures can be represented by a horizontal shift of the ice melting curve as a function of solution a_{w} by Δa_{w} = 0.2703 and 0.2466, respectively. Corresponding hetrogeneous ice nucleation rate coefficients, J_{het}, are (9.6 ± 2.5)×10^{4} and (5.4 ± 1.4)×10^{4} cm^{−2} s^{−1} for LEO and PP containing droplets, respectively, and remain constant along freezing curves characterized by Δa_{w}. Consequently predictions of freezing temperatures and kinetics can be made without knowledge of the solute type when relative humidity and ice nuclei (IN) surface areas are known. The acquired ice nucleation data are applied to evaluate different approaches to fit and reproduce experimentally derived frozen fractions. In addition, we apply a basic formulation of classical nucleation theory (α(T)model) to calculate contact angles and frozen fractions. Contact angles calculated for each ice nucleus as a function of temperature, α(T)model, reproduce exactly experimentally derived frozen fractions without involving freefit parameters. However, assigning the IN a single contact angle for the entire population (singleα model) is not suited to represent the frozen fractions. Application of αPDF, active sites, and deterministic model approaches to measured frozen fractions yield similar good representations. Furthermore, when using a single parameterization of αPDF or active sites distribution to fit all individual a_{w} immersion freezing data simultaneously, frozen fraction curves are not reproduced. This implies that these fitting formulations cannot be applied to immersion freezing of aqueous solutions, and suggests that derived fit parameters do not represent independent particle properties. Thus, from fitting frozen fractions only, the underlying ice nucleation mechanism and nature of the ice nucleating sites cannot be inferred. In contrast to using fitted functions obtained to represent experimental conditions only, we suggest to use experimentally derived J_{het} as a function of temperature and a_{w} that can be applied to conditions outside of those probed in laboratory. This is because J_{het}(T) is independent of time and IN surface areas in contrast to the fit parameters obtained by representation of experimentally derived frozen fractions. Citation: Rigg, Y. J., Alpert, P. A., and Knopf, D. A.: Immersion freezing of water and aqueous ammonium sulfate droplets initiated by humiclike substances as a function of water activity, Atmos. Chem. Phys., 13, 66036622, doi:10.5194/acp1366032013, 2013. 
