Heat flow is an important constraint on planetary formation and evolution. It has been suggested that Martian obliquity cycles might cause periodic collapses in atmospheric pressure, leading to corresponding decreases in regolith thermal conductivity (which is controlled by gas in the pore spaces). Geothermal heat would then build up in the subsurface, potentially affecting present-day heat flow - and thus the measurements made by a heat-flow probe such as the InSight HP$^{3}$ instrument. To gauge the order of magnitude of this effect, we model the diffusion of a putative heat pulse caused by thermal conductivity changes with a simple numerical scheme and compare it to the heat-flow perturbations caused by other effects. We find that an atmospheric collapse to 300 Pa in the last 40 kyr would lead to a present-day heat flow that is up to $2-8\%$ larger than the average geothermal background. Considering the InSight mission with expected $5-15\%$ error bars on the HP$^{3}$ measurement, this perturbation would only be significant in the best-case scenario of full instrument deployment, completed measurement campaign, and a well-modelled surface configuration. The prospects for detecting long-term climate perturbations via spacecraft heat-flow experiments remain challenging.

This paper reviews habitability conditions for a terrestrial planet from the point of view of geosciences. It addresses how interactions between the interior of a planet or a moon and its atmosphere and surface (including hydrosphere and biosphere) can affect habitability of the celestial body. It does not consider in detail the role of the central star but focusses more on surface conditions capable of sustaining life. We deal with fundamental issues of planetary habitability, i.e. the environmental conditions capable of sustaining life, and the above-mentioned interactions can affect the habitability of the celestial body. We address some hotly debated questions including: - How do core and mantle affect the evolution and habitability of planets? - What are the consequences of mantle overturn on the evolution of the interior and atmosphere? - What is the role of the global carbon and water cycles? - What influence do comet and asteroid impacts exert on the evolution of the planet? - How does life interact with the evolution of the Earth's geosphere and atmosphere? - How can knowledge of the solar system geophysics and habitability be applied to exoplanets? In addition, we address the identification of preserved life tracers in the context of the interaction of life with planetary evolution.

Aims: The long-term carbon cycle for planets with a surface entirely covered by oceans works differently from that of the present-day Earth because inefficient erosion leads to a strong dependence of the weathering rate on the rate of volcanism. In this paper, we investigate the long-term carbon cycle for these planets throughout their evolution. Methods: We built box models of the long-term carbon cycle based on CO$_2$ degassing, seafloor-weathering, metamorphic decarbonation, and ingassing and coupled them with thermal evolution models of stagnant lid and plate tectonics planets. Results: The assumed relationship between the seafloor-weathering rate and the atmospheric CO$_2$ or the surface temperature strongly influences the climate evolution for both tectonic regimes. For a plate tectonics planet, the atmospheric CO$_2$ partial pressure is characterised by an equilibrium between ingassing and degassing and depends on the temperature gradient in subduction zones affecting the stability of carbonates. For a stagnant lid planet, partial melting and degassing are always accompanied by decarbonation, such that the combined carbon content of the crust and atmosphere increases with time. Whereas the initial mantle temperature for plate tectonics planets only affects the early evolution, it influences the evolution of the surface temperature of stagnant lid planets for much longer. Conclusions: For both tectonic regimes, mantle cooling results in a decreasing atmospheric CO$_2$ partial pressure. For a plate tectonics planet this is caused by an increasing fraction of subduction zones that avoid crustal decarbonation, and for stagnant lid planets this is caused by an increasing decarbonation depth. This mechanism may partly compensate for the increase of the surface temperature due to increasing solar luminosity with time and hereby contribute to keep planets habitable in the long-term.

