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 more » bin scheme, with low error, if an equivalent volume spherical radius is provided. « less
; ;
Award ID(s):
Publication Date:
Journal Name:
Journal of the Atmospheric Sciences
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 amongmore »ice crystal habit, nucleation parameterizations, and resultant cold-season mesoscale precipitation within detailed bulk microphysical models allowing adaptive habit.

    « 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 andmore »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.« less
  3. 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 ofmore »ρ 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.« less
  4. There are few measurements of the vapor growth of small ice crystals at temperatures below -30°C. Presented here are mass-growth measurements of heterogeneously and homogeneously frozen ice particles grown within an electrodynamic levitation diffusion chamber at temperatures between -44 and -30°C and supersaturations ( s i ) between 3 and 29%. These growth data are analyzed with two methods devised to estimate the deposition coefficient ( α) without the direct use of s i . Measurements of s i are typically uncertain, which has called past estimates of α into question. We find that the deposition coefficient ranges from 0.002 to unity and is scattered with temperature, as shown in prior measurements. The data collectively also show a relationship between α and s i , with α rising (falling) with increasing s i for homogeneously (heterogeneously) frozen ice. Analysis of the normalized mass growth rates reveals that heterogeneously-frozen crystals grow near the maximum rate at low s i , but show increasingly inhibited (low α) growth at high s i . Additionally, 7 of the 17 homogeneously frozen crystals cannot be modeled with faceted growth theory or constant α. These cases require the growth mode to transition from efficient tomore »inefficient in time, leading to a large decline in α. Such transitions may be, in part, responsible for the inconsistency in prior measurements of α.« less
  5. Abstract

    Bin and bulk schemes are the two primary methods to parameterize cloud microphysical processes. This study attempts to reveal how their structural differences (size‐resolved vs. moment‐resolved) manifest in terms of cloud and precipitation properties. We use a bulk scheme, the Arbitrary Moment Predictor (AMP), which uses process parameterizations identical to those in a bin scheme but predicts only moments of the size distribution like a bulk scheme. As such, differences between simulations using AMP's bin scheme and simulations using AMP itself must come from their structural differences. In one‐dimensional kinematic simulations, the overall difference between AMP (bulk) and bin schemes is found to be small. Full‐microphysics AMP and bin simulations have similar mean liquid water path (mean percent difference <4%), but AMP simulates significantly lower mean precipitation rate (−35%) than the bin scheme due to slower precipitation onset. Individual processes are also tested. Condensation is represented almost perfectly with AMP, and only small AMP‐bin differences emerge due to nucleation, evaporation, and sedimentation. Collision‐coalescence is the single biggest reason for AMP‐bin divergence. Closer inspection shows that this divergence is primarily a result of autoconversion and not of accretion. In full microphysics simulations, lowering the diameter threshold separating cloud and rain categorymore »in AMP fromtoreduces the largest AMP‐bin difference to ∼10%, making the effect of structural differences between AMP (and perhaps triple‐moment bulk schemes generally) and bin even smaller than the parameterization differences between the two bin schemes.

    « less