Mixed-phase clouds (clouds that consist of both cloud droplets and ice crystals) are frequently present in the Earth’s atmosphere and influence the Earth’s energy budget through their radiative properties, which are highly dependent on the cloud water phase. In this study, the phase partitioning of cloud water is compared among six global climate models (GCMs) and with Cloud and Aerosol Lidar with Orthogonal Polarization retrievals. It is found that the GCMs predict vastly different distributions of cloud phase for a given temperature, and none of them are capable of reproducing the spatial distribution or magnitude of the observed phase partitioning. While some GCMs produced liquid water paths comparable to satellite observations, they all failed to preserve sufficient liquid water at mixed-phase cloud temperatures. Our results suggest that validating GCMs using only the vertically integrated water contents could lead to amplified differences in cloud radiative feedback. The sensitivity of the simulated cloud phase in GCMs to the choice of heterogeneous ice nucleation parameterization is also investigated. The response to a change in ice nucleation is quite different for each GCM, and the implementation of the same ice nucleation parameterization in all models does not reduce the spread in simulated phase among GCMs. The results suggest that processes subsequent to ice nucleation are at least as important in determining phase and should be the focus of future studies aimed at understanding and reducing differences among the models.