Intense glaciation during the end of Cryogenian time (~635 million years ago) marks the coldest climate state in Earth history - a time when glacial deposits accumulated at low, tropical paleolatitudes. The leading idea to explain these deposits, the snowball Earth hypothesis, predicts globally frozen surface conditions and subfreezing temperatures, with global climate models placing surface temperatures in the tropics between -20 °C and -60 °C. However, precise paleosurface temperatures based upon geologic constraints have remained elusive and the global severity of the glaciation undetermined. Here we make new geologic observations of tropical periglacial, aeolian and fluvial sedimentary structures formed during the end-Cryogenian, Marinoan glaciation in South Australia; these observations allow us to constrain ancient surface temperatures. We find periglacial sand wedges and associated deformation suggest that ground temperatures were sufficiently warm to allow for ductile deformation of a sandy regolith. The wide range of deformation structures likely indicate the presence of a paleoactive layer that penetrated 2-4 m below the ground surface. These observations, paired with a model of ground temperature forced by solar insolation, constrain the local mean annual surface temperature to within a few degrees of freezing. This temperature constraint matches well with our observations of fluvial deposits, which require temperatures sufficiently warm for surface runoff. Although this estimate coincides with one of the coldest near sea-level tropical temperatures in Earth history, if these structures represent peak Marinaon glacial conditions, they do not support the persistent deep freeze of the snowball Earth hypothesis. Rather, surface temperatures near 0 °C allow for regions of seasonal surface melting, atmosphere-ocean coupling and possible tropical refugia for early metazoans. If instead these structures formed during glacial onset or deglaciation, then they have implications for the timescale and character for the transition into or out of a snowball state.

The Tonian-Cryogenian Tambien Group of northern Ethiopia is a mixed carbonate-siliciclastic sequence that culminates in glacial deposits associated with the first of the Cryogenian glaciations-the Sturtian "Snowball Earth. " Tambien Group deposition occurred atop arc volcanics and volcaniclastics of the Tsaliet Group. New U-Pb isotope dilution-thermal ionization mass spectrometry (ID-TIMS) dates demonstrate that the transition between the Tsaliet and Tambien Groups occurred at ca. 820 Ma in western exposures and ca. 795 Ma in eastern exposures, which is consistent with west to east arc migration and deposition in an evolving back-arc basin. The presence of intercalated tuffs suitable for high-precision geochronology within the Tambien Group enable temporal constraints on stratigraphic data sets of the interval preceding, and leading into, the Sturtian glaciation. Recently discovered exposures of Sturtian glacial deposits and underlying Tambien Group strata in the Samre Fold-Thrust Belt present the opportunity to further utilize this unique association of tuffs and carbonate lithofacies. U-Pb ID-TIMS ages from zircons indicate that Tambien Group carbonates were deposited from ca. 820 Ma until 0-2 m.y. before the onset of the Sturtian glaciation, making the group host to a relatively complete carbonate stratigraphy leading into this glaciation. New δ13C and 87Sr/86Sr data and U-Pb ID-TIMS ages from the Tambien Group are used in conjunction with previously published isotopic and geochronologic data to construct newly time-calibrated composite Tonian carbon and strontium isotope curves. Tambien Group δ13C data and U-Pb ID-TIMS ages reveal that a pre-Sturtian sharp negative δ13C excursion (referred to as the Islay anomaly in the literature) precedes the Sturtian glaciation by-18 m.y., is synchronous in at least two separate basins, and is followed by a prolonged interval of positive δ13C values. The composite Tonian 87Sr/86Sr curve shows that, following an extended interval of low and relatively invariant values, inferred seawater 87Sr/86Sr rose ca. 880-770 Ma, then subsequently decreased leading up to the ca. 717 Ma initiation of the Sturtian glaciation. These data, when combined with a simple global weathering model and analyses of the timing and paleolatitude of large igneous province eruptions and arc accretion events, suggest that the 87Sr/86Sr increase was influenced by increased subaerial weathering of radiogenic lithologies as Rodinia rifted apart at low latitudes. The following 87Sr/86Sr decrease is consistent with enhanced subaerial weathering of arc lithologies accreting in the tropics over tens of millions of years, lowering pCO2 and contributing to the initiation of the Sturtian glaciation.

