Thickness and structure of the martian crust from InSight seismic data

Liquid water in the Martian mid-crust

Large volumes of liquid water transiently existed on the surface of Mars more than 3 billion years ago. Much of this water is hypothesized to have been sequestered in the subsurface or lost to space. We use rock physics models and Bayesian inversion to identify combinations of lithology, liquid water saturation, porosity, and pore shape consistent with the constrained mid-crust (∼11.5 to 20 km depths) seismic velocities and gravity near the InSight lander. A mid-crust composed of fractured igneous rocks saturated with liquid water best explains the existing data. Our results have implications for understanding Mars' water cycle, determining the fates of past surface water, searching for past or extant life, and assessing in situ resource utilization for future missions.

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.

Synergies Between Venus & Exoplanetary Observations: Venus and Its Extrasolar Siblings

Here we examine how our knowledge of present day Venus can inform terrestrial exoplanetary science and how exoplanetary science can inform our study of Venus. In a superficial way the contrasts in knowledge appear stark. We have been looking at Venus for millennia and studying it via telescopic observations for centuries. Spacecraft observations began with Mariner 2 in 1962 when we confirmed that Venus was a hothouse planet, rather than the tropical paradise science fiction pictured. As long as our level of exploration and understanding of Venus remains far below that of Mars, major questions will endure. On the other hand, exoplanetary science has grown leaps and bounds since the discovery of Pegasus 51b in 1995, not too long after the golden years of Venus spacecraft missions came to an end with the Magellan Mission in 1994. Multi-million to billion dollar/euro exoplanet focused spacecraft missions such as JWST, and its successors will be flown in the coming decades. At the same time, excitement about Venus exploration is blooming again with a number of confirmed and proposed missions in the coming decades from India, Russia, Japan, the European Space Agency (ESA) and the National Aeronautics and Space Administration (NASA). Here we review what is known and what we may discover tomorrow in complementary studies of Venus and its exoplanetary cousins.

Observational constraints on the process and products of Martian serpentinization

The alteration of olivine-rich rocks to serpentine minerals, (hydr)oxides, and aqueous hydrogen through serpentinization is long thought to have influenced the distribution of habitable environments on early Mars and the evolution of the early Martian hydrosphere and atmosphere. Nevertheless, the planetary importance of Martian serpentinization has remained a matter of debate. To constrain the process and products of Martian serpentinization, we studied serpentinized iron-rich olivines from the 1.1-billion-year Duluth Complex. These data indicate that serpentinized iron-rich olivine would have been accompanied by a fivefold increase in hydrogen production relative to serpentinized terrestrial mantle peridotites. In contrast to previous expectations, this style of serpentinization yields hisingerite as the dominant iron serpentine mineral at comparatively low temperature and pH, consistent with meteorite mineralogy and in situ rover data. The widespread occurrence of oxidized iron-bearing phyllosilicates in highly magnetized regions of the Martian crust supports the hypothesis that serpentinization was more pervasive on early Mars than currently estimated.

Constraints on the martian crust away from the InSight landing site

The most distant marsquake recorded so far by the InSight seismometer occurred at an epicentral distance of 146.3 ± 6.9o, close to the western end of Valles Marineris. On the seismogram of this event, we have identified seismic wave precursors, i.e., underside reflections off a subsurface discontinuity halfway between the marsquake and the instrument, which directly constrain the crustal structure away (about 4100-4500 km) from the InSight landing site. Here we show that the Martian crust at the bounce point between the lander and the marsquake is characterized by a discontinuity at about 20 km depth, similar to the second (deeper) intra-crustal interface seen beneath the InSight landing site. We propose that this 20-km interface, first discovered beneath the lander, is not a local geological structure but likely a regional or global feature, and is consistent with a transition from porous to non-porous Martian crustal materials.

A Positive Feedback Between Crustal Thickness and Melt Extraction for the Origin of the Martian Dichotomy

