Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ice cloud models”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 145 records · Page 8

Primary and secondary ice production: interactions and their relative importance

A discrepancy of up to 5 orders of magnitude between ice crystal and ice nucleating particle (INP) number concentrations was found in the measurements, indicating the potentially important role of secondary ice production (SIP) in the clouds. However, the interactions between primary and SIP processes and their relative importance remain unexplored. In this study, we implemented five different ice nucleation schemes as well as physical representations of SIP processes (i.e., droplet shattering during rain freezing, ice-ice collisional break-up, and rime splintering) in the Community Earth System Model version 2 (CESM2). We ran CESM2 in the single column mode for model comparisons with the DOE Atmospheric Radiation Measurement (ARM) Mixed-Phase Arctic Cloud Experiment (M-PACE) observations. We found that the model experiments with aerosol-aware ice nucleation schemes and SIP processes yield the best simulation results for the M-PACE single-layer mixed-phase clouds. We further investigated the relative importance of ice nucleation and SIP to ice number and cloud phase as well as interactions between ice nucleation and SIP in the M-PACE single-layer mixed-phase clouds. Our results show that SIP contributes 80 % to the total ice formation and transforms ~30 % of pure liquid-phase clouds simulated in the model experiments without considering SIP into mixed-phase clouds. The SIP is not only a result of ice crystals produced from ice nucleation, but also competes with the ice nucleation by reducing the number concentrations of cloud droplets and cloud-borne dust INPs. Conversely, strong ice nucleation also suppresses SIP by glaciating mixed-phase clouds and thereby reducing the amount of precipitation particles (rain and graupel).

54 ENVIRONMENTAL SCIENCES↗

Total Temperature Measurements in Icing Cloud Flows Using a Rearward Facing Probe

This paper reports on temperature and humidity measurements from a series of ice-crystal icing tunnel experiments conducted in June 2018 at the Propulsion Systems Laboratory at the NASA Glenn Research Center. The tests were fundamental in nature and were aimed at investigating the icing processes on a two-dimensional NACA0012 airfoil subjected to artificially generated icing clouds. Prior to the tests on the airfoil, a suite of instruments, including total temperature and humidity probes, were used to characterize the thermodynamic flow and icing cloud conditions of the facility. Two different total temperature probes were used in these tests which included a custom designed rearward facing probe and a commercial self-heating total temperature probe. The rearward facing probe, the main total temperature probe, was designed to reduce and mitigate the contaminating effects of icing and ingestion of ice crystals and water droplets at the probe's inlet. The probe also serves as an air-sample inlet for a light absorption based humidity measurement. The paper includes a section which discusses total temperature and humidity measurement considerations, and another section which provides an analysis of the main probe's performance characteristics. A computational fluid dynamic model of the flow around the probe was also conducted to gain insight into the trajectory of the flow entering the probe inlet. The experiments included a series of tests in which the relative humidity of the facility flow was swept through with increasingly larger values. The data showed that the rearward facing probe can reasonably capture the flow's total temperature and humidity under mild to moderate icing conditions but produces anomalous results under more intense icing conditions. The experimental data was also compared to an in-house developed thermodynamic model which takes into account the interaction of the main flow with the icing cloud. Comparison to the thermodynamic model showed that the rearward facing probe measured the predicted trends.

Agui, Juan H.↗

Attendee Response to Guiding Questions of the September 2014 DOE Workshop on North Slope of Alaska Priorities

We advocate for several modeling approaches suitable to prove understanding of and capability to model aerosol behavior in nucleating droplets and ice and to model regional cloud systems at the mesoscale. The full complement of data that is ideal for particular study types is increasingly well understood; the ASR ISDAC Science Plan provides an excellent example of listing the data required for each study type targeted by that campaign. Whereas discussion here focused on the microscale and mesoscale, other studies should obviously focus on the regional, seasonal and global scale. When considering ice microphysics alone, laser focus on unconstrained model components could be usefully intensified (e.g., collision-coalescence kernels). Similarly, a focus on boundary layer structure would be deeply valuable to a wide range of model types. The expense of observations motivates efficient planning for observational data sets to support specific study approaches; the concept of closure is often useful (e.g., the surface energy budget).

nucleation↗

Semi-Analytic Functions to Calculate the Deposition Coefficients for Ice Crystal Vapor Growth in Bin and Bulk Microphysical Models

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.

54 ENVIRONMENTAL SCIENCES↗

