Cloud Condensation on Insoluble Nuclei

• Cloud condensation on insoluble nuclei is a process that describes the nucleation, activation, and diffusion-limited growth of cloud droplets on non-soluble aerosol particles like mineral dust.
• The model integrates heterogeneous nucleation theory with detailed aerosol size distributions and activation thresholds to predict condensation behavior under varying environmental conditions.
• Observational validations and parameterizations link microphysical processes to improved weather and climate simulations, enhancing precipitation forecasts and radiative transfer modeling across planetary atmospheres.

A cloud condensation model on insoluble nuclei describes the nucleation, activation, and growth of cloud droplets or ice particles on the surfaces of non-soluble aerosol particles such as mineral dust, in contrast to classical Köhler-theory frameworks that emphasize aqueous and soluble seeds. These models capture the physics of heterogeneous nucleation, account for the kinetic and thermodynamic barriers unique to insoluble substrates, and underpin parameterizations relevant for rainfall initiation, visibility, radiative transfer, and climate in diverse planetary atmospheres, including Earth, Mars, and substellar objects (Arielly et al., 25 Jan 2026, Shaposhnikov et al., 2022, Gao et al., 2018, Svensmark et al., 2012).

1. Physical and Mathematical Basis

Cloud droplet activation on insoluble nuclei fundamentally differs from growth on soluble seeds. Instead of the critical supersaturation dictated by Köhler curves, nucleation is viewed as the stochastic formation of a water (or condensable) molecular “shell” on the dust grain. The condensation rate is subject to an energy barrier $E_{\rm cond}(T)$ and characterized by a Boltzmann-condensation probability, absent bulk dissolution.

For water vapor condensation, the key equations (Arielly et al., 25 Jan 2026, Gao et al., 2018) include:

• Water-vapor number density: $n_\nu=P_\nu/(k_B T)$, with $P_\nu=\phi P_{H_2O}(T)$ using empirical expressions for $P_{H_2O}(T)$ (e.g., Buck, Magnuss).
• Activation probability: $$P_{\rm cond}=\exp\left[-\frac{E_{\rm cond}(T)}{k_B T}(T_{\rm dew}-T)\right]$$.
• Total attachment probability for an $N_{\rm mol}$ molecular shell: $$P_T\approx N_{\rm mol}\exp\left[-\frac{E_{\rm cond}(T)}{k_B T}(T_{\rm dew}-T)\right]$$.
• Droplet number density: $n_d = n_\nu P_T$.

The nucleation rate for insoluble nuclei, as in Martian and substellar applications (Shaposhnikov et al., 2022, Gao et al., 2018), generally follows classical nucleation theory adapted for heterogeneous surfaces. For a nucleus of radius $r_{\rm nuc}$, the formation rate $J_{\rm het}$ incorporates the effect of contact angle, surface energy, flux of vapor molecules, and the Zeldovich factor:

$$J_{\rm het} = 4\pi^2 r_{\rm nuc}^2 a_c^2 \Phi c_{\rm surf} Z \exp\left[- \frac{F f(\mu, r_{\rm nuc}/a_c)}{k_B T}\right]$$

where $a_c$ is the critical cluster radius, $F$ is the cluster formation energy, and $f(\mu, x)$ accounts for the geometry.

2. Aerosol Size Distribution and Activation Threshold

The cloud condensation process is highly sensitive to the size distribution of insoluble aerosols. Models typically assume a log-normal or bimodal distribution for mineral dust or ice nuclei (Arielly et al., 25 Jan 2026, Shaposhnikov et al., 2022):

$$n_a(r_a) = \frac{n_{\rm tot}}{r_a\sqrt{2\pi}\sigma} \exp\left[-\frac{(\ln(r_a/T_a))^2}{2 \sigma^2}\right]$$

In the MPI-MGCM Martian hydrological scheme, bimodality is expressed as a superposition of small (effective radius $\sim0.05~\mu$m) and large (variable, $\sim0.1$–$0.5$ $\mu$m) log-normal modes. Vertical binning of this distribution is essential for accurate microphysical evolution and subsequent feedbacks (Shaposhnikov et al., 2022).

The minimum effective size for a dust particle to act as a CCN is determined by fitting the shell-formation model to backscatter data, with $r_{\rm nuc}$ estimates clustering at $46.7\pm1.8$ nm in arid-earth field studies (Arielly et al., 25 Jan 2026).

3. Condensation Microphysics and Growth

