“…PlaSim has been used before to study snowball events on Earth‐like planets (e.g., Boschi et al, ; Linsenmeier et al, ; Lucarini et al, ) and is able to reproduce snowball phenomena found using other models by Pierrehumbert (), Abbot, Voigt, et al (), and Abbot (), such as increased variability in troposphere lapse rate, reduced tropopause height, weakened extratropical winter lapse rates, and warming of snowball states by cloud forcing. Models disagree on the extent of sea ice coverage during snowball events (e.g., Kirschvink, ; Lewis et al, ; Rodehacke et al, ; Yang et al, ), but in PlaSim, snowballs are characterized by complete sea ice cover, and exit from a snowball event begins when sea ice starts to melt, triggering an albedo feedback runaway. Sea ice thickness is limited to 9 m, but we do not restrict the depth of land ice and snow.…”