The Equilibrium Temperature of Planets in Elliptical Orbits

TL;DR

This paper analytically demonstrates that the average equilibrium temperature of planets in elliptical orbits decreases slightly with eccentricity, impacting habitability assessments and the understanding of planetary climates.

Contribution

It provides new analytic solutions for orbital averages, revealing that equilibrium temperature decreases with eccentricity, challenging the common assumption of a direct increase.

Findings

Average equilibrium temperature decreases with eccentricity

Habitable zone limits are wider for elliptical orbits

Implications for climate and habitability models

Abstract

There exists a positive correlation between orbital eccentricity and the average stellar flux that planets receive from their parent star. Often, though, it is assumed that the average equilibrium temperature would correspondingly increase with eccentricity. Here we test this assumption by calculating and comparing analytic solutions for both the spatial and temporal averages of orbital distance, stellar flux, and equilibrium temperature. Our solutions show that the average equilibrium temperature of a planet, with a constant albedo, slowly decreases with eccentricity until converging to a value 90% that of a circular orbit. This might be the case for many types of planets (e.g., hot-jupiters); however, the actual equilibrium and surface temperature of planets also depend on orbital variations of albedo and greenhouse. Our results also have implications in understanding the climate, habitability and the occurrence of potential Earth-like planets. For instance, it helps explain why the limits of the habitable zone for planets in highly elliptical orbits are wider than expected from the mean flux approximation, as shown by climate models.

Figures

Tables

Table 1. Table 1: Calculated temporal averages of orbital distance ⟨ r ⟩ delimited-⟨⟩ 𝑟 \langle r\rangle , stellar flux ⟨ F ⟩ delimited-⟨⟩ 𝐹 \langle F\rangle , and equilibrium temperature ⟨ T e ​ q ⟩ delimited-⟨⟩ subscript 𝑇 𝑒 𝑞 \langle T_{eq}\rangle for some Solar System bodies.

	Literature Valuesa				Temporal Averagesb
Name	$a$ (AU)	$e$	$A$	$T_{eq}$ (K)	$⟨r⟩$ (AU)	$⟨F⟩$	$⟨T_{eq}⟩$ (K)
Mercury	0.38709927	0.20563593	0.068	439.6	0.395	6.819	439
Venus	0.72333566	0.00677672	0.77	226.6	0.723	1.911	227
Earth	1.00000261	0.01671123	0.306	254.0	1.000	1.000	254
Moon	1.00000261	0.01671123	0.11	270.4	1.000	1.000	271
Mars	1.52371034	0.0933941	0.250	209.8	1.530	0.433	210
Pluto	39.48211675	0.2488273	0.4	37.5	40.704	0.001	39

Table 2. Table 2: Estimated average equilibrium temperature ⟨ T e ​ q ⟩ delimited-⟨⟩ subscript 𝑇 𝑒 𝑞 \langle T_{eq}\rangle and effective thermal distance r T subscript 𝑟 𝑇 r_{T} for a modeled ocean planet with increasing eccentricity, but under a constant average stellar flux equal to Earth.

Planet Modela				Estimated Values
$e$	$a$ (AU)	$A$	$⟨T_s⟩$ (K)	$⟨T_{eq}⟩$b (K)	$r_T$c (AU)
0.0	1.000	0.250	267	259	1.000
0.2	1.003	0.295	263	253	1.008
0.4	1.045	0.350	255	242	1.067

Equations

$$F_{\star }=\frac{L_{\star }}{4\pi d^2}=\frac{\sigma R_{\star }^2{T}_{\star }^4}{d^2},$$

$$T_{eq}={\left[\frac{(1-A)F_{\star }}{4\beta ϵ\sigma }\right]}^{\frac{1}{4}}={\left[\frac{(1-A)L_{\star }}{16\beta ϵ\sigma \pi d^2}\right]}^{\frac{1}{4}}=T_{\star }\sqrt{\frac{R_{\star }}{2d}}{\left(\frac{1-A}{\beta ϵ}\right)}^{\frac{1}{4}},$$

