Search NASA⌕ Search

SEARCH · Search NASA

Results for “cloud phase”

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 379 records · Page 21

Above-Cloud Precipitable Water Retrievals using the MODIS 0.94 micron Band with Applications for Multi-Layer Cloud Detection

In order to perform satellite retrievals of cloud properties, it is important to account for the effect of the above-cloud atmosphere on the observations. The solar bands used in the operational MODIS Terra and Aqua cloud optical and microphysical algorithms (visible, NIR, and SWIR spectral windows) are primarily affected by water vapor, and to a lesser extent by well-mixed gases. For water vapor, the above-cloud column amount, or precipitable water, provides adequate information for an atmospheric correction; details of the vertical vapor distribution are not typically necessary for the level of correction required. Cloud-top pressure has a secondary effect due to pressure broadening influences. For well- mixed gases, cloud-top pressure is also required for estimates of above-cloud abundances. We present a method for obtaining above-cloud precipitable water over dark Ocean surfaces using the MODIS 0.94 pm vapor absorption band. The retrieval includes an iterative procedure for establishing cloud-top temperature and pressure, and is useful for both single layer water and ice clouds. Knowledge of cloud thermodynamic phase is fundamental in retrieving cloud optical and microphysical properties. However, in cases of optically thin cirrus overlapping lower water clouds, the concept of a single unique phase is ill- defined and depends, at least, on the spectral region of interest. We will present a method for multi-layer and multi-phase cloud detection which uses above-cloud precipitable water retrievals along with several existing MODIS operational cloud products (cloud-top pressure derived from a C02 slicing algorithm, IR and SWIR phase retrievals). Results are catagorized by whether the radiative signature in the MODIS solar bands is primarily that of a water cloud with ice cloud contamination, or visa-versa. Examples in polar and mid-latitude regions will be shown.

Platnick, S.↗

Improvement of aerosol activation/ice nucleation in a source-oriented WRF-Chem model to study a winter Storm in California

The source-oriented Weather Research and Forecasting chemistry (SOWC) model is modified to improve mixed-phase cloud parameterization by including aerosol aging processes (external mixture) in cloud activation/nucleation. The improved SOWC model is applied to investigate how aerosol mixing state influences cloud and ice formation and atmospheric optical properties during a winter storm. SOWC tracks 6-dimensional chemical variables (X, Z, Y, Size Bins, Source Types, Species) through an explicit simulation of atmospheric chemistry and physics. A new module is implemented into the enhanced SOWC model to simulate microphysics processes with all source-oriented hydrometeors (cloud, ice, rain, snow and graupel) by using the modified Morrison two-moment microphysics scheme. In this study, all aerosol source types can activate to form cloud droplets based on the Köhler theory, but the dust source type is the only source of ice nuclei (IN). Here, the Goddard shortwave and longwave radiation schemes are modified to use a new geometric optics method to estimate cloud optical properties by considering the chemistry components and the physical shape of each hydrometeor.

54 ENVIRONMENTAL SCIENCES↗

Contact Freezing of Water Droplets by Crystalline Organic Acids

The ability of water to freeze into ice crystals in mixed-phase clouds affects physical properties, including particle size, precipitation rates, and radiative properties. The presence of an insoluble particle at the surface of water droplets can promote ice nucleation at temperatures higher than that of pure water, even in the absence of a collision. However, contact freezing remains an underexplored mode of ice nucleation. Here, we present a study of atmospherically relevant organic acids and their role as effective ice-nucleating particles (INP) in contact mode using a Raman-microscope-equipped environmental chamber. We determined contact freezing temperatures induced by solid crystals of docosanol, adipic acid, cis-pinonic acid, fumaric acid, 4-hydroxybenzoic acid, palmitic acid, phthalic acid, sebacic acid, stearic acid, terephthalic acid, and vanillic acid. All solids except fumaric acid promoted contact freezing of water droplets at significantly higher temperatures than pure water in the chamber (−15.0 to −18.5 °C vs −21.3 °C). Physical and chemical properties were identified which correlate with greater effectiveness of INPs in the contact mode, including crystal lattice mismatch with ice, carbon number, and insolubility in water. In conclusion, we suggest that the presence of these organic solids in atmospheric aerosols may promote atmospheric ice nucleation at warm temperatures.

Clouds↗

On the Forward Modeling of Radar Doppler Spectrum Width From LES: Implications for Model Evaluation