Evidence for Secondary Ice Production in Southern Ocean Maritime Boundary Layer Clouds

Abstract Maritime boundary‐layer clouds over the Southern Ocean (SO) have a large shortwave radiative effect. Yet, climate models have difficulties in representing these clouds and, especially, their phase in this observationally sparse region. This study aims to increase the knowledge of SO cloud phase by presenting in‐situ cloud microphysical observations from the Southern Ocean Clouds, Radiation, Aerosol, Transport Experimental Study (SOCRATES). We investigate the occurrence of ice in summertime marine stratocumulus and cumulus clouds in the temperature range between 6 and −25°C. Our observations show that in ice‐containing clouds, maximum ice number concentrations of up to several hundreds per liter were found. The observed ice crystal concentrations were on average one to two orders of magnitude higher than the simultaneously measured ice nucleating particle (INP) concentrations in the temperature range below −10°C and up to five orders of magnitude higher than estimated INP concentrations in the temperature range above −10°C. These results highlight the importance of secondary ice production (SIP) in SO summertime marine boundary‐layer clouds. Evidence for rime splintering was found in the Hallett‐Mossop (HM) temperature range but the exact SIP mechanism active at lower temperatures remains unclear. Finally, instrument simulators were used to assess simulated co‐located cloud ice concentrations and the role of modeled HM rime‐splintering. We found that CAM6 is deficient in simulating number concentrations across the HM temperature range with little sensitivity to the model HM process, which is inconsistent with the aforementioned observational evidence of highly active SIP processes in SO low‐level clouds.

54 ENVIRONMENTAL SCIENCES↗

Cloud-atmosphere impacts on the central Arctic surface energy budget (Final Report)

This project investigated how clouds and atmosphere-sea ice interactions shape the central Arctic surface energy budget using observations from the MOSAiC expedition and coordinated model experiments. Analyses focused on three themes. First, a new method applied to thermistor string data revealed that snow thermal conductivity is systematically higher than assumed in most models, strongly influencing conductive heat flux and sea ice growth, especially under clear skies. Second, a full annual cycle of surface energy budget observations, combined with detailed cloud microphysical retrievals, demonstrated how cloud regimes regulate radiative, turbulent, and conductive fluxes, driving seasonal contrasts in sea ice energy balance: in winter, radiative forcing elicits strong turbulent and conductive responses, while in summer excess energy is primarily partitioned into surface melt. Third, these observations were used to evaluate and improve regional and global forecast models. While some models captured key Arctic cloud-radiation states, most exhibited persistent deficiencies in representing liquid-containing clouds, snow-on-sea-ice processes, and flux responses to radiative forcing. Together, these results provide unprecedented benchmarks for understanding coupled Arctic processes, improve parameterizations of snow, sea ice, and cloud interactions, and highlight critical pathways by which atmospheric variability drives sea ice change in a rapidly evolving Arctic system. Lastly, the project was a resounding success, resulting in 38 peer-reviewed publications and dozens of presentations, while also supporting an early career scientist.

54 ENVIRONMENTAL SCIENCES↗

The Radiative Effects of Martian Water Ice Clouds on the Local Atmospheric Temperature Profile

Mars Pathfinder made numerous discoveries, one of which was a deep temperature inversion that extended from about 15 km down to 8 km above the surface. It has been suggested by Haberle et al. (1999. J. Geophys. Res. 104, 8957-8974.) that radiative cooling by a water ice cloud may generate such an inversion. Clouds can strongly affect the local air temperature due to their ability to radiate efficiently in the infrared and due to the low air mass of the martian atmosphere, which allows the temperature to change during the relatively short lifetime of a cloud. We utilize a time-dependent microphysical aerosol model coupled to a radiative--convective model to explore the effects water ice clouds have on the local martian temperature profile. We constrain the dust and water vapor abundance using data from the Viking Missions and Mars Pathfinder. Water t ice clouds with visible optical depths of r > 0.1 form readily in these simulations. These clouds alter the local air temperature directly, through infrared cooling, and indirectly, by redistributing atmospheric dust. With this model we are able to reproduce the temperature inversions observed by Mars Pathfinder and Mars Global t Surveyor 2000 Academic Press

Colaprete, Anthony↗

An Improved Bulk Microphysical Scheme for Studying Precipitation Processes: Comparisons with Other Schemes