$$F=\frac{L}{r^2},$$

$$T_{eq}=T_o{\left[\frac{(1-A)L}{\beta ϵr^2}\right]}^{\frac{1}{4}},$$

$$g=\frac{G}{\sigma T_s^4}=1-{\left(\frac{T_{eq}}{T_s}\right)}^4,$$

$$T_s=T_o{\left[\frac{(1-A)L}{\beta ϵ(1-g)r^2}\right]}^{\frac{1}{4}}.$$

$$⟨T_{eq}⟩\ne T_o{\left[\frac{(1-A)⟨F⟩}{\beta ϵ}\right]}^{\frac{1}{4}},$$

$$⟨T_{eq}⟩\ne T_o{\left[\frac{(1-A)L}{\beta ϵ⟨r{⟩}^2}\right]}^{\frac{1}{4}},$$

$$\overline{r}=\frac{1}{2\pi }{\int }_0^{2\pi }r\mathrm{d}f=a\sqrt{1-e^2},$$

$$\overline{F}=\frac{1}{2\pi }{\int }_0^{2\pi }F\mathrm{d}f=\frac{L(2+e^2)}{2a^2(1-e^2{)}^2},$$

$${\overline{T}}_{eq}=\frac{1}{2\pi }{\int }_0^{2\pi }{T}_{eq}\mathrm{d}f$$

$$\approx T_o{\left[\frac{(1-A)L}{\beta ϵa^2}\right]}^{\frac{1}{4}}\left[1+\frac{7}{16}{e}^2+\frac{337}{1024}{e}^4+O(e^6)\right],$$

$$⟨r⟩=\frac{1}{T}{\int }_0^{T}r\mathrm{d}t=a\left(1+\frac{e^2}{2}\right),$$

$$⟨F⟩=\frac{1}{T}{\int }_0^{T}F\mathrm{d}t=\frac{L}{a^2\sqrt{1-e^2}},$$

$$⟨T_{eq}⟩=\frac{1}{T}{\int }_0^T{T}_{eq}\mathrm{d}t$$

$$\approx T_o{\left[\frac{(1-A)L}{\beta ϵa^2}\right]}^{\frac{1}{4}}\left[1-\frac{1}{16}{e}^2-\frac{15}{1024}{e}^4+O(e^6)\right],$$

$$r_F=a(1-e^2{)}^{1/4}.$$

$$r_T$$

$$\approx a\left[1+\frac{1}{8}{e}^2+\frac{21}{512}{e}^4+O(e^6)\right].$$

$$⟨T_{eq}⟩=T_o{\left[\frac{(1-A)L}{\beta ϵa^2}\right]}^{\frac{1}{4}}\frac{2\sqrt{1+e}}{\pi }\mathbf{E}\left(\sqrt{\frac{2e}{1+e}}\right),$$

$$r_T=a{\left[\frac{2\sqrt{1+e}}{\pi }\mathbf{E}\left(\sqrt{\frac{2e}{1+e}}\right)\right]}^{-2},$$

Full text

stars: planetary systems, planets and satellites: fundamental parameters, planetary habitability, equilibrium temperature, habitable zone, eccentricity.

††slugcomment: Méndez, A. & Rivera-Valentín, E. G. (2017) ApJL, 837, L1

1 Introduction

The climate and potential habitability of planets depend on a complex interaction of many stellar and planetary properties (Schulze-Makuch et al., 2011). For example, Earth experiences a small annual change of 4 K in its average global surface temperature111Berkeley Earth: http://berkeleyearth.org/data/ mostly due to its obliquity and distribution of land and ocean areas and not its low eccentricity (i.e., $e>0.0167$); however, high eccentricities may have a strong effect on the climate of exoplanets or exomoons (e.g., stability of liquid water between seasonal extremes). This is a very important effect to consider since 75% of all exoplanets with known orbits have orbital eccentricities larger than Earth222NASA Exoplanet Archive: http://exoplanetarchive.ipac.caltech.edu/.