Large-eddy simulations of an observed single-layer Arctic mixed-phase cloud are analyzed to study the value of forward modeling of profiling millimeter wave cloud radar Doppler spectral width for model evaluation. Individual broadening terms and their uncertainties are quantified for the observed spectral width and compared to modeled broadening terms. Modeled turbulent broadening is narrower than the observed values when the turbulent kinetic energy dissipation rate from the subgrid scale model is used in the forward model. The total dissipation rates, estimated with the subgrid scale dissipation rates and the numerical dissipation rates, agree much better with both the retrieved dissipation rates and those inferred from the power spectra of the simulated vertical air velocity. The comparison of the microphysical broadening provides another evaluative measure of the ice properties in the simulation. To accurately retrieve dissipation rates as well as each broadening term from the observations, we suggest a few modifications to previously presented techniques. First, we show that the inertial subrange spectrum filtered with the radar sampling volume is a better underlying model than the unfiltered -5/3 law for the retrieval of the dissipation rate from the power spectra of the mean Doppler velocity. Second, we demonstrate that it is important to filter out turbulence and remove the layer-mean reflectivity-weighted mean fall speed from the observed mean Doppler velocity to avoid overestimation of shear broadening. Finally, we provide a method to quantify the uncertainty in the retrieved dissipation rates, which eventually propagates to the uncertainty in the microphysical broadening.

54 ENVIRONMENTAL SCIENCES↗

Coupling Warm Rain With an Eddy Diffusivity/Mass Flux Parameterization: 1. Model Description and Validation

A new version of the stochastic multiplume Jet Propulsion Laboratory Eddy-Diffusivity/Mass-Flux (JPL-EDMF) parameterization which consistently couples the simplified Khairoutdinov and Kogan (2000), https://doi.org/10.1175/1520-0493(2000)128<0229:ANCPPI>2.0.CO;2, warm phase cloud microphysical parameterization with the parameterization of cloud macrophysical and subgrid scale dynamical processes is described. The new parameterization combines the EDMF approach with an assumed shape of a joint probability density function of thermodynamic and kinematic variables which provide the basis for the computation of all parameterized processes. As far as we are aware this is the first attempt to consistently couple all of these parameterized processes in the EDMF framework. This paper is part one of a two paper series. Here, the JPL-EDMF parameterization is described and benchmark simulations of precipitating stratocumulus and cumulus convection are performed in a single-column-model framework. The parameterization results compare favorably to the reference large-eddy-simulation results. In the second part (Smalley et al., 2022, https://doi.org/10.1029/2021MS002729) the JPL-EDMF parameterization is validated for a wide range of observation-based scenarios covering the continuous transition from subtropical stratocumulus to cumulus convection derived from global reanalysis, and parameterization uncertainties are studied in detail.

54 ENVIRONMENTAL SCIENCES↗

Simulating Southern Ocean Aerosol and Ice Nucleating Particles in the Community Earth System Model Version 2

Southern Ocean (SO) low-level mixed phase clouds have been a long-standing challenge for Earth system models to accurately represent. While improvements to the Community Earth System Model version 2 (CESM2) resulted in increased supercooled liquid in SO clouds and improved model radiative biases, simulated SO clouds in CESM2 now contain too little ice. Previous observational studies have indicated that marine particles are major contributor to SO low-level cloud heterogeneous ice nucleation, a process that initiates a number of cloud processes that govern cloud radiative properties. In this study, we utilize detailed aerosol and ice nucleating particle (INP) measurements from two recent measurement campaigns to assess simulated aerosol abundance, number size distributions, and composition and INP parameterizations for use in CESM2. Our results indicate that CESM2 has a positive bias in simulated surface-level total aerosol surface area at latitudes north of 58°S. Measured INP populations were dominated by marine INPs and we present evidence of refractory INPs present over the SO assumed here to be mineral dust INPs. Results highlight a critical need to assess simulated mineral dust number and size distributions in CESM2 in order to adequately represent SO INP populations and their response to long-term changes in atmospheric transport patterns and land use change. Furthermore, we also discuss important cautions and limitations in applying a commonly used mineral dust INP parameterization to remote regions like the pristine SO.

58 GEOSCIENCES↗

Preserving Tracer Correlations in Moment-Based Atmospheric Transport Models