A North/South difference in crustal thickness is likely at the origin of the Martian dichotomy in topography. Recent crustal thickness maps were obtained by inversion of topography and gravity data seismically anchored at the InSight station. On average, the Martian crust is 51-71 km thick with a southern crust thicker by 18-28 km than the northern one. The origin of this crustal dichotomy is still debated although the hypothesis of a large impact is at present very popular. Here, we propose a new mechanism for the formation of this dichotomy that involves a positive feedback between crustal growth and mantle melting. As the crust is enriched in heat-producing elements, the lid of a one-plate planet is hotter and thinner where the crust is thicker, inducing a larger amount of partial melt below the lid and hence a larger rate of melt extraction and crustal growth. We first demonstrate analytically that larger wavelength perturbations, that is, hemispherical perturbations, grow faster because smaller wavelengths are more attenuated by thermal diffusion. We then use a parameterized thermal evolution model with a well-mixed mantle topped by two different lids characterized by their thermal structures and thicknesses to study the growth of the crust in the two hemispheres. Our results demonstrate that this positive feedback can generate a significant crustal dichotomy.

Late delivery of exotic chromium to the crust of Mars by water-rich carbonaceous asteroids

The terrestrial planets endured a phase of bombardment following their accretion, but the nature of this late accreted material is debated, preventing a full understanding of the origin of inner solar system volatiles. We report the discovery of nucleosynthetic chromium isotope variability (μ⁵⁴Cr) in Martian meteorites that represent mantle-derived magmas intruded in the Martian crust. The μ⁵⁴Cr variability, ranging from -33.1 ± 5.4 to +6.8 ± 1.5 parts per million, correlates with magma chemistry such that samples having assimilated crustal material define a positive μ⁵⁴Cr endmember. This compositional endmember represents the primordial crust modified by impacting outer solar system bodies of carbonaceous composition. Late delivery of this volatile-rich material to Mars provided an exotic water inventory corresponding to a global water layer >300 meters deep, in addition to the primordial water reservoir from mantle outgassing. This carbonaceous material may also have delivered a source of biologically relevant molecules to early Mars.

Surface waves and crustal structure on Mars

We detected surface waves from two meteorite impacts on Mars. By measuring group velocity dispersion along the impact-lander path, we obtained a direct constraint on crustal structure away from the InSight lander. The crust north of the equatorial dichotomy had a shear wave velocity of approximately 3.2 kilometers per second in the 5- to 30-kilometer depth range, with little depth variation. This implies a higher crustal density than inferred beneath the lander, suggesting either compositional differences or reduced porosity in the volcanic areas traversed by the surface waves. The lower velocities and the crustal layering observed beneath the landing site down to a 10-kilometer depth are not a global feature. Structural variations revealed by surface waves hold implications for models of the formation and thickness of the martian crust.

Largest recent impact craters on Mars: Orbital imaging and surface seismic co-investigation

Two >130-meter-diameter impact craters formed on Mars during the later half of 2021. These are the two largest fresh impact craters discovered by the Mars Reconnaissance Orbiter since operations started 16 years ago. The impacts created two of the largest seismic events (magnitudes greater than 4) recorded by InSight during its 3-year mission. The combination of orbital imagery and seismic ground motion enables the investigation of subsurface and atmospheric energy partitioning of the impact process on a planet with a thin atmosphere and the first direct test of martian deep-interior seismic models with known event distances. The impact at 35°N excavated blocks of water ice, which is the lowest latitude at which ice has been directly observed on Mars.

A seismic meteor strike on Mars

A meteor impact and its subsequent seismic waves reveal the crustal structure of Mars

Seismic detection of a deep mantle discontinuity within Mars by InSight

The depth and sharpness of a midmantle seismic discontinuity, associated with the phase transition from mineral olivine to its higher-pressure polymorphs, provide essential clues to understanding the temperature and composition of Martian mantle. Using data from NASA’s InSight mission, we examined five marsquakes located 3,400 to 4,400 km away from the InSight lander and observed triplications of the P and S waves that resulted from the interaction with a seismic discontinuity produced by the postolivine transition. Our observations indicate that the Martian mantle is more iron rich than Earth,and both planets have a similar potential temperature. Our geodynamic modeling further constrains the mantle composition and surface heat flow and indicates that the mantle was cold in the early Noachian.

