We derive synthetic electrical conductivity, seismic velocity, and density distributions from the results of martian mantle convection models affected by basin-forming meteorite impacts. The electrical conductivity features an intermediate minimum in the strongly depleted topmost mantle, sandwiched between higher conductivities in the lower crust and a smooth increase towards almost constant high values at depths greater than 400 km. The bulk sound speed increases mostly smoothly throughout the mantle, with only one marked change at the appearance of $\beta$-olivine near 1100 km depth. An assessment of the detectability of the subsurface traces of an impact suggests that its signature would be visible in both observables at least if efficient melt extraction from the shock-molten target occurs, but it will not always be particularly conspicuous even for large basins; observations with extensive spatial and temporal coverage would improve their detectability. Electromagnetic sounding may be a more promising method for investigating the properties of impact structures than seismology and, as far as the mantle is concerned, than gravity. It could also help constraining the water content of the mantle in general. By comparing electromagnetic sounding data of an impact structure with model predictions, it might also be possible to answer the open question of the efficiency of impact-generated melt extraction.

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 earliest stage of the evolution of a fully assembled planet is profoundly affected by a number of basin-forming impacts large enough to change the dynamics of its deeper interior. These impacts are in some cases quite closely spaced and follow one another in short time intervals, so that their effects interact and result in behavior that may differ from a simple sum of the effects of two individual and isolated impacts. We use two-dimensional models of mantle convection in a Mars-like planet and a simple parameterized representation of the principal effects of impacts to study some of the dynamical effects and interactions of multiple large impacts. In models of only two impacts, we confirm that the dynamical effects of the impacts reinforce each other the closer they are in space and time but that the effects do not always correspond to straightforward superpositions of those of single, isolated impacts. In models with multiple (4-8) impacts with variable sizes, distances, and frequencies, the global response of the mantle is as variable as the impact sequences in the short term, but in the long term the different evolutionary paths converge for several indicator variables such as the mean flow velocity, temperature, or heat flow. Nonetheless, beyond a certain impact frequency and energy, lithospheric instabilities triggered by large impacts occur on a global scale, reinvigorate mantle dynamics for long time spans, and entail a late stage of melt production in addition to the initial melting stage that is not observed in one- or two-impact models. After one or several very large impacts, some lithospheric material may founder and sink to the core-mantle boundary, and if enough of it accumulates there, it enhances the heat flux out of the core for several hundred millions of years, with possible effects on dynamo activity.

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.

Impactors of different types and sizes can produce a final crater of the same diameter on a planet under certain conditions. We derive the condition for such "isocrater impacts" from scaling laws, as well as relations that describe how the different impactors affect the interior of the target planet; these relations are also valid for impacts that are too small to affect the mantle. The analysis reveals that in a given isocrater impact, asteroidal impactors produce anomalies in the interior of smaller spatial extent than cometary or similar impactors. The differences in the interior could be useful for characterizing the projectile that formed a given crater on the basis of geophysical observations and potentially offer a possibility to help constrain the demographics of the ancient impactor population. A series of numerical models of basin-forming impacts on Mercury, Venus, the Moon, and Mars illustrates the dynamical effects of the different impactor types on different planets. It shows that the signature of large impacts may be preserved to the present in Mars, the Moon, and Mercury, where convection is less vigorous and much of the anomaly merges with the growing lid. On the other hand, their signature will long have been destroyed in Venus, whose vigorous convection and recurring lithospheric instabilities obliterate larger coherent anomalies.

We ran several series of two-dimensional numerical mantle convection simulations representing in idealized form the thermochemical evolution of a Mars-like planet. In order to study the importance of compositional buoyancy of melting mantle, the models were set up in pairs of one including all thermal and compositional contributions to buoyancy and one accounting only for the thermal contributions. In several of the model pairs, single large impacts were introduced as causes of additional strong local anomalies, and their evolution in the framework of the convecting mantle was tracked. The models confirm that the additional buoyancy provided by the depletion of the mantle by regular melting can establish a global stable stratification of the convecting mantle and throttle crust production. Furthermore, the compositional buoyancy is essential in the stabilization and preservation of local compositional anomalies directly beneath the lithosphere and offers a possible explanation for the existence of distinct, long-lived reservoirs in the martian mantle. The detection of such anomalies by geophysical means is probably difficult, however; they are expected to be detected by gravimetry rather than by seismic or heat flow measurements. The results further suggest that the crustal thickness can be locally overestimated by up to ~20 km if impact-induced density anomalies in the mantle are neglected.

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 evolution and escape of the martian atmosphere and the planet's water inventory can be separated into an early and late evolutionary epoch. The first epoch started from the planet's origin and lasted $\sim$500 Myr. Because of the high EUV flux of the young Sun and Mars' low gravity it was accompanied by hydrodynamic blow-off of hydrogen and strong thermal escape rates of dragged heavier species such as O and C atoms. After the main part of the protoatmosphere was lost, impact-related volatiles and mantle outgassing may have resulted in accumulation of a secondary CO$_2$ atmosphere of a few tens to a few hundred mbar around $\sim$4--4.3 Gyr ago. The evolution of the atmospheric surface pressure and water inventory of such a secondary atmosphere during the second epoch which lasted from the end of the Noachian until today was most likely determined by a complex interplay of various nonthermal atmospheric escape processes, impacts, carbonate precipitation, and serpentinization during the Hesperian and Amazonian epochs which led to the present day surface pressure.

We compare estimates of atmospheric precipitation during the Martian Noachian-Hesperian boundary 3.8 Gyr ago as calculated in a radiative-convective column model of the atmosphere with runoff values estimated from a geomorphological analysis of dendritic valley network discharge rates. In the atmospheric model, we assume CO2-H2O-N2 atmospheres with surface pressures varying from 20 mb to 3 bar with input solar luminosity reduced to 75% the modern value. Results from the valley network analysis are of the order of a few mm d-1 liquid water precipitation (1.5-10.6 mm d-1, with a median of 3.1 mm d-1). Atmospheric model results are much lower, from about 0.001-1 mm d-1 of snowfall (depending on CO2 partial pressure). Hence, the atmospheric model predicts a significantly lower amount of precipitated water than estimated from the geomorphological analysis. Furthermore, global mean surface temperatures are below freezing, i.e. runoff is most likely not directly linked to precipitation. Therefore, our results strongly favor a cold early Mars with episodic snowmelt as a source for runoff. Our approach is challenged by mostly unconstrained parameters, e.g. greenhouse gas abundance, global meteorology (for example, clouds) and planetary parameters such as obliquity- which affect the atmospheric result - as as well as by inherent problems in estimating discharge and runoff on ancient Mars, such as a lack of knowledge on infiltration and evaporation rates and on flooding timescales, which affect the geomorphological data. Nevertheless, our work represents a first step in combining and interpreting quantitative tools applied in early Mars atmospheric and geomorphological studies.

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 [...]

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.