Cloud microphysical processes play an important role in non-hydrostatic high-resolution simulations. Over the past decade both research and operational numerical weather prediction models have started using more complex cloud microphysical schemes that were originally developed for high-resolution cloud-resolving models. An improved bulk microphysical parameterization (adopted from the Goddard microphysics schemes) has recently implemented into the Weather Research and Forecasting (WRF) model. This bulk microphysical scheme has three different options --- 2ICE (cloud ice & snow), 3ICE-graupel (cloud ice, snow & graupel) and 3ICE-hail (cloud ice, snow & hail). High-resolution model simulations are conducted to examine the impact of microphysical schemes on two different weather events (a midlatitude linear convective system and an Atlantic hurricane). In addition, this bulk microphysical parameterization is compared with WIRF's three other bulk microphysical schemes.

Tao, W. K.↗

Influence of Ice Particle Surface Roughening on the Global Cloud Radiative Effect

Ice clouds influence the climate system by changing the radiation budget and large-scale circulation. Therefore, climate models need to have an accurate representation of ice clouds and their radiative effects. In this paper, new broadband parameterizations for ice cloud bulk scattering properties are developed for severely roughened ice particles. The parameterizations are based on a general habit mixture that includes nine habits (droxtals, hollow/solid columns, plates, solid/hollow bullet rosettes, aggregate of solid columns, and small/large aggregates of plates). The scattering properties for these individual habits incorporate recent advances in light-scattering computations. The influence of ice particle surface roughness on the ice cloud radiative effect is determined through simulations with the Fu-Liou and the GCM version of the Rapid Radiative Transfer Model (RRTMG) codes and the National Center for Atmospheric Research Community Atmosphere Model (CAM, version 5.1). The differences in shortwave (SW) and longwave (LW) radiative effect at both the top of the atmosphere and the surface are determined for smooth and severely roughened ice particles. While the influence of particle roughening on the single-scattering properties is negligible in the LW, the results indicate that ice crystal roughness can change the SW forcing locally by more than 10 W m(exp −2) over a range of effective diameters. The global-averaged SW cloud radiative effect due to ice particle surface roughness is estimated to be roughly 1-2 W m(exp −2). The CAM results indicate that ice particle roughening can result in a large regional SW radiative effect and a small but nonnegligible increase in the global LW cloud radiative effect.

clouds↗

On the Susceptibility of Cold Tropical Cirrus to Ice Nuclei Abundance

Numerical simulations of cirrus formation in the tropical tropopause layer (TTL) are used to evaluate the impact heterogeneous ice nuclei (IN) abundance on cold cloud microphysical properties. The model includes homogeneous and heterogeneous ice nucleation, deposition growthsublimation, and sedimentation. Reanalysis temperature and wind fields are used to force the simulations, with addition of high-frequency wave temperature variability based on in situ measurements. The model results are constrained by recent in situ observations of TTL cirrus and relative humidity, as well as satellite measurements of cirrus bulk properties and occurrence frequencies. Temperature variability driven by high-frequency waves has a dominant influence on TTL cirrus microphysical properties and occurrence frequencies, and inclusion of these waves is required to produce agreement between the simulated and observed cloud properties as well as the abundance of TTL cirrus. With a composition-independent supersaturation threshold for homogeneous freezing of aqueous aerosols, the model produces excessive ice concentrations compared with in situ observations. Inclusion of relatively numerous heterogeneous nuclei ( 100L) in the simulations improves the agreement with observed ice concentrations. However, when IN contribute significantly to TTL cirrus ice nucleation, the occurrence frequency of large supersaturations with respect to ice is less than indicated by in situ measurements. The simulated TTL cirrus extinction statistics agree with observations (within uncertainties) from both in situ and remote-sensing measurements. We find that the sensitivity of TTL cirrus extinction and ice water content statistics to heterogeneous ice nuclei abundance is relatively weak. Likewise, the simulated occurrence frequencies of TTL cirrus are insensitive to ice nuclei abundance, both in terms of cloud frequency height distribution and regional distribution throughout the tropics.

occurences↗

Uncertainty quantification and optimization of precipitating hydrometeor parameters for winter precipitation in a cloud microphysics scheme

