The sensitivities of in cloud and cloud top phase distributions to primary ice formation in ICON-LEM

Mixed phase clouds remain a leading source of uncertainty in our attempt to quantify cloud-climate and aerosol-cloud climate interactions. Nevertheless, recent advances in parametrizing the primary ice formation process, high resolution cloud modelling, and retrievals of cloud phase distributions from satellite data offer an excellent opportunity to conduct closure studies on the sensitivity of the cloud phase to microphysical and dynamical processes. Particularly, the reliability of satellite data to resolve the phase at the top of the cloud provides a promising benchmark to compare model output to. We run large eddy simulations with the new ICOsahedral Non-hydrostatic atmosphere model (ICON) to place bounds on the sensitivity of in cloud and cloud top phase to the primary ice formation process. State of the art primary ice formation parametrizations in the form of the cumulative ice active site density ns are implemented in idealized deep convective cloud simulations. We exploit the ability of ICON-LEM to switch between a two moment microphysics scheme and the newly developed Predicted Particle Properties (P3) scheme by running our simulations in both configurations for comparison. To quantify the sensitivity of cloud phase to primary ice formation, cloud ice content is evaluated against order of magnitude changes in ns at variable convective strengths. Furthermore, we assess differences between in cloud and cloud top phase distributions as well as the potential impact of updraft velocity on the suppression of the Wegener-Bergeron-Findeisen process. The study aims to evaluate our practical understanding of primary ice formation in the context of predicting the structure and evolution of mixed phase clouds.