Sea-ice thermodynamics can determine waterbelt scenarios for Snowball Earth

Abstract. Snowball Earth refers to multiple periods in the Neoproterozoic during which geological evidence indicates that the Earth was largely covered in ice. A Snowball Earth results from a runaway ice–albedo feedback, but there is an ongoing debate about how the feedback stopped: with fully ice-covered oceans or with a narrow strip of open water around the Equator. The latter states are called waterbelt states and are an attractive explanation for Snowball Earth events because they provide a refugium for the survival of photosynthetic aquatic life, while still explaining Neoproterozoic geology. Waterbelt states can be stabilized by bare sea ice in the subtropical desert regions, which lowers the surface albedo and stops the runaway ice–albedo feedback. However, the choice of sea-ice model in climate simulations significantly impacts snow cover on ice and, consequently, surface albedo. Here, we investigate the robustness of waterbelt states with respect to the thermodynamical representation of sea ice. We compare two thermodynamical sea-ice models, an idealized zero-layer Semtner model, in which sea ice is always in equilibrium with the atmosphere and ocean, and a three-layer Winton model that is more sophisticated and takes into account the heat capacity of ice. We deploy the global icosahedral non-hydrostatic atmospheric (ICON-A) model in an idealized aquaplanet setup and calculate a comprehensive set of simulations to determine the extent of the waterbelt hysteresis. We find that the thermodynamic representation of sea ice strongly influences snow cover on sea ice over the range of all simulated climate states. Including heat capacity by using the three-layer Winton model increases snow cover and enhances the ice–albedo feedback. The waterbelt hysteresis found for the zero-layer model disappears in the three-layer model, and no stable waterbelt states are found. This questions the relevance of a subtropical bare sea-ice region for waterbelt states and might help explain drastically varying model results on waterbelt states in the literature.