In order to understand the climate on terrestrial planets orbiting nearby Sun-like stars, one would like to know their thermal inertia. We use a global climate model to simulate the thermal phase variations of Earth-analogs and test whether these data could distinguish between planets with different heat storage and heat transport characteristics. In particular, we consider a temperate climate with polar ice caps (like modern Earth), and a snowball state where the oceans are globally covered in ice. We first quantitatively study the periodic radiative forcing from, and climatic response to, rotation, obliquity, and eccentricity. Orbital eccentricity and seasonal changes in albedo cause variations in the global-mean absorbed flux. The responses of the two climates to these global seasons indicate that the temperate planet has 3 times the bulk heat capacity of the snowball planet due to the presence of liquid water oceans. The temperate obliquity seasons are weaker than one would expect based on thermal inertia alone; this is due to cross-equatorial oceanic and atmospheric energy transport. Thermal inertia and cross-equatorial heat transport have qualitatively different effects on obliquity seasons, insofar as heat transport tends to reduce seasonal amplitude without inducing a phase lag. For an Earth-like planet, however, this effect is masked by the mixing of signals from low thermal inertia regions (sea ice and land) with that from high thermal inertia regions (oceans), which also produces a damped response with small phase lag. We then simulate thermal lightcurves as they would appear to a high-contrast imaging mission (TPF-I/Darwin) and consider the inverse problem of estimating thermal inertia based solely on time-resolved photometry. [Abridged]

Identifying liquid water on the surface of planets is a high priority, as this traditionally defines habitability. One proposed signature of oceans is specular reflection ("glint"), which increases the apparent albedo of a planet at crescent phases. We post-process a global climate model of an Earth-like planet to simulate reflected lightcurves. Significantly, we obtain glint-like phase variations even though we do not include specular reflection in our model. This false positive is the product of two generic properties: 1) for modest obliquities, a planet's poles receive less orbit-averaged stellar flux than its equator, so the poles are more likely to be covered in highly reflective snow and ice, and 2) we show that reflected light from a modest-obliquity planet at crescent phases probes higher latitudes than at gibbous phases, therefore a planet's apparent albedo will naturally increase at crescent phase. We suggest that this "latitude-albedo effect" will operate even for large obliquities: in that case the equator receives less orbit-averaged flux than the poles, and the equator is preferentially sampled at crescent phase. Using rotational and orbital color variations to map the surfaces of directly imaged planets and estimate their obliquity will therefore be a necessary pre-condition for properly interpreting their reflected phase variations. The latitude-albedo effect is a particularly convincing glint false positive for zero-obliquity planets, and such worlds are not amenable to latitudinal mapping. This effect severely limits the utility of specular reflection for detecting oceans on exoplanets.