Interior structures and tidal heating in the TRAPPIST-1 planets
Key Words.:astrobiology – methods: numerical – planets and satellites: general – planets and satellites: interiors
Context:With seven planets, the TRAPPIST-1 system has among the largest number of exoplanets discovered in a single system so far. The system is of astrobiological interest, because three of its planets orbit in the habitable zone of the ultracool M dwarf.
Aims:We aim to determine interior structures for each planet and estimate the temperatures of their rock mantles due to a balance between tidal heating and convective heat transport to assess their habitability. We also aim to determine the precision in mass and radius necessary to determine the planets’ compositions.
Methods:Assuming the planets are composed of uniform-density noncompressible materials (iron, rock, HO), we determine possible compositional models and interior structures for each planet. We also construct a tidal heat generation model using a single uniform viscosity and rigidity based on each planet’s composition.
Results:The compositions for planets b, c, d, and e remain uncertain given the error bars on mass and radius. With the exception of TRAPPIST-1c, all have densities low enough to indicate the presence of significant HO. Planets b and c experience enough heating from planetary tides to maintain magma oceans in their rock mantles; planet c may have surface eruptions of silicate magma, potentially detectable with next-generation instrumentation. Tidal heat fluxes on planets d, e, and f are twenty times higher than Earth’s mean heat flow.
Conclusions:Planets d and e are the most likely to be habitable. Planet d avoids the runaway greenhouse state if its albedo is 0.3. Determining the planet’s masses within 0.1 – 0.5 Earth masses would confirm or rule out the presence of HO and/or iron. Understanding the geodynamics of ice-rich planets f, g, and h requires more sophisticated modeling that can self-consistently balance heat production and transport in both rock and ice layers.
The recent discovery of seven roughly Earth-sized planets orbiting the low-mass star TRAPPIST-1 (Gillon et al., 2017) has vaulted this system to the forefront of exoplanetary characterization. The planets orbit the star with semi-major axes AU, and orbital periods of a few Earth days. Given their proximity to the star, and the star’s low mass and luminosity, the surface of each planet has a moderate temperature, ranging from to 400 K (Gillon et al., 2017). Given the planets’ mean densities, these low temperatures suggest that some might have solid surfaces composed of HO ice and/or rock. The planets also have non-zero orbital eccentricities (Gillon et al., 2017; Wang et al., 2017), so that the tidal forces on the bodies vary with time, resulting in heating of their interiors by tidal dissipation (Luger et al., 2017). Tidal heating is an important energy source in the satellites of the outer planets in our solar system (e.g., McEwen et al., 2004; Greeley et al., 2004; Spencer et al., 2009) and could significantly enhance the habitable volume in the galaxy by warming exoplanets and their satellites (exomoons) (Dobos and Turner, 2015).
Here, “tidal dissipation” (or tidal heating) refers to the process of the dissipation of orbital energy in the interior of a solid body (a moon or a planet). For a single secondary object, that is, a lone moon orbiting a planet, or a lone planet orbiting a star, tidal dissipation rapidly decreases the orbital eccentricity of the secondary, and then ceases when (Murray and Dermott, 1999).
If multiple objects orbit the primary, and the orbiting objects occupy a mean motion resonance, periodic gravitational perturbations will help maintain non-zero eccentricities. This is the case, for example, for the satellites of Jupiter and Saturn. Jupiter’s innermost moon, Io, has a tidal heat flux about 1–2 W/m (Spencer et al., 2000; Veeder et al., 2004). The heat flux is so large that the moon has become volcanically active, erupting silicate magmas and sulfur-rich compounds onto its surface (McEwen et al., 2004). Tidal heating can also be an important heat source in bodies with a significant surface layer of HO (whether it be solid ice or liquid water). Saturn’s icy moon Enceladus, which has a mean radius of 252 km (Porco et al., 2006) and is thus only asteroid-sized, has plumes of water-rich material erupting from tidally induced cracks at its south pole (Porco et al., 2006), and a localized heat flow between 0.1 and 0.25 W/m, roughly three times that of the Earth (Spencer et al., 2009). Tidal deformation and the geodynamical consequences of tidal heating have also created the bizarre geology on the surfaces of Jupiter’s moons Europa (Greeley et al., 2004) and Ganymede (Pappalardo et al., 2004).
The TRAPPIST-1 planets occupy a mean motion resonance chain with orbital period ratios slightly larger than 24 : 15 : 9 : 6 : 4 : 3 : 2 (Unterborn et al., 2017; Luger et al., 2017). The resonances maintain the orbital eccentricities of the planets resulting in long-term tidal heating (Luger et al., 2017). Luger et al. (2017) find that assuming Earth-like tidal dissipation rates, TRAPPIST-1b may have a tidal heat flux larger than Io () with peaks larger than due to the periodic change of its orbital eccentricity. The other six planets in the system would experience weaker tidal heating, with heat fluxes ranging between values similar to the Earth and Io for planets c–e, and less than Earth for planets f–h.
Here, we constructed simple compositional and thermal models of the interiors of each of the TRAPPIST-1 planets. Assuming that the planets are composed of uniform-density (incompressible) rock, iron, and HO, we explored the range of possible compositional and interior models permitted by the error bars in mass and radius for each planet. Because the error bars on mass and radius are so large, more sophisticated models including, for example, compression and thermal evolution, are unlikely to give robust and meaningful results. Our calculations allowed us to determine the precision with which the mass of each planet must be measured to determine whether a given planet harbors a significant reservoir of HO. We computed the amount of tidal dissipation in each object based on a uniform viscosity and rigidity, which depend on the volume fractions of ice and rock in each planet. We also determined the expected temperature in the rock mantles of planets b, c, d, e, and f by setting the tidal heat flux equal to the amount of heat that could be removed by solid-state convection. Some of the planets have globally averaged heat fluxes greater than the heat flux that can trigger a runaway greenhouse effect, which causes irreversible water evaporation and loss (Kasting et al., 1993).
2 Interior structures
2.1 Prior work
Table 1 summarizes the values of mass, radius, mean density, effective temperature, and orbital period for each of the seven TRAPPIST-1 planets. The vast majority of the values originate from Wang et al. (2017). The masses range from almost 0.1 to 1.63 , where kg is the mass of the Earth. The radii are generally close to the radius of the Earth ( km).
Several studies have constrained interior states for the TRAPPIST-1 planets using PREM (Preliminary Earth Reference Model), a density and seismic velocity profile for the Earth (Dziewonski and Anderson, 1981). Zeng et al. (2016) use an equation of state for silicate material based on PREM to constrain iron contents and core sizes in rocky exoplanets as a function of their size and temperature. These results have been used by Gillon et al. (2017) and Wang et al. (2017) to constrain the rock/metal fractions in the TRAPPIST-1 planets. Wang et al. (2017) combine the original discovery data for the TRAPPIST-1 system with 73.6 days of K2 observations to obtain updated masses and radii for the planets, which serve as the primary source of mass and radius information for our study. Wang et al. (2017) also use equations of state for rock and metal based on PREM to construct interior models of the planets, demonstrating that several of the bodies have densities low enough to contain a significant amount of water in their interiors.
Another method is to use equations of state for silicate minerals (e.g., MgSiO) to calculate density and pressure as a function of depth, constrained by the mass and radius of the planet (Unterborn et al., 2017). This detailed analysis has been performed for TRAPPIST-1b and 1c. In these models, phase transitions from one polymorph of MgSiO to another are considered, resulting in a mantle which contains layers of enstatite, pyroxene, and garnet, before a final transition to Bridgmanite at a depth of km. Each of these phases have much lower densities than compressed Bridgmanite, leading to a low-density upper mantle, and permitting a relatively thin HO layer.
Here, we have calculated the simplest possible interior structures for each of the TRAPPIST-1 planets assuming that the rock, metal, and ice layers in their interiors have uniform densities. We calculated the suite of structures that fit the mean mass and radius, and explore the range of possible interior structures based on the uncertainties of planetary masses and radii. Given the large error bars on the masses and radii of the planets, some of which permit planetary mass of zero, simple uniform-density models provide interior structures that are just as valid as those taking into account the pressure and temperature changes with depth inside the bodies (see e.g., Unterborn et al., 2017).
Table 1 shows the mean densities for the planets that we calculated based on their maximum and minimum permitted masses and radii; the large error bars on these quantities yield large error bars on the mean densities, as well. Despite the uncertainties, the densities of the TRAPPIST-1 planets show quite a bit of variation, ranging from values close to the density of water ice under compression, to densities hinting at rocky, and possibly iron-rich compositions.
|Planet||Mass (M||Radius||(kg/m)||T (K)||Orbital Period (days)||Eccentricity|
Our choices of densities for rock, iron, and ice include the effects of compression and phase changes that will occur at depth inside Earth-sized planets. For iron, we assumed a constant density of kg/m, consistent with the density of Earth’s inner core (Dziewonski and Anderson, 1981; Unterborn et al., 2017). For rock in planets b, c, d, and g we assumed a constant density of kg/m, appropriate for compressed Bridgmanite (MgSiO) (McCammon et al., 2016; Unterborn et al., 2017). For HO at low pressures, we assumed a density of kg/m, which accounts for the possibility of both a solid ice I shell and a liquid water ocean. For HO at high pressures, above 209 MPa, we assumed a density of kg/m, which accounts for the possible presence of high-pressure ice polymorphs (hpp) II through VII (Hobbs, 1974), each of which could be present in distinct layers within the planet. When high-pressure ice phases are present, it is likely that ices VI and VII will dominate because they are stable over a broad range of temperatures at the tens of GPa central pressures appropriate for TRAPPIST-1 bodies.
For planets e, f, and h, the minimum permitted planetary densities based on the masses and radii from Wang et al. (2017) are smaller than the density of high-pressure HO polymorphs. In the context of our assumptions, a planet with would correspond to a planet composed of pure HO, which we consider unlikely based on planet and satellite formation models (e.g., Canup and Ward 2002; Morbidelli et al. 2012; Johansen et al. 2015). In Table 1, we reported a minimum mass corresponding to the mass of a planet with the maximum allowed radius, and .
When determining the interior structures in the low-density limit for planets e and f, and for all of the interior structures for planet h, we adopted a density for rock kg/m, the density of Prinn-Fegley rock. Prinn-Fegley rock has been used as a representative composition for the rocky component of the outer planet satellites (Mueller and McKinnon, 1988; Barr and Canup, 2008). This could represent a more hydrated, but compressed rock core.
With the exception of TRAPPIST-1c, all of the seven planets have densities low enough to indicate the presence of HO (Wang et al., 2017; Unterborn et al., 2017). For these bodies, we began by estimating the maximum depth of the ice I/liquid water layer: , where is the surface gravity, and MPa is the maximum depth at which ice I is possible (Hobbs, 1974). Subtracting the mass of ice I, we determined the mean density of the remaining planetary material, which allowed us to constrain the relative sizes of the mantles of high pressure ice polymorphs and rock, and the radius of the iron core. Interior models applied by Wang et al. (2017) suggest that TRAPPIST-1b may contain approximately 25% water by mass (with a large uncertainty between about 0 and 60%). Using this constraint of 25% mass fraction, we were able to obtain a unique interior structure for TRAPPIST-1b that satisfied the mass and radius constraints.
3 Thermal model
Here, we have computed the heat production by tides and heat lost by conduction and convection to constrain tidal heating rates and rock mantle temperatures in each of TRAPPIST-1 exoplanets. The most common method of estimating tidal heating is to assume that the interior of the secondary object (in this case, the planet) has a Maxwell viscoelastic rheology (Love, 1906; Peale and Cassen, 1978), characterized by a shear modulus, , and viscosity, . More complex rheologies (for example, Andrade or a standard linear solid) can also be used, and yield slightly different behaviors (e.g., Efroimsky, 2012). However, given the large uncertainties in the planets’ compositions, we chose the Maxwell model because it is a simple model with simple behavior and relatively few parameters.
In a Maxwell viscoelastic material, heating is maximized in material whose Maxwell time, is comparable to the orbital period of the secondary, (Love, 1906). For ice I, the low-density phase of HO that floats on liquid water, the Maxwell time is comparable to the orbital period of Europa, 85 hours, if the ice is close to its melting point (Barr and Showman, 2009).
In a tidally heated body, the interior temperature depends upon the balance and feedback between internal heat generation by, for example, radiogenic and tidal sources, and heat transported by conduction and possibly solid-state convection (Schubert et al., 2001; Tackley, 2001; Hussmann et al., 2002; Moore, 2006). A body whose interior, or portion of its interior, has a viscosity and rigidity close to the values for which , will undergo intense tidal heating. A body with a slightly higher viscosity, that is, one that is cooler, can be warmed by tidal heating until . The body can experience inner melting, which in silicate bodies leads to volcanism (Peale, 2003), and in ice/rock bodies, leads to the maintenance of liquid water oceans (Ojakangas and Stevenson, 1989). Thus, the tidal heating rate in a solid body and its interior structure (namely, composition, temperature, and phase as a function of depth), depend on one another, and the body may reach an equilibrium structure where heat generation is equal to heat loss.
Figure 1 shows a schematic illustration of the methods we used to calculate the internal temperature inside each of the TRAPPIST-1 planets based on their different interior structures. Given values for the volume fractions of each of the planets’ constituents (), we constructed a volume-averaged viscosity () and rigidity (). Based on these average values, we calculated the response of the planet to the tidal forcing, expressed by the imaginary part of the Love number, . This value was used to calculate the surface heat flux from tidal dissipation . The heat generated from tidal dissipation is likely removed from the planets’ deep interiors via solid-state convection, the dominant heat transfer mechanism in the terrestrial planets and ice/rock satellites in our Solar System. To determine the temperature in the planets’ deep interiors, we balanced and the convective heat flux (), and found the value of temperature () for which . If is above the solidus temperature for rock, the planet may possess a magma ocean in its interior. The globally averaged energy flux on the planet is the sum of the energy flux from the mean stellar irradiation () and tidal heating,
where is taken at the top of the planet’s atmosphere. The stellar radiation can be expressed as
where is the stellar luminosity, is the Bond albedo of the planet, and are the semi-major axis and the eccentricity of the planet’s orbit, respectively. Greenhouse and cloud feedbacks were neglected in the calculation of . The albedo of each of the planets was set to 0.3, an Earth-like value.
3.1 Heat generation
Here, we have focused on tidal heating as the dominant source of energy for heating the interiors of the TRAPPIST-1 planets. This is because the magnitude of tidal heating expected in rock- or ice-dominated bodies at the orbital periods associated with the TRAPPIST-1 planets, a few to ten watts per meter squared, is several orders of magnitude higher than the heat flows associated with radiogenic heating, which is typically of the order of tens of milliwatts per meter squared (Henning et al., 2009; Henning and Hurford, 2014). Residual heat of accretion is also a possible heat source in any solid planet. The energy per unit mass released during accretion, is sufficient to melt planets as large as the TRAPPIST-1 bodies, and is still a significant contribution to the Earth’s heat flux. Even if the TRAPPIST-1 bodies still retained 20% of their accretional heat, the resulting heat flux would still be much smaller than the tidal heat flux.
The tidal heat flux was calculated using a viscoelastic model for a homogeneous body described by Dobos and Turner (2015), which was originally developed by Henning et al. (2009). Although the planet was assumed to be homogeneous, we mimicked the effect of multiple materials by determining the effective viscosity and rigidity for each layers and weighting them by their volume fractions. To approximate the material properties of a planet composed of several different materials, we calculated the volume fraction of each material contained in the planet: for a given material, , the volume fraction , where is the total volume of the planet. We calculated values of the volume fractions of ice I (), high-pressure ice polymorphs (), rock (), and iron () contained in each planet, using the mass and radius of the planet as constraints, assuming uniform densities for each of the materials. A single uniform viscosity and rigidity for the planet was approximated by
where is the viscosity of the material. We used a similar relationship to construct a single value of the shear modulus () that approximates the behavior of the entire planet:
The total amount of energy produced by tidal dissipation in a synchronously rotating body (Segatz et al., 1988),
where is the orbital frequency of the planet, is the planet’s orbital period, is the gravitational constant, is the planet’s eccentricity, and is its radius. The quantity is the imaginary part of the Love number, which describes how the planet’s gravitational potential is disturbed by its tidal distortion. The value of depends on the structure and rheology in the body (Segatz et al., 1988).
There are a variety of different ways of calculating for a solid planet. In some studies, is expressed as , where is a constant parameter describing the fraction of orbital energy dissipated per tidal cycle (see, e.g., Murray and Dermott 1999 for discussion, and Papaloizou et al. 2017 for example application of this approach to the TRAPPIST-1 system). The value of for solid planetary bodies is not well-constrained, and is generally thought to be between and 200 (Goldreich and Soter, 1966; Murray and Dermott, 1999).
A constant- approach, however, does not take into account the feedback between the thermal state of a planet and its response to tidal forces. A warm planet will deform more in response to tidal forces than a cold, rigid planet. Thus, it is more accurate to calculate based on the viscosity and rigidity of the planet’s interior. A common approach is to assume that the planet is composed of a single material (for example, rock). Given the density, shear modulus, and viscosity of the planetary material, and an assumption about the rheology (for example Maxwell, Andrade, etc., Efroimsky (2012)), it is possible to derive an analytic expression for (e.g., Harrison 1963; Henning et al. 2009; Dobos and Turner 2015, see Eq. 6 here).
A more sophisticated approach is to assume that the planet is composed of, perhaps two or three materials, each of which has a different density, shear modulus, and viscosity. In this case, analytic, although quite cumbersome expressions can still be obtained for (Harrison, 1963; Remus et al., 2012; Beuthe, 2013; Remus et al., 2015). Many of these studies rely upon simplifying assumptions to obtain closed-form analytic solutions (e.g., Remus et al. 2012, 2015); in most cases, these simplifying assumptions may not apply to the TRAPPIST-1 planets. Another approach is to calculate numerically, by representing the planet by ten or more layers, each of which has a density, shear modulus, and viscosity dictated by the composition of the layer, but also thermal state of the material (e.g., Sabadini et al. 1982; Moore 2006; Roberts and Nimmo 2008).
Given the large uncertainties in the masses and radii of the TRAPPIST-1 planets, calculation of interior thermal states using more sophisticated models (for example, one in which the heat generation and transport in each different compositional layer inside the planet is treated separately) would be premature. Therefore, we used a Maxwell viscoelastic rheology to calculate the value of Im():
The globally averaged tidal heat flux was calculated from
It is worth noting that this tidal heating model is valid only for small orbital eccentricities, that is, for . At higher eccentricities, additional higher-order terms must be included in the tidal potential, which can increase the amount of heating (Mignard, 1980). Fortunately, this constraint is fulfilled for each of the seven planets, according to the data of Gillon et al. (2017) and Wang et al. (2017). In our model, we did not consider gravitational interactions of the planets with each other, only the star-planet (two-body) tidal interaction was considered in each case.
3.2 Heat loss
The tidal heat generated inside the planet was assumed to be transported to the surface by solid-state convection, then conducted across a surface boundary layer where it is radiated into space. The heat flux across a convecting mantle,
where is the thermal conductivity, is the difference in temperature between the hot, convecting mantle () and is the temperature at the surface, and is the thickness of the convecting layer. The parameter , the Nusselt number, describes the relative efficiency of convective heat transport versus conductive heat transport. In an internally heated planet with a Newtonian rheology, is related to the thermal, physical, and rheological properties of the mantle (Solomatov and Moresi, 2000),
where , is the activation energy in the flow law for the convecting material and J/mol/K is the universal gas constant. The parameter is the Rayleigh number, which describes the vigor of convection,
where the coefficient of thermal expansion , thermal diffusivity , , and is the viscosity evaluated at the mantle temperature. Using Eq. 10 to evaluate the Rayleigh number in Eq. 9, we obtained an expression for the convective heat flux which did not depend on the thickness of the convecting layer, or the temperature at the surface of the planet (Solomatov and Moresi, 2000; Barr, 2008),
We assumed that the majority of tidal heat will be produced in, and transported within, the rock mantle of the planet. Thus, the parameters in Eq. 11 were evaluated using the properties of rock: kJ/mol (see Sect. 3.4.1), W/m/K, 1/K, and J/kg (Solomatov and Moresi, 2000).
To determine the equilibrium (or equilibria) between and , we evaluated Eq. 7 and Eq. 11 for a variety of values of , the temperature in the rock mantle of each planet. Two equilibria are possible between the tidal heat flux and convective cooling flux which correspond to two equilibrium temperatures of , which we referred to as (Moore, 2006; Dobos and Turner, 2015). The left panel of Fig. 2 illustrates schematically how the tidal and convective heat fluxes vary as a function of temperature. One equilibrium can exist at a temperature well below the solidus, corresponding to a balance between tidal heat generation and transport in a purely solid planetary mantle (point (A) in the left panel, and also see the middle panel of Fig. 2). In this case, the amount of tidal heat generated is sufficiently low to be removed by convection, even if the convective flow is relatively sluggish. This equilibrium is unstable – if the tidal heat generation increases, convection cannot remove the additional tidal heat (Moore, 2003; Dobos and Turner, 2015). Depending on the tidal forcing frequency and the rheology of the planet, a second equilibrium can be achieved in which the planet is partially molten (point (B) in the left panel, and the right panel of Fig. 2). Tidal heat generation decreases as a function of temperature and melt fraction above the solidus because partially molten rock is less dissipative than warm, solid rock. However, the convective heat flux increases sharply as the presence of melt decreases the mantle viscosity. This high-temperature equilibrium is stable (Moore, 2003; Dobos and Turner, 2015) – as the tidal heating rate increases, the viscosity of the mantle decreases, permitting more efficient convective heat flow and resulting in cooling of the mantle.
3.3 Runaway greenhouse limit
To characterize the habitability of the planets, we determined whether the planet might enter a runaway greenhouse state. When the globally averaged heat flux on a planet exceeds a critical value, then surface water slowly evaporates. The water molecules in the upper atmosphere disintegrate to hydrogen and oxygen due to photodissociation by ultraviolet photons, and then hydrogen molecules can easily escape to space. Hence, entering the runaway greenhouse phase could imply that water molecules would not be able to re-form later, even if the temperature decreases (Kasting et al., 1993).
To determine the runaway greenhouse limit, we used the formulation of Pierrehumbert (2010, Chapter 4). If the globally averaged energy flux on the planet from stellar irradiation and tidal heat () exceeded the minimum energy flux on top of a water-rich atmosphere then it triggers the runaway greenhouse effect:
where , is the Stefan-Boltzmann constant, J/kg is the latent heat, J/(kg K) is the gas constant for water vapor, is the gray absorption coefficient, Pa is the reference pressure for absorption, , Pa and K.
3.4 Viscoelastic material properties
For rock, we used the viscosity and rigidity model used for exomoon interiors described in Dobos and Turner (2015) which is based on the works of Henning et al. (2009), Moore (2003), and Fischer and Spohn (1990). The behavior of the rock changes as a function of temperature, as the temperature crosses the solidus point ( K), across a “breakdown” point ( K) at which the volume fraction of solid rock crystals is equal to the volume fraction of melted rock. At this point, the behavior of the crystal and melt mixture begins to behave more like a liquid (Renner et al., 2000). When the rock is completely melted, (above the liquidus (, where K) its viscosity and shear modulus are quite small. For temperatures below the solidus, GPa, and the viscosity is described by,
For temperatures above the solidus, but below the breakdown temperature (), the shear modulus depends on temperature (Fischer and Spohn, 1990),
where K and . In this regime, the viscosity of rock depends on the volume fraction of melt present (Moore, 2003), , which varies linearly with temperature between the liquidus and solidus,
where is the melt fraction coefficient (, Moore, 2003). Above the breakdown temperature, but still below the liquidus, the shear modulus, now controlled by the behavior of the volumetrically dominant silicate melt, is a constant at Pa. The viscosity follows the Einstein-Roscoe relationship (Moore, 2003),
Above the liquidus, Pa s, and the viscosity continues to depend on temperature, obeying Eq. (16) with .
3.4.2 Ice I and liquid water
where the molar volume m/mol, diffusion coefficient , and activation energy kJ/mol. Grain sizes for ice in natural systems on Earth are limited to a few millimeters, due to the presence of silicate particles which inhibit grain growth (Durand et al., 2006). Temperature and strain rate conditions in terrestrial glaciers are somewhat similar to those expected in the outer planet satellites in our Solar System, leading to the suggestion that mm in the ice I mantles of those bodies (Barr and McKinnon, 2007). In the absence of other information, it is reasonable to assume that the process would work in a similar fashion on other planets: we adopted mm as a nominal value. Similar to rock, the shear modulus for ice I is weakly temperature-dependent (Gammon et al., 1983),
where Pa, K.
The melting temperature of ice I is K. Above the melting point, , and Pa s which we used for the liquid water layers.
3.4.3 High-pressure ice polymorphs (HPPs)
Experimental data for the ductile behavior of the high-pressure ice polymorphs are more limited than experimental studies for ice I. It is known that each of the high-pressure polymorphs has a distinct flow law, and many of them are strongly non-Newtonian, so that the viscosity depends strongly on stress (Durham et al., 1997b). Similar behavior is observed in rock far from the melting point (e.g., Karato et al. 1995). We used a flow law that approximates the laboratory-determined flow behavior for ices VI and VII from Durham et al. (1997b, a). Our tidal model assumed a Maxwell behavior, in which viscosity does not depend on stress, so we assumed a constant stress of 0.3 MPa. The approximate hpp flow law,
where and kJ/mol, gives a viscosity at the approximate melting point of ices VI and VII, K, of Pa s, somewhat similar to nominal near-melting-point viscosities for ice I. Data for the shear modulus of ices VI and VII are also limited, but at temperatures close to the melting point, Pa (Shimizu et al., 1996), where this value represents an average over the range for which experimental exist, which roughly corresponds to the pressure conditions expected within the TRAPPIST-1 planets’ interiors. We assumed was independent of temperature below the melting point. Above the melting point, high-pressure ices take on the properties of liquid water: , and Pa s.
The viscosity and rigidity of iron were assumed to be constant, and equal to the values deduced for the inner core of the Earth. The rigidity of iron, Pa, based on PREM (Dziewonski and Anderson, 1981). The viscosity of iron in Earth’s inner core is thought to lie between Pa s (Koot and Dumberry, 2011); we used Pa s.
To test the validity of our approach, we first used our methods to calculate the equilibrium interior temperature and tidal heat flux from Jupiter’s moon Io, for which the tidal heat flux has been measured by spacecraft to be W/m (Spencer et al., 2000; Veeder et al., 2004). We used the following parameters for Io: mass kg, rock volume fraction , iron fraction (Schubert et al., 2004), orbital period days, eccentricity , and km. With these values, our model yielded a tidal heat flux W/m, comparable in magnitude to that observed, and an equilibrium rock mantle temperature K, above the rock solidus, and consistent with the inference of a partially molten mantle beneath Io’s surface (Khurana et al., 2011).
Because the uncertainties in the measured planetary radius and mass are quite large (in some cases, they permit zero mass), it is important to provide limits to our calculations, as well. Here, we have reported a range of interior structures for the lowest and highest densities that can be estimated from the minimum and maximum radii and masses (see Table 1). In the cases of planets e, f, and h, we have modified some values in order to avoid the possibility of zero mass and the lack of a rocky core. For these planets, when determining the interior structure for the low-density case (and also for the other density values in the case of planet h), we assumed the rock to be Prinn-Fegley rock, as described in Sect. 2.2. These calculations with the lowest and highest densities effectively “bracket” the range of possible interior structures.
With these new structures, we recalculated the tidal fluxes, as well. Also, the uncertainties of the orbital eccentricities propagate into the tidal heating calculation. For this reason, we recalculated the tidal fluxes using the lowest and highest eccentricity values given in Table 1. Hence, we got nine different values for the tidal flux (one mean value for each density and the upper and lower limits provided by the minimum and maximum eccentricities).
When calculating the incident stellar flux, the largest uncertainty comes from the unknown albedo of the planets. For these calculations, we used the albedo value of the Earth, . Because the true albedos of the planets are not known, we repeated the calculations for 0.1 and 0.5 planetary albedos, too, and used the difference as error bars for the stellar (and consequently for the global) flux of the planets.
For the runaway greenhouse flux, the uncertainty came from the change of the surface gravity which depends on the mass and radius of the planet which are highly uncertain. We used the maximum and minimum values of stellar flux obtained from the minimum and maximum density calculations, and indicate them as uncertainties of the global flux of the planets. The uncertainties of the tidal flux were neglected here, because those are two orders of magnitude lower compared to the uncertainties in the stellar flux.
4 Results and discussion
Table 2 summarizes the volume fraction of each constituent material in each of the TRAPPIST-1 planets. Because of the large uncertainties in the measured mass and radius of the planets, we first calculated three end-member structures, corresponding to the highest-density planet, a structure for the mean-density planet, in which the planet has the best-fit mass and radius from Wang et al. (2017), and the low-density planet corresponding to the lowest permitted density. For planet h, the mean density value was very close to the minimum that could be obtained with assuming a low-density rock core, thus, we reported only the mean and maximum density cases.
Given the large error bars on mass and radius for each planet, it was not always possible to determine the bulk composition or the number of different compositional layers inside the planet – for planets b, c, d, and e, water-free structures were permitted if the planets had densities close to their upper limits. Likewise, each of these planets could be composed of rock and HO, completely devoid of iron, if their true densities were close to their lower limits. For planets b, c, d, and e, we also computed the limit where the structure qualitatively changed, for example, one of the layers disappeared or a new one appeared, consistent with each planet’s possible values of mass and radius.
The values in Table 2 served as inputs in our calculations of heat generation by tides and loss by interior processes. Table 3 summarizes the values of equilibrium mantle temperature, tidal heat flux, and global energy flux at the top of the atmosphere from our geophysical models. In addition to the uncertainties provided by the densities (mass and radius values), the orbital eccentricity altered the tidal heating flux on the planet, too. Calculating with the three density and three eccentricity values for each planet, we obtained nine different results for the equilibrium temperatures and the tidal heating fluxes. (For planet h it resulted only in six different values.) For simplicity, we only presented the one case calculated with the mean density and mean eccentricity value, and the upper and lower uncertainties from among all other cases. The equilibrium values of mantle temperature obtained in our calculations, as well as the tidal heat flux, , and threshold heat flux for a runaway greenhouse, , are summarized in Table 3.
4.1 Planet b
The top row of Fig. 3 illustrates our proposed interior structures for TRAPPIST-1b. The mass of planet b is comparable to the planet Venus, but with a radius of , its mean density of kg/m is significantly lower than the mean densities of either Venus or Earth. Given our assumed densities for rock and iron, the highest permitted density for TRAPPIST-1b, kg/m, permitted the planet to be composed of these two materials, without a significant reservoir of HO (see left panel of Fig. 3). However, this represents an extreme end-member scenario. The absolute lower limit of mean density, kg/m, permitted an interior structure with no iron, solely composed of rock and HO (see right panel of Fig. 3). Given the current estimates of mass and density, the true structure of TRAPPIST-1b lies somewhere between these end-members (for example, Fig. 3, middle), with an iron core, rock mantle, mantle of high pressure polymorphs, and a surface layer of HO. This is consistent with more detailed interior models that take into account compression of rock and metal at depth, which suggest that TRAPPIST-1b must contain more than 6 to 8% HO (Unterborn et al., 2017), likely about HO by mass (Wang et al., 2017). It seems likely that the planet does have some liquid water, even if it is just a relatively thin layer at the surface (Unterborn et al., 2017). Due to its proximity to the star, the effective surface temperature on TRAPPIST-1b, K. Thus, the planet must be in a runaway greenhouse state.
The top, left panel of Fig. 4 illustrates how the interior structure and composition of TRAPPIST-1b varies as a function of planetary mass and radius. The figure depicts the range of masses and radii for which TRAPPIST-1b is composed of ice, rock, and metal in a four-layered interior structure similar to the middle panel in the top row of Fig. 3. At smaller masses and larger radii, TRAPPIST-1b is composed of ice and rock; for larger masses and smaller radii, the planet is made of rock and iron. If the error bars on the mass for TRAPPIST-1b were reduced to , it would be possible to determine which of the three compositions best matched the planet’s true composition.
We found that the equilibrium mantle temperature for TRAPPIST-1b was above the rock solidus temperature, suggesting that the rock component of TRAPPIST-1b can be partially molten. If TRAPPIST-1b has a thick mantle of solid hpp ice, buoyant melt will rise to the top of the mantle, where it can erupt at the rock/ice boundary. This will drive melting at the base of the high-pressure ice polymorph mantle, driving ice melting and the rise of buoyant water through the hpp ice mantle (Barr et al., 2001).
The surface heat flux from tidal heat for TRAPPIST-1b, W/m, is twice that of Jupiter’s moon Io (Spencer et al., 2000; Veeder et al., 2004), which is one of the most volcanically active planetary bodies in our Solar System (Schubert et al., 2004). Thus, we expect TRAPPIST-1b to be volcanically active, but any activity would be concentrated at the ice/rock boundary, or at the bottom of a thick surface ocean. This raises the possibility of habitable regions inside TRAPPIST-1b in the vicinity of, for example, hydrothermal vents at the ice/rock boundary, as has been suggested for Jupiter’s moon Europa (Vance et al., 2007). However, the global energy flux on TRAPPIST-1b is about three times larger than the runaway greenhouse limit for the planet, indicating that the surface temperature on the planet may be too high to support life.
4.2 Planet c
TRAPPIST-1c has the highest mean density among the planets in the system, with an estimated kg/m. The simplest permitted structure with such a high density is a large iron core, with a rock mantle, and no HO layer at the surface (see middle panel in the second row of Fig. 3). This structure is consistent with other interior models, which show that the planet is roughly 50% iron by mass (Wang et al., 2017), and less than 6 to 8 % HO by mass (Unterborn et al., 2017). Figure 3 also illustrates the two extreme end-member interior structures that were possible given the large error bars on the mean density for TRAPPIST-1c. Taking the highest possible density for TRAPPIST-1c, the simplest interior structure was essentially a massive iron core with a thin layer of rock at the surface. The other extreme possibility was a rock-rich planet with no metal core and a mantle of HO (see right panel of Fig. 3).
The top, right panel of Fig. 4 illustrates the pairs of mass and radius for TRAPPIST-1c where a rock/iron composition is permitted, and the small region of parameter space where TRAPPIST-1c could be composed of rock and HO. Further improvements to the mass measurement of TRAPPIST-1c could help distinguish between these models; if we knew the mass to within , we could determine whether TRAPPIST-1c might harbor a layer of HO near its surface. Consistent with others’ interior models, we found that HO-bearing structures for TRAPPIST-1c are possible only if its density is close to the lower limit permitted by the error bars.
The equilibrium mantle temperature in the rock layer inside TRAPPIST-1c was K, above the solidus temperature for rock, and a surface tidal heat flux of 1.3 W/m, comparable to Jupiter’s moon Io. Because it is unlikely that TRAPPIST-1c has a surface layer of HO, it is possible that TRAPPIST-1c may have surface eruptions of silicate magma. The high heat flow may be large enough to be observed in infrared wavelengths. The heat flow on Io is known to be spatially and temporally variable as new volcanic vents form, or eruptions occur from previously dormant calderas (e.g., Veeder et al. 2012). The surface of Io also harbors numerous lava lakes whose surfaces periodically rupture and founder, resulting in momentary increases in the global heat flow (Veeder et al., 2012). Thus, as instrumentation improves, searching for excess thermal emission from the TRAPPIST-1 system may provide an avenue for the discovery of volcanism in the system.
In the mean density case, the iron core had a radius of km which is larger than the iron core of the Earth (radius km). This implies that planet c might have a strong magnetic field protecting the surface from stellar wind and flare erosions. However, due to the close-in orbit, the global energy flux is almost twice as large as the runaway greenhouse limit, which probably results in an uninhabitable surface.
4.3 Planet d
Planet d is much smaller than its inner two neighbors, with a mass , and radius . The third row of Fig. 3 illustrates the range of possible interior structures permitted by the mass and radius estimates. In the high-density limit, the planet could be composed solely of rock and iron, with an iron core of radius 2820 km. In the low-density limit, the planet could have no iron and be composed of rock and HO (right panel of Fig. 3). If planet d has any significant water content, it would be concentrated near the surface of the planet and would form a liquid water ocean at the surface because K.
Figure 4 illustrates the range of possible compositions for TRAPPIST-1d. For the vast majority of masses and radii, a significant HO component was required to match the observations. If the mass of the planet is close to its maximum value, and radius close to the minimum permitted value, structures with no HO can match the observations. Determining the mass of TRAPPIST-1d to would shed light on whether TRAPPIST-1d has significant water content.
The equilibrium interior temperature in the rocky component of TRAPPIST-1d K was extremely close to the solidus temperature for rock, indicating that the mantle could be partially molten, but certainly low-viscosity and convective. Because TRAPPIST-1d likely has a thick layer of HO on the surface, any melt extraction from the mantle would result in melting of the solid ice mantle. With an effective surface temperature K, TRAPPIST-1d is likely to be covered by a global water ocean. With a global heat flux from tidal heating of 160 mW/m, TRAPPIST-1d experiences a tidal heat flow twenty times the mean heat flow of the Earth, which provides an abundance of geothermal energy to drive chemistry in its surface ocean.
We only found an equilibrium temperature for the high-density case for high orbital eccentricities. In the cases of lower eccentricity, the tidal heating rate was not high enough to drive convection in the body, hence no equilibrium temperature can be found. In other words, this indicated that the planet may not experience significant tidal heating, and/or that a more sophisticated tidal model may be needed for these cases to determine the balance between tidal heating and convection. The error bars for the eccentricity of planet d permitted the possibility that the planet’s orbit is almost circular, implying insignificant tidal heating.
According to the calculations of Kopparapu et al. (2017), planet d is in the runaway greenhouse state. In contrast, we found that the global energy flux is lower than the runaway greenhouse limit, if the planetary albedo is . However, Kopparapu et al. (2017) uses a climate model for Earth-like bodies, but our model Pierrehumbert (2010) took into account the mass of the planet (through its surface gravity), which was significantly smaller than the Earth (). Because the planet was very close to the runaway greenhouse limit (if not exceeded it), liquid water on its surface was only possible if the planet was tidally locked to the star. Because of the synchronization, a high albedo cloud above the substellar point can reflect most incoming radiation (Yang et al., 2013; Kopparapu et al., 2016).
4.4 Planets e, f, and g
As a group, TRAPPIST-1e, f, and g have somewhat similar masses, and radii, yielding similar mean densities. Each of the three planets is less massive than Earth but with a radius comparable to that of the Earth, yielding mean densities in between those for ice and rock, similar to the mean densities of outer planet satellites and dwarf planets in our Solar System, including Pluto (Stern et al., 2015), Enceladus (Porco et al., 2006), and Ganymede (Pappalardo et al., 2004).
The mass of planet e is quite uncertain, : the error bars encompass . As a result, the mean density of planet e is also not well-constrained: kg/m (Wang et al., 2017). In the context of our assumed densities for ice, rock, and iron, the minimum density permitted for planet e was kg/m, the density of the high-pressure ice polymorphs. With this density, the minimum permitted mass for planet e is , with a radius of . Assuming a low-density core composed of Prinn-Fegley rock with kg/m, TRAPPIST-1e could harbor a small rock core, even though its minimum density was comparable to (see Table 1). The fourth row of Fig. 3 illustrates the end-member interior structures, ranging from a water-free planet to one composed of almost 100% HO (right). With an effective surface temperature K (Gillon et al., 2017) below the melting point of water ice at low pressures (273 K), it is possible that TRAPPIST-1e has a solid ice surface. The relative sizes of the liquid water ocean and polymorph layers depended on the temperature in the interior of the planet, but the maximum thickness of the ice I layer can be constrained by the depth to the 209 MPa phase transition among ice I, liquid water, and ice III (Hobbs, 1974).
The bottom, right panel of Fig. 4 illustrates the pairs of mass and radius for which TRAPPIST-1e had an HO-rich composition. If the mass of TRAPPIST-1e is close to the maximum permitted value, the planet could contain only rock and iron, but for the vast majority of mass/radius pairs, a thick ice mantle is predicted. To rule out a water-free structure, one would have to determine the mass of TRAPPIST-1e to within .
Planets f and g have masses and radii similar to planet e, and thus, similar mean densities and interior structures. Both planets are less dense than TRAPPIST-1e, and so it seems likely that both of the planets have a significant amount of HO. Similar to planet e, the mass of TRAPPIST-1f is not well-constrained, and the error bars permit the planet to have a low-density structure with a mean density equal to that of the high-pressure ice polymorphs. Because there is a thin layer of ice I or water at the surface, the low-density structure could have a very small Prinn-Fegley rock core.
For planets e and f, we found that their rock mantles would have equilibrium temperatures near the solidus of rock, similar to planet d. Also their global heat fluxes from tidal heating were to 180 mW/m similar to planet d. Both planets f and g are far enough from their parent star to have K (Gillon et al., 2017), suggesting that their surfaces could be solid HO ice. Given that the tidal heat flow is about twenty times stronger than Earth, it seems likely that both planets could harbor liquid water oceans, perched between the surface layer of ice I and underlain by a mantle of hpp ices, similar to the jovian moons Ganymede and Callisto (Schubert et al., 2004).
We were not able to determine an equilibrium temperature for planet g, nor for the low-density case of planet f. This could be because tidal forces from the star are too small to significantly deform the rock cores of these planets, implying that tidal dissipation in the rock mantles of each planet is negligible. Tidal dissipation could generate heat in the ice mantles of the planet, but a more sophisticated model for would be needed to determine magnitude of the heat and its partitioning between the hpp ice mantle and ice I shell.
4.5 Planet h
The outermost planet, TRAPPIST-1h, has a mass of only , and a radius of (Wang et al., 2017). Given our assumed densities of rock and ice, TRAPPIST-1h, with a mean density kg/m, should contain only HO, and no rock or metal.
Although it was permitted in the context of our assumptions, we thought it unlikely that TRAPPIST-1h would be composed of pure water ice. The vast majority of ice-rich bodies in the outer solar system contain some fraction of rock (one possible exception being Saturn’s small moon Tethys, which may not be primordial (Salmon and Canup, 2017)). In our solar system, almost all of the outer planet satellites and dwarf planets contain some rock, and from a formation standpoint, it seems unlikely that TRAPPIST-1h is completely devoid of rock. For this reason we determined interior structures using a low-density mineral assemblage commonly used as a proxy for the rocky component of outer Solar System satellites: Prinn-Fegley rock, for the planet, and a planetary mean density of 1300 kg/m.
We found that a model TRAPPIST-1h with a density could harbor a core of rock 1030 km in radius, with the remainder of the planet composed of HO in various phases (see right panel in the last row of Fig. 3). If TRAPPIST-1h had a density close to its maximum permitted value, kg/m, it could be composed of almost pure rock, with only km of HO near the surface (see left panel in the last row of Fig. 3). With its low effective surface temperature, K (Gillon et al., 2017), planet h undoubtedly has a solid ice surface. Similar to our null result for TRAPPIST-1g, we were not able to determine an equilibrium temperature dictated by the balance of tidal heat generation and transport in TRAPPIST-1h. Further knowledge of the composition and interior structure of TRAPPIST-1g and h, as well as a more sophisticated tidal model, will be required to assess the role of tidal heating in its interior evolution.
5 Implications for habitability
TRAPPIST-1b may contain a large amount of HO, and is likely volcanically active, but its activity is probably concentrated at the ice/rock boundary deep inside the planet. Due to its proximity to the star, the planet is very likely to be in a runaway greenhouse state (Kopparapu et al., 2017; Bolmont et al., 2017; Bourrier et al., 2017). TRAPPIST-1c is not likely to have much water or ice, and could have surface eruptions of silicate magma in the style of Io. Detection of excess thermal emission from the TRAPPIST-1 system could provide indirect evidence for surface volcanism on these bodies.
Using 1-dimensional and 3-dimensional climate models, Gillon et al. (2017) deduce that the surface temperatures of planets e, f, and g are suitable for harboring water oceans on their surface. However, using the runaway greenhouse limit of Pierrehumbert (2010), we found that planet d, too, may be habitable, if its albedo is . Planet d might be covered by a global water ocean that can provide a favorable environment for the appearance of life. This conclusion supports the finding of Vinson and Hansen (2017), who calculate the surface temperature of the TRAPPIST-1 planets due to stellar irradiation and tidal heating due to the circularization of the orbits, and found that planets d, e, and f might be habitable.
Planets e, f, g, and h are likely to have icy surfaces with liquid oceans underneath. According to Checlair et al. (2017), however, if the planets are tidally locked and have an active carbon cycle, then they should not have solid ice surfaces. The three-dimensional climate models of Turbet et al. (2017) show that TRAPPIST-1e is able to sustain its water content in liquid phase regardless of its atmospheric composition, provided that the planet is synchronously locked and abundant in water. In the atmospheres of planets f and g, a few bars of carbon dioxide is needed to facilitate the melting of ice on their surfaces. Most studies find that planet e is most likely to be habitable among the TRAPPIST-1 planets. For example, according to Wolf (2017), even a thin atmosphere could be enough to sustain temperatures of K locally around the substellar point, melting % of the surface ice. Vinson and Hansen (2017) also find that including tidal heating to the energy budget, the surface temperature of planet e is around 280 K.
According to the tidal heating calculations of Papaloizou et al. (2017), planet e, and in most of the investigated cases, also planet f, could be both in the conservative habitable zone and in the tidal habitable zone. The tidal habitable zone indicates an intermediate level of tidal heating flux on the planetary surface (between 0.04 and 2 W/m) which could drive plate tectonics required for the carbonate-silicate cycle that stabilizes the atmospheric temperature of the planet (Barnes et al., 2009). Papaloizou et al. (2017) also find that planet d orbits in the optimistic habitable zone, and its tidal heating is sufficient to be in the tidal habitable zone, too.
Bourrier et al. (2017) find that during the assumed 8 Gyr lifetime of the system, planets b through f may have lost, together, more than 20 Earth oceans worth of water due to hydrodynamic escape. However, the estimation is highly uncertain because of the wide possible ranges in planetary masses, and their values should be considered upper limits for hydrodynamical water loss. Magnetic interaction with the host star, however, can result in additional loss of water (Dong et al., 2017).
In any case, on planets b, e, f, g, and h, life might appear in the tidally heated (subsurface) ocean close to hydrothermal vents. Although our tidal heating model did not constrain equilibrium temperatures for planets g and h, they may still experience tidal heating. Our model might not be suitable for low heating rates: tidal heating can also be balanced by conduction alone if the planet’s mantle is too viscous to convect. We have also neglected planet–planet tidal interactions which may significantly contribute to the tidal heating in each planet, because they orbit in a close mean motion resonant chain. In addition, using tidal evolution models, Barnes (2017) find that planets in the habitable zone of late type M dwarfs are very likely to enter a 1:1 spin–orbit resonance in less than 1 Gyr. If the planet’s rotation is fast, then the orbital eccentricity is expected to grow slowly, which will lead to somewhat stronger tidal heating.
Lingam and Loeb (2017) investigate the probability of habitability of three TRAPPIST-1 planets, e, f, and g, orbiting inside the circumstellar habitable zone. By calculating the likelihood functions for habitability based on surface temperature and atmospheric loss of the planets (due to hydrodynamic escape and stellar winds), they find that for planet e, which gave the most promising results, the likelihood function is at least one to two orders of magnitude lower than in the case of Earth.
Any form of life is threatened however, by strong stellar activity. Analysis of K2 observations of TRAPPIST-1 showed frequent flares which present a serious threat for planetary atmospheres (Vida et al., 2017; Garraffo et al., 2017). Strong magnetic fields are required to protect the atmosphere from erosion, and life from high energy radiation.
The seven planets of the TRAPPIST-1 system have been the focus of intense observation and debate since the discovery of the system. Several of the planets are roughly Earth-sized, and their mean densities are not dissimilar from the terrestrial planets and large icy satellites in our Solar System. Moreover, the effective surface temperatures of each of the planets are relatively modest, suggesting solid or water-covered surfaces.
Each of the seven TRAPPIST-1 planets are quite close to the parent star, as indicated by their orbital periods which range from 1.5 to about 20 days. The planets also occupy a mean-motion resonance (Gillon et al., 2017) and have non-zero orbital eccentricities. This implies that the planets experience tidal heating, which can be a significant energy source.
Here, we have constructed interior models for the planets using constant densities for ice, rock, and metal. The masses and radii of each of the planets are not well-constrained and have large uncertainties; for this reason, we applied our model using the mean density of each planet and also the extrema obtained from the error bars of observations. Within the range of estimated masses and radii for planets b, c, d, and e, each planet could be composed of solely iron and rock (our “high-density” cases). If the true densities of planets b, d, and e are close to the mean values implied by observations, each of these planets could have a layer of HO near the surface. The large density of planet c points toward a rock/iron composition, HO is permitted only if the true density is closer to the lower limit of the estimated value.
Planets b and c likely harbor magma oceans in their interiors. The surface heat fluxes from tidal heat range from hundreds of milliwatts per meter squared (about twenty times higher than the Earth’s heat flow) to watts per meter squared, comparable to the heat flow estimated for Jupiter’s volcanically active moon, Io. Further work is required to reconcile the seemingly water-rich composition of planet b with its proximity to the TRAPPIST-1 star. The majority of permitted compositions for planet c indicate that the body is water-poor, raising the possibility that melt extracted from the magma ocean could erupt onto the surface. Further observations of the TRAPPIST-1 system may be able to detect the thermal or chemical signature of this activity, if present.
Of the seven planets in the TRAPPIST-1 system, we consider planets d and e to be the most likely candidates for habitable environments. Planet d is likely to be covered by a global water ocean, and can be habitable if its albedo is 0.3. Planetary albedos may be estimated from photometric measurements during occultations (see e.g. Rowe et al., 2008) which may give an additional constraint on the habitability of TRAPPIST-1d. As discussed in Section 5, planet e is likely to have liquid water on its surface, too, according to climate models (Turbet et al., 2017; Wolf, 2017).
Planets f, g, and h are undoubtedly HO-rich, with thick mantles of ice, possible subsurface liquid water oceans, and solid surfaces composed of ice I, the low-density phase of ice present in terrestrial glaciers and ice sheets. It is worth noting that the present estimates of the mean density of planet h permit a 100% HO composition, but this seems unlikely given that even the most ice-rich bodies in our Solar System (for example, the Kuiper Belt Objects) have a significant amount of rock in their interiors (Barr and Schwamb, 2016). The geodynamical evolution of the outer planets of TRAPPIST-1 should be studied in more detail with a layered tidal and thermal model that can self-consistently determine the amount of tidal heating and heat transport by convection in both the rock and ice mantles.
Given the uncertainties in the compositions of each of the planets, we have employed a relatively simple geophysical model to determine the balance between tidal heat production and transport in each planet. Our approach provides more realistic values for than a constant- approach because we take into account the feedback between the temperature of the planet. We used a tidal model that assumes each planet is a homogeneous sphere, but mimic the effect of the presence of multiple materials by determining an effective viscosity and rigidity for the planet based on the viscosities and rigidities of its component materials. We focused our efforts on determining the temperature of the rock mantle in each planet by balancing heat produced with the amount of heat transported by solid-state convection in the rock layer of each planet. When updated masses and radii are available, we will be pursuing more detailed calculations of tidal heating using multi-layered models and numerical techniques.
Further observations of the TRAPPIST-1 system are required to provide better constraints on the composition and interior structure of each planet. When improved compositional constraints are obtained, more sophisticated modeling techniques that calculate tidal dissipation and energy transport in a multi-layered planet can provide improved estimates of the tidal heat flux and mantle temperature. Our knowledge of the compositions of the four interior planets, chiefly, their water content, can significantly be improved if their masses can be determined to within to , which will require further observations of the TRAPPIST-1 system.
We thank the anonymous referee for a thorough review of the manuscript and helpful comments. We also thank Dr. László Molnár for initiating this collaboration, and Dr. Jessica Irving, Dr. Tess Caswell, and Dr. Ákos Kereszturi for helpful conversation. VD has been supported by the Hungarian National Research, Development and Innovation Office (NKFIH) grants K119993, K-115709 and GINOP-2.3.2-15-2016-00003.
- Barnes (2017) Barnes, R. (2017). Tidal Locking of Habitable Exoplanets. Celestial Mechanics and Dynamical Astronomy, 129:509–536.
- Barnes et al. (2009) Barnes, R., Jackson, B., Greenberg, R., and Raymond, S. N. (2009). Tidal Limits to Planetary Habitability. ApJ, 700:L30–L33.
- Barr (2008) Barr, A. C. (2008). Mobile lid convection beneath Enceladus’ south polar terrain. Journal of Geophysical Research (Planets), 113:7009.
- Barr and Canup (2008) Barr, A. C. and Canup, R. M. (2008). Constraints on gas giant satellite formation from the interior states of partially differentiated satellites. Icarus, 198(1):163–177.
- Barr and McKinnon (2007) Barr, A. C. and McKinnon, W. B. (2007). Convection in ice i shells and mantles with self-consistent grain size. Journal of Geophysical Research: Planets, 112(E2).
- Barr et al. (2001) Barr, A. C., Pappalardo, R. T., and Stevenson, D. J. (2001). Rise of Deep Melt into Ganymede’s Ocean and Implications for Astrobiology. In Lunar and Planetary Science Conference, volume 32 of Lunar and Planetary Science Conference.
- Barr and Schwamb (2016) Barr, A. C. and Schwamb, M. E. (2016). Interpreting the densities of the kuiper belt’s dwarf planets. Monthly Notices of the Royal Astronomical Society, 460(2):1542–1548.
- Barr and Showman (2009) Barr, A. C. and Showman, A. P. (2009). Heat Transfer in Europa’s Icy Shell, page 405.
- Beuthe (2013) Beuthe, M. (2013). Spatial patterns of tidal heating. Icarus, 223:308–329.
- Bolmont et al. (2017) Bolmont, E., Selsis, F., Owen, J. E., Ribas, I., Raymond, S. N., Leconte, J., and Gillon, M. (2017). Water loss from terrestrial planets orbiting ultracool dwarfs: implications for the planets of TRAPPIST-1. MNRAS, 464:3728–3741.
- Bourrier et al. (2017) Bourrier, V., de Wit, J., Bolmont, E., Stamenković, V., Wheatley, P. J., Burgasser, A. J., Delrez, L., Demory, B.-O., Ehrenreich, D., Gillon, M., Jehin, E., Leconte, J., Lederer, S. M., Lewis, N., Triaud, A. H. M. J., and Van Grootel, V. (2017). Temporal Evolution of the High-energy Irradiation and Water Content of TRAPPIST-1 Exoplanets. AJ, 154:121.
- Canup and Ward (2002) Canup, R. M. and Ward, W. R. (2002). Formation of the Galilean Satellites: Conditions of Accretion. Astron. J, 124:3404–3423.
- Checlair et al. (2017) Checlair, J., Menou, K., and Abbot, D. S. (2017). No snowball on habitable tidally locked planets. ApJ, 845:132.
- Dobos and Turner (2015) Dobos, V. and Turner, E. L. (2015). Viscoelastic Models of Tidally Heated Exomoons. ApJ, 804:41.
- Dong et al. (2017) Dong, C., Jin, M., Lingam, M., Airapetian, V. S., Ma, Y., and van der Holst, B. (2017). Atmospheric escape from the TRAPPIST-1 planets and implications for habitability. ArXiv e-prints, 1705.05535v2.
- Durand et al. (2006) Durand, G., Weiss, J., Lipenkov, V., Barnola, J., Krinner, G., Parrenin, F., Delmonte, B., Ritz, C., Duval, P., Röthlisberger, R., et al. (2006). Effect of impurities on grain growth in cold ice sheets. Journal of Geophysical Research: Earth Surface, 111(F1).
- Durham et al. (1997a) Durham, W. B., Kirby, S. H., and Stern, L. A. (1997a). Correction to “Creep of water ices at planetary conditions: A compilation,” by William B. Durham, Stephen H. Kirby, and Laura A. Stern. Journal of Geophysical Research, 102(E12):28725–28725.
- Durham et al. (1997b) Durham, W. B., Kirby, S. H., and Stern, L. A. (1997b). Creep of water ices at planetary conditions: A compilation. Journal of Geophysical Research, 102(E7):16293–16302.
- Dziewonski and Anderson (1981) Dziewonski, A. M. and Anderson, D. L. (1981). Preliminary reference Earth model. Physics of the Earth and Planetary Interiors, 25:297–356.
- Efroimsky (2012) Efroimsky, M. (2012). Tidal dissipation compared to seismic dissipation: In small bodies, Earths, and super-Earths. ApJ, 746:150.
- Fischer and Spohn (1990) Fischer, H.-J. and Spohn, T. (1990). Thermal-orbital histories of viscoelastic models of Io (J1). Icarus, 83:39–65.
- Gammon et al. (1983) Gammon, P. H., Kiefte, H., Clouter, M. J., and Denner, W. W. (1983). Elastic constants of artificial and natural ice samples by brillouin spectroscopy. Journal of Glaciology, 29(103):433–460.
- Garraffo et al. (2017) Garraffo, C., Drake, J. J., Cohen, O., Alvarado-Gomez, J. D., and Moschou, S. P. (2017). The Threatening Environment of the TRAPPIST-1 Planets. ApJ, 843(L33).
- Gillon et al. (2017) Gillon, M., Triaud, A. H. M. J., Demory, B.-O., Jehin, E., Agol, E., Deck, K. M., Lederer, S. M., de Wit, J., Burdanov, A., Ingalls, J. G., Bolmont, E., Leconte, J., Raymond, S. N., Selsis, F., Turbet, M., Barkaoui, K., Burgasser, A., Burleigh, M. R., Carey, S. J., Chaushev, A., Copperwheat, C. M., Delrez, L., Fernandes, C. S., Holdsworth, D. L., Kotze, E. J., Van Grootel, V., Almleaky, Y., Benkhaldoun, Z., Magain, P., and Queloz, D. (2017). Seven temperate terrestrial planets around the nearby ultracool dwarf star TRAPPIST-1. Nature, 542:456–460.
- Goldreich and Soter (1966) Goldreich, P. and Soter, S. (1966). Q in the Solar System. Icarus, 5:375–389.
- Goldsby and Kohlstedt (2001) Goldsby, D. and Kohlstedt, D. (2001). Superplastic deformation of ice: Experimental observations. Journal of Geophysical Research: Solid Earth, 106(B6):11017–11030.
- Greeley et al. (2004) Greeley, R., Chyba, C. F., Head, III, J. W., McCord, T. B., McKinnon, W. B., Pappalardo, R. T., and Figueredo, P. H. (2004). Geology of Europa. In Bagenal, F., Dowling, T. E., and McKinnon, W. B., editors, Jupiter. The Planet, Satellites and Magnetosphere, pages 329–362.
- Harrison (1963) Harrison, J. C. (1963). An Analysis of the Lunar Tides. J. Geophys. Res., 68:4269–4280.
- Henning and Hurford (2014) Henning, W. G. and Hurford, T. (2014). Tidal Heating in Multilayered Terrestrial Exoplanets. ApJ, 789:30.
- Henning et al. (2009) Henning, W. G., O’Connell, R. J., and Sasselov, D. D. (2009). Tidally Heated Terrestrial Exoplanets: Viscoelastic Response Models. ApJ, 707:1000–1015.
- Hobbs (1974) Hobbs, P. V. (1974). Ice physics.
- Hussmann et al. (2002) Hussmann, H., Spohn, T., and Wieczerkowski, K. (2002). Thermal Equilibrium States of Europa’s Ice Shell: Implications for Internal Ocean Thickness and Surface Heat Flow. Icarus, 156:143–151.
- Johansen et al. (2015) Johansen, A., Mac Low, M.-M., Lacerda, P., and Bizzarro, M. (2015). Growth of asteroids, planetary embryos, and Kuiper belt objects by chondrule accretion. Science Advances, 1:1500109.
- Karato et al. (1995) Karato, S.-I., Zhang, S., and Wenk, H.-R. (1995). Superplasticity in Earth’s Lower Mantle: Evidence from Seismic Anisotropy and Rock Physics. Science, 270:458–461.
- Kasting et al. (1993) Kasting, J. F., Whitmire, D. P., and Reynolds, R. T. (1993). Habitable Zones around Main Sequence Stars. Icarus, 101:108–128.
- Khurana et al. (2011) Khurana, K. K., Jia, X., Kivelson, M. G., Nimmo, F., Schubert, G., and Russell, C. T. (2011). Evidence of a global magma ocean in ioâs interior. Science, 332(6034):1186–1189.
- Koot and Dumberry (2011) Koot, L. and Dumberry, M. (2011). Viscosity of the Earth’s inner core: Constraints from nutation observations. Earth and Planetary Science Letters, 308:343–349.
- Kopparapu et al. (2017) Kopparapu, R. k., Wolf, E. T., Arney, G., Batalha, N. E., Haqq-Misra, J., Grimm, S. L., and Heng, K. (2017). Habitable Moist Atmospheres on Terrestrial Planets near the Inner Edge of the Habitable Zone around M Dwarfs. ApJ, 845:5.
- Kopparapu et al. (2016) Kopparapu, R. k., Wolf, E. T., Haqq-Misra, J., Yang, J., Kasting, J. F., Meadows, V., Terrien, R., and Mahadevan, S. (2016). The Inner Edge of the Habitable Zone for Synchronously Rotating Planets around Low-mass Stars Using General Circulation Models. ApJ, 819:84.
- Lingam and Loeb (2017) Lingam, M. and Loeb, A. (2017). Physical constraints on the likelihood of life on exoplanets. ArXiv e-prints, 1707.02996.
- Love (1906) Love, A. E. H. (1906). A treatise on the mathematical theory of elasticity. Cambridge University Press.
- Luger et al. (2017) Luger, R., Sestovic, M., Kruse, E., Grimm, S. L., Demory, B.-O., Agol, E., Bolmont, E., Fabrycky, D., Fernandes, C. S., Van Grootel, V., Burgasser, A., Gillon, M., Ingalls, J. G., Jehin, E., Raymond, S. N., Selsis, F., Triaud, A. H. M. J., Barclay, T., Barentsen, G., Howell, S. B., Delrez, L., de Wit, J., Foreman-Mackey, D., Holdsworth, D. L., Leconte, J., Lederer, S., Turbet, M., Almleaky, Y., Benkhaldoun, Z., Magain, P., Morris, B. M., Heng, K., and Queloz, D. (2017). A seven-planet resonant chain in TRAPPIST-1. Nature Astronomy, 1:0129.
- McCammon et al. (2016) McCammon, C., Caracas, R., Glazyrin, K., Potapkin, V., Kantor, A., Sinmyo, R., Prescher, C., Kupenko, I., Chumakov, A., and Dubrovinsky, L. (2016). Sound velocities of bridgmanite from density of states determined by nuclear inelastic scattering and first-principles calculations. Progress in Earth and Planetary Science, 3:10.
- McEwen et al. (2004) McEwen, A. S., Keszthelyi, L. P., Lopes, R., Schenk, P. M., and Spencer, J. R. (2004). The lithosphere and surface of Io, pages 307–328.
- Mignard (1980) Mignard, F. (1980). The evolution of the lunar orbit revisited. II. Moon and Planets, 23:185–201.
- Moore (2003) Moore, W. B. (2003). Tidal heating and convection in Io. Journal of Geophysical Research (Planets), 108:5096.
- Moore (2006) Moore, W. B. (2006). Thermal equilibrium in Europa’s ice shell. Icarus, 180:141–146.
- Morbidelli et al. (2012) Morbidelli, A., Lunine, J. I., O’Brien, D. P., Raymond, S. N., and Walsh, K. J. (2012). Building terrestrial planets. Ann. Rev. Earth. Planet. Sci., 40:251–275.
- Mueller and McKinnon (1988) Mueller, S. and McKinnon, W. B. (1988). Three-layered models of ganymede and callisto: Compositions, structures, and aspects of evolution. Icarus, 76(3):437–464.
- Murray and Dermott (1999) Murray, C. D. and Dermott, S. F. (1999). Solar System Dynamics. Cambridge University Press.
- Ojakangas and Stevenson (1989) Ojakangas, G. W. and Stevenson, D. J. (1989). Thermal state of an ice shell on Europa. Icarus, 81:220–241.
- Papaloizou et al. (2017) Papaloizou, J. C. B., Szuszkiewicz, E., and Terquem, C. (2017). The TRAPPIST-1 system: Orbital evolution, tidal dissipation, formation and habitability. MNRAS, page stx2980.
- Pappalardo et al. (2004) Pappalardo, R. T., Collins, G. C., Head, III, J. W., Helfenstein, P., McCord, T. B., Moore, J. M., Prockter, L. M., Schenk, P. M., and Spencer, J. R. (2004). Geology of Ganymede, pages 363–396.
- Peale (2003) Peale, S. J. (2003). Tidally Induced Volcanism. Celestial Mechanics and Dynamical Astronomy, 87:129–155.
- Peale and Cassen (1978) Peale, S. J. and Cassen, P. (1978). Contribution of tidal dissipation to lunar thermal history. Icarus, 36:245–269.
- Pierrehumbert (2010) Pierrehumbert, R. T. (2010). Principles of Planetary Climate.
- Porco et al. (2006) Porco, C. C., Helfenstein, P., Thomas, P. C., Ingersoll, A. P., Wisdom, J., West, R., Neukum, G., Denk, T., Wagner, R., Roatsch, T., Kieffer, S., Turtle, E., McEwen, A., Johnson, T. V., Rathbun, J., Veverka, J., Wilson, D., Perry, J., Spitale, J., Brahic, A., Burns, J. A., Del Genio, A. D., Dones, L., Murray, C. D., and Squyres, S. (2006). Cassini Observes the Active South Pole of Enceladus. Science, 311:1393–1401.
- Remus et al. (2012) Remus, F., Mathis, S., Zahn, J.-P., and Lainey, V. (2012). Anelastic tidal dissipation in multi-layer planets. A&A, 541:A165.
- Remus et al. (2015) Remus, F., Mathis, S., Zahn, J.-P., and Lainey, V. (2015). The surface signature of the tidal dissipation of the core in a two-layer planet. A&A, 573:A23.
- Renner et al. (2000) Renner, J., Evans, B., and Hirth, G. (2000). On the rheologically critical melt fraction. Earth and Planetary Science Letters, 181(4):585–594.
- Roberts and Nimmo (2008) Roberts, J. H. and Nimmo, F. (2008). Tidal heating and the long-term stability of a subsurface ocean on Enceladus. Icarus, 194:675–689.
- Rowe et al. (2008) Rowe, J. F., Matthews, J. M., Seager, S., Miller-Ricci, E., Sasselov, D., Kuschnig, R., Guenther, D. B., Moffat, A. F. J., Rucinski, S. M., Walker, G. A. H., and Weiss, W. W. (2008). The Very Low Albedo of an Extrasolar Planet: MOST Space-based Photometry of HD 209458. ApJ, 689:1345–1353.
- Sabadini et al. (1982) Sabadini, R., Boschi, E., and Yuen, D. A. (1982). Polar wandering and the forced responses of a rotating, multilayered, viscoelastic planet. J. Geophys. Res., 87:2885–2903.
- Salmon and Canup (2017) Salmon, J. and Canup, R. M. (2017). Accretion of saturn’s inner mid-sized moons from a massive primordial ice ring. 836:109.
- Schubert et al. (2004) Schubert, G., Anderson, J. D., Spohn, T., and McKinnon, W. B. (2004). Interior Composition, Structure and Dynamics of the Galilean Satellites, pages 281–306.
- Schubert et al. (2001) Schubert, G., Turcotte, D. L., and Olson, P. (2001). Mantle Convection in the Earth and Planets.
- Segatz et al. (1988) Segatz, M., Spohn, T., Ross, M. N., and Schubert, G. (1988). Tidal dissipation, surface heat flow, and figure of viscoelastic models of Io. Icarus, 75:187–206.
- Shimizu et al. (1996) Shimizu, H., Nabetani, T., Nishiba, T., and Sasaki, S. (1996). High-pressure Elastic properties of the VI and VII phase of Ice in dense HO and DO. Physical Review B, 53(10):6107–6110.
- Solomatov and Moresi (2000) Solomatov, V. and Moresi, L.-N. (2000). Scaling of time-dependent stagnant lid convection: Application to small-scale convection on earth and other terrestrial planets. Journal of Geophysical Research: Solid Earth, 105(B9):21795–21817.
- Spencer et al. (2009) Spencer, J. R., Barr, A. C., Esposito, L. W., Helfenstein, P., Ingersoll, A. P., Jaumann, R., McKay, C. P., Nimmo, F., and Waite, J. H. (2009). Enceladus: An Active Cryovolcanic Satellite, page 683.
- Spencer et al. (2000) Spencer, J. R., Rathbun, J. A., Travis, L. D., Tamppari, L. K., Barnard, L., Martin, T. Z., and McEwen, A. S. (2000). Io’s Thermal Emission from the Galileo Photopolarimeter- Radiometer. Science, 288:1198–1201.
- Stern et al. (2015) Stern, S. A., and co-authors (2015). The Pluto system: Initial results from its exploration by New Horizons. Science, 350:aad1815.
- Tackley (2001) Tackley, P. J. (2001). Convection in Io’s asthenosphere: Redistribution of nonuniform tidal heating by mean flows. J. Geophys. Res., 106:32971–32982.
- Turbet et al. (2017) Turbet, M., Bolmont, E., Leconte, J., Forget, F., Selsis, F., Tobie, G., Caldas, A., Naar, J., and Gillon, M. (2017). Climate diversity on cool planets around cool stars with a versatile 3-D Global Climate Model: the case of TRAPPIST-1 planets. ArXiv e-prints, 1707.06927.
- Unterborn et al. (2017) Unterborn, C. T., Desch, S. J., Hinkel, N., and Lorenzo, A. (2017). Constraining the Compositions of the TRAPPIST-1 Planets to Trace Snow Lines and Migration in M Dwarf Disks. ArXiv e-prints, 1706.02689.
- Vance et al. (2007) Vance, S., Harnmeijer, J., Kimura, J., Hussmann, H., DeMartin, B., and Brown, J. M. (2007). Hydrothermal systems in small ocean planets. Astrobiology, 7(6):987–1005.
- Veeder et al. (2012) Veeder, G. J., Davies, A. G., Matson, D. L., Johnson, T. V., Williams, D. A., and Radebaugh, J. (2012). Io: Volcanic thermal sources and global heat flow. Icarus, 219:701–722.
- Veeder et al. (2004) Veeder, G. J., Matson, D. L., Johnson, T. V., Davies, A. G., and Blaney, D. L. (2004). The polar contribution to the heat flow of Io. Icarus, 169:264–270.
- Vida et al. (2017) Vida, K., Kővári, Z., Pál, A., Oláh, K., and Kriskovics, L. (2017). Frequent Flaring in the TRAPPIST-1 System–Unsuited for Life? ApJ, 841:124.
- Vinson and Hansen (2017) Vinson, A. M. and Hansen, B. M. S. (2017). On the Spin States of Habitable Zone Exoplanets Around M Dwarfs: The Effect of a Near-Resonant Companion. MNRAS, 472(3):3217–3229.
- Wang et al. (2017) Wang, S., Wu, D.-H., Barclay, T., and Laughlin, G. P. (2017). Updated Masses for the TRAPPIST-1 Planets. ArXiv e-prints, 1704.04290.
- Wolf (2017) Wolf, E. T. (2017). Assessing the Habitability of the TRAPPIST-1 System Using a 3D Climate Model. ApJ, 839:L1.
- Yang et al. (2013) Yang, J., Cowan, N. B., and Abbot, D. S. (2013). Stabilizing Cloud Feedback Dramatically Expands the Habitable Zone of Tidally Locked Planets. ApJ, 771:L45.
- Zeng et al. (2016) Zeng, L., Sasselov, D. D., and Jacobsen, S. B. (2016). Mass-Radius Relation for Rocky Planets Based on PREM. ApJ, 819:127.