The evolution of the interior of stagnant-lid bodies is comparatively easier to model and predict with respect to the Earth's due to the absence of the large uncertainties associated with the physics of plate tectonics, its onset time and efficiency over the planet's history. Yet, the observational record for these bodies is both scarcer and sparser with respect to the Earth's. It is restricted to a limited number of samples and a variety of remote-sensing measurements of billions-of-years-old surfaces whose actual age is difficult to determine precisely. Combining these observations into a coherent picture of the thermal and convective evolution of the planetary interior represents thus a major challenge. In this chapter, we review key processes and (mostly geophysical) observational constraints that can be used to infer the global characteristics of mantle convection and thermal evolution of the interior of Mercury, the Moon and Mars, the three major terrestrial bodies where a stagnant lid has likely been present throughout most of their history.

We explore the application of machine learning based on mixture density neural networks (MDNs) to the interior characterization of low-mass exoplanets up to 25 Earth masses constrained by mass, radius, and fluid Love number $k_2$. We create a dataset of 900$\:$000 synthetic planets, consisting of an iron-rich core, a silicate mantle, a high-pressure ice shell, and a gaseous H/He envelope, to train a MDN using planetary mass and radius as inputs to the network. For this layered structure, we show that the MDN is able to infer the distribution of possible thicknesses of each planetary layer from mass and radius of the planet. This approach obviates the time-consuming task of calculating such distributions with a dedicated set of forward models for each individual planet. While gas-rich planets may be characterized by compositional gradients rather than distinct layers, the method presented here can be easily extended to any interior structure model. The fluid Love number $k_2$ bears constraints on the mass distribution in the planets' interior and will be measured for an increasing number of exoplanets in the future. Adding $k_2$as an input to the MDN significantly decreases the degeneracy of the possible interior structures.

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.

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.

The magma ocean (MO) is a crucial stage in the build-up of terrestrial planets. Its solidification and the accompanying outgassing of volatiles set the conditions for important processes occurring later or even simultaneously, such as solid-state mantle convection and atmospheric escape. To constrain the duration of a global-scale Earth MO we have built and applied a 1D interior model coupled alternatively with a grey H2O/CO2 atmosphere or with a pure H2O atmosphere treated with a line-by-line model described in a companion paper by Katyal et al. (2019). We study in detail the effects of several factors affecting the MO lifetime, such as the initial abundance of H2O and CO2, the convection regime, the viscosity, the mantle melting temperature, and the longwave radiation absorption from the atmosphere. In this specifically multi-variable system we assess the impact of each factor with respect to a reference setting commonly assumed in the literature. We find that the MO stage can last from a few thousand to several million years. By coupling the interior model with the line-by-line atmosphere model, we identify the conditions that determine whether the planet experiences a transient magma ocean or it ceases to cool and maintains a continuous magma ocean. We find a dependence of this distinction simultaneously on the mass of the outgassed H2O atmosphere and on the MO surface melting temperature. We discuss their combined impact on the MO's lifetime in addition to the known dependence on albedo, orbital distance and stellar luminosity and we note observational degeneracies that arise thereby for target exoplanets.

The evolution of Earth's early atmosphere and the emergence of habitable conditions on our planet are intricately coupled with the development and duration of the magma ocean phase during the early Hadean period (4 to 4.5 Ga). In this paper, we deal with the evolution of the steam atmosphere during the magma ocean period. We obtain the outgoing longwave radiation using a line-by-line radiative transfer code GARLIC. Our study suggests that an atmosphere consisting of pure H$_{2}$O, built as a result of outgassing extends the magma ocean lifetime to several million years. The thermal emission as a function of solidification timescale of magma ocean is shown. We study the effect of thermal dissociation of H$_{2}$O at higher temperatures by applying atmospheric chemical equilibrium which results in the formation of H$_{2}$ and O$_{2}$ during the early phase of the magma ocean. A 1-6\% reduction in the OLR is seen. We also obtain the effective height of the atmosphere by calculating the transmission spectra for the whole duration of the magma ocean. An atmosphere of depth ~100 km is seen for pure water atmospheres. The effect of thermal dissociation on the effective height of the atmosphere is also shown. Due to the difference in the absorption behavior at different altitudes, the spectral features of H$_{2}$ and O$_{2}$ are seen at different altitudes of the atmosphere. Therefore, these species along with H$_{2}$O have a significant contribution to the transmission spectra and could be useful for placing observational constraints upon magma ocean exoplanets.

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.

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.

Planetary embryos form protoplanets via mutual collisions, which can lead to the development of magma oceans. During their solidification, large amounts of the mantles' volatile contents may be outgassed. The resulting H$_2$O/CO$_2$ dominated steam atmospheres may be lost efficiently via hydrodynamic escape due to the low gravity and the high stellar EUV luminosities. Protoplanets forming later from such degassed building blocks could therefore be drier than previously expected. We model the outgassing and subsequent hydrodynamic escape of steam atmospheres from such embryos. The efficient outflow of H drags along heavier species (O, CO$_2$, noble gases). The full range of possible EUV evolution tracks of a solar-mass star is taken into account to investigate the escape from Mars-sized embryos at different orbital distances. The envelopes are typically lost within a few to a few tens of Myr. Furthermore, we study the influence on protoplanetary evolution, exemplified by Venus. We investigate different early evolution scenarios and constrain realistic cases by comparing modeled noble gas isotope ratios with observations. Starting from solar values, consistent isotope ratios (Ne, Ar) can be found for different solar EUV histories, as well as assumptions about the initial atmosphere (either pure steam or a mixture with accreted H). Our results generally favor an early accretion scenario with a small amount of accreted H and a low-activity Sun, because in other cases too much CO$_2$ is lost during evolution, which is inconsistent with Venus' present atmosphere. Important issues are likely the time at which the initial steam atmosphere is outgassed and/or the amount of CO$_2$ which may still be delivered at later evolutionary stages. A late accretion scenario can only reproduce present isotope ratios for a highly active young Sun, but then very massive steam atmospheres would be required.

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 impact heat accumulated during the late stage of planetary accretion can melt a significant part or even the entire mantle of a terrestrial body, giving rise to a global magma ocean. [...] Assuming fractional crystallization of the magma ocean, dense cumulates are produced close to the surface, largely due to iron enrichment in the evolving magma ocean liquid (Elkins-Tanton et al., 2003). A gravitationally unstable mantle thus forms, which is prone to overturn. We investigate the cumulate overturn and its influence on the thermal evolution of Mars using mantle convection simulations in 2D cylindrical geometry. We present a suite of simulations using different initial conditions and a strongly temperature-dependent viscosity. We assume that all radiogenic heat sources have been enriched during the freezing-phase of the magma ocean in the uppermost 50 km and that the initial steam-atmosphere created by the degassing of the freezing magma ocean was rapidly lost, implying that the surface temperature is set to present-day values. In this case, a stagnant lid forms rapidly on top of the convective interior preventing the uppermost dense cumulates to sink, even when allowing for a plastic yielding mechanism. Below this dense stagnant lid, the mantle chemical gradient settles to a stable configuration. The convection pattern is dominated by small-scale structures, which are difficult to reconcile with the large-scale volcanic features observed over Mars' surface and partial melting ceases in less than 900 Ma. Assuming that the stagnant lid can break because of additional mechanisms and allowing the uppermost dense layer to overturn, a stable density gradient is obtained, with the densest material and the entire amount of heat sources lying above the CMB. This stratification leads to a strong overheating of the lowermost mantle [...]