Secondary ice production (SIP) can increase ice crystal number concentration (ICNC) by several orders of magnitude, particularly in clean clouds with low concentrations of ice-nucleating particles (INPs). The most common SIP process in models is rime splintering (RS) also called as the Hallett-Mossop process. The generally adopted RS-formulation gives 350 splinters per milligram of rimed ice at the temperature of 268 K. We used large-eddy simulations to examine if rime splintering could explain the high ICNC observed during the ACLOUD (Arctic CLoud Observations Using airborne measurements during polar Day) campaign where cloud temperatures close to 268 K are favourable for rime splintering. With the default model setup, the splinter production rate had to be multiplied by a factor ten to close the gap between modelled and observed ICNCs. Similar changes have been made in other modelling studies. The factor of ten multiplier helped to trigger SIP so that it became a self-sustaining process, fully independent of the primary freezing initiated by INPs. Our simulations reached realistic steady-state ICNCs and maintained stable mixed-phase clouds through the 24 h simulation time. Additional sensitivity tests showed that the efficiency of SIP depends strongly on model parametrizations (e.g., fall velocity–mass–dimension parametrizations and those describing the dependency of SIP on temperature and particle size and habit) and air temperature, so that simulations with a modified setup were able to reproduce the observed ICNCs without the factor of ten multiplier.
Research Council of Finland322532359342HORIZON EUROPE Climate, Energy and Mobility101137639101137680Introduction
Shallow mixed-phase clouds (MPCs) are common over high-latitude marine regions . Their important role in the formation of precipitation and the radiation budget make them highly sensitive elements in global climate and weather models . Clouds reflect most of the incoming short-wave solar radiation, but they also cause long-wave warming, which can cause sea-ice melting over high-latitudes. This is one of the main drivers of the Arctic amplification .
Global climate models struggle representing MPCs mainly because of low spatial and temporal resolution. MPC are inherently unstable thermodynamic systems, highly sensitive to turbulent surface fluxes and aerosol perturbations in the cloud condensation nuclei (CCN) and ice-nucleating particles (INP) number concentrations . Another challenge in modelling MPC comes from their susceptibility to experience secondary ice production (SIP), which can produce ice crystal number concentrations (ICNC) up to four orders of magnitude higher than INP number concentrations . Field campaigns such as the Mixed-Phase Arctic Cloud Experiment (M-PACE) , the Ny-Ålesund AeroSol Cloud Experiment (NASCENT) and the Aerosol-Cloud Coupling And Climate Interactions in the Arctic (ACCACIA) offer robust observational evidence of the occurrence of SIP in Arctic clouds, not only due to large differences between ICNC and INP number concentrations, but also due to the presence of fragments of frozen drops, needles and sheaths as well as broken dendrite branches in images obtained from in-cloud sampling systems .
The three most studied SIP processes relevant for mixed-phase clouds are: (1) rime splintering also known as the Hallett-Mossop process, (2) droplet shattering during freezing and (3) ice-ice collisional breakup. Fragile ice crystals like dendrites may break up mechanically when colliding with another large ice particle . The process is most effective at temperatures around -15 °C. When the surface of a large drizzle droplet freezes, e.g. by contact with an ice particle, the resulting pressure increase within the droplet may cause the ice surface to break, which releases small ice fragments, or the droplet may eject an ice particle . Recent observations e.g., have shown that the process can take place at temperatures well above -10 °C but then the droplet needs to be large to produce significant number of secondary ice particles. The most studied SIP process is rime splintering where fragile heavily rimed ice particles release ice splinters when colliding with large drizzle droplets, although the exact mechanism is not well known . The process is most effective at temperatures close to -5 °C.
In this study we performed ten meter-resolution large-eddy simulations (LESs) using the University of California Los Angeles Large Eddy Simulation model combined with the two-moment cloud microphysics scheme of (UCLALES-SB) and with the Sectional Aerosol module for Large Scale Applications cloud microphysics scheme (UCLALES-SALSA) to investigate the interplay between primary and secondary ice production processes, which can determine the phase and longevity of Arctic mixed-phase clouds. Our LES study is based on observations from the ACLOUD (Arctic CLoud Observations Using airborne measurements during polar Day) campaign carried out at north-west of Svalbard (Norway) in May–June 2017 . Due to the observed conditions showing the absence of large drizzle drops and ice particles and cloud top temperatures close to -5 °C, we expected that rime splintering is the dominating SIP process. Our first simulations confirmed this expectation, but these also showed that rime splintering is not strong enough to produce significant ice concentrations. First, we will examine if we can artificially increase rime splintering to reproduce the observed ice concentration. Then we will examine the impacts of meteorological and modelling uncertainties on secondary ice production. We will also examine if the results are sensitive on microphysics by comparing two-moment and sectional representations. Overall, our study aims to quantify the potential of Hallett-Mossop process in representing secondary ice production in such warm mixed-phase clouds.
MethodsThe ACLOUD campaign
Current LES simulations are based on observations from the Arctic CLoud Observations Using airborne measurements during polar Day (ACLOUD) campaign . The campaign included airborne observations from the Arctic boundary layer and clouds to understand their roles in Arctic amplification. Low-level clouds were frequently observed during a warm and moist period from 30 May to 12 June 2017 . Here we focus on three flights with the Polar 6 aircraft conducted on the 2, 4, and 5 June when mixed-phase clouds were observed. These observations and the data analysis are described in detail by , so a summary focused on the model simulations is given here. During the three flights, there was a uniform non-precipitating cloud deck above pack ice in a region North of Svalbard (the current focus region is 8.5–12.0° E and 81.1–81.4° N). Our simulations will be focused on the 2 June flight, which is the one with the highest observed ice crystal concentrations and cloud top temperatures, but we will use the two other flights to assess the impact of variability of meteorological conditions.
As reported by , cloud liquid (LWP) and ice (IWP) water paths were in the range of 48–82 gm-2 and 4.1–9.5 gm-2, respectively, for the three flights. Cloud base heights were between 100 and 200 m while the cloud top height was consistently about 440 m. The maximum cloud droplet number concentrations (CDNCs) from average vertical profiles Fig. 6 varied between 66 and 152 cm-3. These correlate with the above-cloud aerosol concentrations ranging from 125 to 173 cm-3. Measured droplet size distributions showed the absence of large drizzle droplets. Namely, the largest liquid droplets were about 30 µm in diameter (Supplement in ). The observed cloud top temperatures ranged from -6.7 to -4.6 °C. Figure a in the Appendix shows the observed temperatures along with the idealized initial profiles for the model simulations (see Sect. ).
Non-spherical ice particles in the diameter range from 9 to 1550 µm were detected using three different instruments as explained by . Particle shattering due to collisions with the probes had an impact on concentrations at the lower part of the cloud, and if such shattering was observed then ice particles smaller than 200 µm were excluded. The maximum ice crystal number concentrations at the upper part of the cloud are about 10 L-1 while the concentrations for particles larger than 200 µm at the lower part of the cloud are about 1 L-1Fig. 7. The expected range of ice concentration is thus 1–10 L-1, which is in line with other observations from that region e.g.,. The ice particle shape observations showed that most particles were single crystals including needles and columns. In addition, significant fraction (38.5 %) of the observed ice crystals were rimed.
Ice-nucleating particles (INPs) are needed to initiate primary ice formation at the observed cloud temperatures . There were no airborne INP measurements, but some ship-based measurements are available although only at -22.5 °C temperature . Even at such a low temperature, measured daily (2–5 June 2017) INP concentrations are in the order of 0.1 L-1. These measurements rule out possible pollution episodes and indicate that low INP concentrations can be expected for the region of interest. The best literature estimates for INP concentration at the cloud top temperature of about -5 °C is in the order of 1×10-3L-1e.g.,. This is significantly less than the observed cloud ice concentration of about 1 L-1Fig. 7. The three orders of magnitude difference between the INP and ice concentrations indicates that there is at least one active SIP process.
UCLALES-SALSA
UCLALES-SALSA is the LES model used in this study. SALSA refers to the sectional aerosol-cloud microphysics, which was added to UCLALES as an additional module. UCLALES with the default “SB” two-moment bulk microphysics by is a commonly used LES model especially for liquid clouds e.g.,. In this study we also use the more recently updated two-moment bulk ice microphysics , which is also used in large scale models . For clarity, this model version is referred to as UCLALES-SB. Both model versions share the same LES framework including radiative transfer and surface interactions and only their microphysics differ.
Computationally light but simplified SB microphysics allows conducting hundreds of simulations, which is useful for tasks like sensitivity tests where the impacts of model parameters on predictions is quantified. SALSA microphysics allows explicit modelling of aerosol-cloud-ice processes but this comes with a significant computational cost. By comparing predictions from both SB and SALSA microphysics, we can see the impact of the level of microphysical details.
SB microphysics
Liquid clouds in UCLALES-SB are diagnostic, which means that cloud water mixing ratio is diagnosed by using the saturation adjustment approach and a fixed cloud droplet number concentration is specified as a model input. The two-moment rain microphysics by describe both total mass and number concentrations while size distribution is assumed to follow a fixed gamma-distribution. Rain drop formation is based on a autoconversion parametrization and then the droplets can grow by condensation of water vapour and by collecting cloud droplets and smaller rain drops. The liquid-cloud scheme was extended for mixed-phase and ice clouds by , , and . The solid particle types include ice, snow, graupel, and hail, which have both mass and number as prognostic variables. The two-moment scheme accounts for various interactions between liquid and solid particles. The details are given in the original publications, so only a brief description is given here. The ice category represents small ice crystals formed by ice nucleation that are growing mainly by deposition of water vapour. In the absence of prognostic aerosols and thus INPs, primary ice nucleation is parametrized, so that the in-cloud ice crystal number concentration depends only on temperature. Collisions of ice with cloud droplets and larger rain drops leads to rimed ice, and depending on the resulting particle size, those particles are described by the snow, graupel and hail categories. Further riming and accretion lead to even larger particles and the resulting type is determined based on the size and type of colliding particles.
In this study we use hydrometeor parametrizations (fall velocity–mass–dimension parametrizations, parameters of the gamma-distribution, and mass limits) from with the exception that fall velocity–mass–dimension parametrization for ice which is from . This change was made because the mass–dimension parametrization of ice from produces exceptionally high dimensions compared with those from any other parametrization used in our simulations. used mass–dimension parametrizations from to calculate the ice water path (IWP) from the measured particle size. This parametrization, which happens to be the same as the snow parametrization, gives the same mass for 1.2 mm particles as the parametrization by . For particles smaller than 1.2 mm, the parametrization gives higher mass than the parametrization by . Because simulated particles are typically smaller than 1.2 mm, the parametrization gives smaller dimension than that from . This will be examined in Sect. .
The only SIP process included is rime splintering (Eq. ). Splinter production rates (dNi/dt, s-1) are parametrized as product of constant 350 mg-1 giving the number of splinters per milligram of rime, temperature-dependent efficiency term f(T), and water mass riming rate dmrime/dt (kgs-1) e.g.. The efficiency is linear between the minimum (zero at 265 K), optimal (one at 268 K) and the maximum (zero at 270 K) temperatures. We assume that splinters are small and therefore assign them to the ice category.
dNi/dt=350mg-1f(T)dmrime/dt
This parametrization includes rime mass (dmrime/dt) from collisions between any liquid droplet and solid ice particle. Notably, any limits for droplet diameter such as 25 µm minimum e.g., are excluded as it would require calculating incomplete gamma functions. Also, the size limits are more important for parametrizations where the number of splinters depends on the number of droplets collected . Some studies have also limitations for particle types that can produce splinters, for example, , , and exclude ice particles while include those. In this case, ice happens to be the dominant frozen particle type, so excluding it would essentially prevent rime splintering.
UCLALES-SB has mass concentration and the average maximum dimension thresholds for all collisions including riming, which have no physical meaning but presumably were used to reduce computational costs. For example, the default minimum total ice water mixing ratio and dimension are 10−5kgm-3 and 150 µm, respectively, for ice-cloud collisions (the corresponding limits for cloud are 10−6kgm-3 and 10 µm). A dimension of 150 µm is not a real limitation for the currently used mass–dimension parametrization of ice. However, 10−5kgm-3 is a high value considering the low primary ice concentration and the slow depositional growth rates. Thus, for all simulations here, we set the solid particle concentration limits to 10−9kgm-3, which is the same as the threshold concentration for rain. made the same conclusion on concentration limits when simulating cumulus clouds over the Southern Ocean. Likewise, reduced thresholds so that rime splintering could happen in their simulations with microphysics based on the Ny-Ålesund Aerosol Cloud Experiment (NASCENT). Similar adjustments have been made by , , and .
SALSA microphysics
Cloud microphysics in UCLALES-SALSA is treated using a sectional (bin) approach where aerosol, cloud, rain, and ice are described using several size sections (bins) for which microphysical processes are calculated. The liquid and ice cloud microphysics are originally described by and , respectively. The bins keep track of chemical composition (mass of solutes and water) and the number of aerosol particles and hydrometeors. Aerosol and cloud droplet size bins are based on the dry particle size, which includes solutes but not water and assumes a spherical particle shape. Rain droplet size bins are based on the wet size that accounts for the droplet volume including solutes and water. For ice we also use liquid water-equivalent wet size bins which are independent of the assumed ice particle shape. The wet size is basically the same as the size of a liquid droplet resulting in from ice being melted. Here, the aerosol is described with 12 logarithmically distributed dry size bins from 10 nm to 3 µm. Cloud bins are based on aerosol bins, but they do not include the first three bins (nucleation mode) as these particles are too small to activate. Seven logarithmically distributed rain bins range from 50 to 2000 µm and ten ice bins from 10 to 2000 µm (water-equivalent wet diameter).
Water is allowed to partition between vapour, liquid and ice phases based on equilibrium conditions at the droplet (so-called κ-Köhler; ) and ice particle surfaces. Water vapour flux is diffusion-limited and related to ambient saturation ratio, thus this non-equilibrium approach allows the prediction of supersaturation and cloud activation without any additional parametrizations. Here cloud activation means that when the wet size of an aerosol bin exceeds the critical droplet size, it is moved (partially or completely) to a corresponding cloud bin. Rain drop formation can be based on either autoconversion-like bulk parametrization from SB microphysics e.g., or counting the cloud-cloud collision where the resulting droplet size exceeds a threshold often set to 20 µm . Independent of the origin, rain drops will grow mainly by colliding with smaller cloud and rain droplets and eventually precipitate if conditions are suitable. Because liquid precipitation was not observed and simulated rain water paths were negligible, we will use the simple autoconversion-like bulk parametrization.
For this study, we implemented the same rime splintering parametrization as used in the SB microphysics (Eq. ). SALSA microphysics has certain concentration limits for all processes including riming, but the limits represent numerical accuracy of the model. Additional size limits are available for calculating the riming rate for the rime splintering process. The limit was set to 10 µm for cloud droplets, rain drops and ice particles. This means that the smallest cloud droplet bins can be excluded while all rain drops and ice particles are typically larger than 10 µm. When the size and particle type limits for riming and rime splintering are essentially removed, both SALSA and SB microphysics have similar chances to produce secondary ice particles.
In addition to the rime splintering, we implemented the and parametrizations for droplet shattering (DS) and ice-ice collisional breakup (IIBR), respectively. The latter is based on parametrization developed by but has revised parameters. Previously, implemented these processes into another SALSA version, so here we implement these into our SALSA version. Some IIBR parameters depend on rime mass fraction, but this is not predicted in our SALSA version. Therefore, we use parameters for unrimed particles, which have rime mass fraction less than 0.5. parametrization expects that only those collisions where ice crystals are not sticking together can produce secondary ice. We set the ice-ice sticking efficiency to 0.2, which is similar to the rime mass fraction dependent values used by and . Here DS includes both cloud and rain droplets, just as rime splintering. The parametrization for Mode 1 (small ice, large drop) droplet shattering is limited to droplet diameters larger than 50 µm and temperatures below 270.15 K while Mode 2 (large ice, small drop) is limited to droplet diameters larger than 150 µm . Size or temperature limits are not applied to IIBR.
LES setup
Adjustable model parameters are given in Table in the Appendix. The LES domain covers a horizontal area of 10 km× 10 km and extends vertically up to 1 km. Horizontal resolution is 100 m and vertical resolution is 10 m below 600 m. Above 600 m, vertical resolution increases by 3 % for each vertical level. Simulations have a maximum time step of 1 s and the total simulation time is 24 h where the first hour is with liquid clouds only (spin-up). The spin-up is used to allow the development of turbulence in a liquid cloud before particle sedimentation and rain and ice microphysics are fully included. For short-wave radiation, the solar zenith angle is fixed to 60° to match with the observations made during 1–2 h around the local noon in early June. In addition, following , sea surface albedo is set to 0.5, which represents partial ice cover over pack ice. Long-wave emissions are based on the surface temperature, which is set to be the same as that of the initial atmospheric profile at the lowest model level. For pack ice we set the surface roughness to 0.04 msee, e.g.,. Latent and sensible heat fluxes are set to 15 and 0 Wm-2, respectively, based on initial tests where the fluxes were simulated. Large scale subsidence is described with a constant divergence of 5 × 10−6s-1. This relatively large value was selected to limit the increase in cloud top height and LWP driven by radiative cooling.
Statistics are calculated every 2 min, and this is the output frequency of domain mean statistics. Horizontally and time-averaged profiles, and vertically integrated instantaneous column outputs are saved every 10 min.
For SB microphysics, we set the base case CDNC to 80 × 106kg-1≈ 100 × 106m-3, which is in the range of the observed maximum values from 66 × 106 to 152 × 106m-3. For SALSA we assume ammonium sulfate aerosol (hygroscopicity parameter κ= 0.6) so that the shape of the initial unimodal log-normal size distribution (Dg=106nm and σ=1.81) is based on observations (see Fig. b) and the total aerosol concentration is set to 150 × 106kg-1. With this initial aerosol, simulated CDNC will be similar with the fixed value of 80 × 106kg-1 used by the SB microphysics. For SALSA ice, we use the same fall velocity–mass–dimension parametrizations from as with UCLALES-SB.
Because we are focusing on SIP, we use a simple primary ice formation approach where the in-cloud INP concentration is given as an input value. In practice, this means that cloud droplets freeze until the total ice concentration reaches the given INP concentration. As explained in Sect. , the expected INP concentration can be as low as 10−3L-1 or about 1 kg-1 (based on literature) but should not be larger than 0.1 L-1 or about 100 kg-1 (observed at -22.5 °C). So, with these LES simulations, we will test different INP concentrations including 1, 10 and 100 kg-1. These INP concentrations are from one to three orders of magnitude lower than the observed ice concentration of at least 1 L-1 or about 1000 kg-1. We will take this as our minimum target value, thus in all following simulations we aim to reach ice concentration of 1000 kg-1. We will also conduct a simulation where INP concentrations is set to 1000 kg-1 and SIP is switched off. This represents a modelling approach where the observed ice concentrations can be reached even without SIP by using unrealistically high INP concentrations.
The initial temperature and humidity profiles were reconstructed based on the observed cloud extent, LWP, and cloud top temperature while also noting that these can change during the simulations. For example, ice formation and precipitation decrease LWP. The default profiles have specified surface temperature, relative humidity (RH) and pressure set to 1027 hPa, which allow the calculation of liquid water potential temperature (θL) and total water mixing ratio (rt) at the surface. These are assumed to be constant throughout a well-mixed boundary layer. A linear water vapour mixing ratio and potential temperature jumps (-0.9gkg-1 and 8 K, respectively) are assumed for the inversion layer (from 380 to 470 m), and above that the change is based on fixed gradients of -0.001gkg-1m-1 for rt and 0.01 Km-1 for θL. Figure a shows the observations along with the absolute temperatures calculated from these profiles (based on the saturation adjustment method). Additionally, Fig. b shows the observed and average wind components, which are used in all LES simulations.
The two additional initial profiles (cool and moist) were generated to test the impact of observational variability of the meteorological parameters. The cool profiles were generated by decreasing θL by 2 K and by decreasing rt so that LWP is not changing. The moist profiles were generated by increasing rt by 0.08 gkg-1. In this case, the latent heating within the cloud layer increases the absolute temperature compared with that of the default profile. Figure shows the initial temperature and moisture profiles as well as absolute temperature and RH based on the saturation adjustment method.
Initial (a) liquid water potential temperature (θL, K, solid lines) and (b) total water mixing ratio (rt, gkg-1) profiles for the LES simulations including the base case and the cool and moist cases. Panel (a) shows also the absolute temperature (T, K, dashed lines) and (c) shows RH (%), which were calculated from θL and rt by using the saturation adjustment method.
Results
The first test simulations with UCLALES-SALSA showed that the SIP processes were not able to produce significant amounts of ice even with the highest INP concentration of 100 kg-1. Panel (a) in Fig. shows the domain mean ice crystal number concentration (ICNC) for grid cells containing ice from two simulations: one without SIP and one with the three SIP processes switched on. Clearly, the ice concentrations are practically identical. Panel (b) shows the primary freezing rates and contributions from each SIP process (when SIP is on). Rime splintering (RS) produces about an order of magnitude less ice than primary freezing, but RS is still at least three orders of magnitude more efficient than ice-ice collisional breakup (IIBR) and droplet shattering (DS). Thus, in the following simulations we will focus on RS, but leave the other SIP processes on in all SALSA simulations. SB microphysics includes just rime splintering, and it is equally inefficient in the corresponding simulations. Next we will use the computationally fast UCLALES-SB for exploring suitable adjustments for RS SIP.
Simulated (a) domain mean ice crystal number concentration (ICNC) and (b) primary and secondary ice production rates for two SALSA simulations with INP concentration 100 kg-1 and SIP switched on (solid lines) and off (dashed lines).
Base case
showed that in addition to removing size and concentration limits, rime splintering had to be artificially increased to have an impact on ice concentration. So, to exceed the observed minimum ice crystal number concentration of about 1000 kg-1, we will first artificially increase secondary ice production simply by multiplying rime splintering rate (Eq. ) by a constant factor. The impacts of other possible adjustment will be examined in Sect. by means of sensitivity tests. Figure shows UCLALES-SB simulations with different INP concentrations (1, 10, 100, and 1000 kg-1) when rime splintering (RS) rates are multiplied by a constant factor (0, 1, 5, and 10). Panel (a) shows the domain mean ice crystal number concentration (ICNC) for grid cells containing ice and panels (b) and (c) show the horizontally averaged ice (IWP) and liquid (LWP) water paths, respectively. It should be noted that the SB microphysics includes ice, snow, graupel and hail categories, but only ice is shown here and in the following figures. This is because ice concentration is typically about two orders of magnitude higher than the concentration of any other solid particle type.
Simulated (a) ice crystal number concentration (ICNC), (b) ice water path (IWP) and (c) liquid water path (LWP) for the cases with different INP concentrations (1, 10, 100, and 1000 kg-1) and multipliers for the rime splintering (RS) secondary ice production rate (0, 1, 5, and 10) from UCLALES-SB. The thick gray line shows the target ICNC of 1000 kg-1.
The target (observed) ICNC of 1000 kg-1 is reached when INP concentration is set to that value and rime splintering is switched off by setting the multiplier to zero (the one dashed line). Simulations with more realistic INP concentrations of 1, 10 and 100 kg-1 and RS switched on (the unit multipliers) produce ice concentrations that are practically the same as the INP concentration, which means that RS secondary ice production is insignificant. Increasing RS by a factor 5 does help, but it requires about 15 h until the target ice concentration is reached. This is a long time compared to, for example, diurnal temperature variations which could trigger or prevent secondary ice production. Increasing the rate by a factor 10 means that SIP starts almost immediately after the spin-up and the target ice concentration is reached within 9 h. Interestingly, the factor of ten increase is enough for INP concentrations ranging from 1 to 100 kg-1 to reach the same steady-state ice concentration of about 2000 kg-1. Basically, this means that a strong enough SIP becomes self-sustaining, so it no longer needs or depends on the primary ice formation. The same behaviour is seen in the simulations with a factor of 5 increase, but there is a significant time delay and the steady-state ice concentration is lower (about 1000 kg-1). Clearly, the more efficient SIP is able to maintain a higher steady-state ICNC and IWP, which is closely related to ICNC.
When ice concentration is in the order of 1000 kg-1 (with or without SIP), the cloud starts to precipitate ice, and the continuous ice production and removal causes the decrease in LWP. This decrease in liquid cloud water reduces SIP but has no direct impact on the INP concentration. Thus, cloud properties are different depending on how ice formation is modelled: with high INP concentration (≈ ICNC) without SIP or low INP concentration (≪ ICNC) and SIP producing most ice particles. Naturally, SIP accounts for the feedback between ice production and cloud water and ice removal mechanisms such as precipitation.
From now on, the factor of 10 increase in secondary ice production is considered as the base case setting for SB microphysics. Interestingly, to match with the observed ice concentrations, needed to apply the same factor of ten adjustment to their rime splintering parametrization. noted that when only the rime splintering process is accounted for, ice production had to be increased by about a factor of 10–20 to obtain a good agreement with the observed ice concentrations.
Comparison of cloud microphysical schemes
Figure shows additional domain mean statistics from the four UCLALES-SB simulations (solid lines) described above and from the corresponding UCLALES-SALSA simulations (dashed lines). The reference simulations have INP concentration set to 1000 kg-1 and SIP is switched off (RS × 0). The three SIP simulations with INP concentrations set to 1, 10 and 100 kg-1 have the rime splintering (RS) ice production rate multiplied by ten (RS × 10). Time in this figure and in all SALSA simulations is limited to 15 h, because this is enough for both SB (see Fig. ) and SALSA simulations to reach a steady-state.
Simulated (a) ice crystal number concentration (ICNC), (b) ice water path (IWP), (c) liquid water path (LWP), (d) RS secondary ice production rate, and (e) primary cloud droplet freezing rate for the cases with different INP concentrations (1, 10, 100, and 1000 kg-1) and multipliers for the rime splintering (RS) secondary ice production rate (0 or 10). The thick gray line shows the target ICNC of 1000 kg-1.
The first thing that Fig. shows is that initially SB has higher SIP rates and ice concentrations, but SALSA has higher steady-state SIP rates so that the final ICNC is about 3000 kg-1 while this for SB this is 2000 kg-1. Overall, however, SB and SALSA predictions are similar when considering the differences in complexity between the microphysical models and the case where SIP increases ice concentrations by orders of magnitude. In this case, the complexity influences the time needed to run a 24 h simulation: about 1 h 20 min with the two-moment SB and 29 h with the sectional SALSA (parallel run with 100 CPU cores), i.e., there is a factor of 20 difference in computational costs. Clearly, the efficiency of SB makes it useful for conducting large numbers of test simulations while SALSA can provide additional details about the process.
Another thing that Fig. confirms is that RS SIP rate (panel d) exceeds the primary ice production rate (panel e) within two to eight hours depending on the INP concentration. When ice concentration becomes large enough (about 1000 kg-1), contribution from primary freezing becomes negligible. Thus, SIP maintains a feedback loop where INPs are not needed any more. This was confirmed by a test simulation, where switching off the primary freezing after 6 h had negligible impact.
The third thing that Fig. reveals is that the SIP rate correlates linearly with ice crystal number concentration and ice water path (IWP). Indeed, calculating spatial and temporal correlations between SIP rate and various model outputs reveals that the highest absolute Pearson's correlation coefficients are seen for ice number concentration, water vapour deposition rate, and IWP. Table shows Pearson's correlation coefficients for these variables (and LWP as a reference) calculated for both SB and SALSA simulations where the INP concentration is set to 100 kg-1 and SIP rate is multiplied by a factor of 10. Temporal correlation is calculated for the domain mean time series outputs and spatial correlation is calculated for snapshots of column-averaged or integrated 2D model outputs taken from the 7th hour. These three independent, i.e., not directly related to the SIP rate like riming rate, variables clearly stand out. ICNC, water vapour deposition rate, and IWP represent the 1st, 2nd and 3rd moment of the ice size distribution, respectively. The most obvious explanation is that SIP requires cloud droplet–ice collisions. Because cloud droplets and LWP are more evenly distributed, ice crystals are more important for the spatial correlation. The negative temporal correlation between SIP rate and LWP is related to the fact that ice is produced at the expense of liquid water, which is also apparent from Fig. . found a positive spatial correlation between observed vertical air velocity and SIP rates, but this is not that clear in our simulations, because the correlation coefficients (0.36 for SB and 0.25 for SALSA) are smaller than those for LWP. In fact, it looks like higher vertical velocities mean higher LWPs, which support ice production.
Spatial and temporal correlation between SIP rate and the three independent variables with the highest Pearson's correlation coefficients (Deposition rate, IWP, and ICNC) and additionally also LWP. Spatial correlations are calculated for the column-averaged model outputs taken from hour 7 while temporal correlation is calculated for the domain mean output time series. All simulations have INP concentration set to 100 kg-1 and SIP rate multiplied by a factor of 10.
As an example, Fig. a shows the temporal correlation between domain mean SIP rate and IWP for all SB and SALSA simulations where RS SIP is enabled. Figure b shows a snapshot of vertically integrated SIP rate contours over IWP colourmap (SB simulation, INP = 100 kg-1, RS × 10, time = 7 h). Clearly, the SIP rate contours match with the regions with high IWP.
Simulated (a) domain mean RS SIP rate as a function of IWP for the six simulations from the first 15 simulation hours and (b) instantaneous (t= 7 h) SIP rate contours (the change in column ICNC due to SIP, m-2s-1) and IWP colourmap for the SB simulations with INP = 100 kg-1 and SIP rate multiplied by a factor of 10.
Vertical distributions
Because the time to reach the target ice concentration of 1000 kg-1 depends on simulation (see Fig. ), the horizontally averaged profiles of the cloud parameters (Fig. ) are selected from the time step when ICNC reaches 1000 kg-1 for the first time (here we take the concentration at the altitude of 355 m). Because the profiles are fairly similar for simulations with INP concentrations of 1, 10 and 100 kg-1 and SIP rate multiplied by a factor of 10, we show only the first one in addition to the simulation with INP concentration of 1000 kg-1 without SIP. Panel (c) shows that the fixed CDNC for the SB microphysics matches well with the prognostic CDNC from the SALSA simulations. This is the case because total aerosol number was adjusted for this purpose. The most obvious difference is that SIP produces ICNC profiles that have a maximum within the cloud and lower values below while the profiles without SIP are essentially uniform and even increasing with decreasing altitude. The ACLOUD observations cannot show if the profiles should be uniform or not. This is because the observations are subjected to the typical detection limits (cannot see the smallest particles) as well as particle shattering effects , which impacts depend on altitude.
Profiles of (a) ICNC, (b) ice water mixing ratio, (c) CDNC, and (d) liquid water mixing ratio from the selected simulations. The time (hh:mm) for each simulation is selected so that ICNC first exceeds 1000 kg-1 at the altitude of 355 m.
Figure shows additional statistics about vertical distributions. Panel (a) shows the minimum absolute temperatures. These and especially their minimum values (cloud top temperatures) are similar for all simulations. Panel (b) show that the parametrized primary ice nucleation takes place mostly at the cloud top, but for the no-SIP simulation (1000 kg-1, RS × 0) the rates are significant for the whole cloud layer. SIP rates are distributed more evenly based on the cloud temperature and liquid water mixing ratio (Fig. d). The rime splintering process takes place between 265 and 270 K and the maximum efficiency is at 268 K, which is lower than the minimum cloud top temperatures. In addition, the maximum rate is seen at about 400 m altitude (high liquid water mixing ratios) where the temperatures are about 269 K. This indicates that a cooler temperature profile where the minimum temperature is slightly below 268 K would increase SIP rates. This will be examined in the next section. Due to the different distributions of the primary and secondary ice production, ice crystals in the SIP runs are larger (panel d) at the altitudes below the SIP region. There is also a difference between SB and SALSA microphysics so that the mean diameter is larger in SB simulations.
Profiles of (a) minimum absolute temperature, (b) primary freezing rate, (c) RS secondary ice production rate, and (d) ice crystal dimension from the selected simulations. The time (hh:mm) for each simulation is selected so that ICNC first exceeds 1000 kg-1 at the altitude of 355 m.
Sensitivity tests
Here we conduct sensitivity tests based on both observational and model variables that are most influential for SIP (see Table for model settings). We used simple trial and error method to determine a multiplier for the SIP rate so that the simulated ice concentration reaches the observed ice concentration of about 1000 kg-1, but a lower ice concentration was accepted in some cases where the higher SIP multiplier would have resulted in ice concentrations well above 1000 kg-1. These simulations are made with the computationally light SB microphysics, because as shown above, SALSA produces qualitatively similar results. We will focus on the case with the highest INP concentration of 100 kg-1 thus the base case simulation has SIP rate multiplied by ten (RS × 10). Here we limit simulation time to 10 h.
The observational variables that will be examined include CDNC, cloud temperature and liquid water path (LWP). Cloud temperature has a direct impact on the rime splintering process while CDNC and LWP have an impact on cloud dynamics. Figure shows how these observational uncertainties influence RS secondary ice production. When simulations are initialized with the humid total water mixing ratio profile (see Fig. ), LWP is about 30 % higher, but this has a relatively small impact on ice concentration (the Moist simulation has the same RS multiplier as in the base case). The cool profile (θL reduced by 2 K), on the other hand, has a clear impact on SIP. Just a factor of two multiplier is enough to start significant ice production in the Cool simulation. Initially, the ice concentration increases rapidly but soon the increased precipitation removal starts to limit ice production. Reducing CDNC from 80 × 106kg-1 by 50 % to 40×106kg-1 means that cloud droplets are larger, which means higher fall velocities, so also the riming rate increases. As a result, reducing the multiplier by 40 % from 10 to 6 is enough for the CDNC/2 simulation to reach and overpass the target ICNC of 1000 kg-1 around the 8th hour.
Sensitivity tests based on the observed variability of moisture (Moist), temperature (Cool) and cloud droplet number concentration (CDNC/2), mass-dimension-velocity (m–D–v) parametrizations from (SB06 ice and snow m–D–v), and temperature efficiency (f(T) in Eq. ). The last test is with cool profiles and SB06 ice m–D–v parametrization. In each simulation the RS SIP rate is multiplied by a factor so that ice concentration increases to about 1000 kg-1. INP concentration is 100 kg-1 in all simulations.
Modelling uncertainties are also significant and not only related to the rime splintering parametrization (the number of fragments per accumulated mass of rime, temperature limits, and possible size limits). From the many adjustable model parameters, mass-dimension-velocity (m–D–v) parametrizations seem to have the largest impact on SIP. In Fig. , we show simulations where the current ice parametrization is replaced by ice and snow parametrizations from , SB06. The SB06 ice parametrization represents an extreme parametrization regarding ice crystal size, which is increased by almost 100 % (from about 300 to 600 µm). Also, the fall velocity changes, but this has negligible impact on SIP. With this parametrization, SIP becomes more efficient so that a factor of four increase for SIP rate is enough (SB06 ice m–D–v). The SB06 snow parametrization includes the same m–v parametrization as is used for the current ice, so we only change the m–D parametrization which is the same as used by (from ) for calculating ice mass from the observed ice crystal shapes. This parametrization drastically reduces particle radius from about 300 to 100 µm. This reduced SIP so that it must be multiplied by a factor of 3.5 in addition to the original 10 (SB06 snow m–D–v). The importance of the m–D parametrization can be understood by the fact that collision kernel, which is used for calculating riming rate, is related to the square of the dimension while fall velocity has only linear dependency.
The currently used triangular temperature efficiency curve (linear from zero at 265 K≈-8 °C to one at 268 K≈-5 °C and back to zero at 270 K≈-3 °C) is less efficient compared with some other alternatives. For example, used piecewise constant efficiency curve so that it is one for -6 °C <T<-4 °C and 0.5 for the other temperatures between -8 °C <T<-2 °C . have unit efficiency between -8 °C <T<-3 °C and 0.01 elsewhere . has parabolic efficiency for temperature range from -2 to -8 °C. As a temperature efficiency test, we modified the efficiency so that it is one between -8 and -3 °C and zero elsewhere. This increases SIP so that only a factor of three increase is needed for the rime splintering (Temperature efficiency).
Overall, this sensitivity study suggests that with the cooler temperature profiles, slightly lower CDNC, the ice m–D–v parametrization from SB06, and the more efficient temperature dependency, the LES can reach the observed ice concentration of about 1000 kg-1 without modifying the rime splintering parametrization, and indeed this is the case. This is shown by the last sensitivity test (Cool, SB06 ice m–D–v) in Fig. . Here we have cooler temperature profile and use SB06 ice m–D–v parametrization, but some other combinations of the adjustments would have the same effect.
Conclusions
Here we used observations by to initialize LES simulations that aimed at reproducing the high ice concentrations observed in a relatively warm mixed-phase cloud deck where secondary ice production (SIP) was expected to dominate over the primary freezing initiated by INPs. With cloud top temperatures of about -5 °C, the dominant SIP process was rime splintering also known as the Hallett-Mossop process. With the default microphysical setup the model was not able to produce secondary ice, even after giving up from the commonly applied size and particle type limitations, so we artificially increased the rime splintering SIP rate by a constant factor. A factor of ten increase was well enough for the base case so that SIP was able to first rapidly increase the ice crystal number concentration (ICNC) from the primary ice concentration as low as 1 kg-1 to above the observed minimum value of about 1000 kg-1, and then maintain that over several hours. Basically this means that a strong enough SIP can become self-sustaining and thus be independent on the primary freezing. Interestingly, the factor of ten increase worked well for the two cloud microphysics models used in this study: the detailed sectional SALSA and the fast two-moment SB . The factor of ten happens to be the same enhancement as used in some previous studies .
With the artificially increased rime splintering SIP rate, SALSA produced steady-state ice concentrations of about 3000 kg-1 while this for SB was 2000 kg-1. Although statistically different, the values are surprisingly similar considering the differences between microphysical models and the case where SIP increases ice concentrations by orders of magnitude. Although SB and SALSA produce qualitatively similar results, their computational costs differ a lot. Namely, the computational costs of sectional SALSA microphysics are about 20 times higher than those of the bulk SB microphysics. Here, the efficiency of SB made it useful for conducting large numbers of sensitivity test simulations.
An alternative for artificially adjusting SIP rate is adjusting temperature efficiency or other model parametrizations (mass–dimension–fall velocity) or setup (temperature) to increase SIP. The triangular temperature efficiency curve (linear from zero at 265 K≈-8 °C to one at 268 K≈-5 °C and back to zero at 270 K≈-3 °C) used in the current rime splintering parametrization is less efficient compared with some other alternatives. For example, used piecewise constant efficiency curve so that it is one for -6 °C <T<-4 °C and 0.5 for the other temperatures between -8 °C <T<-2 °C . had unit efficiency between -8 °C <T<-3 °C and 0.01 elsewhere . Using any one of these would increase SIP rates. Moisture content has a much smaller effect while cloud droplet number concentration (CDNC) can have a significant effect especially when CDNC is low so that cloud droplets become larger. Mass–dimension–fall velocity parametrizations are important for the riming process, so different parametrizations can either increase or decrease SIP rates. At least in this case, suitable combination of those was able to initiate self-sustaining secondary ice production without using any artificial multiplier for the rime splintering rate. On the other hand, defining size or particle type limits (e.g., excluding ice particles like , , and ) may completely prevent SIP for certain cloud types, especially shallow clouds that have relatively small ice particles and narrow range of in-cloud temperatures. Ideally, such conditions should be replaced by smooth probability terms that reduce ice production in the case of unfavourable conditions. Overall, our results support the previous findings about the high sensitivity of SIP on various model setups and environmental conditions, which is a challenge for modelling.
For the shallow clouds in this study, the other potential SIP mechanisms are droplet shattering (DS) and ice-ice collisional breakup (IIBR). All SALSA simulations account for DS and IIBR based on parametrizations described in detail by , but these SIP processes were at least two orders of magnitude less efficient in producing secondary ice when compared with rime splintering (RS). Clearly, the current shallow cloud with relatively warm temperatures and small droplet sizes is more suitable for RS than for DS or IIBR. It is also clear that other SIP processes (at least droplet shattering and ice-ice collisional breakup) should be accounted for when conditions are more suitable for them. We focused on these processes in our previous study on cumulus congestus where cloud temperatures are lower and thus more favourable for IIBR and DS.
Regardless of the exact mechanism, it is essential to account for SIP rather than fix INP or ice crystal number directly. The simulations showed that vertical ice distributions and ice crystal sizes are different depending on how they were simulated. Moreover, using SIP allows the negative feedback between precipitation and ice production, which allows the development of stable mixed-phase clouds. Increase in the ice concentration will deplete liquid water, which in turn reduces SIP rates and stabilizes the cloud phase partitioning. With high fixed ice concentrations, the clouds are more likely to dissipate or glaciate, which is an issue seen in many large-scale models. Although the simulated vertical ice profiles were different with and without SIP, it was not possible to see if the observations match better with either one of those. However, this could be possible in future studies.
A recent study has questioned the existence of the rime splintering process. Our study cannot confirm that the process is real, but at least the simulated ice concentrations match well with the observed concentrations. The currently (and commonly) used parametrization is simple enough for models that have simplified microphysics (e.g., large-scale models) so it is useful at least for now.
Simulation settings
Figure a shows the observed temperatures from three research flights during the warm period and Fig. b shows wind speed components from flight 11 . The solid black lines indicate the default LES initialization based on observations from flight 11 (2 June 2017). The temperature profiles are adiabatic below cloud top and adjusted to match with the observed cloud top temperature and LWP. Adiabatic temperature and cloud water profiles are reconstructed based on given surface temperature and RH. To account for the radiative cooling, the minimum temperature seen at the cloud top is set to be slightly warmer than the observed minimum temperature of -4.56 °C. The initial rapid cooling decreases simulated minimum temperatures to -4.5 °C and then the cooling continues at a slower rate of -0.03 °C h−1. LWP is calculated from the cloud water profiles, and to account for the increasing LWP seen in simulations with low ice concentration, the initial value is lower than the observed LWP. The cooler initial profile is reconstructed by reducing surface temperature by 2 K and adjusting surface RH so that LWP is the same as in the default case. The third profile is generated by increasing surface RH so that LWP increases. Overall, these setups cover the observed cloud top temperature and LWP ranges. The initial wind profiles were calculated from the observations as a weighted mean. The weight for altitude x (based on the LES grid) for wind velocity observation at altitude zi is wi(x)=exp(-((x-zi)/(50m))2).
Figure a shows the observed total aerosol number concentration and Fig. b shows the average ambient aerosol size distribution . The average aerosol size distribution includes observations from time period 12:12:40–12:18:58 when sampling ambient aerosol (marked with the larger dots in panel a). The log-normal fit in panel (b) was used to initialize aerosol size distribution for SALSA except that the total aerosol number concentration was set to 150 × 106kg-1 (195 cm-3 when air density is 1.3 kgm-3) so that the simulated CDNC matched with the observed value of about 80 × 106kg-1.
Table shows the model parameters and settings for all LES simulations. These are also described in the main text.
Observed (a) temperature profiles from three research flights during the warm period, and (b) wind profiles for flight 11 . The solid black lines indicate the default LES initialization based on observations from flight 11 (2 June 2017). The cool and moist temperature profiles are used for additional sensitivity tests described in the main text.
Observed (a) total aerosol number concentration time series and (b) the average ambient aerosol size distribution. The data is from 2 June 2017, flight . The black and grey colours in (a) indicate time periods when measuring ambient aerosol and cloud particle residuals, respectively. The ambient aerosol size distribution in (b) is averaged from time period 12:12:40–12:18:58, which is marked with the larger dots in (a). Panel (b) also shows a log-normal fit to the data covering the SALSA aerosol bins (bin limits indicated by the grey vertical lines).
Model parameters and other simulation settings.
ParameterCommon defaultsHorizontal domainNx=Ny=100, dx=dy=100mVertical domainNz=85, dz=10m below 600 m and stretched by 1.03 above thatTimeMaximum step = 1 s, total = 86 400 s, spin-up = 3600 sOutputsStatistics every 120 s, averages every 600 sMean windsu=9.0ms-1, v=6.0ms-1Θ00270 KFour-stream radiationSZA = 60°, background = mid-latitude winter atmosphere (kmlw.lay), α=0.5, SST = 272 KLarge-scale divergence5×10-6s-1Surfacez0=0.04m, SHF = 0 Wm-2, LHF = 15 Wm-2Primary iceCloud droplets freeze when Si>0 and ICNC < INP concentrationINP concentrations1, 10, 100, and 1000 kg-1RS temperature efficiencyf(T) is linear between Tmin=265K, Topt=268K, and Tmax=270KIce v=av⋅mbvav=27.7, bv=0.216Ice D=aD⋅mbDaD=0.835, bD=0.390SB defaultsFixed CDNC80 × 106kg-1Riming limitsDroplets: rc>10-6kgm-3, Dc>10µmIce: ri>10-9kgm-3, Di>0µmSALSA defaultsInitial aerosolLog-normal size distribution: Ntot=150×106kg-1, Dg=0.106µm, σ=1.81Composed of sulfate: M=132.14gmol-1, ρ=1770kgm-3, ν=2.49Riming limitsDroplets: Nc>10-3m-3, Dc>10µmIce: Ni>10-6m-3, Di>10µmAerosol bins12 logarithmically spaced bins between 10 and 3000 nmRain bins7 logarithmically spaced bins between 50 and 2000 µmIce bins10 logarithmically spaced bins between 10 and 2000 µmSensitivity testsMoistMoist initial profilesCoolCool initial profiles and SST = 269.985 KCDNC/2CDNC =40×106kg-1SB06 ice m–D–vav=317, bv=0.363, aD=0.217, bD=0.302SB06 snow m–D–vav=27.7, bv=0.216, aD=8.156, bD=0.526RS temperature efficiencyf(T)=1 between 265 and 270 KCode and data availability
Brief description of the simulations, source code of UCLALES-SALSA, and the simulation data used in this publication are available from 10.5281/zenodo.18184323.
Author contributions
TR designed and conducted UCLALES-SALSA simulations. TR, SC, MP, and SR have contributed to developing the UCLALES-SALSA model. EJ provided the observational data used in this study. TR prepared the manuscript with contributions from all co-authors.
Competing interests
The contact author has declared that none of the authors has any competing interests.
Disclaimer
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims made in the text, published maps, institutional affiliations, or any other geographical representation in this paper. The authors bear the ultimate responsibility for providing appropriate place names. Views expressed in the text are those of the authors and do not necessarily reflect the views of the publisher.
Acknowledgements
The authors wish to acknowledge CSC – IT Center for Science, Finland, for computational resources.
Financial support
This research has been supported by the Research Council of Finland (decision numbers 322532 and 359342) and by the European Union's Horizon Europe CleanCloud (grant agreement no. 101137639) and CERTAINTY (no. 101137680) projects.
Review statement
This paper was edited by Luisa Ickes and reviewed by two anonymous referees.
ReferencesAckerman, A. S., vanZanten, M. C., Stevens, B., Savic-Jovcic, V., Bretherton, C. S., Chlond, A., Golaz, J.-C., Jiang, H., Khairoutdinov, M., Krueger, S. K., Lewellen, D. C., Lock, A., Moeng, C.-H., Nakamura, K., Petters, M. D., Snider, J. R., Weinbrecht, S., and Zulauf, M.: Large-Eddy Simulations of a Drizzling, Stratocumulus-Topped Marine Boundary Layer, Mon. Weather Rev., 137, 1083–1110, 10.1175/2008MWR2582.1, 2009.Ahola, J., Korhonen, H., Tonttila, J., Romakkaniemi, S., Kokkola, H., and Raatikainen, T.: Modelling mixed-phase clouds with the large-eddy model UCLALES–SALSA, Atmos. Chem. Phys., 20, 11639–11654, 10.5194/acp-20-11639-2020, 2020.Atlas, R. L., Bretherton, C. S., Blossey, P. N., Gettelman, A., Bardeen, C., Lin, P., and Ming, Y.: How Well Do Large-Eddy Simulations and Global Climate Models Represent Observed Boundary Layer Structures and Low Clouds Over the Summertime Southern Ocean?, J. Adv. Model. Earth Sy., 12, e2020MS002205, 10.1029/2020MS002205, 2020. Blahak, U.: Towards a better representation of high density ice particles in a state-of-the-art two-moment bulk microphysical scheme, 15th International Conference on Clouds and Precipitation, 7–11 July 2008, Cancun, Mexico, International Commission on Clouds and Precipitation, 2008.Brown, P. R. A. and Francis, P. N.: Improved Measurements of the Ice Water Content in Cirrus Using a Total-Water Probe, J. Atmos. Ocean. Tech., 12, 410–414, 10.1175/1520-0426(1995)012<0410:IMOTIW>2.0.CO;2, 1995.Calderón, S. M., Hyttinen, N., Kokkola, H., Raatikainen, T., Lawson, R. P., and Romakkaniemi, S.: Secondary ice formation in cumulus congestus clouds: insights from observations and aerosol-aware large-eddy simulations, Atmos. Chem. Phys., 25, 14479–14500, 10.5194/acp-25-14479-2025, 2025.Cesana, G. and Storelvmo, T.: Improving climate projections by understanding how cloud phase affects radiation, J. Geophys. Res.-Atmos., 122, 4594–4599, 10.1002/2017JD026927, 2017.Cotton, W. R., Tripoli, G. J., Rauber, R. M., and Mulvihill, E. A.: Numerical Simulation of the Effects of Varying Ice Crystal Nucleation Rates and Aggregation Processes on Orographic Snowfall, J. Appl. Meteorol. Clim., 25, 1658–1680, 10.1175/1520-0450(1986)025<1658:NSOTEO>2.0.CO;2, 1986.Eirund, G. K., Possner, A., and Lohmann, U.: Response of Arctic mixed-phase clouds to aerosol perturbations under different surface forcings, Atmos. Chem. Phys., 19, 9847–9864, 10.5194/acp-19-9847-2019, 2019.Ehrlich, A., Wendisch, M., Lüpkes, C., Buschmann, M., Bozem, H., Chechin, D., Clemen, H.-C., Dupuy, R., Eppers, O., Hartmann, J., Herber, A., Jäkel, E., Järvinen, E., Jourdan, O., Kästner, U., Kliesch, L.-L., Köllner, F., Mech, M., Mertes, S., Neuber, R., Ruiz-Donoso, E., Schnaiter, M., Schneider, J., Stapf, J., and Zanatta, M.: A comprehensive in situ and remote sensing data set from the Arctic CLoud Observations Using airborne measurements during polar Day (ACLOUD) campaign, Earth Syst. Sci. Data, 11, 1853–1881, 10.5194/essd-11-1853-2019, 2019.Ferrier, B. S.: A Double-Moment Multiple-Phase Four-Class Bulk Ice Scheme. Part I: Description, J. Atmos. Sci., 51, 249–280, 10.1175/1520-0469(1994)051<0249:ADMMPF>2.0.CO;2, 1994.Field, P. R., Lawson, R. P., Brown, P. R. A., Lloyd, G., Westbrook, C., Moisseev, D., Miltenberger, A., Nenes, A., Blyth, A., Choularton, T., Connolly, P., Buehl, J., Crosier, J., Cui, Z., Dearden, C., DeMott, P., Flossmann, A., Heymsfield, A., Huang, Y., Kalesse, H., Kanji, Z. A., Korolev, A., Kirchgaessner, A., Lasher-Trapp, S., Leisner, T., McFarquhar, G., Phillips, V., Stith, J., and Sullivan, S.: Secondary Ice Production: Current State of the Science and Recommendations for the Future, Meteorol. Monogr., 58, 7.1–7.20, 10.1175/AMSMONOGRAPHS-D-16-0014.1, 2017.Gierens, R., Kneifel, S., Shupe, M. D., Ebell, K., Maturilli, M., and Löhnert, U.: Low-level mixed-phase clouds in a complex Arctic environment, Atmos. Chem. Phys., 20, 3459–3481, 10.5194/acp-20-3459-2020, 2020.Grzegorczyk, P., Wobrock, W., Canzi, A., Niquet, L., Tridon, F., and Planche, C.: Investigating secondary ice production in a deep convective cloud with a 3D bin microphysics model: Part I – Sensitivity study of microphysical processes representations, Atmos. Res., 313, 107774, 10.1016/j.atmosres.2024.107774, 2025.Hallett, J. and Mossop, S. C.: Production of secondary ice particles during the riming process, Nature, 249, 26–28, 10.1038/249026a0, 1974.Hartmann, J., Lüpkes, C., and Chechin, D.: 1Hz resolution aircraft measurements of wind and temperature during the ACLOUD campaign in 2017, PANGAEA, 10.1594/PANGAEA.902849, 2019.Hohenegger, C., Korn, P., Linardakis, L., Redler, R., Schnur, R., Adamidis, P., Bao, J., Bastin, S., Behravesh, M., Bergemann, M., Biercamp, J., Bockelmann, H., Brokopf, R., Brüggemann, N., Casaroli, L., Chegini, F., Datseris, G., Esch, M., George, G., Giorgetta, M., Gutjahr, O., Haak, H., Hanke, M., Ilyina, T., Jahns, T., Jungclaus, J., Kern, M., Klocke, D., Kluft, L., Kölling, T., Kornblueh, L., Kosukhin, S., Kroll, C., Lee, J., Mauritsen, T., Mehlmann, C., Mieslinger, T., Naumann, A. K., Paccini, L., Peinado, A., Praturi, D. S., Putrasahan, D., Rast, S., Riddick, T., Roeber, N., Schmidt, H., Schulzweida, U., Schütte, F., Segura, H., Shevchenko, R., Singh, V., Specht, M., Stephan, C. C., von Storch, J.-S., Vogel, R., Wengel, C., Winkler, M., Ziemen, F., Marotzke, J., and Stevens, B.: ICON-Sapphire: simulating the components of the Earth system and their interactions at kilometer and subkilometer scales, Geosci. Model Dev., 16, 779–811, 10.5194/gmd-16-779-2023, 2023.Huang, Y., Blyth, A. M., Brown, P. R. A., Choularton, T. W., Connolly, P., Gadian, A. M., Jones, H., Latham, J., Cui, Z., and Carslaw, K.: The development of ice in a cumulus cloud over southwest England, New J. Phys., 10, 105021, 10.1088/1367-2630/10/10/105021, 2008.Järvinen, E., Nehlert, F., Xu, G., Waitz, F., Mioche, G., Dupuy, R., Jourdan, O., and Schnaiter, M.: Investigating the vertical extent and short-wave radiative effects of the ice phase in Arctic summertime low-level clouds, Atmos. Chem. Phys., 23, 7611–7633, 10.5194/acp-23-7611-2023, 2023.Kanji, Z. A., Ladino, L. A., Wex, H., Boose, Y., Burkert-Kohn, M., Cziczo, D. J., and Krämer, M.: Overview of Ice Nucleating Particles, Meteorol. Monogr., 58, 1.1–1.33, 10.1175/AMSMONOGRAPHS-D-16-0006.1, 2017.Keinert, A., Spannagel, D., Leisner, T., and Kiselev, A.: Secondary Ice Production upon Freezing of Freely Falling Drizzle Droplets, J. Atmos. Sci., 77, 2959–2967, 10.1175/JAS-D-20-0081.1, 2020.Korolev, A., McFarquhar, G., Field, P. R., Franklin, C., Lawson, P., Wang, Z., Williams, E., Abel, S. J., Axisa, D., Borrmann, S., Crosier, J., Fugal, J., Krämer, M., Lohmann, U., Schlenczek, O., Schnaiter, M., and Wendisch, M.: Mixed-Phase Clouds: Progress and Challenges, Meteorol. Monogr., 58, 5.1–5.50, 10.1175/AMSMONOGRAPHS-D-17-0001.1, 2017.Kudzotsa, I., Phillips, V. T. J., Dobbie, S., Formenton, M., Sun, J., Allen, G., Bansemer, A., Spracklen, D., and Pringle, K.: Aerosol indirect effects on glaciated clouds. Part I: Model description, Q. J. Roy. Meteor. Soc., 142, 1958–1969, 10.1002/qj.2791, 2016.Li, G., Wieder, J., Pasquier, J. T., Henneberger, J., and Kanji, Z. A.: Predicting atmospheric background number concentration of ice-nucleating particles in the Arctic, Atmos. Chem. Phys., 22, 14441–14454, 10.5194/acp-22-14441-2022, 2022.Listowski, C., Delanoë, J., Kirchgaessner, A., Lachlan-Cope, T., and King, J.: Antarctic clouds, supercooled liquid water and mixed phase, investigated with DARDAR: geographical and seasonal variations, Atmos. Chem. Phys., 19, 6771–6808, 10.5194/acp-19-6771-2019, 2019.Luke, E. P., Yang, F., Kollias, P., Vogelmann, A. M., and Maahn, M.: New insights into ice multiplication using remote-sensing observations of slightly supercooled mixed-phase clouds in the Arctic, P. Natl. Acad. Sci. USA, 118, e2021387118, 10.1073/pnas.2021387118, 2021.McFarquhar, G. M. and Cober, S. G.: Single-Scattering Properties of Mixed-Phase Arctic Clouds at Solar Wavelengths: Impacts on Radiative Transfer, J. Climate, 17, 3799–3813, 10.1175/1520-0442(2004)017<3799:SPOMAC>2.0.CO;2, 2004.Mertes, S., Kästner, U., and Macke, A.: Airborne in-situ measurements of the aerosol absorption coefficient, aerosol particle number concentration and size distribution of cloud particle residuals and ambient aerosol particles during flight P6_206_ACLOUD_2017_1706021001, PANGAEA, 10.1594/PANGAEA.900414, 2019.Mioche, G., Jourdan, O., Ceccaldi, M., and Delanoë, J.: Variability of mixed-phase clouds in the Arctic with a focus on the Svalbard region: a study based on spaceborne active remote sensing, Atmos. Chem. Phys., 15, 2445–2461, 10.5194/acp-15-2445-2015, 2015.Mioche, G., Jourdan, O., Delanoë, J., Gourbeyre, C., Febvre, G., Dupuy, R., Monier, M., Szczap, F., Schwarzenboeck, A., and Gayet, J.-F.: Vertical distribution of microphysical properties of Arctic springtime low-level mixed-phase clouds over the Greenland and Norwegian seas, Atmos. Chem. Phys., 17, 12845–12869, 10.5194/acp-17-12845-2017, 2017.Morrison, H., Thompson, G., and Tatarskii, V.: Impact of Cloud Microphysics on the Development of Trailing Stratiform Precipitation in a Simulated Squall Line: Comparison of One- and Two-Moment Schemes, Mon. Weather Rev., 137, 991–1007, 10.1175/2008MWR2556.1, 2009.Murray, B. J., Carslaw, K. S., and Field, P. R.: Opinion: Cloud-phase climate feedback and the importance of ice-nucleating particles, Atmos. Chem. Phys., 21, 665–679, 10.5194/acp-21-665-2021, 2021.Noppel, H., Blahak, U., Seifert, A., and Beheng, K. D.: Simulations of a hailstorm and the impact of CCN using an advanced two-moment cloud microphysical scheme, Atmos. Res., 96, 286–301, 10.1016/j.atmosres.2009.09.008, 2010.Pasquier, J. T., Henneberger, J., Ramelli, F., Lauber, A., David, R. O., Wieder, J., Carlsen, T., Gierens, R., Maturilli, M., and Lohmann, U.: Conditions favorable for secondary ice production in Arctic mixed-phase clouds, Atmos. Chem. Phys., 22, 15579–15601, 10.5194/acp-22-15579-2022, 2022.Petters, M. D. and Kreidenweis, S. M.: A single parameter representation of hygroscopic growth and cloud condensation nucleus activity, Atmos. Chem. Phys., 7, 1961–1971, 10.5194/acp-7-1961-2007, 2007.Phillips, V. T. J., Yano, J.-I., and Khain, A.: Ice Multiplication by Breakup in Ice–Ice Collisions. Part I: Theoretical Formulation, J. Atmos. Sci., 74, 1705–1719, 10.1175/JAS-D-16-0224.1, 2017.Phillips, V. T. J., Patade, S., Gutierrez, J., and Bansemer, A.: Secondary Ice Production by Fragmentation of Freezing Drops: Formulation and Theory, J. Atmos. Sci., 75, 3031–3070, 10.1175/JAS-D-17-0190.1, 2018.Prenni, A. J., Harrington, J. Y., Tjernström, M., DeMott, P. J., Avramov, A., Long, C. N., Kreidenweis, S. M., Olsson, P. Q., and Verlinde, J.: Can Ice-Nucleating Aerosols Affect Arctic Seasonal Climate?, B. Am. Meteorol. Soc., 88, 541–550, 10.1175/BAMS-88-4-541, 2007.Previdi, M., Smith, K. L., and Polvani, L. M.: Arctic amplification of climate change: a review of underlying mechanisms, Environ. Res. Lett., 16, 093003, 10.1088/1748-9326/ac1c29, 2021.Raatikainen, T.: Code and data for “Can rime splintering explain the ice production in Arctic mixed-phase clouds?”, Zenodo [data set], 10.5281/zenodo.18184323, 2026.Rangno, A. L. and Hobbs, P. V.: Ice particles in stratiform clouds in the Arctic and possible mechanisms for the production of high ice concentrations, J. Geophys. Res.-Atmos., 106, 15065–15075, 10.1029/2000JD900286, 2001.Reisner, J., Rasmussen, R. M., and Bruintjes, R. T.: Explicit forecasting of supercooled liquid water in winter storms using the MM5 mesoscale model, Q. J. Roy. Meteor. Soc., 124, 1071–1107, 10.1002/qj.49712454804, 1998.Schäfer, B., David, R. O., Georgakaki, P., Pasquier, J. T., Sotiropoulou, G., and Storelvmo, T.: Simulations of primary and secondary ice production during an Arctic mixed-phase cloud case from the Ny-Ålesund Aerosol Cloud Experiment (NASCENT) campaign, Atmos. Chem. Phys., 24, 7179–7202, 10.5194/acp-24-7179-2024, 2024.Seidel, J. S., Kiselev, A. A., Keinert, A., Stratmann, F., Leisner, T., and Hartmann, S.: Secondary ice production – no evidence of efficient rime-splintering mechanism, Atmos. Chem. Phys., 24, 5247–5263, 10.5194/acp-24-5247-2024, 2024.Seifert, A.: On the Parameterization of Evaporation of Raindrops as Simulated by a One-Dimensional Rainshaft Model, J. Atmos. Sci., 65, 3608–3619, 10.1175/2008JAS2586.1, 2008.Seifert, A. and Beheng, K. D.: A double-moment parameterization for simulating autoconversion, accretion and selfcollection, Atmos. Res., 59–60, 265–281, 10.1016/S0169-8095(01)00126-0, 2001.Seifert, A. and Beheng, K. D.: A two-moment cloud microphysics parameterization for mixed-phase clouds. Part 1: Model description, Meteorol. Atmos. Phys., 92, 45–66, 10.1007/s00703-005-0112-4, 2006.Seifert, A. and Heus, T.: Large-eddy simulation of organized precipitating trade wind cumulus clouds, Atmos. Chem. Phys., 13, 5631–5645, 10.5194/acp-13-5631-2013, 2013.Seifert, A., Köhler, C., and Beheng, K. D.: Aerosol-cloud-precipitation effects over Germany as simulated by a convective-scale numerical weather prediction model, Atmos. Chem. Phys., 12, 709–725, 10.5194/acp-12-709-2012, 2012.Seifert, A., Blahak, U., and Buhr, R.: On the analytic approximation of bulk collision rates of non-spherical hydrometeors, Geosci. Model Dev., 7, 463–478, 10.5194/gmd-7-463-2014, 2014.Sotiropoulou, G., Sullivan, S., Savre, J., Lloyd, G., Lachlan-Cope, T., Ekman, A. M. L., and Nenes, A.: The impact of secondary ice production on Arctic stratocumulus, Atmos. Chem. Phys., 20, 1301–1316, 10.5194/acp-20-1301-2020, 2020.Sotiropoulou, G., Vignon, É., Young, G., Morrison, H., O'Shea, S. J., Lachlan-Cope, T., Berne, A., and Nenes, A.: Secondary ice production in summer clouds over the Antarctic coast: an underappreciated process in atmospheric models, Atmos. Chem. Phys., 21, 755–771, 10.5194/acp-21-755-2021, 2021.Stevens, B. and Seifert, A.: Understanding macrophysical outcomes of microphysical choices in simulations of shallow cumulus convection, J. Meteorol. Soc. Jpn. Ser. II, 86A, 143–162, 10.2151/jmsj.86A.143, 2008.Stevens, B., Moeng, C.-H., and Sullivan, P. P.: Large-Eddy Simulations of Radiatively Driven Convection: Sensitivities to the Representation of Small Scales, J. Atmos. Sci., 56, 3963–3984, 10.1175/1520-0469(1999)056<3963:LESORD>2.0.CO;2, 1999.Stevens, B., Moeng, C.-H., Ackerman, A. S., Bretherton, C. S., Chlond, A., de Roode, S., Edwards, J., Golaz, J.-C., Jiang, H., Khairoutdinov, M., Kirkpatrick, M. P., Lewellen, D. C., Lock, A., Müller, F., Stevens, D. E., Whelan, E., and Zhu, P.: Evaluation of Large-Eddy Simulations via Observations of Nocturnal Marine Stratocumulus, Mon. Weather Rev., 133, 1443–1462, 10.1175/MWR2930.1, 2005.Sullivan, S. C., Hoose, C., and Nenes, A.: Investigating the contribution of secondary ice production to in-cloud ice crystal numbers, J. Geophys. Res.-Atmos., 122, 9391–9412, 10.1002/2017JD026546, 2017.Sullivan, S. C., Barthlott, C., Crosier, J., Zhukov, I., Nenes, A., and Hoose, C.: The effect of secondary ice production parameterization on the simulation of a cold frontal rainband, Atmos. Chem. Phys., 18, 16461–16480, 10.5194/acp-18-16461-2018, 2018a.Sullivan, S. C., Hoose, C., Kiselev, A., Leisner, T., and Nenes, A.: Initiation of secondary ice production in clouds, Atmos. Chem. Phys., 18, 1593–1610, 10.5194/acp-18-1593-2018, 2018b.Takahashi, T., Nagao, Y., and Kushiyama, Y.: Possible High Ice Particle Production during Graupel–Graupel Collisions, J. Atmos. Sci., 52, 4523–4527, 10.1175/1520-0469(1995)052<4523:PHIPPD>2.0.CO;2, 1995.Tonttila, J., Maalick, Z., Raatikainen, T., Kokkola, H., Kühn, T., and Romakkaniemi, S.: UCLALES–SALSA v1.0: a large-eddy model with interactive sectional microphysics for aerosol, clouds and precipitation, Geosci. Model Dev., 10, 169–188, 10.5194/gmd-10-169-2017, 2017.Tonttila, J., Afzalifar, A., Kokkola, H., Raatikainen, T., Korhonen, H., and Romakkaniemi, S.: Precipitation enhancement in stratocumulus clouds through airborne seeding: sensitivity analysis by UCLALES-SALSA, Atmos. Chem. Phys., 21, 1035–1048, 10.5194/acp-21-1035-2021, 2021.van der Dussen, J. J., de Roode, S. R., Ackerman, A. S., Blossey, P. N., Bretherton, C. S., Kurowski, M. J., Lock, A. P., Neggers, R. A. J., Sandu, I., and Siebesma, A. P.: The GASS/EUCLIPSE model intercomparison of the stratocumulus transition as observed during ASTEX: LES results, J. Adv. Model. Earth Sy., 5, 483–499, 10.1002/jame.20033, 2013.Weiss, A. I., King, J., Lachlan-Cope, T., and Ladkin, R.: On the effective aerodynamic and scalar roughness length of Weddell Sea ice, J. Geophys. Res.-Atmos., 116, 10.1029/2011JD015949, 2011.Wendisch, M., Macke, A., Ehrlich, A., Lüpkes, C., Mech, M., Chechin, D., Dethloff, K., Velasco, C. B., Bozem, H., Brückner, M., Clemen, H.-C., Crewell, S., Donth, T., Dupuy, R., Ebell, K., Egerer, U., Engelmann, R., Engler, C., Eppers, O., Gehrmann, M., Gong, X., Gottschalk, M., Gourbeyre, C., Griesche, H., Hartmann, J., Hartmann, M., Heinold, B., Herber, A., Herrmann, H., Heygster, G., Hoor, P., Jafariserajehlou, S., Jäkel, E., Järvinen, E., Jourdan, O., Kästner, U., Kecorius, S., Knudsen, E. M., Köllner, F., Kretzschmar, J., Lelli, L., Leroy, D., Maturilli, M., Mei, L., Mertes, S., Mioche, G., Neuber, R., Nicolaus, M., Nomokonova, T., Notholt, J., Palm, M., van Pinxteren, M., Quaas, J., Richter, P., Ruiz-Donoso, E., Schäfer, M., Schmieder, K., Schnaiter, M., Schneider, J., Schwarzenböck, A., Seifert, P., Shupe, M. D., Siebert, H., Spreen, G., Stapf, J., Stratmann, F., Vogl, T., Welti, A., Wex, H., Wiedensohler, A., Zanatta, M., and Zeppenfeld, S.: The Arctic Cloud Puzzle: Using ACLOUD/PASCAL Multiplatform Observations to Unravel the Role of Clouds and Aerosol Particles in Arctic Amplification, B. Am. Meteorol. Soc., 100, 841–871, 10.1175/BAMS-D-18-0072.1, 2019. Young, G., Jones, H. M., Choularton, T. W., Crosier, J., Bower, K. N., Gallagher, M. W., Davies, R. S., Renfrew, I. A., Elvidge, A. D., Darbyshire, E., Marenco, F., Brown, P. R. A., Ricketts, H. M. A., Connolly, P. J., Lloyd, G., Williams, P. I., Allan, J. D., Taylor, J. W., Liu, D., and Flynn, M. J.: Observed microphysical changes in Arctic mixed-phase clouds when transitioning from sea ice to open ocean, Atmos. Chem. Phys., 16, 13945–13967, 10.5194/acp-16-13945-2016, 2016.Young, G., Lachlan-Cope, T., O'Shea, S. J., Dearden, C., Listowski, C., Bower, K. N., Choularton, T. W., and Gallagher, M. W.: Radiative Effects of Secondary Ice Enhancement in Coastal Antarctic Clouds, Geophys. Res. Lett., 46, 2312–2321, 10.1029/2018GL080551, 2019.Zhao, X., Liu, X., Phillips, V. T. J., and Patade, S.: Impacts of secondary ice production on Arctic mixed-phase clouds based on ARM observations and CAM6 single-column model simulations, Atmos. Chem. Phys., 21, 5685–5703, 10.5194/acp-21-5685-2021, 2021.Ziegler, C. L., Ray, P. S., and MacGorman, D. R.: Relations of Kinematics, Microphysics and Electrification in an Isolated Mountain Thunderstorm, J. Atmos. Sci., 43, 2098–2115, 10.1175/1520-0469(1986)043<2098:ROKMAE>2.0.CO;2, 1986.