A linear non-diffusive algorithm for advective transport is developed that greatly improves the detail at which aerosols and clouds can be represented in atmospheric models. Linear advection schemes preserve tracer correlations but the most basic linear scheme is rarely used by atmospheric modelers on account of its excessive numerical diffusion. Higher-order schemes are in widespread use, but these present new problems as nonlinear adjustments are required to avoid occurrences of negative concentrations, spurious oscillations, and other non-physical effects. Generally successful at reducing numerical diffusion during the advection of individual tracers, for example, particle number or mass, the higher-order schemes fail to preserve even the simplest of correlations between interrelated tracers. As a result, important attributes of aerosol and cloud populations including radial moments of particle size distributions, molecular precursors related through chemical equilibria, aerosol mixing state, and distribution of cloud phase are poorly represented. We introduce a new transport scheme, minVAR, that is both non-diffusive and preservative of tracer correlations, thereby combining the best features of the basic and higher-order schemes while enabling new features such as the tracking of sub-grid information at arbitrarily fine scales with high computational efficiency.

54 ENVIRONMENTAL SCIENCES↗

Sea spray aerosol as a unique source of ice nucleating particles

Ice nucleating particles (INPs) are vital for ice initiation in, and precipitation from, mixed-phase clouds. A source of INPs from oceans within sea spray aerosol (SSA) emissions has been suggested in previous studies but remained unconfirmed. Here, we show that INPs are emitted using real wave breaking in a laboratory flume to produce SSA. The number concentrations of INPs from laboratory-generated SSA, when normalized to typical total aerosol number concentrations in the marine boundary layer, agree well with measurements from diverse regions over the oceans. Data in the present study are also in accord with previously published INP measurements made over remote ocean regions. INP number concentrations active within liquid water droplets increase exponentially in number with a decrease in temperature below 0 °C, averaging an order of magnitude increase per 5 °C interval. The plausibility of a strong increase in SSA INP emissions in association with phytoplankton blooms is also shown in laboratory simulations. Nevertheless, INP number concentrations, or active site densities approximated using “dry” geometric SSA surface areas, are a few orders of magnitude lower than corresponding concentrations or site densities in the surface boundary layer over continental regions. These findings have important implications for cloud radiative forcing and precipitation within low-level and midlevel marine clouds unaffected by continental INP sources, such as may occur over the Southern Ocean.

54 ENVIRONMENTAL SCIENCES↗

Simulated Polarimetric Fields of Ice Vapor Growth Using the Adaptive Habit Model. Part I: Large-Eddy Simulations

Abstract The bulk adaptive habit model (AHM) explicitly predicts ice particle aspect ratio, improving the representation of microphysical processes and properties, including ice–liquid-phase partitioning. With the unique ability to predict ice particle shape and density, the AHM is combined with an offline forward operator to produce fields of simulated polarimetric variables. An evaluation of AHM-forward-simulated dual-polarization radar signatures in an idealized Arctic mixed-phase cloud is presented. Interpretations of those signatures are provided through microphysical model output using the large-eddy simulation mode of the Weather Research and Forecasting Model. Vapor-grown ice properties are associated with distinct observable signatures in polarimetric radar variables, with clear sensitivities to the simulated ice particle properties, including ice number, size, and distribution shape. In contrast, the liquid droplet number has little influence on both polarimetric and microphysical variables in the case presented herein. Polarimetric quantities are sensitive to the dominating crystal habit type in a volume, with enhancements for aspect ratios much lower or higher than unity. This synthesis of a microphysical model and a polarimetric forward simulator is a first step in the evaluation of detailed AHM microphysics.

54 ENVIRONMENTAL SCIENCES↗

Ice nucleation ability of loess from the northwestern United States

The heterogeneous nucleation of ice processes involving loess particles that influences the formation of mixed-phase clouds are poorly understood. Here, the ice nucleating ability of wind-blown dust or loess accumulated from the past glaciated area was investigated at three temperatures: -26, -30, and -34 °C and at below and above saturation with respect to liquid water conditions. Total six loess samples from different regions across Columbia Basin province, WA, USA were collected, dry dispersed, size-selected at mobility diameter 200 nm, and investigated for their ice nucleation efficiency. To understand the effect of atmospheric processing during long-range transport on their ice nucleating ability, similar experiments were also performed on acid-treated loess samples. Additionally, the ice nucleating properties of Arizona Test Dust (ATD) were investigated as a surrogate for natural mineral dust particles to test the experimental approach. Results show that treated particles have lower ice nucleation efficiency compared to untreated particles at all temperature and saturation with respect to liquid water conditions. Comparison based on ice-active site density ( N s ) metric indicate that loess particles at saturation with respect to liquid water conditions are marginally more efficient than the mineral and soil dust values reported in the literature, but they have lower efficiencies than the predicted N s efficiency of K-feldspar particles at supercooled temperatures greater than -38 °C.