The effects of orbital eccentricity on planetary habitability is a complex problem that has been explored by one- and three-dimensional climate models (e.g., Williams & Pollard, 2002; Dressing et al., 2010; Dobrovolskis, 2013; Wang et al., 2014; Armstrong et al., 2014; Dobrovolskis, 2015; Linsenmeier et al., 2015; Bolmont et al., 2016; Way & Georgakarakos, 2017). The suitability for planets to maintain surface liquid water can be constrained by their orbital location as defined by the habitable zone (HZ), the circular region around a star in which liquid water could exist on the surface of a rocky planet (Huang, 1959; Dole, 1964; Kasting et al., 1993). Estimates on the occurrence of Earth-like worlds around stars $\eta_{\oplus}$ depend on recognizing those in the HZ (Traub, 2015).

The conservative inner edge of the HZ is given by the runaway of liquid water (i.e., surface temperatures above $\sim$340 K), corresponding to a rapid increase in stratospheric water vapor content, which leads to loss of an Earth-like ocean inventory within four to five billion years. The conservative outer edge is given by the maximum greenhouse necessary to keep temperatures above freezing (i.e., surface temperature above 273 K). These limits assume an Earth-like planet with N${}_{\text{2}}$-CO${}_{\text{2}}$-H${}_{\text{2}}$O atmospheres and supported by a carbonate-silicate cycle (Kasting et al., 1993). Further developments on the limits of the HZ for main-sequence stars are discussed elsewhere (e.g., Pierrehumbert & Gaidos, 2011; Joshi & Haberle, 2012; Yang et al., 2013; Shields et al., 2013; Kopparapu et al., 2013, 2014, 2016; Kitzmann, 2016, 2017). The HZ for circular orbits can be defined either in terms of orbital distance, stellar flux, or equilibrium temperature. Usually, a planetary albedo of 0.3 (similar to Earth or Jupiter) is used for the calculation of the equilibrium temperature as a basis for comparisons.

The equilibrium temperature is useful for comparing the temperature regimes of planets, but is not sufficient to draw conclusions on their climate (Leconte et al., 2013). Surface temperatures are controlled by a complex interaction of albedo and greenhouse effects. Planetary albedo can be estimated from equilibrium temperatures determined spectroscopically or from planet/star flux ratios (Cowan & Agol, 2011; Cowan et al., 2012). In practice, it might be easier to go the other way around and estimate equilibrium temperatures from planetary albedo using thermal and optical phase curve photometry, which requires less light gathering power than spectroscopy (Mallama, 2009; von Paris et al., 2016).

The equilibrium temperature $T_{eq}$ of the planet results from a balance between the incident stellar flux $F_{\star}$ from the host star and that absorbed by the planet (Selsis et al., 2007), where $F_{\star}$ and $T_{eq}$ are given from the Stefan-Boltzmann law

[TABLE]

[TABLE]

where $L_{\star}$, $R_{\star}$, $T_{\star}$ are the luminosity, radius, and the effective temperature of the star, respectively, $d$ is the distance of the planet from the star, $\epsilon$ is the broadband thermal emissivity (usually $\epsilon\approx 1$), and $\sigma$ is the Stefan-Boltzmann constant. The planetary albedo or bond albedo $A$ is the fraction of incident radiation, over all wavelengths, which is scattered by the combined effect of the surface and atmosphere of the planet. The factor $\beta$ is the fraction of the planet surface that re-radiates the absorbed flux, $\beta=1$ for fast rotators and $\beta\approx 0.5$ for tidally locked planets without oceans or atmospheres (Kaltenegger & Sasselov, 2011). Equation 2 ignores other energy sources such as formation energy, tidal deformation, and radiogenic decay. For simplicity, these equations can be expressed using solar and terrestrial units where the distance between the planet and the star, $r=d/a_{\earth}$, are in astronomical units (AU), and the stellar flux, $L=L_{\star}/L_{\sun}$, in solar units such that $F$ and $T_{eq}$ become

[TABLE]

[TABLE]

