Search

The prevailing paradigm for the climatic effects of large asteroid or comet impacts is a reduction in sunlight and significant short-term cooling caused by atmospheric aerosol loading. Here we show, using global climate model experiments, that the large increases in stratospheric water vapor that can occur upon impact with the ocean cause radiative forcings of over +20 W m−2 in the case of 10 km sized bolides. The result of such a positive forcing is rapid climatic warming, increased upper ocean stratification, and potentially disruption of upper ocean ecosystems. Since two thirds of the world's surface is ocean, we suggest that some bolide impacts may actually warm climate overall. For impacts producing both stratospheric water vapor and aerosol loading, radiative forcing by water vapor can reduce or even cancel out aerosol-induced cooling, potentially causing 1–2 decades of increased temperatures in both the upper ocean and on the land surface. Such a response, which depends on the ratio of aerosol to water vapor radiative forcing, is distinct from many previous scenarios for the climatic effects of large bolide impacts, which mostly account for cooling from aerosol loading. Finally, we discuss how water vapor forcing from bolide impacts may have contributed to two well-known phenomena: extinction across the Cretaceous/Paleogene boundary and the deglaciation of the Neoproterozoic snowball Earth.

The Snowball Earth bifurcation, or runaway ice-albedo feedback, is defined for particular boundary conditions by a critical CO2 and a critical sea-ice cover (SI), both of which are essential for evaluating hypotheses related to Neoproterozoic glaciations. Previous work has shown that the Snowball Earth bifurcation, denoted as (CO2, SI)*, differs greatly among climate models. Here, we study the effect of bare sea-ice albedo, sea-ice dynamics and ocean heat transport on (CO2, SI)* in the atmosphere–ocean general circulation model ECHAM5/MPI-OM with Marinoan (~ 635 Ma) continents and solar insolation (94% of modern). In its standard setup, ECHAM5/MPI-OM initiates a~Snowball Earth much more easily than other climate models at (CO2, SI)* ≈ (500 ppm, 55%). Replacing the model's standard bare sea-ice albedo of 0.75 by a much lower value of 0.45, we find (CO2, SI)* ≈ (204 ppm, 70%). This is consistent with previous work and results from net evaporation and local melting near the sea-ice margin. When we additionally disable sea-ice dynamics, we find that the Snowball Earth bifurcation can be pushed even closer to the equator and occurs at a hundred times lower CO2: (CO2, SI)* ≈ (2 ppm, 85%). Therefore, the simulation of sea-ice dynamics in ECHAM5/MPI-OM is a dominant determinant of its high critical CO2 for Snowball initiation relative to other models. Ocean heat transport has no effect on the critical sea-ice cover and only slightly decreases the critical CO2. For disabled sea-ice dynamics, the state with 85% sea-ice cover is stabilized by the Jormungand mechanism and shares characteristics with the Jormungand climate states. However, there is no indication of the Jormungand bifurcation and hysteresis in ECHAM5/MPI-OM. The state with 85% sea-ice cover therefore is a soft Snowball state rather than a true Jormungand state. Overall, our results demonstrate that differences in sea-ice dynamics schemes can be at least as important as differences in sea-ice albedo for causing the spread in climate models' estimates of the Snowball Earth bifurcation. A detailed understanding of Snowball Earth initiation therefore requires future research on sea-ice dynamics to determine which model's simulation is most realistic.

This work was supported by the NERC under Grant number NE/L001764/1. AJB is funded by NERC support of the Isotope Community Support Facility at SUERC. We are grateful to S. Mojzsis and I. Mukherjee for critical reviews, and A. Spencer and K. Mackay for help with sampling