Geophysical evidence for an enriched molten silicate layer above Mars's core.

The detection of deep reflected S waves on Mars inferred a core size of 1,830 ± 40 km (ref. 1), requiring light-element contents that are incompatible with experimental petrological constraints. This estimate assumes a compositionally homogeneous Martian mantle, at odds with recent measurements of anomalously slow propagating P waves diffracted along the core-mantle boundary2. An alternative hypothesis is that Mars's mantle is heterogeneous as a consequence of an early magma ocean that solidified to form a basal layer enriched in iron and heat-producing elements. Such enrichment results in the formation of a molten silicate layer above the core, overlain by a partially molten layer3. Here we show that this structure is compatible with all geophysical data, notably (1) deep reflected and diffracted mantle seismic phases, (2) weak shear attenuation at seismic frequency and (3) Mars's dissipative nature at Phobos tides. The core size in this scenario is 1,650 ± 20 km, implying a density of 6.5 g cm-3, 5-8% larger than previous seismic estimates, and can be explained by fewer, and less abundant, alloying light elements than previously required, in amounts compatible with experimental and cosmochemical constraints. Finally, the layered mantle structure requires external sources to generate the magnetic signatures recorded in Mars's crust.

Marsquake Locations and 1-D Seismic Models for Mars From InSight Data.

We present inversions for the structure of Mars using the first Martian seismic record collected by the InSight lander. We identified and used arrival times of direct, multiples, and depth phases of body waves, for 17 marsquakes to constrain the quake locations and the one-dimensional average interior structure of Mars. We found the marsquake hypocenters to be shallower than 40 km depth, most of them being located in the Cerberus Fossae graben system, which could be a source of marsquakes. Our results show a significant velocity jump between the upper and the lower part of the crust, interpreted as the transition between intrusive and extrusive rocks. The lower crust makes up a significant fraction of the crust, with seismic velocities compatible with those of mafic to ultramafic rocks. Additional constraints on the crustal thickness from previous seismic analyses, combined with modeling relying on gravity and topography measurements, yield constraints on the present-day thermochemical state of Mars and on its long-term history. Our most constrained inversion results indicate a present-day surface heat flux of 22 ± 1 mW/m2, a relatively hot mantle (potential temperature: 1740 ± 90 K) and a thick lithosphere (540 ± 120 km), associated with a lithospheric thermal gradient of 1.9 ± 0.3 K/km. These results are compatible with recent seismic studies using a reduced data set and different inversion approaches, confirming that Mars' potential mantle temperature was initially relatively cold (1780 ± 50 K) compared to that of its present-day state, and that its crust contains 10-12 times more heat-producing elements than the primitive mantle.

A Reconstruction Algorithm for Temporally Aliased Seismic Signals Recorded by the InSight Mars Lander.

In December 2018, the NASA InSight lander successfully placed a seismometer on the surface of Mars. Alongside, a hammering device was deployed at the landing site that penetrated into the ground to attempt the first measurements of the planetary heat flow of Mars. The hammering of the heat probe generated repeated seismic signals that were registered by the seismometer and can potentially be used to image the shallow subsurface just below the lander. However, the broad frequency content of the seismic signals generated by the hammering extends beyond the Nyquist frequency governed by the seismometer's sampling rate of 100 samples per second. Here, we propose an algorithm to reconstruct the seismic signals beyond the classical sampling limits. We exploit the structure in the data due to thousands of repeated, only gradually varying hammering signals as the heat probe slowly penetrates into the ground. In addition, we make use of the fact that repeated hammering signals are sub-sampled differently due to the unsynchronized timing between the hammer strikes and the seismometer recordings. This allows us to reconstruct signals beyond the classical Nyquist frequency limit by enforcing a sparsity constraint on the signal in a modified Radon transform domain. In addition, the proposed method reduces uncorrelated noise in the recorded data. Using both synthetic data and actual data recorded on Mars, we show how the proposed algorithm can be used to reconstruct the high-frequency hammering signal at very high resolution.

Upper mantle structure of Mars from InSight seismic data.

For 2 years, the InSight lander has been recording seismic data on Mars that are vital to constrain the structure and thermochemical state of the planet. We used observations of direct (P and S) and surface-reflected (PP, PPP, SS, and SSS) body-wave phases from eight low-frequency marsquakes to constrain the interior structure to a depth of 800 kilometers. We found a structure compatible with a low-velocity zone associated with a thermal lithosphere much thicker than on Earth that is possibly related to a weak S-wave shadow zone at teleseismic distances. By combining the seismic constraints with geodynamic models, we predict that, relative to the primitive mantle, the crust is more enriched in heat-producing elements by a factor of 13 to 20. This enrichment is greater than suggested by gamma-ray surface mapping and has a moderate-to-elevated surface heat flow.