Carbonate mud represents one of the most important geochemical archives for reconstructing ancient climatic, environmental, and evolutionary change from the rock record. Mud also represents a major sink in the global carbon cycle. Yet, there remains no consensus about how and where carbonate mud is formed. Here, we present stable isotope and trace-element data from carbonate constituents in the Bahamas, including ooids, corals, foraminifera, and algae. We use geochemical fingerprinting to demonstrate that carbonate mud cannot be sourced from the abrasion and mixture of any combination of these macroscopic grains. Instead, an inverse Bayesian mixing model requires the presence of an additional aragonite source. We posit that this source represents a direct seawater precipitate. We use geological and geochemical data to show that "whitings" are unlikely to be the dominant source of this precipitate and, instead, present a model for mud precipitation on the bank margins that can explain the geographical distribution, clumped-isotope thermometry, and stable isotope signature of carbonate mud. Next, we address the enigma of why mud and ooids are so abundant in the Bahamas, yet so rare in the rest of the world: Mediterranean outflow feeds the Bahamas with the most alkaline waters in the modern ocean (>99.7th-percentile). Such high alkalinity appears to be a prerequisite for the nonskeletal carbonate factory because, when Mediterranean outflow was reduced in the Miocene, Bahamian carbonate export ceased for 3-million-years. Finally, we show how shutting off and turning on the shallow carbonate factory can send ripples through the global climate system.

Tides and Earth-Moon system evolution are coupled over geological time. Tidal energy dissipation on Earth slows Earth's rotation rate, increases obliquity, lunar orbit semi-major axis and eccentricity, and decreases lunar inclination. Tidal and core-mantle boundary dissipation within the Moon decrease inclination, eccentricity and semi-major axis. Here we integrate the Earth-Moon system backwards for 4.5 Ga with orbital dynamics and explicit ocean tide models that are "high-level" (i.e., not idealized). To account for uncertain plate tectonic histories, we employ Monte Carlo simulations, with tidal energy dissipation rates (normalized relative to astronomical forcing parameters) randomly selected from ocean tide simulations with modern ocean basin geometry and with 55, 116, and 252 Ma reconstructed basin paleogeometries. The normalized dissipation rates depend upon basin geometry and Earth's rotation rate. Faster Earth rotation generally yields lower normalized dissipation rates. The Monte Carlo results provide a spread of possible early values for the Earth-Moon system parameters. Of consequence for ocean circulation and climate, absolute (un-normalized) ocean tidal energy dissipation rates on the early Earth may have exceeded today's rate due to a closer Moon. Prior to ∼3Ga , evolution of inclination and eccentricity is dominated by tidal and core-mantle boundary dissipation within the Moon, which yield high lunar orbit inclinations in the early Earth-Moon system. A drawback for our results is that the semi-major axis does not collapse to near-zero values at 4.5 Ga, as indicated by most lunar formation models. Additional processes, missing from our current efforts, are discussed as topics for future investigation.

Geological evidence indicates that grounded ice sheets reached sea level at all latitudes during two long-lived Cryogenian (58 and ≥5 My) glaciations. Combined uranium-lead and rhenium-osmium dating suggests that the older (Sturtian) glacial onset and both terminations were globally synchronous. Geochemical data imply that CO₂ was 10² PAL (present atmospheric level) at the younger termination, consistent with a global ice cover. Sturtian glaciation followed breakup of a tropical supercontinent, and its onset coincided with the equatorial emplacement of a large igneous province. Modeling shows that the small thermal inertia of a globally frozen surface reverses the annual mean tropical atmospheric circulation, producing an equatorial desert and net snow and frost accumulation elsewhere. Oceanic ice thickens, forming a sea glacier that flows gravitationally toward the equator, sustained by the hydrologic cycle and by basal freezing and melting. Tropical ice sheets flow faster as CO₂ rises but lose mass and become sensitive to orbital changes. Equatorial dust accumulation engenders supraglacial oligotrophic meltwater ecosystems, favorable for cyanobacteria and certain eukaryotes. Meltwater flushing through cracks enables organic burial and submarine deposition of airborne volcanic ash. The subglacial ocean is turbulent and well mixed, in response to geothermal heating and heat loss through the ice cover, increasing with latitude. Terminal carbonate deposits, unique to Cryogenian glaciations, are products of intense weathering and ocean stratification. Whole-ocean warming and collapsing peripheral bulges allow marine coastal flooding to continue long after ice-sheet disappearance. The evolutionary legacy of Snowball Earth is perceptible in fossils and living organisms.