where $T_{o}=[L_{\sun}/(16\pi\sigma a_{\earth}^{2})]^{\frac{1}{4}}$ = 278.5 K (i.e., the equilibrium temperature of Earth for zero albedo). The derivations presented in this paper use solar and terrestrial units from now on as shown in equations 3 and 4.

Surface temperatures, though, depend on the equilibrium temperature and greenhouse effect of the planets. The normalized greenhouse effect $g$ can be used to connect the equilibrium and surface temperatures as

[TABLE]

where $G$ is the greenhouse effect or forcing (W/m2), and $T_{s}$ is the surface temperature (Raval & Ramanathan, 1989). The normalized greenhouse is very convenient since it summarizes a complex planetary property in a unitless quantity, as the planetary albedo does. Together equations 4 and 5 combine to give the surface temperature of a planet as

[TABLE]

The final surface temperature of a planet is controlled by its heat distribution $\beta$, emissivity $\epsilon$, bond albedo $A$, and greenhouse $g$; all properties that depend on combined surface and atmospheric properties. It is also important to clarify that the average of equation 6 is in fact the effective surface temperature $\langle T_{se}\rangle$ and not necessarily the actual surface temperature $\langle T_{s}\rangle$, depending on the spatial and temporal scales considered (i.e., global, daily, or annual averages). For example, the annual global average surface temperature of Mars is $\langle T_{s}\rangle\approx$ 202 K, but its effective surface temperature is $\langle T_{se}\rangle\approx$ 214 K (Haberle, 2013). This is an important distinction necessary to calculate greenhouse effects.

Equations 3 and 4 are more complicated for elliptical orbits. It is often assumed that the average equilibrium temperature $\langle T_{eq}\rangle$ of a planet in an elliptic orbit can be simply calculated from its average stellar flux $\langle F\rangle$ or distance $\langle r\rangle$. Leconte et al. (2013) recognized that the temporal average of the equilibrium temperature is lower than the equilibrium temperature computed from the temporally averaged flux, although they did not calculate it. Thus, the expressions

[TABLE]

[TABLE]

are incorrect interpretations that have been used extensively in the literature. For example, Kaltenegger & Sasselov (2011) (equation 6 in their paper) and Cowan & Agol (2011) (equation 9 in their paper) suggested a formulation equivalent to equation 7. Also Hinkel & Kane (2013) (table 2 in their paper) used this interpretation to calculate the equilibrium temperature for planets in elliptical orbits. Formulations similar to equations 7 or 8 might introduce large errors in estimates of the average equilibrium temperature of planets, especially for large eccentricities.

In this paper, we computed analytical solutions for the average orbital distance, stellar flux, and equilibrium temperature of planets in elliptical orbits. Orbital averages were computed with respect to spatial (i.e., geometric interpretation) and temporal (i.e., physical interpretation) coordinates for comparison purposes. In §2 and §3, we show correct formulations using both spatial and temporal averages for elliptical orbits, respectively. To our knowledge, we determined a new analytical solution for the average equilibrium temperature of planets in elliptical orbits. In §4, we apply our result to define an effective orbit for the HZ. The implications of our results are discussed in §5. Our study is motivated by planets in the HZ, but it is applicable to any planet in elliptical orbit with nearly constant albedo (e.g., Hot-Jupiters).

2 Spatial Averages for Elliptical Orbits

Spatial averages for orbital distance $\overline{r}$, stellar flux $\overline{F}$, and equilibrium temperature $\overline{T}_{eq}$ do not have a practical observational application, but were computed here for comparison purposes. They could be calculated with respect to the eccentric anomaly $E$ or true anomaly $f$ (orbital longitude). Equations 3 and 4 were integrated over the true anomaly (making the substitution $r=a(1-e^{2})/(1+e\cos{f})$) for a full orbit to obtain the corresponding averages,

[TABLE]

[TABLE]

[TABLE]

