A planet’s crust bears witness to the history of planetary formation and evolution, but for Mars, no absolute measurement of crustal thickness has been available. Here, we determine the structure of the crust beneath the InSight landing site on Mars using both marsquake recordings and the ambient wavefield. By analyzing seismic phases that are reflected and converted at subsurface interfaces, we find that the observations are consistent with models with at least two and possibly three interfaces. If the second interface is the boundary of the crust, the thickness is 20 ± 5 kilometers, whereas if the third interface is the boundary, the thickness is 39 ± 8 kilometers. Global maps of gravity and topography allow extrapolation of this point measurement to the whole planet, showing that the average thickness of the martian crust lies between 24 and 72 kilometers. Independent bulk composition and geodynamic constraints show that the thicker model is consistent with the abundances of crustal heat-producing elements observed for the shallow surface, whereas the thinner model requires greater concentration at depth.

Clues to a planet’s geologic history are contained in its interior structure, particularly its core. We detected reflections of seismic waves from the core-mantle boundary of Mars using InSight seismic data and inverted these together with geodetic data to constrain the radius of the liquid metal core to 1830 ± 40 kilometers. The large core implies a martian mantle mineralogically similar to the terrestrial upper mantle and transition zone but differing from Earth by not having a bridgmanite-dominated lower mantle. We inferred a mean core density of 5.7 to 6.3 grams per cubic centimeter, which requires a substantial complement of light elements dissolved in the iron-nickel core. The seismic core shadow as seen from InSight’s location covers half the surface of Mars, including the majority of potentially active regions—e.g., Tharsis—possibly limiting the number of detectable marsquakes.

ol 185 was a typical sol on Mars (a Mars sol is 24 h 39.5 min long, and we number sols starting from landing). The ground acceleration spectrogram recorded by the very broadband (VBB) instrument of SEIS 1-3 (Seismic Experiment for Interior Structure; Fig. 1a) is dominated by the noise produced by the weakly turbulent night-time winds and by the powerful, thermally driven convective turbulence during the day 4. Around 17:00 local mean solar time (lmst), the wind fluctuations die out quite suddenly and the planet remains very quiet into the early night hours. Several distinctive features can be seen every sol on Mars. Lander vibrations activated by the wind appear as horizontal thin lines with frequency varying daily as a result of temperature variations of the lander; almost invisible during quiet hours, they are not excited by seismic events (for example, the lander mode at 4 Hz in Fig. 1a). We also observe a pronounced ambient resonance at 2.4 Hz, strongest on the vertical component, with no clear link to wind strength but excited by all the seismic vibrations at that frequency. The relative excitations of the 2.4 Hz and 4 Hz modes serve as discriminants for the origin of ground vibrations recorded by SEIS, allowing us to distinguish between local vibrations induced by atmospheric or lander activity and more distant sources of ground vibrations. On Sol 185, two weak events can also be spotted in the quiet hours of the early evening, one with a broadband frequency content and a second 80 min later, centred on the 2.4 Hz resonance band (Fig. 1a).

Aims. We aim to present a generalized Bayesian inference method for constraining interiors of super Earths and sub-Neptunes. Our methodology succeeds in quantifying the degeneracy and correlation of structural parameters for high dimensional parameter spaces. Specifically, we identify what constraints can be placed on composition and thickness of core, mantle, ice, ocean, and atmospheric layers given observations of mass, radius, and bulk refractory abundance constraints (Fe, Mg, Si) from observations of the host star's photospheric composition. Methods. We employed a full probabilistic Bayesian inference analysis that formally accounts for observational and model uncertainties. Using a Markov chain Monte Carlo technique, we computed joint and marginal posterior probability distributions for all structural parameters of interest. We included state-of-the-art structural models based on self-consistent thermodynamics of core, mantle, high-pressure ice, and liquid water. Furthermore, we tested and compared two different atmospheric models that are tailored for modeling thick and thin atmospheres, respectively. Results. First, we validate our method against Neptune. Second, we apply it to synthetic exoplanets of fixed mass and determine the effect on interior structure and composition when (1) radius; (2) atmospheric model; (3) data uncertainties; (4) semi-major axes; (5) atmospheric composition (i.e., a priori assumption of enriched envelopes versus pure H/He envelopes); and (6) prior distributions are varied. Conclusions. Our main conclusions are: (1) given available data, the range of possible interior structures is large; quantification of the degeneracy of possible interiors is therefore indispensable for meaningful planet characterization. (2) Our method predicts models that agree with independent estimates of Neptune's interior. (3) Increasing the precision in mass and radius leads to much improved constraints on ice mass fraction, size of rocky interior, but little improvement in the composition of the gas layer, whereas an increase in the precision of stellar abundances enables to better constrain mantle composition and relative core size; (4) for thick atmospheres, the choice of atmospheric model can have significant influence on interior predictions, including the rocky and icy interior. The preferred atmospheric model is determined by envelope mass. This study provides a methodology for rigorously analyzing general interior structures of exoplanets which may help to understand how exoplanet interior types are distributed among star systems. This study is relevant in the interpretation of future data from missions such as TESS, CHEOPS, and PLATO.

We invert the Martian tidal response and mean mass and moment of inertia for chemical composition, thermal state, and interior structure. The inversion combines phase equilibrium computations with a laboratory‐based viscoelastic dissipation model. The rheological model, which is based on measurements of anhydrous and melt‐free olivine, is both temperature and grain size sensitive and imposes strong constraints on interior structure. The bottom of the lithosphere, defined as the location where the conductive geotherm meets the mantle adiabat, occurs deep within the upper mantle (∼200–400 km depth) resulting in apparent upper mantle low‐velocity zones. Assuming an Fe‐FeS core, our results indicate (1) a mantle with a Mg# (molar Mg/Mg+Fe) of ∼0.75 in agreement with earlier geochemical estimates based on analysis of Martian meteorites; (2) absence of bridgmanite‐ and ferropericlase‐dominated basal layer; (3) core compositions (15–18.5 wt% S), core radii (1,730–1,840 km), and core‐mantle boundary temperatures (1620–1690°C) that, together with the eutectic‐like core compositions, suggest that the core is liquid; and (4) bulk Martian compositions with a Fe/Si (weight ratio) of 1.66–1.81. We show that the inversion results can be used in tandem with geodynamic simulations to identify plausible geodynamic scenarios and parameters. Specifically, we find that the inversion results are largely reproducible by stagnant lid convection models for a range of initial viscosities (∼1018–1020 Pa s) and radioactive element partitioning between crust and mantle around 0.01–0.1. The geodynamic models predict a mean surface heat flow between 15 and 25 mW/m2.