Milankovitch Cycles for a Circumstellar Earth-analog within $\alpha$ Centauri-like Binaries

An Earth-analog orbiting within the habitable zone of $\alpha$ Centauri B was shown to undergo large variations in its obliquity, or axial tilt, which affects the planetary climate by altering the radiative flux for a given latitude (Quarles et al. 2019). We examine the potential implications of these obliquity variations for climate through Milankovitch cycles using an energy balance model with ice sheets. Similar to previous studies, the largest amplitude obliquity variations from spin-orbit resonances induce snowball states within the habitable zone, while moderate variations can allow for persistent ice caps or an ice belt. Particular outcomes for the global ice distribution can depend on the planetary orbit, obliquity, spin precession, binary orbit, and which star the Earth-analog orbits. An Earth-analog with an inclined orbits relative to the binary orbit can periodically transition through several global ice distribution states and risk runaway glaciation when periods of ice caps and an ice belt overlap. When determining the potential habitability for planets in stellar binaries, more care must be taken due to the orbital and spin dynamics.