54 ENVIRONMENTAL SCIENCES↗

Development of Cloud and Precipitation Property Retrieval Algorithms and Measurement Simulators from ASR Data

What has made the ASR program unique is the amount of information that is available. The suite of recently deployed instruments significantly expands the scope of the program (Mather and Voyles, 2013). The breadth of this information allows us to pose sophisticated process-level questions. Our ASR project, now entering its third year, has been about developing algorithms that use this information in ways that fully exploit the new capacity of the ARM data streams. Using optimal estimation (OE) and Markov Chain Monte Carlo (MCMC) inversion techniques, we have developed methodologies that allow us to use multiple radar frequency Doppler spectra along with lidar and passive constraints where data streams can be added or subtracted efficiently and algorithms can be reformulated for various combinations of hydrometeors by exchanging sets of empirical coefficients. These methodologies have been applied to boundary layer clouds, mixed phase snow cloud systems, and cirrus.

54 ENVIRONMENTAL SCIENCES↗

Marine Ice Nuclei Collections – MAGIC (MAGIC-IN) Final Campaign Summary

This campaign augmented measurements obtained via deployment of the Atmospheric Radiation Measurement (ARM) Climate Research Facility’s ARM Mobile Facility (AMF) in the Marine ARM GPCI1 Investigation of Clouds (MAGIC) field campaign. The measurements, comprised of shipboard aerosol collections obtained during the five legs of the summer 2013 cruises, were sent for offline processing to measure ice nucleating particle (INP) number concentrations. The forty-three sample periods each represented, nominally, 24-hour segments during outbound and inbound transits of the Horizon Spirit. The samples were collected at locations between Los Angeles and Hawaii. Eight samples have been analyzed for immersion freezing temperature spectra thus far, using funding from other grants. Remaining samples are being frozen until support for further processing is obtained. Future analyses will investigate the inorganic/organic proportions of ice nuclei, in addition to determining the genetic composition of the overall biological community associated with INPs. Resulting correlations will be compared with other archived aerosol quantities, meteorological and ocean data (e.g., temperature, wind speed, sea surface temperature, etc…) and satellite ocean color products. These findings will ultimately aid in parameterizing oceanic (e.g., sea spray) INP emissions in regional and global scale models, when illustrating aerosol connections to cloud phases and properties. Independent future analyses of frozen filter samples, as proposed by collaborating investigators at the time of this report, will include single particle analyses of marine boundary layer aerosol compositions and morphology. The MAGIC-IN data are considered representative of the oligotrophic, low Chlorophyll-a (with the exception of near-shore) ocean regions, which exist along the MAGIC transect. Current analyses suggest that INP numbers in the marine boundary layer over this region are typically low, compared to existing measurements over marine areas and those collected in the laboratory as the result of realistic sea spray particle generation. These findings, along with separate studies, confirm the existence of highly variable emission sources for INP from oceans, (though weaker than land-based emissions at modestly cooled temperatures).

54 ENVIRONMENTAL SCIENCES↗

Examining the Ice-Nucleating Particles from the Southern Great Plains Field Campaign Report