Constraining the thermal and compositional state of the mantle is crucial for deciphering the formation and evolution of Mars. Mineral physics predicts that Mars’ deep mantle is demarcated by a seismic discontinuity arising from the pressure-induced phase transformation of the mineral olivine to its higher-pressure polymorphs, making the depth of this boundary sensitive to both mantle temperature and composition. Here, we report on the seismic detection of a midmantle discontinuity using the data collected by NASA’s InSight Mission to Mars that matches the expected depth and sharpness of the postolivine transition. In five teleseismic events, we observed triplicated P and S waves and constrained the depth of this discontinuity to be 1,006 ± 40 km by modeling the triplicated waveforms. From this depth range, we infer a mantle potential temperature of 1,605 ± 100 K, a result consistent with a crust that is 10 to 15 times more enriched in heat-producing elements than the underlying mantle. Our waveform fits to the data indicate a broad gradient across the boundary, implying that the Martian mantle is more enriched in iron compared to Earth. Through modeling of thermochemical evolution of Mars, we observe that only two out of the five proposed composition models are compatible with the observed boundary depth. Our geodynamic simulations suggest that the Martian mantle was relatively cold 4.5 Gyr ago (1,720 to 1,860 K) and are consistent with a present-day surface heat flow of 21 to 24 mW/m².

In Situ Regolith Seismic Velocity Measurement at the InSight Landing Site on Mars

Interior exploration using Seismic Investigations, Geodesy and Heat Transport's (InSight) seismometer package Seismic Experiment for Interior Structure (SEIS) was placed on the surface of Mars at about 1.2 m distance from the thermal properties instrument Heat flow and Physical Properties Package (HP3) that includes a self-hammering probe. Recording the hammering noise with SEIS provided a unique opportunity to estimate the seismic wave velocities of the shallow regolith at the landing site. However, the value of studying the seismic signals of the hammering was only realized after critical hardware decisions were already taken. Furthermore, the design and nominal operation of both SEIS and HP3 are nonideal for such high-resolution seismic measurements. Therefore, a series of adaptations had to be implemented to operate the self-hammering probe as a controlled seismic source and SEIS as a high-frequency seismic receiver including the design of a high-precision timing and an innovative high-frequency sampling workflow. By interpreting the first-arriving seismic waves as a P-wave and identifying first-arriving S-waves by polarization analysis, we determined effective P- and S-wave velocities ofv P = 11 9 - 21 + 45  m/s andv S = 6 3 - 7 + 11  m/s, respectively, from around 2,000 hammer stroke recordings. These velocities likely represent bulk estimates for the uppermost several 10s of cm of regolith. An analysis of the P-wave incidence angles provided an independent v P /v S ratio estimate of1.8 4 - 0.35 + 0.89 that compares well with the traveltime based estimate of1.8 6 - 0.25 + 0.42 . The low seismic velocities are consistent with those observed for low-density unconsolidated sands and are in agreement with estimates obtained by other methods.

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.

Early crustal processes revealed by the ejection site of the oldest martian meteorite

The formation and differentiation of the crust of Mars in the first tens of millions of years after its accretion can only be deciphered from incredibly limited records. The martian breccia NWA 7034 and its paired stones is one of them. This meteorite contains the oldest martian igneous material ever dated: ~4.5 Ga old. However, its source and geological context have so far remained unknown. Here, we show that the meteorite was ejected 5-10 Ma ago from the north-east of the Terra Cimmeria-Sirenum province, in the southern hemisphere of Mars. More specifically, the breccia belongs to the ejecta deposits of the Khujirt crater formed 1.5 Ga ago, and it was ejected as a result of the formation of the Karratha crater 5-10 Ma ago. Our findings demonstrate that the Terra Cimmeria-Sirenum province is a relic of the differentiated primordial martian crust, formed shortly after the accretion of the planet, and that it constitutes a unique record of early crustal processes. This province is an ideal landing site for future missions aiming to unravel the first tens of millions of years of the history of Mars and, by extension, of all terrestrial planets, including the Earth.

The History of Water in Martian Magmas From Thorium Maps