where $\mathbf{E}$ is the complete elliptic integral of the second kind (Weisstein, 2016; GSL, 2016) and $a$ is the semi-major axis. Equation 9 is sometimes used as the average distance of planets in elliptical orbits. This is a geometric average and the more physically meaningful quantity, the temporal average (see §3), should be used instead. The equations for stellar flux (equation 10) and equilibrium temperature (equation 11) are new results of this paper. Both increase with eccentricity while the average distance decreases (figure 1); however, these results are of little practical value since they are geometric averages.

3 Temporal Averages for Elliptical Orbits

Temporal averages for orbital distance $\langle r\rangle$, stellar flux $\langle F\rangle$, and equilibrium temperature $\langle T_{eq}\rangle$ were computed. They could be calculated with respect to time $t$ or the mean anomaly $M$. Equations 3 and 4 were integrated over time (making the substitutions $r=a(1-e\cos{E})$, $M=E-e\sin{E}$, and $M=(2\pi/T)t$) for a full orbital period $T$ to obtain the corresponding averages

[TABLE]

[TABLE]

[TABLE]

where $\mathbf{E}$ is the complete elliptic integral of the second kind (Weisstein, 2016; GSL, 2016) and $a$ is the semi-major axis. Equation 13 is the correct physical average distance of planets, but equation 9 is sometimes used instead. Equation 14 is a well known expression for the average stellar flux of elliptical orbits (Williams & Pollard, 2002). Equation 15 is a new result of this study. These equations show that the average distance ($a\leq\langle r\rangle<\frac{3}{2}a$) and stellar flux ($F|_{e=0}\leq\langle F\rangle<\infty$) increases with eccentricity while the average equilibrium temperature ($T_{eq}|_{e=0}\leq\langle T_{eq}\rangle<\frac{2\sqrt{2}}{\pi}T_{eq}|_{e=0}$) slowly decreases until converging to $\sim 90\%$ of the equilibrium temperature for circular orbits (figure 2).

4 The Habitable Zone

The HZ can be defined either in terms of distance, stellar flux, or equilibrium temperature for circular orbits, using equations 3 and 4 to convert from one another. This is not trivial for elliptical orbits because these quantities diverge differently with eccentricity. The usual approach is to compare the average stellar flux of the planet in an elliptical orbit (equation 14) with the stellar flux limits of the HZ, the so called ‘mean flux approximation’ (Bolmont et al., 2016). Equivalently, Barnes et al. (2008) suggested an HZ definition for elliptical orbits based on the mean flux approximation, but in terms of orbital distance. Their eccentric habitable zone (EHZ) is the range of orbits for which a planet receives as much flux over one orbit as a planet on a circular orbit in the HZ. This consists of comparing the semi-major axis of the planet with the HZ limits for circular orbits scaled by $(1-e^{2})^{-1/4}$.

Here we prefer to avoid redefining the limits of the HZ, which adds a complication for multiplanetary systems (i.e., each planet having a different HZ). Instead, we suggest to define some effective distance for planets in elliptical orbits, not necessarily equal to their average orbital distance as given by equations 9 or 13. Based on the mean flux approximation, the effective flux distance $r_{F}$, the equivalent circular orbit with the same average stellar flux $\langle F\rangle$ as the elliptic orbit is given by

[TABLE]

Alternatively, using equation 15 we define a new effective thermal distance $r_{T}$, the equivalent circular orbit with the same average equilibrium temperature $\langle T_{eq}\rangle$ as the elliptic orbit, given by

[TABLE]

The effective distance based on the mean flux approximation (equation 17) and our mean thermal approximation (equation 18) are not equivalent, one decreasing while the other increasing with eccentricity, respectively. The difference between these two solutions are less than 5% for eccentricities below 0.7, so our solution is a small correction that is more significant for high eccentric orbits. Nevertheless, we argue that the mean thermal approximation provides a better characterization of the limits of the HZ than the mean flux approximation for elliptical orbits since the HZ is constrained by the thermal regime of the planet. Therefore, we suggest that equation 18 is better to determine whether planets in elliptical orbits are within the HZ limits333VPL’s HZ Calculator: http://depts.washington.edu/naivpl/content/hz-calculator.

5 Discussion

