Causes of ice age intensification across the Mid-Pleistocene Transition

During the Mid-Pleistocene Transition (MPT; 1,200-800 kya), Earth's orbitally paced ice age cycles intensified, lengthened from ∼40,000 (∼40 ky) to ∼100 ky, and became distinctly asymmetrical. Testing hypotheses that implicate changing atmospheric CO₂ levels as a driver of the MPT has proven difficult with available observations. Here, we use orbitally resolved, boron isotope CO₂ data to show that the glacial to interglacial CO₂ difference increased from ∼43 to ∼75 μatm across the MPT, mainly because of lower glacial CO₂ levels. Through carbon cycle modeling, we attribute this decline primarily to the initiation of substantive dust-borne iron fertilization of the Southern Ocean during peak glacial stages. We also observe a twofold steepening of the relationship between sea level and CO₂-related climate forcing that is suggestive of a change in the dynamics that govern ice sheet stability, such as that expected from the removal of subglacial regolith or interhemispheric ice sheet phase-locking. We argue that neither ice sheet dynamics nor CO₂ change in isolation can explain the MPT. Instead, we infer that the MPT was initiated by a change in ice sheet dynamics and that longer and deeper post-MPT ice ages were sustained by carbon cycle feedbacks related to dust fertilization of the Southern Ocean as a consequence of larger ice sheets.