Water inventories in Martian magmas are poorly constrained. Meteorite-based estimates range widely, from 102 to >104 ppm H2O, and are likely variably influenced by degassing. Orbital measurements of H primarily reflect water cycled and stored in the regolith. Like water, Th behaves incompatibly during mantle melting, but unlike water Th is not prone to degassing and is relatively immobile during aqueous alteration at low temperature. We employ Th as a proxy for original, mantle-derived H2O in Martian magmas. We use regional maps of Th from Mars Odyssey to assess variations in magmatic water across major volcanic provinces and through time. We infer that Hesperian and Amazonian magmas had ∼100-3,000 ppm H2O, in the lower range of previous estimates. The implied cumulative outgassing since the Hesperian, equivalent to a global H2O layer ∼1-40 m deep, agrees with Mars' present-day surface and near-surface water inventory and estimates of sequestration and loss rates.

Repetitive marsquakes in Martian upper mantle

Marsquakes excite seismic wavefield, allowing the Martian interior structures to be probed. However, the Martian seismic data recorded by InSight have a low signal-to-noise ratio, making the identification of marsquakes challenging. Here we use the Matched Filter technique and Benford’s Law to detect hitherto undetected events. Based on nine marsquake templates, we report 47 newly detected events, >90% of which are associated with the two high-quality events located beneath Cerberus Fossae. They occurred at all times of the Martian day, thus excluding the tidal modulation (e.g., Phobos) as their cause. We attribute the newly discovered, low-frequency, repetitive events to magma movement associated with volcanic activity in the upper mantle beneath Cerberus Fossae. The continuous seismicity suggests that Cerberus Fossae is seismically highly active and that the Martian mantle is mobile.

The authors detect 47 hitherto unreported low-frequency marsquakes originating from Cerberus Fossae at all times of the Martian day. The matched filter technique confirms repetitive events implying that the Martian mantle is dynamically active.

Geometry and Segmentation of Cerberus Fossae, Mars: Implications for Marsquake Properties

The NASA InSight mission to Mars successfully landed on 26 November 2018 in Elysium Planitia. It aims to characterize the seismic activity and aid in the understanding of the internal structure of Mars. We focus on the Cerberus Fossae region, a giant fracture network ∼1,200 km long situated east of the InSight landing site where M ∼3 marsquakes were detected during the past 2 years. It is formed of five main fossae located on the southeast of the Elysium Mons volcanic rise. We perform a detailed mapping of the entire system based on high-resolution satellite images and Digital Elevation Models. The refined cartography reveals a range of morphologies associated with dike activity at depth. Width and throw measurements of the fossae are linearly correlated, suggesting a possible tectonic control on the shapes of the fossae. Widths and throws decrease toward the east, indicating the long-term direction of propagation of the dike-induced graben system. They also give insights into the geometry at depth and how the possible faults and fractures are rooted in the crust. The exceptional preservation of the fossae allows us to detect up to four scales of segmentation, each formed by a similar number of 3-4 segments/subsegments. This generic distribution is comparable to continental faults and fractures on Earth. We anticipate higher stress and potential marsquakes within intersegment zones and at graben tips.

Scattering Attenuation of the Martian Interior through Coda Wave Analysis

We investigate the scattering attenuation characteristics of the Martian crust and uppermost mantle to understand the structure of the Martian interior. We examine the energy decay of the spectral envelopes for 21 high-quality Martian seismic events from Sol 128 to Sol 500 of InSight operations. We use the model of Dainty et al. (1974b) to approximate the behavior of energy envelopes resulting from scattered wave propagation through a single diffusive layer over an elastic half-space. Using a grid search, we mapped the layer parameters that fit the observed InSight data envelopes. The single diffusive layer model provided better fits to the observed energy envelopes for High Frequency (HF) and Very High Frequency (VF) than for the Low Frequency (LF) and Broadband (BB) events. This result is consistent with the suggested source depths (Giardini et al., 2020) for these families of events and their expected interaction with a shallow scattering layer. The shapes of the observed data envelopes do not show a consistent pattern with event distance, suggesting that the diffusivity and scattering layer thickness is non-uniform in the vicinity of InSight at Mars. Given the consistency in the envelope shapes between HF and VF events across epicentral distances and the tradeoffs between the parameters that control scattering, the dimensions of the scattering layer remain unconstrained but require that scattering strength decreases with depth and that the rate of decay in scattering strength is fastest near the surface. This is generally consistent with the processes that would form scattering structures in planetary lithospheres.