Our temporal solution shows that the average equilibrium temperature of planets slightly decreases with eccentricity until converging to a value about 10% less than for circular orbits. This result might seem contradictory because the average stellar flux strongly increases with eccentricity; however, the average distance also increases with increasing eccentricity. The net effect is a small decrease in the equilibrium temperature with eccentricity as long as the planetary albedo stays constant throughout the orbit. That might be the case for bodies without atmospheres; however, the bond albedo of planets with condensables (e.g., oceans) in elliptic orbits might change significantly from the hot periastron (e.g., due to the absorption from water vapor) to the cold apastron (e.g., due to the reflectivity of ice).

We were able to reproduce well known expressions for the average distance and average stellar flux (equations 9, 13, and 14) among our analytic solutions. The complete elliptic integral of the second kind were calculated using the GNU Scientific Library function gsl_sf_ellint_Ecomp (GSL, 2016) and the Mathematica function EllipticE (note that the Mathematica implementation of this function requires the square of the argument). Our results were also validated with numerical solutions. The equations 3 and 4 were independently numerically integrated for a full orbit in IDL using the RK4 procedure, which is a fourth-order Runge-Kutta method. Both the analytic and numerical solutions agree as shown in figures 1 and 2.

As discussed in the introduction §1, there are different and incorrect interpretations in the literature of the equilibrium temperature of planets in elliptical orbits (equations 7 and 8). They might introduce large errors especially for large eccentricities (figure 3). The most common interpretation is one based on the average stellar flux (equation 7), also known as the ‘mean flux approximation’. This interpretation gives errors larger than 5% for the average equilibrium temperature of planets with eccentricities higher than 0.48, and over 10% for $e>0.64$. In fact, the interpretation based on the average distance (equation 8) is more consistent with our result (i.e., also decreases with eccentricity). In any case, average values should always be between the corresponding minimum orbital distance at periastron $r_{p}=a(1-e)$ and maximum at apastron $r_{a}=a(1+e)$.

A similar, but incorrect analytic solution to our average equilibrium temperature from equation 15 was proposed by Brasser et al. (2014) (equation 13 in their paper). Their solution suggests that the equilibrium temperature increases instead with eccentricity, but it is not self-consistent and for eccentricities below 0.8 gives values larger than the equilibrium temperature at periastron (i.e., the maximum value). For example, the solution of Brasser et al. (2014) gives 230 K for the average equilibrium temperature for Mars based on its orbital parameters and assuming a 0.25 bond albedo. This average value is outside the possible range of 200 K to 220 K between apastron and periastron, respectively. Our solution gives 210 K, which is consistent with the NASA Ames Mars General Circulation Model (Haberle, 2013). Our calculated values are also in agreement from those in the literature for other Solar System objects (table 1).

Unfortunately, previous studies on the climate of planets in elliptical orbits did not include equilibrium temperature estimates to directly compare with our results. Nevertheless, our results might explain previous conflicting studies, using surface temperature as a proxy. There are different ways to study the effect of eccentricity on the climate of planets (e.g., at constant semi-major axis, stellar flux, etc.). Bolmont et al. (2016) studied the effect of eccentricity on the climate of tidally-locked ocean planets using a Global Climate Model LMDz under constant stellar flux. They showed that the higher the eccentricity of the planet and the higher the luminosity of the star, the less reliable is the mean flux approximation. A similar conclusion was previously obtained by Dressing et al. (2010) for Earth-like planets in high eccentric orbits with a one-dimensional energy balance climate model (EBM). Our equilibrium temperature estimates positively correlate with the surface temperature results of Bolmont et al. (2016), which show that temperature decreases with eccentricity, contrary to the mean flux approximation (table 2). In this case, our results characterize the decreasing thermal response with eccentricity and corresponding increase in effective orbital distance.