The precipitating hydrometeor parameters used in cloud microphysics schemes carry inherent uncertainties. The quantification of these uncertainties, together with parameter optimization, can significantly improve precipitation forecasts. This study investigates the effects of 13 parameters in the Weather Research and Forecasting (WRF) Double-Moment 6-class (WDM6) microphysics scheme, which define the hydrometeor characteristics such as fall velocity–diameter and mass–diameter relationships, as well as the shape parameter of the drop size distribution for precipitating particles such as rain, snow, and graupel on simulated winter precipitation. A comparison between the model's pre-defined parameters and observations from the International Collaborative Experiments for the PyeongChang 2018 Olympic and Paralympic winter games (ICE-POP 2018) field campaign reveals that the fall velocity–diameter relationship for rain, the mass–diameter relationships for snow and graupel, and the shape parameters for all precipitating particles in the WDM6 scheme deviate from the median values observed by the two-dimensional video disdrometer (2DVD). To quantify parameter sensitivities, a perturbed parameter ensemble (PPE) of 256 simulations was conducted within parameter ranges constrained by 2DVD observations for three winter precipitation cases. Bayesian optimization was then applied to identify parameter sets that minimized the root mean square error (RMSE) for each case, achieving reductions of up to 30.2 %. These results demonstrate that ensemble-based uncertainty quantification and parameter optimization can help identify key parameters and provide a pathway to improving precipitation simulation performance. In addition, measurement sites can be strategically selected based on regions that show high sensitivity to variations in hydrometeor characteristic parameters.

Bayesian optimization↗

Diagnosing Aircraft Icing Potential from Satellite Cloud Retrievals

The threat for aircraft icing in clouds is a significant hazard that routinely impacts aviation operations. Accurate diagnoses and forecasts of aircraft icing conditions requires identifying the location and vertical distribution of clouds with super-cooled liquid water (SLW) droplets, as well as the characteristics of the droplet size distribution. Traditional forecasting methods rely on guidance from numerical models and conventional observations, neither of which currently resolve cloud properties adequately on the optimal scales needed for aviation. Satellite imagers provide measurements over large areas with high spatial resolution that can be interpreted to identify the locations and characteristics of clouds, including features associated with adverse weather and storms. This paper describes new techniques for interpreting cloud products derived from satellite data to infer the flight icing threat to aircraft. For unobscured low clouds, the icing threat is determined using empirical relationships developed from correlations between satellite imager retrievals of liquid water path and droplet size with icing conditions reported by pilots (PIREPS). For deep ice over water cloud systems, ice and liquid water content (IWC and LWC) profiles are derived by using the imager cloud properties to constrain climatological information on cloud vertical structure and water phase obtained apriori from radar and lidar observations, and from cloud model analyses. Retrievals of the SLW content embedded within overlapping clouds are mapped to the icing threat using guidance from an airfoil modeling study. Compared to PIREPS and ground-based icing remote sensing datasets, the satellite icing detection and intensity accuracies are approximately 90% and 70%, respectively, and found to be similar for both low level and deep ice over water cloud systems. The satellite-derived icing boundaries capture the reported altitudes over 90% of the time. Satellite analyses corresponding to the time and location of several recent aviation accidents and with icing PIREPS are also presented that reveal skill in identifying severe icing conditions. These results demonstrate the utility of satellite cloud retrievals for quantitatively diagnosing the potential for icing conditions on temporal and spatial scales that should be useful to the aviation community. Plans are being developed to deliver these new satellite products to the GOES-R Proving Ground in the near future so that they can be evaluated in operational applications.

Smith, William L., Jr.↗

GCSS Cirrus Parcel Model Comparison Project