The Tharsis mantle source of depleted shergottites revealed by 90 million impact craters

The only martian rock samples on Earth are meteorites ejected from the surface of Mars by asteroid impacts. The locations and geological contexts of the launch sites are currently unknown. Determining the impact locations is essential to unravel the relations between the evolution of the martian interior and its surface. Here we adapt a Crater Detection Algorithm that compile a database of 90 million impact craters, allowing to determine the potential launch position of these meteorites through the observation of secondary crater fields. We show that Tooting and 09-000015 craters, both located in the Tharsis volcanic province, are the most likely source of the depleted shergottites ejected 1.1 million year ago. This implies that a major thermal anomaly deeply rooted in the mantle under Tharsis was active over most of the geological history of the planet, and has sampled a depleted mantle, that has retained until recently geochemical signatures of Mars' early history.

Improving Constraints on Planetary Interiors With PPs Receiver Functions

Seismological constraints obtained from receiver function (RF) analysis provide important information about the crust and mantle structure. Here, we explore the utility of the free-surface multiple of the P-wave (PP) and the corresponding conversions in RF analysis. Using earthquake records, we demonstrate the efficacy of PPs-RFs before illustrating how they become especially useful when limited data is available in typical planetary missions. Using a transdimensional hierarchical Bayesian deconvolution approach, we compute robust P-to-S (Ps)- and PPs-RFs with InSight recordings of five marsquakes. Our Ps-RF results verify the direct Ps converted phases reported by previous RF analyses with increased coherence and reveal other phases including the primary multiple reverberating within the uppermost layer of the Martian crust. Unlike the Ps-RFs, our PPs-RFs lack an arrival at 7.2 s lag time. Whereas Ps-RFs on Mars could be equally well fit by a two- or three-layer crust, synthetic modeling shows that the disappearance of the 7.2 s phase requires a three-layer crust, and is highly sensitive to velocity and thickness of intra-crustal layers. We show that a three-layer crust is also preferred by S-to-P (Sp)-RFs. While the deepest interface of the three-layer crust represents the crust-mantle interface beneath the InSight landing site, the other two interfaces at shallower depths could represent a sharp transition between either fractured and unfractured materials or thick basaltic flows and pre-existing crustal materials. PPs-RFs can provide complementary constraints and maximize the extraction of information about crustal structure in data-constrained circumstances such as planetary missions.

Potential Pitfalls in the Analysis and Structural Interpretation of Seismic Data from the Mars InSight Mission

The Seismic Experiment for Interior Structure (SEIS) of the InSight mission to Mars, has been providing direct information on Martian interior structure and dynamics of that planet since it landed. Compared to seismic recordings on Earth, ground motion measurements acquired by SEIS on Mars are made under dramatically different ambient noise conditions, but include idiosyncratic signals that arise from coupling between different InSight sensors and spacecraft components. This work is to synthesize what is known about these signal types, illustrate how they can manifest in waveforms and noise correlations, and present pitfalls in structural interpretations based on standard seismic analysis methods. We show that glitches, a type of prominent transient signal, can produce artifacts in ambient noise correlations. Sustained signals that vary in frequency, such as lander modes which are affected by variations in temperature and wind conditions over the course of the Martian Sol, can also contaminate ambient noise results. Therefore, both types of signals have the potential to bias interpretation in terms of subsurface layering. We illustrate that signal processing in the presence of identified nonseismic signals must be informed by an understanding of the underlying physical processes in order for high fidelity waveforms of ground motion to be extracted. While the origins of most idiosyncratic signals are well understood, the 2.4 Hz resonance remains debated and the literature does not contain an explanation of its fine spectral structure. Even though the selection of idiosyncratic signal types discussed in this paper may not be exhaustive, we provide guidance on best practices for enhancing the robustness of structural interpretations.

Seismic detection of the martian core

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.

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.