skip to main content


Title: Semi-Analytic Functions to Calculate the Deposition Coefficients for Ice Crystal Vapor Growth in Bin and Bulk Microphysical Models.
Abstract Numerical cloud models require estimates of the vapor growth rate for ice crystals. Current bulk and bin microphysical parameterizations generally assume that vapor growth is diffusion limited, though some parameterizations include the influence of surface attachment kinetics through a constant deposition coefficient. A parameterization for variable deposition coefficients is provided herein. The parameterization is an explicit function of the ambient ice supersaturation and temperature, and an implicit function of crystal dimensions and pressure. The parameterization is valid for variable surface types including growth by dislocations and growth by step nucleation. Deposition coefficients are predicted for the two primary growth directions of crystals, allowing for the evolution of the primary habits. Comparisons with benchmark calculations of instantaneous mass growth indicate that the parameterization is accurate to within a relative error of 1%. Parcel model simulations using Lagrangian microphysics as a benchmark indicate that the bulk parameterization captures the evolution of mass mixing ratio and fall speed with typical relative errors of less than 10%, whereas the average axis lengths can have errors of up to 20%. The bin model produces greater accuracy with relative errors often less that 10%. The deposition coefficient parameterization can be used in any bulk and bin scheme, with low error, if an equivalent volume spherical radius is provided.  more » « less
Award ID(s):
1824243
NSF-PAR ID:
10250488
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Journal of the Atmospheric Sciences
ISSN:
0022-4928
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Ice crystal habit significantly impacts ice crystal processes such as growth by vapor deposition. Despite this, most bulk microphysical models disregard this natural shape effect and assume ice to grow spherically. This paper focuses on how the evolution of ice crystal shape and choice of ice nucleation parameterization in the adaptive habit microphysics model (AHM) influence the lake-effect storm that occurred during intensive observing period 4 (IOP4) of the Ontario Winter Lake-effect Systems (OWLeS) field campaign. This localized snowstorm produced total accumulated liquid-equivalent precipitation amounts up to 17.92 mm during a 16-h time period, providing a natural laboratory to investigate the ice–liquid partitioning within the cloud, various microphysical process rates, the accumulated precipitation magnitude, and its associated spatial distribution. Two nucleation parameterizations were implemented, and aerosol data from a size-resolved advanced particle microphysics (APM) model were ingested into the AHM for use in parameterizing ice and cloud condensation nuclei. Simulations allowing ice crystals to grow nonspherically produced 1.6%–2.3% greater precipitation while altering the nucleation parameterization changed the type of accumulating hydrometeors. In addition, all simulations were highly sensitive to the domain resolution and the source of initial and boundary conditions. These findings form the foundational understanding of relationships among ice crystal habit, nucleation parameterizations, and resultant cold-season mesoscale precipitation within detailed bulk microphysical models allowing adaptive habit.

     
    more » « less
  2. Abstract. Ice growth from vapor deposition is an important process for the evolution of cirrus clouds, but the physics of depositional ice growth at the low temperatures (<235 K) characteristic of the upper troposphere and lower stratosphere is not well understood. Surface attachment kinetics, generally parameterized as a deposition coefficient αD, control ice crystal habit and also may limit growth rates in certain cases, but significant discrepancies between experimental measurements have not been satisfactorily explained. Experiments on single ice crystals have previously indicated the deposition coefficient is a function of temperature and supersaturation, consistent with growth mechanisms controlled by the crystal's surface characteristics. Here we use observations from cloud chamber experiments in the Aerosol Interactions and Dynamics in theAtmosphere (AIDA) aerosol and cloud chamber to evaluate surface kinetic models in realistic cirrus conditions. These experiments have rapidly changing temperature, pressure, and ice supersaturation such that depositional ice growth may evolve from diffusion limited to surface kinetics limited over the course of a single experiment. In Part 1, we describe the adaptation of a Lagrangian parcel model with the Diffusion Surface Kinetics Ice Crystal Evolution (DiSKICE) model (Zhang and Harrington, 2014) to the AIDA chamber experiments. We compare the observed ice water content and saturation ratios to that derived under varying assumptions for ice surface growth mechanisms for experiments simulating ice clouds between 180 and 235 K and pressures between 150 and 300 hPa. We found that both heterogeneous and homogeneous nucleation experiments at higher temperatures (>205 K) could generally be modeled consistently with either a constant deposition coefficient or the DiSKICE model assuming growth on isometric crystals via abundant surface dislocations. Lower-temperature experiments showed more significant deviations from any depositional growth model, with different ice growth rates for heterogeneous and homogeneous nucleation experiments. 
    more » « less
  3. Abstract

    The drag coefficient, Stanton number and Dalton number are of particular importance for estimating the surface turbulent fluxes of momentum, heat and water vapor using bulk parameterization. Although these bulk transfer coefficients have been extensively studied over the past several decades in marine and large‐lake environments, there are no studies analyzing their variability for smaller lakes. Here, we evaluated these coefficients through directly measured surface fluxes using the eddy‐covariance technique over more than 30 lakes and reservoirs of different sizes and depths. Our analysis showed that the transfer coefficients (adjusted to neutral atmospheric stability) were generally within the range reported in previous studies for large lakes and oceans. All transfer coefficients exhibit a substantial increase at low wind speeds (<3 m s−1), which was found to be associated with the presence of gusts and capillary waves (except Dalton number). Stanton number was found to be on average a factor of 1.3 higher than Dalton number, likely affecting the Bowen ratio method. At high wind speeds, the transfer coefficients remained relatively constant at values of 1.6·10−3, 1.4·10−3, 1.0·10−3, respectively. We found that the variability of the transfer coefficients among the lakes could be associated with lake surface area. In flux parameterizations at lake surfaces, it is recommended to consider variations in the drag coefficient and Stanton number due to wind gustiness and capillary wave roughness while Dalton number could be considered as constant at all wind speeds.

     
    more » « less
  4. Abstract

    Ice habit diagrams published prior to 2009—and many since—do not accurately describe in situ observations of ice shapes as a function of temperature and moisture. Laboratory studies and analysis of field observations by Bailey and Hallett in a series of papers in 2002, 2004, and 2009 corrected several errors from earlier studies, but their work has not been widely disseminated. We present a new, simplified diagram based on Bailey and Hallett’s work that focuses on several ice growth forms arising from the underlying surface processes by which mass is added to a crystal: tabular, columnar, branched, side branched, two types of polycrystalline forms, and a multiple growth regime at low ice supersaturations. To aid interpretation for a variety of applications, versions of the ice growth diagram are presented in terms of relative humidity with respect to water as well as the traditional formats of relative humidity with respect to ice and vapor density excess. These diagrams are intended to be understandable and useful in classroom settings at the sophomore undergraduate level and above. The myriad shapes of pristine snow crystals can be described as the result of either a single growth form or a sequence of growth forms. Overlays of data from upper-air soundings on the ice growth diagrams aid interpretation of expected physical properties and processes in conditions of ice growth.

     
    more » « less
  5. Abstract An electrodynamic levitation thermal-gradient diffusion chamber was used to grow 268 individual, small ice particles (initial radii of 8–26 μ m) from the vapor, at temperatures ranging from −65° to −40°C, and supersaturations up to liquid saturation. Growth limited by attachment kinetics was frequently measured at low supersaturation, as shown in prior work. At high supersaturation, enhanced growth was measured, likely due to the development of branches and hollowed facets. The effects of branching and hollowing on particle growth are often treated with an effective density ρ eff . We fit the measured time series with two different models to estimate size-dependent ρ eff values: the first model decreases ρ eff to an asymptotic deposition density ρ dep , and the second models ρ eff by a power law with exponent P . Both methods produce similar results, though the fits with ρ dep typically have lower relative errors. The fit results do not correspond well with models of isometric or planar single-crystalline growth. While single-crystalline columnar crystals correspond to some of the highest growth rates, a newly constructed geometric model of budding rosette crystals produces the best match with the growth data. The relative frequency of occurrence of ρ dep and P values show a clear dependence on ice supersaturation normalized to liquid saturation. We use these relative frequencies of ρ dep and P to derive two supersaturation-dependent mass–size relationships suitable for cloud modeling applications. 
    more » « less