The Cirrus Parcel Model Comparison Project, a project of GCSS Working Group on Cirrus Cloud Systems (WG2), involves the systematic comparison of current models of ice crystal nucleation and growth for specified, typical, cirrus cloud environments. The goal of this project is to document and understand the factors resulting in significant inter-model differences. The intent is to foment research leading to model improvement and validation. In Phase 1 of the project reported here, simulated cirrus cloud microphysical properties are compared for situations of "warm" (-40 C) and "cold" (-60 C) cirrus subject to updrafts of 4, 20 and 100 cm/s, respectively. Five models participated. These models employ explicit microphysical schemes wherein the size distribution of each class of particles (aerosols and ice crystals) is resolved into bins. Simulations are made including both homogeneous and heterogeneous ice nucleation mechanisms. A single initial aerosol population of sulfuric acid particles is prescribed for all simulations. To isolate the treatment of the homogeneous freezing (of haze drops) nucleation process, the heterogeneous nucleation mechanism is disabled for a second parallel set of simulations. Qualitative agreement is found for the homogeneous-nucleation-only simulations, e.g., the number density of nucleated ice crystals increases with the strength of the prescribed updraft. However, non-negligible quantitative differences are found. Detailed analysis reveals that the homogeneous nucleation formulation, aerosol size, ice crystal growth rate (particularly the deposition coefficient), and water vapor uptake rate are critical components that lead to differences in predicted microphysics. Systematic bias exists between results based on a modified classical theory approach and models using an effective freezing temperature approach to the treatment of nucleation. Each approach is constrained by critical freezing data from laboratory studies, but each includes assumptions that can only be justified by further laboratory data. Consequently, it is not yet clear if the two approaches can be made consistent. Large haze particles may deviate considerably from equilibrium size in moderate to strong updrafts (20-100 cm/s) at -60 C when the commonly invoked equilibrium assumption is lifted. The resulting difference in particle-size-dependent solution concentration of haze particles may significantly affect the ice nucleation rate during the initial nucleation interval. The uptake rate for water vapor excess by ice crystals is another key component regulating the total number of nucleated ice crystals. This rate, the product of ice number concentration and ice crystal diffusional growth rate, which is sensitive to the deposition coefficient when ice particles are small, partially controls the peak nucleation rate achieved in an air parcel and the duration of the active nucleation time period. The effects of heterogeneous nucleation are most pronounced in weak updraft situations. Vapor competition by the nucleated (heterogeneous) ice crystals limits the achieved ice supersaturation and thus suppresses the contribution of homogeneous nucleation. Correspondingly, ice crystal number density is markedly reduced. Definitive laboratory and atmospheric benchmark data are needed for the heterogeneous nucleation process. Inter-model differences are correspondingly greater than in the case of the homogeneous nucleation process acting alone.

Lin, Ruei-Fong↗

An Improved Convection Parameterization with Detailed Aerosol–Cloud Microphysics for a Global Model

Abstract A new microphysical treatment that includes aerosol–cloud interactions and secondary ice production (SIP) mechanisms is implemented in the convection scheme of the Community Atmosphere Model, version 6 (CAM6). The approach is to embed a 1D Lagrangian parcel model in the bulk convective plume of the existing deep convection parameterization. Aerosol activation, growth processes including collision/coalescence, and three processes of SIP mechanisms, two of which are normally overlooked in atmospheric models, are represented in this embedded parcel model. These microphysical processes are treated with a hybrid bin/bulk scheme and a high spatial and temporal resolution for the integration of the embedded parcel in 1D, allowing vertical velocity to determine the microphysical evolution following the in-cloud motion during ascent. Simulations of an observed case (Midlatitude Continental Convective Clouds Experiment) of a mesoscale convective system in Oklahoma, United States, with a single-column model (SCAM) version of CAM, are compared with aircraft in situ and ground-based observations of microphysical properties from the convection and precipitation. Results from the validation show the new microphysical scheme has a good representation of the ice initiation in the bulk convective plume, including the known and empirically quantified pathways of primary and secondary initiation, with benefits for the accuracy of properties of its supercooled cloud liquid. The sensitivity simulations and use of tagging tracers for the validated simulation confirm that the newly included SIP mechanisms are of paramount importance for convective microphysics and can be successfully treated in the global model.

54 ENVIRONMENTAL SCIENCES↗

Preparation for Scaling Studies of Ice-Crystal Icing at the NRC Research Altitude Test Facility

This paper describes experiments conducted at the National Research Council (NRC) of Canadas Research Altitiude Test Facility between March 26 and April 11, 2012. The tests, conducted collaboratively between NASA and NRC, focus on three key aspects in preparation for later scaling work to be conducted with a NACA 0012 airfoil model in the NRC Cascade rig: (1) cloud characterization, (2) scaling model development, and (3) ice-shape profile measurements. Regarding cloud characterization, the experiments focus on particle spectra measurements using two shadowgraphy methods, cloud uniformity via particle scattering from a laser sheet, and characterization of the SEA Multi-Element probe. Overviews of each aspect as well as detailed information on the diagnostic method are presented. Select results from the measurements and interpretation are presented which will help guide future work.

Ice formation↗

The impact of cloud microphysics and ice nucleation on Southern Ocean clouds assessed with single-column modeling and instrument simulators