We are witness to a great and increasing interest in internal structure, composition and evolution of exoplanets. However, direct measurements of exoplanetary mass and radius cannot be uniquely interpreted in terms of interior structure, justifying the need for an additional observable. The second degree fluid Love number, $k_2$, is proportional to the concentration of mass towards the body's center, hence providing valuable additional information about the internal structure. When hydrostatic equilibrium is assumed for the planetary interior, $k_2$ is a direct function of the planetary shape. Previous attempts were made to link the observed tidally and rotationally induced planetary oblateness in photometric light curves to $k_2$ using ellipsoidal shape models. Here, we construct an analytical 3D shape model function of the true planetary mean radius, that properly accounts for tidal and rotational deformations. Measuring the true planetary mean radius is critical when one wishes to compare the measured $k_2$ to interior theoretical expectations. We illustrate the feasibility of our method and show, by applying a Differential Evolution Markov Chain to synthetic data of WASP-121b, that a precision $\leq$ 65 ppm/$\sqrt{2\,min}$ is required to reliably retrieve $k_2$ with present understanding of stellar limb darkening, therefore improving recent results based on ellipsoidal shape models. Any improvement on stellar limb darkening would increase such performance.

The physical properties of cometary nuclei observed today relate to their complex history and help to constrain their formation and evolution. In this article, we review some of the main physical properties of cometary nuclei and focus in particular on the thermal, mechanical, structural and dielectric properties, emphasizing the progress made during the Rosetta mission. Comets have a low density of 480 $\pm$ 220 kg m-3 and a low permittivity of 1.9 - 2.0, consistent with a high porosity of 70 - 80 %, are weak with a very low global tensile strength $<$100 Pa, and have a low bulk thermal inertia of 0 - 60 J K-1 m-2 s-1/2 that allowed them to preserve highly volatiles species (e.g. CO, CO2, CH4, N2) into their interior since their formation. As revealed by 67P/Churyumov-Gerasimenko, the above physical properties vary across the nucleus, spatially at its surface but also with depth. The broad picture is that the bulk of the nucleus consists of a weakly bonded, rather homogeneous material that preserved primordial properties under a thin shell of processed material, and possibly covered by a granular material; this cover might in places reach a thickness of several meters. The properties of the top layer (the first meter) are not representative of that of the bulk nucleus. More globally, strong nucleus heterogeneities at a scale of a few meters are ruled out on 67P small lobe.

The habitability of a planet depends on various factors, such as delivery of water during the formation, the co-evolution of the interior and the atmosphere, as well as the stellar irradiation which changes in time. Since an unknown number of rocky exoplanets may operate in a one-plate convective regime, i.e., without plate tectonics, we aim at understanding under which conditions planets in such a stagnant-lid regime may support habitable surface conditions. Understanding the interaction of the planetary interior and outgassing of volatiles with the atmosphere in combination with the evolution of the host star is crucial to determine the potential habitability. M-dwarf stars in particular possess a high-luminosity pre-main sequence phase which endangers the habitability of planets around them via water loss. We therefore explore the potential of secondary outgassing from the planetary interior to rebuild a water reservoir allowing for habitability at a later stage. We compute the boundaries of the habitable zone around M, K, G, and F-dwarf stars using a 1D cloud-free radiative-convective climate model accounting for the outgassing history of CO2 and H2O from an interior evolution and outgassing model for different interior compositions and stellar luminosity evolutions. The outer edge of the habitable zone strongly depends on the amount of CO2 outgassed from the interior, while the inner edge is mainly determined via the stellar irradiation, as soon as a sufficiently large water reservoir has been outgassed. A build-up of a secondary water reservoir for planets around M-dwarf stars is possible even after severe water loss during the high luminosity pre-main sequence phase as long as some water has been retained within the mantle. Earth-like stagnant-lid planets allow for habitable surface conditions within a continuous habitable zone that is dependent on interior composition.