The growth of a cloud droplet on an insoluble nucleus follows diffusion-limited condensation, subject to the Kelvin effect for small radii:

• Mass flux: $$J_m=4\pi r D_v (\rho_{v,\infty}-\rho_{v,\text{surf}})F(Kn, \alpha)$$.
• Kelvin effect: $$p_{s,v}(r) = p_{s,v}(\infty)\exp[(2\gamma_w M_w)/(R T \rho_w r)]$$.
• Growth rate: $$\frac{dr}{dt} \approx \frac{D_v M_w}{\rho_w R T r}[p_v-p_{s,v}(r)]$$ (Svensmark et al., 2012).

The particle's condensation history is modified by the supply of condensable vapor and the competition with coagulation losses, as well as (in the case of ion-induced nucleation) catalytic processes that can replenish limiting vapor-phase reactants.

On Mars, condensational growth equations further include corrections for free-molecular regime diffusion and thermal transfer resistances, along with explicit sticking coefficients and the Kelvin–Cunningham slip correction for terminal velocities (Shaposhnikov et al., 2022).

4. Model Validation and Observational Constraints

Validation of cloud condensation on insoluble nuclei employs direct comparison of model-predicted optical backscatter with ceilometer or lidar observations and simultaneous radiosonde vertical profiles (Arielly et al., 25 Jan 2026). The transformation involves computing layerwise droplet number densities and translating these into predicted backscatter signals:

$B'(z) = A\,n_d(z)\,\exp[-X\,n_d(z)]$

where calibration constants $A$ and $X$ are adjusted to match observed profiles. In field campaigns, the model yields a high mean Pearson correlation ($\bar\rho=0.85\pm0.05$) between predicted and measured backscatter, validating the approach for determining CCN-activation thresholds and the minimum size of effective nuclei.

Sensitivity analyses demonstrate the robustness of derived $r_{\rm nuc}$ to uncertainties in relative humidity, temperature, and instrumental calibration.

5. Parameterization and Implementation in Weather and Climate Models

Modern parameterizations for numerical weather prediction (NWP) and general circulation models (GCMs) now incorporate activation schemes that replace or supplement classical Köhler solubility-based lookups with explicit insoluble-shell-formation steps in dust-rich environments (Arielly et al., 25 Jan 2026, Shaposhnikov et al., 2022, Gao et al., 2018). The central activation probability is computed as

$$P_T(T,\phi)=g(T,\phi)\exp\left[-4 f_{\rm pack}\frac{r_{\rm nuc}}{r_{\rm H_2O}}g(T,\phi)\right],\quad g(T,\phi)=\frac{k_B (T_{\rm dew}-T)}{E_{\rm cond}(T)}$$

Such schemes allow CCN number concentration to be efficiently predicted from routinely available temperature and humidity data, together with real-time or climatological CCN size input.

In planetary and brown-dwarf contexts, the eddy-diffusion–sedimentation frameworks use a sedimentation efficiency parameter $f_{\rm sed}$ encapsulating the effects of nucleation rates, material properties, mixing, and gravitational settling, all ultimately linked to the microphysical parameters of insoluble nuclei (Gao et al., 2018).

6. Broader Implications and Current Challenges

Direct linkage of dust size distributions to CCN activation reduces false alarms in light-rain forecasts and corrects for systematic biases in heavy-precipitation detection near desert boundaries (Arielly et al., 25 Jan 2026). Low-data-overhead models, requiring only a single parameter for effective CCN size, enable rapid operational forecasting and resource management in arid zones.

Ion-induced chemistry, catalytically producing condensable acids on cluster ions, may further enhance the supply of vapor for CCN growth, especially under elevated ionization (cosmic ray, gamma) conditions. This mechanism explains empirical findings of sustained CCN-number sensitivity at >50 nm to ionization rate, contrary to standard condensation–coagulation models (Svensmark et al., 2012).

Ongoing research seeks to constrain the variability of key parameters such as contact angle, to validate the year-round existence of fine-mode dust in all environments, and to extend these formulations to interactive dust lifting and coagulation schemes, especially for extraterrestrial atmospheres (Shaposhnikov et al., 2022, Gao et al., 2018).

The cloud–condensation–on–insoluble–nuclei paradigm now underpins parameterizations from terrestrial weather models to planetary and substellar atmosphere simulations, emphasizing the critical role of mineral dust and analogous insoluble particles in modulating cloud and precipitation properties across planetary environments.