Exo-Milankovitch Cycles. I. Orbits and Rotation States (The Astronomical Journal, 2018)

The obliquity of the Earth, which controls our seasons, varies by only ~2fdg5 over ~40,000 years, and its eccentricity varies by only ~0.05 over 100,000 years. Nonetheless, these small variations influence Earth’s ice ages. For exoplanets, however, variations can be significantly larger. Previous studies of the habitability of moonless Earth-like exoplanets have found that high obliquities, high eccentricities, and dynamical variations can extend the outer edge of the habitable zone by preventing runaway glaciation (snowball states). We expand upon these studies by exploring the orbital dynamics with a semianalytic model that allows us to map broad regions of parameter space. We find that, in general, the largest drivers of obliquity variations are secular spin–orbit resonances. We show how the obliquity varies in several test cases, including Kepler-62 f, across a wide range of orbital and spin parameters. These obliquity variations, alongside orbital variations, will have a dramatic impact on the climates of such planets.

Continue Reading →

The Effect of Orbital Configuration on the Possible Climates and Habitability of Kepler-62f (Astrobiology, 2016)

As lower-mass stars often host multiple rocky planets, gravitational interactions among planets can have significant effects on climate and habitability over long timescales. Here we explore a specific case, Kepler-62f, a potentially habitable planet in a five-planet system with a K2V host star. N-body integrations reveal the stable range of initial eccentricities for Kepler-62f is 0.00⩽e⩽0.32, absent the effect of additional, undetected planets. We simulate the tidal evolution of Kepler-62f in this range and find that, for certain assumptions, the planet can be locked in a synchronous rotation state. Simulations using LMD Generic GCM indicate that with 3 bars of CO2 in its atmosphere, Kepler-62f would only be warm enough for surface liquid water at the upper limit of this eccentricity range, providing it has a high planetary obliquity (between 60∘ and 90∘). A climate similar to modern-day Earth is possible for the entire range of stable eccentricities if atmospheric CO2 is increased to 5-bar levels. In a low-CO2 case, simulations with CCSM4 and LMD Generic GCM indicate that increases in planetary obliquity and orbital eccentricity coupled with an orbital configuration that places the summer solstice at or near pericenter permit regions of the planet with above-freezing surface temperatures. This may melt ice sheets formed during colder seasons. If Kepler-62f is synchronously rotating and has an ocean, CO2 levels above 3 bars would be required to distribute enough heat to the night side of the planet to avoid atmospheric freeze-out and permit a large enough region of open water at the planet’s substellar point to remain stable. Overall, we find multiple plausible combinations of orbital and atmospheric properties that permit surface liquid water on Kepler-62f.

Continue Reading →

3D Modeling of GJ1214b’s Atmosphere: Vertical Mixing Driven by an Anti-Hadley Circulation (The Astrophysical Journal, 2015)

GJ1214b is a warm sub-Neptune transiting in front of a nearby M dwarf star. Recent observations indicate the presence of high and thick clouds or haze whose presence requires strong atmospheric mixing. In order to understand the transport and distribution of such clouds/haze, we study the atmospheric circulation and the vertical mixing of GJ1214b with a 3D General Circulation Model for cloud-free hydrogen-dominated atmospheres (metallicity of 1, 10, and 100 times the solar value) and for a water-dominated atmosphere. We analyze the effect of the atmospheric metallicity on the thermal structure and zonal winds. We also analyze the zonal mean meridional circulation and show that it corresponds to an anti-Hadley circulation in most of the atmosphere with upwelling at mid-latitude and downwelling at the equator on average.

Continue Reading →

3D Modeling of GJ1214b’s Atmosphere: Formation of Inhomogeneous High Clouds and Observational Implications (the Astrophysical Journal Letters, 2015)

The warm sub-Neptune GJ1214b has a featureless transit spectrum that may be due to the presence of high and thick clouds or haze. Here, we simulate the atmosphere of GJ1214b with a 3D General Circulation Model for cloudy hydrogen-dominated atmospheres, including cloud radiative effects. We show that the atmospheric circulation is strong enough to transport micrometric cloud particles to the upper atmosphere and generally leads to a minimum of cloud at the equator.

Continue Reading →