Is the present-day water-land ratio a necessary outcome of the evolution of plate tectonic planets with a similar age, volume, mass, and total water inventory as the Earth? This would be the case - largely independent of initial conditions - if Earth's present-day continental volume were at a stable unique equilibrium with strong self-regulating mechanisms of continental growth steering the evolution to this state. In this paper, we question this conjecture. Instead we suggest that positive feedbacks in the plate tectonics model of continental production and erosion may dominate and show that such a model can explain the history of continental growth. We investigate the main mechanisms that contribute to the growth of the volume of the continental crust. In particular, we analyze the effect of the oceanic plate speed, depending on the area and thickness of thermally insulating continents, on production and erosion mechanisms. Effects that cause larger continental production rates for larger values of continental volume are positive feedbacks. In contrast, negative feedbacks act to stabilize the continental volume. They are provided by the increase of the rate of surface erosion, subduction erosion, and crustal delamination with the continental volume. We systematically analyze the strengths of positive and negative feedback contributions to the growth of the continental crust. Although the strengths of some feedbacks depend on poorly known parameters, we conclude that a net predominance of positive feedbacks is plausible. We explore the effect of the combined feedback strength on the feasibility of modeling the observed small positive net continental growth rate over the past 2-3 billion years.

A long-term sublimation model to explain how Phaethon could provide the Geminid stream is proposed. We find that it would take $\sim6$ Myr or more for Phaethon to lose all of its internal ice (if ever there was) in its present orbit. Thus, if the asteroid moved from the region of a 5:2 or 8:3 mean motion resonance with Jupiter to its present orbit less than $1$ Myr ago, it may have retained much of its primordial ice. The dust mantle on the sublimating body should have a thickness of at least $15$ m but the mantle could have been less than $1$ m thick $1000$ years ago. We find that the total gas production rate could have been as large as $10^{27}\rm~s^{-1}$ then, and the gas flow could have been capable of lifting dust particles of up to a few centimeters in size. Therefore, gas production during the past millennium could have been sufficient to blow away enough dust particles to explain the entire Geminid stream. For present-day Phaethon, the gas production is comparatively weak. But strong transient gas release with a rate of $\sim4.5\times10^{19}\rm~m^{-2}s^{-1}$ is expected for its south polar region when Phaethon moves from $0^\circ$ to $2^\circ$ mean anomaly near perihelion. Consequently, dust particles with radii of $<\sim260~\mu m$ can be blown away to form a dust tail. In addition, we find that the large surface temperature variation of $>600$ K near perihelion can generate sufficiently large thermal stress to cause fracture of rocks or boulders and provide an efficient mechanism to produce dust particles on the surface. The time scale for this process should be several times longer than the seasonal thermal cycle, thereby dominating the cycle of appearance of the dust tail.

Context. The mass and radius of a planet directly provide its bulk density, which can be interpreted in terms of its overall composition. Any measure of the radial mass distribution provides a first step in constraining the interior structure. The fluid Love number $k_2$ provides such a measure, and estimates of $k_2$ for extrasolar planets are expected to be available in the coming years thanks to improved observational facilities and the ever-extending temporal baseline of extrasolar planet observations. Aims. We derive a method for calculating the Love numbers $k_n$ of any object given its density profile, which is routinely calculated from interior structure codes. Methods. We used the matrix-propagator technique, a method frequently used in the geophysical community. Results. We detail the calculation and apply it to the case of GJ 436b, a classical example of the degeneracy of mass-radius relationships, to illustrate how measurements of $k_2$ can improve our understanding of the interior structure of extrasolar planets. We implemented the method in a code that is fast, freely available, and easy to combine with preexisting interior structure codes. While the linear approach presented here for the calculation of the Love numbers cannot treat the presence of nonlinear effects that may arise under certain dynamical conditions, it is applicable to close-in gaseous extrasolar planets like hot Jupiters, likely the first targets for which $k_2$ will be measured.