Our results are more relevant for planets in very eccentric orbits $e>0.5$ and constant albedo (e.g., hot-jupiters). In any case, they seem better to describe the thermal regime of planets than the mean flux approximation (figure 3). Temperatures could increase with eccentricity as also shown by climate models (Williams & Pollard, 2002; Dressing et al., 2010). Planets with condensable species (e.g., Earth-like) could experience strong variations on their albedo and heat redistribution during their elliptical orbit. Kasting et al. (1993) showed that the planetary albedo not only depends on the surface and atmospheric properties of a planet, but it is also affected by both the stellar flux and the spectrum of the star. For example, it could go below 0.2 for Earth-like planets receiving over 20% the stellar flux that Earth receives around a Sun-like or cooler star. Heat redistribution is also more efficient in longer period tidally-locked planets (Bolmont et al., 2016). Future studies will consider these important dynamical corrections.

6 Conclusion

Here we determined analytic solutions of the spatial and temporal averages of orbital distance, stellar flux, and equilibrium temperature for planets in elliptical orbits and constant albedo. The solutions were validated (see §5) by reproducing known solutions with the same approach, obtaining equal results with numerical solutions, comparing with some Solar System bodies, and producing results consistent with previous climate models for planets in elliptic orbits. In particular, we determined that:

The average equilibrium temperature $\langle T_{eq}\rangle$ of planets in elliptical orbits slowly decreases with eccentricity $e$ until a converging value $\sim 90\%$ ($\frac{2\sqrt{2}}{\pi}$) the equilibrium temperature for circular orbits, assuming a constant albedo $A$ and heat redistribution $\beta$. The temperature is given by

[TABLE]

where $T_{o}$= 278.5 K, $L$ is the star luminosity (solar units), $a$ is the semi-major axis (AU), $\epsilon$ is the emissivity, and $\mathbf{E}$ is the complete elliptic integral of the second kind (Weisstein, 2016; GSL, 2016). 2. 2.

The potential stability of surface liquid water on rocky planets in elliptical orbits is better characterized by a mean thermal approximation rather than the mean flux approximation. Using observational data, a planet is inside the HZ if their effective thermal distance $r_{T}$, given by

[TABLE]

is within the HZ limits for circular orbits as defined by Kasting et al. (1993). In general, planets move farther outside the HZ (i.e., become colder) with increasing eccentricity if all other conditions stay constant; however, orbital variations of albedo, greenhouse, and heat redistribution might accentuate or invert this trend. 3. 3.

All climate models of planets should include at least global calculations of albedo $A$, normalized greenhouse $g$, equilibrium $T_{eq}$, and surface temperatures $T_{s}$ as a function of orbital time, not just orbital longitude. The distinction between the use of the surface temperature or the effective surface temperature should also be clear (Haberle, 2013). These global parameters help to compare different climate models and to translate them to simple parametric functions that are necessary to determine whether planets are in the HZ from observational data.

The results of this paper might be used to reevaluate the effective location in the HZ of many exoplanets with known eccentricity (using equation 21), especially those small enough to be considered potentially habitable. For example, the planet Proxima Cen b is probably in a circular orbit (Anglada-Escudé et al., 2016), but the outer two planets in the red-dwarf star Wolf 1061 are in eccentric orbits crossing the HZ (Wright et al., 2016; Kane et al., 2016). Also, many giant planets with elliptical orbits that cross the HZ might support habitable exomoons (Forgan & Kipping, 2013; Hinkel & Kane, 2013). Results for individual exoplanets are available and updated in the PHL’s Exoplanet Orbital Catalog444Exoplanet Orbital Catalog: http://phl.upr.edu/projects/habitable-exoplanets-catalog/catalog.

This work was supported by the Planetary Habitability Laboratory (PHL) and the Center for Research and Creative Endeavors (CIC) of the University of Puerto Rico at Arecibo (UPR Arecibo). Part of this study was originally presented in a poster at the 2014 STScI Spring Symposium: Habitable Worlds Across Space And Time (April 28 - May 1, 2014). This research has made use of the NASA Exoplanet Archive, which is operated by the California Institute of Technology, under contract with the National Aeronautics and Space Administration under the Exoplanet Exploration Program. We gratefully acknowledge the anonymous referee, whose comments improved the quality of the paper.