The recent U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) user facility field campaign, named Examining the Ice-Nucleating Particles from Southern Great Plains (ExINP-SGP), targeted experimentally characterizing abundance and other properties of ambient ice-nucleating particles (INPs) at the SGP site in Oklahoma (36° 36' 18" N, 97° 29' 6" W) during the period of October 1 to November 14, 2019. This campaign was funded through the DOE Office of Science Early Career Research Program (DE-SC0018979) as one of three major campaigns to be conducted at three ARM observatories, including the Eastern North Atlantic (ENA) station in 2020–2021 and North Slope of Alaska (NSA) site in 2021–2023, besides SGP. Different INP episodes were assessed to develop an SGP-relevant ice nucleation parameterization that would help understand convective and mixed-phase cloud systems typically observed in this region.

54 ENVIRONMENTAL SCIENCES↗

Ice processes in Antarctica: identification via multi-wavelength active and passive measurements and model evaluation (Final Report)

This project is based on new multi-frequency radar observations collected during the AWARE field campaign in 2016 at McMurdo Station by the ARM Mobile facility. Three scientific goals have been achieved. We have developed ice and mixed-phase cloud microphysics retrievals from multi-wavelength radar observations and ancillary observations, generally applicable at ARM sites. These technique generally better characterize the size of the ice crystals, their equivalent water content and the degree of riming better than using single wavelength techniques. We have characterized ice microphysical properties and ice processes via multi-frequency signatures. We have identified a unique mode in the McMurdo dataset appearing at quite cold temperature that we have attributed to riming processes. Triple frequency radar measurements unfortunately were available only for a short period but our work lays the foundation for building climatological and testing their replicability in models. Simultaneously case studies from the AWARE campaign have been exploited to evaluate and improve the NASA GISS ModelE3 climate model, including selection of case studies for combined goals of (i) running side-by-side using large-eddy simulation (LES) and ModelE3 in single-column model (SCM) mode and (ii) fingerprinting cloud processes in multi-wavelength radar observations. The different components of this work have contributed to 15 published papers with two currently under review.

54 ENVIRONMENTAL SCIENCES↗

Analysis of isothermal and cooling-rate-dependent immersion freezing by a unifying stochastic ice nucleation model

Immersion freezing is an important ice nucleation pathway involved in the formation of cirrus and mixed-phase clouds. Laboratory immersion freezing experiments are necessary to determine the range in temperature, T , and relative humidity, RH, at which ice nucleation occurs and to quantify the associated nucleation kinetics. Typically, isothermal (applying a constant temperature) and cooling-rate-dependent immersion freezing experiments are conducted. In these experiments it is usually assumed that the droplets containing ice nucleating particles (INPs) all have the same INP surface area (ISA); however, the validity of this assumption or the impact it may have on analysis and interpretation of the experimental data is rarely questioned. Descriptions of ice active sites and variability of contact angles have been successfully formulated to describe ice nucleation experimental data in previous research; however, we consider the ability of a stochastic freezing model founded on classical nucleation theory to reproduce previous results and to explain experimental uncertainties and data scatter. A stochastic immersion freezing model based on first principles of statistics is presented, which accounts for variable ISA per droplet and uses parameters including the total number of droplets, N tot , and the heterogeneous ice nucleation rate coefficient, J het ( T ). This model is applied to address if (i) a time and ISA-dependent stochastic immersion freezing process can explain laboratory immersion freezing data for different experimental methods and (ii) the assumption that all droplets contain identical ISA is a valid conjecture with subsequent consequences for analysis and interpretation of immersion freezing. The simple stochastic model can reproduce the observed time and surface area dependence in immersion freezing experiments for a variety of methods such as: droplets on a cold-stage exposed to air or surrounded by an oil matrix, wind and acoustically levitated droplets, droplets in a continuous-flow diffusion chamber (CFDC), the Leipzig aerosol cloud interaction simulator (LACIS), and the aerosol interaction and dynamics in the atmosphere (AIDA) cloud chamber. Observed time-dependent isothermal frozen fractions exhibiting non-exponential behavior can be readily explained by this model considering varying ISA. An apparent cooling-rate dependence of J het is explained by assuming identical ISA in each droplet. When accounting for ISA variability, the cooling-rate dependence of ice nucleation kinetics vanishes as expected from classical nucleation theory. Finally, the model simulations allow for a quantitative experimental uncertainty analysis for parameters N tot , T , RH, and the ISA variability. We discuss the implications of our results for experimental analysis and interpretation of the immersion freezing process.

54 ENVIRONMENTAL SCIENCES↗

Physicochemical characterization and source apportionment of Arctic ice-nucleating particles observed in Ny-Ålesund in autumn 2019

Abstract. Ice-nucleating particles (INPs) initiate primary ice formation in Arctic mixed-phase clouds (MPCs), altering cloud radiative properties and modulating precipitation. For atmospheric INPs, the complexity of their spatiotemporal variations, heterogeneous sources, and evolution via intricate atmospheric interactions challenge the understanding of their impact on microphysical processes in Arctic MPCs and induce an uncertain representation in climate models. In this work, we performed a comprehensive analysis of atmospheric aerosols at the Arctic coastal site in Ny-Ålesund (Svalbard, Norway) from October to November 2019, including their ice nucleation ability, physicochemical properties, and potential sources. Overall, INP concentrations (NINP) during the observation season were approximately up to 3 orders of magnitude lower compared to the global average, with several samples showing degradation of NINP after heat treatment, implying the presence of proteinaceous INPs. Particle fluorescence was substantially associated with INP concentrations at warmer ice nucleation temperatures, indicating that in the far-reaching Arctic, aerosols of biogenic origin throughout the snow- and ice-free season may serve as important INP sources. In addition, case studies revealed the links between elevated NINP and heat lability, fluorescence, high wind speeds originating from the ocean, augmented concentration of coarse-mode particles, and abundant organics. Backward trajectory analysis demonstrated a potential connection between high-latitude dust sources and high INP concentrations, while prolonged air mass history over the ice pack was identified for most scant INP cases. The combination of the above analyses demonstrates that the abundance, physicochemical properties, and potential sources of INPs in the Arctic are highly variable despite its remote location.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Assessing the cloud radiative bias at Macquarie Island in the ACCESS-AM2 model

Abstract. As a long-standing problem in climate models, large positive shortwave radiation biases exist at the surface over the Southern Ocean, impacting the accurate simulation of sea surface temperature, atmospheric circulation, and precipitation. Underestimations of low-level cloud fraction and liquid water content are suggested to predominantly contribute to these radiation biases. Most model evaluations for radiation focus on summer and rely on satellite products, which have their own limitations. In this work, we use surface-based observations at Macquarie Island to provide the first long-term, seasonal evaluation of both downwelling surface shortwave and longwave radiation in the Australian Community Climate and Earth System Simulator Atmosphere-only Model version 2 (ACCESS-AM2) over the Southern Ocean. The capacity of the Clouds and the Earth’s Radiant Energy System (CERES) product to simulate radiation is also investigated. We utilize the novel lidar simulator, the Automatic Lidar and Ceilometer Framework (ALCF), and all-sky cloud camera observations of cloud fraction to investigate how radiation biases are influenced by cloud properties. Overall, we find an overestimation of +9.5±33.5 W m−2 for downwelling surface shortwave radiation fluxes and an underestimation of -2.3±13.5 W m−2 for downwelling surface longwave radiation in ACCESS-AM2 in all-sky conditions, with more pronounced shortwave biases of +25.0±48.0 W m−2 occurring in summer. CERES presents an overestimation of +8.0±18.0 W m−2 for the shortwave and an underestimation of -12.1±12.2 W m−2 for the longwave in all-sky conditions. For the cloud radiative effect (CRE) biases, there is an overestimation of +4.8±28.0 W m−2 in ACCESS-AM2 and an underestimation of -7.9±20.9 W m−2 in CERES. An overestimation of downwelling surface shortwave radiation is associated with an underestimated cloud fraction and low-level cloud occurrence. We suggest that modeled cloud phase is also having an impact on the radiation biases. Our results show that the ACCESS-AM2 model and CERES product require further development to reduce these radiation biases not just in shortwave and in all-sky conditions, but also in longwave and in clear-sky conditions.

54 ENVIRONMENTAL SCIENCES↗

The impact of aerosol mixing state on immersion freezing: insights from classical nucleation theory and particle-resolved simulations

Immersion freezing, initiated by ice-nucleating particles (INPs) in supercooled aqueous droplets, plays an important role in the formation of ice crystals within clouds. The efficiency of immersion freezing depends strongly on INP composition and, crucially, on the mixing state – how chemical species are distributed across the particle population. Here, we quantify the impact of aerosol mixing state on immersion freezing using a combined theoretical and particle-resolved modeling approach. We derive analytical expressions for the frozen fraction of internally and externally mixed INP populations based on classical nucleation theory, showing that the frozen fraction is sensitive to whether ice-active species are present in all particles or only in a subset of the population. We introduce a multi-species immersion freezing scheme into the particle-resolved model PartMC, using the water activity-based immersion freezing model (ABIFM) to compute freezing probabilities for mixed-composition particles. To improve computational efficiency, we implement a Binned Tau-Leaping algorithm and demonstrate an order-of-magnitude speedup with minimal accuracy loss. Simulations reproduce the analytical trends in limiting cases and extend the analysis to more general aerosol populations, where mixing state continues to exert a substantial control on frozen fraction. Sensitivity analyses across particle size, species type, and cooling condition reveal that the mixing state effect is most pronounced when small amounts of highly efficient INPs are mixed with less efficient materials. These findings underscore the need to represent aerosol mixing state explicitly in models of heterogeneous ice nucleation to reduce uncertainty in cloud-phase partitioning.

54 ENVIRONMENTAL SCIENCES↗