Abstract. Supercooled liquid clouds are common at higher latitudes (especially over the Southern Ocean) and are critical for constraining climate projections. We take advantage of the Macquarie Island Cloud and Radiation Experiment (MICRE) to perform an analysis of observed and simulated cloud processes over the Southern Ocean in a region and season dominated by supercooled liquid clouds. Using a single-column version of the European Centre for Medium-Range Weather Forecasts (ECMWF) Integrated Forecast System (IFS), we compare two different cloud microphysical schemes to ground-based observations of cloud, precipitation, and radiation over a 2.5-month period (1 January–17 March 2017). Both schemes are able to reproduce aspects of the cloud and radiation observations during MICRE to within the uncertainty of the data when the thermodynamic profile is prescribed with relaxation. There are differences in water mass and representation of reflectivity between the schemes. A sensitivity study of the cloud microphysics schemes, one a bulk one-moment scheme and the other a two-moment scheme with prediction of mass and number, indicates that several key processes create differences between the schemes. Surface radiative fluxes and total water path are highly sensitive to the formation and fall speed of precipitation. The prediction of hydrometeor number with the two-moment scheme yields a better comparison with observed reflectivity and radiative fluxes, despite predicting higher liquid water contents than observed. With the two-moment scheme, we are also able to test the sensitivity of the results to the input of liquid cloud condensation nuclei (CCN) and ice nuclei (IN). The cloud properties and resulting radiative effects are found to be sensitive to the CCN and IN concentrations. More CCN and IN increase liquid and ice water paths, respectively. Thus, both the dynamic environment and aerosols, integrated through the cloud microphysics, are important for properly representing Southern Ocean cloud radiative effects.

58 GEOSCIENCES↗

Cirrus Parcel Model Comparison Project: The Critical Components to Simulate Cirrus Initiation Explicitly - Phase 1

The Cirrus Parcel Model Comparison Project, a project of the GCSS (GEWEX Cloud System Studies) Working Group on Cirrus Cloud Systems, involves the systematic comparison of current models of ice crystal nucleation and growth for specified, typical, cirrus cloud environments. In Phase I of the project reported here, simulated cirrus cloud microphysical properties are compared for situations of "warm" (40 C) and "cold" (-60 C) cirrus, both subject to updrafts of 4, 20 and 100 centimeters per second. Five models participated. The various models employ explicit microphysical schemes wherein the size distribution of each class of particles (aerosols and ice crystals) is resolved into bins or treated separately. Simulations are made including both the homogeneous and heterogeneous ice nucleation mechanisms. A single initial aerosol population of sulfuric acid particles is prescribed for all simulations. To isolate the treatment of the homogeneous freezing (of haze droplets) nucleation process, the heterogeneous nucleation mechanism is disabled for a second parallel set of simulations. Qualitative agreement is found for the homogeneous-nucleation- only simulations, e.g., the number density of nucleated ice crystals increases with the strength of the prescribed updraft. However, significant quantitative differences are found. Detailed analysis reveals that the homogeneous nucleation rate, haze particle solution concentration, and water vapor uptake rate by ice crystal growth (particularly as controlled by the deposition coefficient) are critical components that lead to differences in predicted microphysics. Systematic bias exists between results based on a modified classical theory approach and models using an effective freezing temperature approach to the treatment of nucleation. Each approach is constrained by critical freezing data from laboratory studies, but each includes assumptions that can only be justified by further laboratory research. Consequently, it is not yet clear if the two approaches can be made consistent. Large haze particles may deviate considerably from equilibrium size in moderate to strong updrafts (20-100 centimeters per second) at -60 C when the commonly invoked equilibrium assumption is lifted. The resulting difference in particle-size- dependent solution concentration of haze particles may significantly affect the ice particle formation rate during the initial nucleation interval. The uptake rate for water vapor excess by ice crystals is another key component regulating the total number of nucleated ice crystals. This rate, the product of particle number concentration and ice crystal diffusional growth rate, which is particularly sensitive to the deposition coefficient when ice particles are small, modulates the peak particle formation rate achieved in an air parcel and the duration of the active nucleation time period. The effects of heterogeneous nucleation are most pronounced in weak updraft situations. Vapor competition by the heterogeneously nucleated ice crystals may limit the achieved ice supersaturation and thus suppresses the contribution of homogeneous nucleation. Correspondingly, ice crystal number density is markedly reduced. Definitive laboratory and atmospheric benchmark data are needed for the heterogeneous nucleation process. Inter-model differences are correspondingly greater than in the case of the homogeneous nucleation process acting alone.

Lin, Ruei-Fong↗