Plate tectonics is a fundamental component for the habitability of the Earth. Yet whether it is a recurrent feature of terrestrial bodies orbiting other stars or unique to the Earth is unknown. The stagnant lid may rather be the most common tectonic expression on such bodies. To understand whether a stagnant-lid planet can be habitable, i.e. host liquid water at its surface, we model the thermal evolution of the mantle, volcanic outgassing of H$_2$O and CO$_2$, and resulting climate of an Earth-like planet lacking plate tectonics. We used a 1D model of parameterized convection to simulate the evolution of melt generation and the build-up of an atmosphere of H$_2$O and CO$_2$ over 4.5 Gyr. We then employed a 1D radiative-convective atmosphere model to calculate the global mean atmospheric temperature and the boundaries of the habitable zone (HZ). The evolution of the interior is characterized by the initial production of a large amount of partial melt accompanied by a rapid outgassing of H$_2$O and CO$_2$. At 1 au, the obtained temperatures generally allow for liquid water on the surface nearly over the entire evolution. While the outer edge of the HZ is mostly influenced by the amount of outgassed CO$_2$, the inner edge presents a more complex behaviour that is dependent on the partial pressures of both gases. At 1 au, the stagnant-lid planet considered would be regarded as habitable. The width of the HZ at the end of the evolution, albeit influenced by the amount of outgassed CO$_2$, can vary in a non-monotonic way depending on the extent of the outgassed H$_2$O reservoir. Our results suggest that stagnant-lid planets can be habitable over geological timescales and that joint modelling of interior evolution, volcanic outgassing, and accompanying climate is necessary to robustly characterize planetary habitability.

The Dawn mission confirms predictions that the asteroid 4 Vesta is differentiated with an iron-rich core, a silicate mantle and a basaltic crust, and confirms Vesta as the parent body of the HED meteorites. To better understand its early evolution, we perform numerical calculations of the thermo-chemical evolution adopting new data obtained by the Dawn mission such as mass, bulk density and size of the asteroid. We have expanded the model of Neumann et al. (2012) that includes accretion, compaction, melting and associated changes of material properties and partitioning of 26Al, advective heat transport, and differentiation by porous flow, to include convection and effective cooling in a magma ocean. Depending on the melt fraction, the heat transport by melt segregation is modelled either by porous flow or by convection and heat flux of a magma ocean with a high effective thermal conductivity. We show that partitioning of 26Al and its transport with the silicate melt is crucial for the formation of a magma ocean. Due to the accumulation of 26Al in the sub-surface (for formation times t0<1.5 Ma), a shallow magma ocean with a thickness of 1 to a few tens of km (depending on the silicate melt viscosity) forms. The lifetime of the shallow magma ocean is O(10^4)-O(10^6) years and convection in this layer is accompanied by the extrusion of 26Al at the surface. The interior differentiates from the outside inward with a mantle that is depleted in 26Al and core formation is completed within ~0.3 Ma. The lower mantle experiences melting below 45% suggesting a harzburgitic to dunitic composition. Our results support the formation of eucrites by the extrusion of early partial melt and cumulative eucrites and diogenites may form from the crystallizing shallow magma ocean. Silicate melt is present for up to 150 Ma, and core convects for ~100 Ma, supporting the idea of an early magnetic field.

Radioisotopic ages for meteorites and their components provide constraints on the evolution of small bodies: timescales of accretion, thermal and aqueous metamorphism, differentiation, cooling and impact metamorphism. Realising that the decay heat of short-lived nuclides (e.g. 26Al, 60Fe), was the main heat source driving differentiation and metamorphism, thermal modeling of small bodies is of utmost importance to set individual meteorite age data into the general context of the thermal evolution of their parent bodies, and to derive general conclusions about the nature of planetary building blocks in the early solar system. As a general result, modelling easily explains that iron meteorites are older than chondrites, as early formed planetesimals experienced a higher concentration of short-lived nuclides and more severe heating. However, core formation processes may also extend to 10 Ma after formation of Calcium-Aluminum-rich inclusions (CAIs). A general effect of the porous nature of the starting material is that relatively small bodies (< few km) will also differentiate if they form within 2 Ma after CAIs. A particular interesting feature to be explored is the possibility that some chondrites may derive from the outer undifferentiated layers of asteroids that are differentiated in their interiors. This could explain the presence of remnant magnetization in some chondrites due to a planetary magnetic field.