Geological evidence suggests that marine ice extended to the equator at least twice during the Neoproterozoic (∼750-635 Myr ago) 1, 2 , inspiring the Snowball Earth hypothesis that the Earth was globally ice covered 3, 4 . In a possible Snowball Earth climate, ocean circulation and mixing processes would have determined the melting and freezing rates that determine ice thickness 5,6 , influenced the survival of photosynthetic life 4,5,7-9 , and may provide important constraints for the interpretation of geochemical and sedimentological observations 4,10 . Here we show that in a Snowball Earth, the ocean would have been well mixed and characterized by a dynamic circulation 11 , with vigorous equatorial meridional overturning circulation, zonal equatorial jets, a well-developed eddy field, strong coastal upwelling and convective mixing. This is in contrast to the sluggish ocean often expected in a snowball scenario 3 .As a result of vigorous convective mixing, the ocean temperature, salinity and density were either uniform in the vertical or weakly stratified in a few locations. Our results are based on a coupled ice flow-ocean circulation model driven by a weak geothermal heat flux under a global ∼1 km ice cover. Compared to the modern ocean, the snowball ocean had far larger vertical mixing rates, and comparable horizontal mixing by ocean eddies. The strong circulation and coastal upwelling resulted in melting rates near continents as much as ten times larger than previously estimated 6, 7 . While we cannot resolve the ongoing debate on the existence of a global snowball ice cover 10,12,13 , we discuss the implications to nutrient supply for photosynthetic activity and to banded iron formations. The new insights and constraints on ocean dynamics may help resolve the Snowball Earth controversy when combined with future geochemical and geological
observations.The flow of thick ice over a Snowball ocean ("sea glaciers", characterized by very different dynamics from thinner sea ice 14 ), has received significant attention over the past few years [5][6][7]9,14,15 . Similarly, the role and dynamics of atmospheric circulation and heat transport, CO2 concentration, cloud feedbacks, and continental configuration have been studied [16][17][18] , as has the role of dust over the Snowball ice cover 17,19 . In contrast, despite its importance, the ocean circulation during Snowball events has received little attention. The few studies that used full ocean General Circulation Models (GCMs) concentrated mostly on the ocean role in Snowball initiation and aftermath 20,21 . None of these studied accounted for the combined effects of thick ice cover flow and driving by geothermal heating 11,13,22,23 , yet 11 simulated an ocean under a 200 m thick ice cover with no geothermal heat flux, and calculated a non steady--state solution with near--uniform temperature and salinity, and vanishing Eulerian circulation together with a strong parameterized eddy--induced high latitude circulation cells.To allow simulating the special circumstan...