Europlanet Science Congress 2022
Palacio de Congresos de Granada, Spain
18 – 23 September 2022
Europlanet Science Congress 2022
Palacio de Congresos de Granada, Spain
18 September – 23 September 2022
Planetary Seismology and Geophysics


The investigations of planet's interior are key to understand the origin and evolution of the planet. Seismology has proven to be a very powerful tool for determining the internal structure of Earth, leading to the direct measurement of the thickness of crust, mantle and core, as well as obtaining detailed 3-dimensional structure models.

The first successful extraterrestrial seismic experiment was carried out by the Apollo missions from 1969 to 1972, when a network of seismometers was deployed along with other in situ geophysical instruments. A renaissance of planetary seismology occurred almost 50 years later, when the InSight (Interior Exploration using Seismic Investigations, Geodesy and Heat Transport) lander deployed successfully a seismometer on the surface of Mars. Since its landing, on November 26, 2018, data from this mission have shown that Mars is a seismically active planet, and have placed key constraints on the thickness of the crust, core radius, and the presence of a low velocity zone in the upper mantle. Improved constraints for the interior structure of both the Moon and Mars have since been obtained through the joint analysis of seismic datasets with other geophysical observables, such as gravity, surface topography, magnetic field measurements, and heat flow.

This session invites abstracts concerning the following topics:
- Studies based on the analysis of existing extraterrestrial seismic datasets.
- Analyses of geophysical data that rely on seismic constraints or that offer information for the design of future seismic experiments.
- Works that propose new methodologies and mission concepts for future extraterrestrial seismic experiments.

Co-organized by MITM
Convener: Philippe Lognonné | Co-conveners: Melanie Drilleau, Foivos Karakostas, Mark Panning, Simon C. Stähler, Mark Wieczorek
| Tue, 20 Sep, 15:30–17:00 (CEST)|Room Machado, Wed, 21 Sep, 10:00–11:30 (CEST)|Room Machado
| Attendance Mon, 19 Sep, 18:45–20:15 (CEST) | Display Mon, 19 Sep, 08:30–Wed, 21 Sep, 11:00|Poster area Level 1

Session assets

Discussion on Slack

Orals: Tue, 20 Sep | Room Machado

Chairpersons: Philippe Lognonné, Mark Panning, Mark Wieczorek
Bruce Banerdt and the InSight Science Team

After an auspicious beginning, with the Ranger, Apollo and Viking missions carrying multiple seismometers into space, planetary seismology lay dormant for 40 years, save for a cruelly-extinguished hope on the Mars’96 mission. InSight put an end to this seismological drought with its landing on Mars in November, 2018, and subsequent deployment of the SEIS instrument two months later (Banerdt et al., 2020; Lognonné et al., 2019, 2020).

After a slow start (uncannily, no seismic events were detected for the first two months after deployment), the subsequent three-plus years  have been an unqualified triumph for planetary seismology. Roughly 2700 seismic events have been detected and cataloged (Clinton et al., 2021; Ceylan et al., 2022), ranging from more than a thousand tiny SF events, thought to be locally-sourced thermal events, to a massive magnitude 5 marsquake more than 2000 km away, but still with a signal-to-noise ratio over 160,000. This catalog includes more than 30 Low Frequency and Broad-Band marsquake signals of quality A or B, meaning they have clear body-wave phase arrivals and, in the case of quality A, identifiable polarization. These have been clearly identified as teleseismic events, which taken together describe a set of ray paths that have probed the interior of the planet from the surface down to its core.

The analysis of this seismic data set has already transformed our knowledge of the interior structure of Mars. We have for the first time been able to directly measure the thickness of the crust beneath InSight (Knapmeyer-Endrun, et al., 2021; Kim et al., 2021; Durán et al., 2022; Drilleau et al., 2022) and the size of the martian core (Stähler et al., 2021). Analysis of the data has further allowed us to place tight constraints on the density of the core and crust, and the thermal profile of the mantle (Stähler et al., 2021; Khan et al., 2021, 2022; Wieczorek et al., 2022). In addition, SEIS’s observations have proved invaluable for studying the near-surface structure beneath InSight on scales from tens of centimeters to hundreds of meters (e.g., Lognonné et al., 2020; Banerdt et al., 2020; Kenda et al., 2020; Hobiger et al., 2021; Murdoch et al., 2021; Compaire et al., 2021, 2022).

InSight’s seismometer has also been able to study the current tectonics of Mars, providing a measure of the size distribution of seismic events (and thus a window into the rate of moment release) and a map of their geographic distribution (in particular, the strong concentration of Eastern hemisphere activity around Cerberus Fossae; Giardini et al., 2020; Stähler et al., 2022) and even revealing the source mechanisms responsible for some marsquakes, from which the orientation of stresses in the lithosphere can be inferred (Brinkman et al., 2021).

InSight and SEIS have proven conclusively the value of planetary seismology in general, and a single high-quality seismometer in particular, in understanding the solid planets of our solar system. The results described in this talk are only a subset of the work that has been done, and was derived from just the first round of analysis of SEIS data. As the seismic community becomes more familiar with the data set and it becomes more widely used, we can expect new insights into the geophysics of Mars for decades to come.



Banerdt, W. B., and Smrekar, S. E., et al., 2020. Initial results from the InSight mission on Mars. Nat. Geosci. 13, 183–189.

Brinkman, N., et al., 2021. First focal mechanisms of marsquakes, J. Geophys. Res. Planets, e2020JE006546.

Ceylan, S., et al., 2022. The marsquake catalogue from InSight, sols 0-1011, Phys. Earth Planet. Int., in press.

Clinton, J., et al., 2021. The Marsquake Catalogue from InSight, Sols 0-478, Phys. Earth Planet. Int. 310, 106597.

Compaire, N., et al., 2021. Autocorrelation of the ground vibrations recorded by the SEIS-InSight seismometer on Mars, J. Geophys. Res. Planets 126, e2020JE006498.

Compaire, N., et al., 2022. Seasonal variations of subsurface seismic velocities monitored by the SEIS-InSight seismometer on Mars, Geophys. J. Int. 229:2, 776–799.

Drilleau, M., et al., 2022. Marsquake location and 1-D seismic models for Mars from InSight data, J. Geophys. Res. Planets, in press.

Durán, C., et al., 2022. Seismology on Mars: An analysis of direct, reflected, and converted seismic body waves with implications for interior structure, Phys. Earth Planet. Int., in press.

Hobiger, M., et al., 2021. The shallow structure of Mars from inversion of high-frequency ambient seismic vibrations Rayleigh wave ellipticity at the InSight landing site, Nat. Comm. 12, 6756,.

Kenda, B., et al., 2020. Subsurface structure at the InSight landing site from compliance measurements by seismic and meteorological experiments, J. Geophys. Res. Planets 125, 6, e2020JE006387.

Khan, A., et al., 2021. Upper mantle structure of Mars from InSight seismic data, Science, 373, 434-438.

Khan, A., et al., 2022. Geophysical and cosmochemical constraints on the bulk composition of Mars , Earth Planet. Sci. Lett. 578.

Kim, D., et al., 2021. Improving constraints on planetary interiors with PPs receiver functions, J. Geophys. Res. Planets 126, e2021JE006983.

Knapmeyer-Endrun, B., et al., 2021. Thickness and structure of the martian crust from InSight seismic data, Science, 373, 438-443.

Lognonné, P., et al., 2020. Constraints on the shallow elastic and anelastic structure of Mars from InSight seismic data. Nat. Geosci. 13, 213–220.

Murdoch, N., et al., 2021. Constraining Martian regolith and vortex parameters from combined seismic and meteorological measurements, J. Geophys. Res. Planets 126, e2020JE006410.

Stähler, S. C., et al., 2021. Seismic detection of the martian core, Science, 373, 443-448.

Stähler, S., et al., 2022. Marsquakes indicate dike-induced tectonics in Cerberus Fossae, Mars, Nature Geoscience, in press.

Wieczorek, M., et al., 2022. InSight constraints on the global character of the Martian crust, J. Geophys. Res. Planets, e2022JE007298.

How to cite: Banerdt, B. and the InSight Science Team: InSight’s Contributions to Planetary Seismology and Geophysics, Europlanet Science Congress 2022, Granada, Spain, 18–23 Sep 2022, EPSC2022-773,, 2022.

Attilio Rivoldini, Sébastien Le Maistre, Alfonso Caldiero, Marie Yseboodt, Rose-Marie Baland, Mikael Beuthe, Tim Van Hoolst, Veronique Dehant, William Folkner, Dustin Buccino, Daniel Kahan, Jean-Charles Marty, Daniele Antonangeli, James Badro, Melanie Drilleau, Alex Konopliv, Marie-Julie Peters, Ana-Catalina Plesa, Henri Samuel, and Nicola Tosi and the InSight/RISE team

We report the results of more than 2 years of monitoring the rotation of Mars with the RISE instrument on InSight. Small periodic variations of the spin axis orientation, called nutations, can be extracted from the Doppler data with enough precision to identify the influence of the Martian fluid core. For the first time for a planetary body other than the Earth, we can measure the period of the Free Core Nutation (FCN), which is a rotational normal mode arising from the misalignment of the rotation axes of the core and mantle. In this way, we confirm the liquid state of the core and estimate its moment of inertia as well as its size. 

The FCN period depends on the dynamical flattening of the core and on its ability to deform. Since the shape and gravity field of Mars deviate significantly from those of a uniformly rotating fluid body, deviations from that state can also be expected for the core. Models accounting for the dynamical shape of Mars can thus be tested by comparing core shape predictions to nutation constraints. The observed FCN period can be accounted for by interior models having a very thick lithosphere loaded by degree-two mass anomalies at the bottom. 

The combination of nutation data and interior structure modeling allows us to deduce the radius of the core and to constrain its density, and thus, to address the nature and abundance of light elements alloyed to iron. The inferred core radius agrees with previous estimates based on geodesy and seismic data. The large fraction of light elements required to match the core density implies that its liquidus is significantly lower than the expected core temperature, making the presence of an inner core highly unlikely. Besides, the existence of an inner core would lead to an additional rotational normal mode the signature of which has not been detected in the RISE data. 

How to cite: Rivoldini, A., Le Maistre, S., Caldiero, A., Yseboodt, M., Baland, R.-M., Beuthe, M., Van Hoolst, T., Dehant, V., Folkner, W., Buccino, D., Kahan, D., Marty, J.-C., Antonangeli, D., Badro, J., Drilleau, M., Konopliv, A., Peters, M.-J., Plesa, A.-C., Samuel, H., and Tosi, N. and the InSight/RISE team: A view into the deep interior of Mars from nutation measured by InSight RISE, Europlanet Science Congress 2022, Granada, Spain, 18–23 Sep 2022, EPSC2022-1101,, 2022.

Anna Horleston, John Clinton, Savas Ceylan, Taichi Kawamura, Simon C. Stähler, Constantinos Charalambous, Nikolaj L. Dahmen, Cecilia Duran, Doyeon Kim, Matthieu Plasman, Géraldine Zenhäusern, Fabian Euchner, Martin Knapmeyer, Domenico Giardini, Philippe Lognonné, William T. Pike, Mark Panning, Suzanne Smrekar, and William B. Banerdt

Introduction: NASA’s InSight lander has been operating on the surface of Mars for over 1200 Martian solar days (sols). Despite decreasing power, the lander has acquired near continuous seismic data throughout the mission to date. The Marsquake Service (MQS) performs daily analysis of the data and catalogues the seismicity. The catalogue currently contains over 1300 local, regional and teleseismic marsquakes and more than 1300 events that are likely associated with thermal cracking close to the lander (Table 1). Here we will present the latest update on the marsquake catalogue.

Martian Seismic Data: The Martian seismic dataset is often heavily contaminated by atmospheric effects, especially during the daytime and in the Northern hemisphere autumn and winter (Figure 1)[1,2]. This environmental noise can obscure the seismicity which is often at low amplitudes. During spring and summer, the atmospheric conditions are much more favourable from shortly before sunset until approximately sunrise and much of the MQS catalogue is detected during these hours. Other transient effects on the data that can interfere with event detection are discussed in [3]. One notable feature of the seismic noise on Mars is a resonance at 2.4 Hz which is observed during quiet periods and is also excited by seismic events that contain energy at that frequency. This resonance is likely caused by sub-surface structure [4] and greatly enhances the detection rate of higher frequency events.

Figure 1: Daily spectrograms of VBB vertical component stacked from sol 72 to sol 1233 with catalogued marsquakes overlain. White lines are where data were not recorded. The black dotted line represents the equivalent sol in mission year 1 to sol 1233. The black dashed line marks 1 year since SEIS was fully deployed.


Marsquake Types and Qualities: MQS categorises marsquakes based on their frequency content and their location quality [2,5]. A full description of the categories and qualities can be found in [2], here we list the general characteristics:

Low Frequency (LF): All energy is below 2.4 Hz and most of the time below 1 Hz.

Broadband (BB): Energy is predominantly at low frequencies but can include energy at and above 2.4 Hz.

2.4 Hz: Energy is contained in a narrow frequency band around 2.4 Hz, reaching a maximum of 4 Hz. These events are thought to be low amplitude versions of High Frequency events.

High Frequency (HF): Energy is predominantly at 2.4 Hz but also extends above 4 Hz. There may also be energy at low frequencies but the dominant energy is at higher frequencies (as opposed to BB events).

Very High Frequency (VF): Energy is at high frequencies, sometimes reaching up to 30 Hz. Energy is strongest on the horizontal components, especially at higher frequencies.

Super High Frequency (SF): very short duration (~10s of seconds) events with mainly horizontal energy between 5-10 Hz. These events are highly temperature dependent and are likely due to thermal cracks observed in the near surface close to the lander [6].

For all except SF events a location quality is assigned based on the ability to identify phase picks and locate the event as follows:

Quality A: Multiple clear phases are identified, and clear polarisation can be defined. The source is well located with distance and backazimuth.

Quality B: Either the distance or the backazimuth is clear but not both. A full location cannot be provided. When clear phases are identified distance is calculated.

Quality C: Clear seismic signal is observed but phases are not readily assigned, nor polarisation identified. A distance may be assigned but will be uncertain.

Quality D: The signal has low signal to noise either due to a low amplitude signal or high contamination from environmental noise. No location is given and there is a small risk that the event may not be tectonic in origin.

The Marsquake Catalogue: The marsquake catalogue is released every three months with a three-month delay. The latest release is V10 [7]. The catalogue includes not only the type and quality of the events but also general information such as location, origin, magnitude and any secondary phases that have been picked. Table 1 shows the total seismicity catalogued up to sol 1233.


Table 1

Event Type


















2.4 Hz

























Noteworthy New Events: The V10 catalogue [7] and data release [8] contain some significant new events including 6 new Quality A events. Most notable is S1094b, BB, QA. This event, with MwMa=4.0 is one of the largest BB events ever detected and originates in a region where no previous seismicity had been seen with a distance of ~60° and a backazimuth of 40°.

We will also introduce the largest event ever recorded on Mars, S1222a.

Summary: The InSight lander has provided an incredible seismic dataset and the MQS has meticulously catalogued the seismicity. Now well into its extended mission the seismicity from the second Martian year of operations has been exceptional and MQS will continue to add to the catalogue.


References: [1] Lognonné et al., 2020, Nature Geo., 13(3), 213-220. DOI: 10.1038/s41561-020-0536-y, [2] Clinton et al., 2021, Physics of the Earth and Planetary Interiors, 310, 106595, DOI: 10.1016/j.pepi.2020.106595, [3] Ceylan et al., 2021, Physics of the Earth and Planetary Interiors, 310, 106597, DOI: 10.1016/j.pepi.2020.106597, [4] Hobiger et al., 2021, Nature Comms, 12(1), 1-13, DOI: 10.1038/s41467-021-26957-7, [5] Giardini et al., 2020, Nature Geo., 13(3) 205-212, DOI: 10.1038/s41561-020-0539-8, [6] Dahmen et al., 2020, J. Geophys. Res., 126, p.e.2020JE006599, DOI: 10.1029/2020JE006599, [7] InSight Marsquake Service (2022), Mars Seismic Catalogue, InSight Mission; V10 2022-04-01, ETHZ, IPGP, JPL, ICL, Univ. Bristol,, [8] InSight Mars SEIS Data Service, 2022, SEIS raw data, InSight Mission, IPGP, JPL, CNES, ETHZ, ICL, MPS, ISAE-Supaero, LPG, MFSC, DOI: 10.18715/SEIS.INSIGHT.XB_2016

How to cite: Horleston, A., Clinton, J., Ceylan, S., Kawamura, T., Stähler, S. C., Charalambous, C., Dahmen, N. L., Duran, C., Kim, D., Plasman, M., Zenhäusern, G., Euchner, F., Knapmeyer, M., Giardini, D., Lognonné, P., Pike, W. T., Panning, M., Smrekar, S., and Banerdt, W. B.: The Seismicity of Mars as Recorded by InSight's Marsquake Service, Europlanet Science Congress 2022, Granada, Spain, 18–23 Sep 2022, EPSC2022-880,, 2022.

Alexander Stott, Raphael Garcia, Armand Chédozeau, Aymeric Spiga, Naomi Murdoch, Baptiste Pinot, David Mimoun, Constantinos Charalambous, Anna Horleston, Scott King, Taichi Kawamura, Nikolaj Dahmen, Salma Barkaoui, Philippe Lognonné, and William Banerdt


The seismic data collected by the SEIS (seismic experiment for the interior structure of Mars) [1] package on the NASA InSight mission provide a way to obtain the structure and properties of Mars [2,3]. As the seismometer is installed on the surface, it is also sensitive to atmospherically generated signals from the pressure and wind [4]. These can constrain the ability to observe seismic events. To that end, in Stott et al. 2022 [5], we recently proposed machine learning models to predict the seismic energy from the pressure, wind and temperature data recorded by the dedicated sensors in the APSS (auxiliary payload sensor suite) and TWINS (temperature and wind for InSight) sensor packages. This yields a tool to analyse and discover marsquakes.

A machine learning model for predicting seismic noise

The inputs used for the machine learning model are:

  • The pressure and the envelope of the pressure
  • The wind speed and direction from two sensor booms
  • The air temperature from each sensor boom

These inputs are fed into the machine learning models to predict the seismic velocity energy for the vertical (Z) and horizontal (North and East) axes in the low frequency (0.4-1Hz) and high frequency (2.2-2.6Hz) regions, as these cover the predominant excitations in both categories of marsquakes in the catalogue [6].

Two methods were proposed, the first using gaussian process regression [7] and the second using a neural network [8]. The neural network was found to have superior prediction performance. However, the gaussian process regression model could be used to better dig into the physical relationships between wind, pressure and the ground motion owing to the form and the probabilistic nature of the method.

Marsquake analysis

Figure 1 shows the prediction and observed data for a period containing a marsquake. Before and after the marsquake the observed and predicted data match well, but during the marsquake they diverge indicating the event not to be generated by atmospheric effects. This event is classified as a 2.4Hz event. This class of event typically only exhibits excitation around the 2.4Hz resonance [6,9,10]. However, the prediction in the 0.4-1Hz bandwidth indicates there is some separate energy. Examining the machine learning prediction across the marsquake catalogue shows low frequency energy occurs consistently in the higher frequency classes of events for greater event magnitudes.

Detecting new events and seismic detectivity

The proposed machine learning models identify marsquakes as separate from the background atmospherically generated noise. This can be used to discover new marsquakes as outliers. An algorithm based on this concept was proposed and is demonstrated in Figure 2. Using this approach, a further 32 good candidate events are found along with another 34 suggestions.


Machine learning is a nascent technology in planetary science, whereby it provides a powerful tool to fit datasets using mathematical complexity. Here, we present a model to fit the atmospheric data to the seismic data from the NASA InSight mission. In doing so we

  • Provide insights on the relationship between ground motion and the wind, pressure and temperature
  • Provide a tool to analyse marsquakes
  • Enable the detection of new marsquakes

In this way the machine learning models offer new information about marsquakes and also enable the study of their detectivity to help place constraints on the seismicity of Mars.


[1] Lognonné et al. SSR 2019

[2] Giardini et al. Nat. Geosci 2020

[3] Stähler et al. Science 2022

[4] Charalambous et al. JGR: Planets 2021

[5] Stott et al. Submitted 2022

[6] Clinton et al. PEPI 2021

[7] Williams and Rasmussen, MIT Press, 2006

[8] Goodfellow et al. MIT Press 2016

[9] Compaire et al. JGR: Planets 2021

[10] Dahmen et al. BSSA 2021

Figure 1: The S0534a 2.4Hz event with the machine learning model prediction. The top panel shows the observed seismic energy compared to the prediction where the energy during the event is seen to diverge. The second and third panels show the SNR (signal to noise ratio) calculated between the observed and predicted signal for the HF (2.2-2.6Hz) and LF (0.4-1Hz) bandwidths respectively. The bottom panel shows the velocity spectrogram during the event.

Figure 2: Figure showing the detection of a new event with the algorithm in [5]. The top panel shows the prediction from the machine learning model and the observed data. The middle panel shows the detrended residual and the smoothed MAE (mean absolute error) highlighting the outlying points indicating an event. The bottom panel shows the histogram (over a 2 hour window) of the MAE where the points belonging to the event are seen to be outliers in the distribution.

How to cite: Stott, A., Garcia, R., Chédozeau, A., Spiga, A., Murdoch, N., Pinot, B., Mimoun, D., Charalambous, C., Horleston, A., King, S., Kawamura, T., Dahmen, N., Barkaoui, S., Lognonné, P., and Banerdt, W.: Using machine learning to separate atmospherically generated noise from marsquakes, Europlanet Science Congress 2022, Granada, Spain, 18–23 Sep 2022, EPSC2022-833,, 2022.

Raphael F. Garcia and the InSight and CTX science teams

The InSight (Interior Exploration using Seismic Investigations, Geodesy and Heat Transport) lander is providing an unprecedented set of high frequency records of ground deformations, pressure and wind at the surface of Mars. Seismic and acoustic waves from impacts have been observed by SEIS seismometer through the ground movements that they generate, and their sources as impacts have been confirmed with orbital images of newly appearing craters. Three impacts have been identified and located, and one additional seismic event is consistent with a date-constrained new impact. Several other events are currently being investigated as well. Arrival times and polarization of seismic and acoustic waves were used to estimate impact locations. These estimated locations were subsequently confirmed by orbital imaging of associated craters and temporal matches to the events times, with previous orbital images showing no craters. Crater dimensions and estimates of meteoroid trajectories from images allow us to understand and model the recorded seismograms. The precise source locations provided by impacts as compared to tectonic sources provide direct constraints on the structure of the martian interior. First arrival seismic waves confirm the previously determined crustal models, and the dispersion of trapped acoustic waves confirms the current models of sound speed and wind in the atmosphere. In addition, these observations provide the first ground-truth for distance-amplitude scaling relationships between impacts and the mechanical waves for Mars. They confirm the relationship between the seismic moment of impacts and the vertical impactor momentum, and demonstrate the capability of planetary seismology to constrain impact rates and internal structure of terrestrial planetary objects throughout the solar system.

How to cite: Garcia, R. F. and the InSight and CTX science teams: First seismo-acoustic location and orbital imaging of recent impacts on Mars: constraints on atmosphere and interior structure and impact processes, Europlanet Science Congress 2022, Granada, Spain, 18–23 Sep 2022, EPSC2022-10,, 2022.

Cecilia Durán, Amir Khan, Savas Ceylan, and Domenico Giardini

The InSight (Interior Exploration using Seismic Investigations, Geodesy and Heat Transport) mission has been collecting high-quality seismic data on Mars for over 1200 Martian solar days (sols), providing the first direct observations of its interior structure. In this study, we report on a complete analysis of low-frequency (LF) seismic evens (main energy below 1 Hz), that have proven essential for understanding the deep interior of Mars [1], [3]-[5]. As Martian seismograms are generally characterized by a low signal-to-noise ratio and strong scattering, data selection and processing is central for the identification and consistent picking of body-wave arrivals. To this end, we employ complementary approaches, including: 1) narrow-band-filtered time-domain envelopes; 2) polarised waveforms and their time-domain envelopes; 3) polarisation analysis; and 4) waveform matching. An example of the application of the first two approaches to identify body-wave arrivals is shown in Fig 1. The careful application of this 4-step processing scheme to each LF event has resulted in a significant increase in body wave phase picks, encompassing direct, reflected, and converted body waves, relative to earlier analyses [1]. With our method, all phases are obtained through iterative refinement, which relies on inverting an initial set of phase picks for models of interior structure, event distance, and depth, while predicting travel times for additional seismic phases. Based on the predictions, we re-pick (every pick is subject to our 4-step processing scheme), thereby enlarging our differential travel-time data set (all picks are relative to the main P- or PP-wave arrival), and subsequently re-invert for an updated set of interior structure models, event distances, and depths. Through application of iterative refinement to new seismic data, we continuously update information on radial velocity models of the crust, mantle, and core, as well as source parameters, and phase picks. The recent identification of the most distant marsquakes recorded by the InSight mission and the first core-diffracted P-wave [2] allows constraints on the P-wave velocity structure of the lower mantle (see distant event in Fig 2) that hitherto could not be constrained because of lack of lower-mantle-traversing P-waves. In addition, on sol 1222, InSight detected the largest event ever recorded on Mars, providing a unique opportunity to understand the interior of Mars through the identification of waves that could not be identified in earlier events because of their low signal-to-noise ratio. Our iterative refinement procedure demonstrates that through continued analysis of new events, our knowledge and understanding of the interior of Mars will steadily improve.

Figure 1: Vertical-component time-domain envelopes and polarised traces for event S0235b across different frequency bands (filter banks). Grey envelopes depict unpolarised time-domain envelopes. Polarised traces and their envelopes are shown in black (lines) and colored envelopes, respectively. Initial seismic phase picks for S0235b and their uncertainties are shown as vertical dotted lines and orange bars, respectively. . The vertical grey bar in panel (A) represents a glitch. Panel (B) shows a zoomed-in version of the P-wave phases in (A). The vertical dashed line in (B) denoted by X, indicates another possible seismic wave arrival. Amplitudes of traces and envelopes are normalized by their maxima and scaled for better visualisation.


Figure 2: Body-wave (P, S, PP, SS, PPP, SSS, ScS and Pdiff) ray path geometry for 15 low-frequency events. Events are indicated by green stars and InSight is represented by a black triangle. The colorbar indicates ray path density in the form of the number of rays that pass through a certain area based on inverted models [1], which accounts for the diffuseness of ray paths. The radial sensitivity computed as the integrated ray path density with epicentral distance is shown at the bottom.



[1] C. Durán, A. Khan, S. Ceylan, et al., “Seismology on Mars: An analysis of direct, reflected, and converted seismic body waves with implications for interior structure,” Physics of the Earth and Planetary Interiors, vol. 325, p. 106 851, 2022, issn: 0031-9201. doi: 10.1016/j.pepi.2022.106851.

[2] A. C. Horleston, J. F. Clinton, S. Ceylan, et al., “The Far Side of Mars: Two Distant Marsquakes Detected by InSight,” The Seismic Record, vol. 2, no. 2, pp. 88–99, Apr. 2022, issn: 2694-4006. doi: 10.1785/0320220007.

[3] A. Khan, S. Ceylan, M. van Driel, et al., “Upper mantle structure of Mars from InSight seismic data,” Science, vol. 373, no. 6553, pp. 434–438, 2021, issn: 0036-8075. doi: 10.1126/science.abf2966.

[4] A. Khan, P. Sossi, C. Liebske, A. Rivoldini, and D. Giardini, “Geophysical and cosmochemical evidence for a volatile-rich Mars,” Earth and Planetary Science Letters, vol. 578, p. 117 330, 2022, issn: 0012-821X. doi: 10.1016/ j.epsl.2021.117330.

[5] S. C. Stähler, A. Khan, W. B. Banerdt, et al., “Seismic detection of the martian core,” Science, vol. 373, no. 6553, pp. 443–448, 2021, issn: 0036-8075. doi: 10.1126/science.abi7730.

How to cite: Durán, C., Khan, A., Ceylan, S., and Giardini, D.: Seismology on Mars: Analysis and identification of body-waves with implications for Martian structure, Europlanet Science Congress 2022, Granada, Spain, 18–23 Sep 2022, EPSC2022-969,, 2022.

Henri Samuel, Mélanie Drilleau, Raphael Garcia, Quancheng Huang, Attilio Rivoldini, Philippe Lognonné, and Bruce Banerdt

Introduction: The identification of deep reflected phases in the InSight mission seismic recording as core reflected phases have led to the first seismic detection of the Martian core [1]. These results indicated that the core size of Mars spans the higher end of InSight pre-mission estimates, implying a large fraction of Sulfur in the core together with smaller fractions of O, C, and H.

The presence of well-separated silicate and metallic envelopes indicates that Mars experienced an early global magma ocean stage whose crystallisation likely led to the formation of a compositionally distinct layer at the bottom of the mantle [2]. Such a layer is expected to be heavily enriched in heat-producing elements and in iron, leading to long-term stability with little mixing between the layer and the overlying mantle. The presence of this enriched basal layer often yields the development of a molten silicate layer above the core that could act as a deep seismic reflector [3].

Objectives and approach:  To test the compatibility of Martian mantle layering with seismic data, we conducted Monte Carlo Markov chain inversions in which the long-term thermo-chemical history of Mars’ main envelopes is embedded into the forward problem (Fig. 1). Contrary to more classical approaches that directly invert for seismic velocities and density along a radial planetary profile, we explore a different and smaller model space that consists in parameters that govern the thermo-chemical evolution of Mars: the planetary initial thermal state, the core size, the mantle rheology, or the crustal enrichment in heat-producing element with respect to the mantle. For each model we evolve the system for 4.5 Gyr.


Figure 1: Schematic representation of a classical [5] (left) and a geodynamic parameterization [4] (used here) to infer Mars’ seismic structure via Monte Carlo Markov chain inversion of seismic data.

This results in thermal structures at present-day that are converted into seismic structures [4]. This approach allows for a more coherent and more constrained profile than in classical inversions (Fig. 1), due to the more informative prior considered, and allows reconstructing the long-term history of the planet. Our inversion approach also considers an enriched silicate layer above the core-mantle boundary and we therefore invert for the layer thickness and for its thermal conductivity. We used the most recent travel time dataset that contains considerably more shallow and deep phases compared to previous studies.

Altogether, this allows us to test the hypothesis of the presence of a molten layer at the top of Mars' core, along with the associated consequences on the interpretation of seismic, geodetic, and geochemical data.


Results: Our inversion results indicate that the presence of a basal mantle layer is at least as compatible with seismic data than a homogeneous Martian structure (Fig. 2).  

Figure 2: Data fit for the inversion output considering a basal mantle layer (left) or a homogeneous mantle (right). Warm and cold colors indicate respectively high and small probability densities.


Figure 3: Inversion results without (top) and with (bottom) a mantle stratification. Prior (grey) and posterior (blue) distribution for several inverted parameters. The plain and dashed red lines indicate the posterior mean value and 1-sigma range.

Fig. 3 shows that the basal mantle layer yields a thinner, crust less enriched in heat-producing element than in the case of a homogeneous mantle. It also shows that a basal mantle layer is also compatible with geodetic constraints on k2 values: 0.174±0.008 [7], in good agreement with theoretical predictions [3].


Conclusions: Our new inversion results indicate that the presence of a basal mantle layer is compatible with seismic data. The basal mantle layer stores a significant fraction of heat-producing elements and depletes the rest of the mantle. This reduces shallow mantle heating and melting leading to thinner crusts that are more compatible with receiver functions and gravity data constraints [8]. Finally, the presence of the enriched basal mantle layer is also compatible with the Love number values for Mars [7]. We will present more advanced inversion results using an augmented database that will contain additional phases, including phases that are most sensitive to the deep mantle structure (e.g., Pdiff).


[1] Stähler, S., et al., Science 373, 443–448 (2021)

[2] Elkins-Tanton, L. et al., JGR., doi:10.1029/2005JE002480 (2003)

[3] Samuel, H. et al., JGR, doi:10.1029/2020JE006613 (2021).

[4] Drilleau, M. et al., G. J. Int., 226, 1615-1644 (2021)

[5] Drilleau M., et al., G. J. Int., 195, 1165-1183  (2013)

[6] Khan, A., et al., Science, 373, 434-438 (2021)

[7] Konopliv, A. et al., JGR, 47, doi:  10.1029/2020GL090568 (2020)

[8] Knapmeyer-Endrun, B. et al., Science, 373 (2021)



How to cite: Samuel, H., Drilleau, M., Garcia, R., Huang, Q., Rivoldini, A., Lognonné, P., and Banerdt, B.: Testing the Presence Deep Martian Mantle Layering in the light of InSight Seismic Data, Europlanet Science Congress 2022, Granada, Spain, 18–23 Sep 2022, EPSC2022-297,, 2022.

Orals: Wed, 21 Sep | Room Machado

Chairpersons: Melanie Drilleau, Foivos Karakostas, Simon C. Stähler
Martin Knapmeyer, Savas Ceylan, Constantinos Charalambous, John Clinton, Nikolaj Dahmen, Cecilia Duran, Anna Horleston, Taichi Kawamura, Doyeon Kim, Matthieu Plasman, Simon Stähler, Géraldine Zenhäusern, Renee C. Weber, Domenico Giardini, Mark Panning, Philippe Lognonné, and Bruce Banerdt

On the 1222nd sol of the InSight mission (or May the Fourth of 2022), a seismic event was detected that turned out to be the largest marsquake recorded so far. At a moment magnitude of 4.7, event S1222a released as much seismic moment as all seismic events previously catalogued by InSight together, and greatly surpasses the second largest event, S0976a.

The frequency of occurrence of earthquakes follows a power law with an exponent (the b value) close to 1 over a wide range of magnitudes. The class of Low Frequency marsquakes, to which S1222a belongs, shows a similar behaviour. At low magnitudes, the slope of the cumulative distribution suggests that the InSight marsquake catalog for Low Frequency events is representative for events with moment magnitude exceeding 3 (Figure 1). Up to and including the occurrence of event S0976a one could however guess that events larger than magnitude approx. 3.6 might be less frequent than predicted by this power law. Event S1222a mends this apparent decrease, and the b value derived from all events is 1 within the formal uncertainty. Hence the resulting distribution follows a power law for events exceeding magnitude 3.

Figure 1 Size frequency distribution of Low Frequency and Broadband events. Red: status before occurrence of S1222a, blue: afterwards. The dashed lines represent the Maximum Likelihood solutions for fitting a Power Law distribution, with b-values as indicated in the legend.

Since S0976a occurred at an epicentral distance of about 140° (Horleston et al., 2022) at night time(02:26 LMST), but was recorded with amplitudes high enough to be visible during typical day time noise, we can infer that the catalog is complete for events the size of S0976a or larger. Since the slope of the distribution appears to be stable over the entire magnitude range from 3 up to that of S1222a, we conclude that the catalogue of Low Frequency events is not only representative, but complete for events larger than 3.

The magnitude of S1222a is slightly to the right of the value predicted by the power law, but, compared to other deviations, not excessively so. Also, the distribution does not show an indication of increasing slope any more. With about 30 events of magnitude 3 during the 3.3 years of InSight observations so far, this suggests that events larger than S1222a can be expected with respective recurrence rates: A magnitude 5 event appears likely about once in 11 years, i.e. on a decadal scale.

We use the approach of Knapmeyer et al. (2018) to estimate the seismic moment rate of Mars, and also the corner magnitude above which the size-frequency distribution becomes much steeper and events larger than the corner magnitude extremely unlikely (our analysis uses the Tapered Gutenberg-Richter distribution with an exponential roll-off above the corner frequency). The simple idea of Knapmeyer et al. (2018) was that, since a few large events release most of the seismic moment, the moment of the largest event ever observed, scaled with the duration covered by the catalog and a factor depending on b value and corner moment, provides a good estimation of the moment rate. Including not only the one largest, but the n largest events (with n between 1 and 10 or so) provides an even better estimation. In the 2018 study we have shown that the approach yields a surprisingly good estimation of the Earth's moment rate after a few months of registration, rather than after the many decades necessary to witness a magnitude 9 event.

We employ here an estimation based on the 5 largest events in the catalog. Knapmeyer et al. (2018) have shown that this low number may already provide a reasonable rate estimate. At the same time we attempt to avoid using events too small to be globally detectable. Future analyses will show if the use of 5 events is actually the best choice. The inclusion of event S1222a shifts the estimated moment rate slightly upwards, as well as the estimated corner moment. Figure 2 shows that both parameters are slightly below those of the WeakMany model of Knapmeyer et al., (2006). Compared to the catalog from before its occurrence, the parameter modifications due to S1222a are hardly significant (and not shown here). At least as important is the reduction of the uncertainty region in the moment-rate / corner-moment parameter space. We cannot yet exclude that the moment rate of Mars is as low as that of the Moon, but if it were, it would be rather unlikely to observe the event sequence that InSight observed.

Figure 2 Estimation of Moment Rate and Corner Moment. For each Pixel, 10000 synthetic catalogs were evaluated to compare their parameters with the five largest events from the InSight catalog (KS5 estimation). Colour indicates the probability that the five largest events drawn from a tapered Gutenberg Richter distribution in 3.3 years reproduce the KS5 estimation obtained from the observed events, to within an 80 %  interdecile around the median. The vertical line indicates the moment rate for the Moon (observed, from Shallow Moonquakes), markers indicate the WeakMany and Medium models of Knapmeyer et al. (2006). The moment rate of the Earth corresponds to an equivalent magnitude of about 8.5 (Knapmeyer et al., 2018). Horizontal lines indicate the magnitude of the largest event and its uncertainty. The maximum of the distribution is marked by a small white circle.



Horleston, A. C., et al. (2022). The Far Side of Mars: Two Distant Marsquakes Detected by InSight, The Seismic Record. 2(2), 88–99, doi: 10.1785/0320220007

Knapmeyer, M., Oberst, J., Hauber, E., Wählisch, M., Deuchler, C., Wagner, R. (2006). Working models for spatial distribution and level of Mars' seismicity. Journal of Geophysical Research, vol. 111, E11006, doi:10.1029/2006JE002708

Knapmeyer, M., et al. (2018). Estimation of the Seismic Moment Rate from an Incomplete Seismicity Catalog, in the Context of the InSight Mission to Mars, Bull. Seis. Soc. Am., vol. 109, No. 3, 1125-1147, doi: 10.1785/0120180258

Taylor et al., (2013). Estimates of seismic activity on the Cerberus Fossae region of Mars, Journal of Geophysical Research, vol. 118, 2570-2581, doi:10.1002/2013JE004469

How to cite: Knapmeyer, M., Ceylan, S., Charalambous, C., Clinton, J., Dahmen, N., Duran, C., Horleston, A., Kawamura, T., Kim, D., Plasman, M., Stähler, S., Zenhäusern, G., Weber, R. C., Giardini, D., Panning, M., Lognonné, P., and Banerdt, B.: Seismic Moment Rate of Mars after Event S1222a, Europlanet Science Congress 2022, Granada, Spain, 18–23 Sep 2022, EPSC2022-845,, 2022.

Spectral analyses of the Largest Marsquake S1222a: Implications for the attenuation, the source parameters and the site effect
Taichi Kawamura, Sabrina Menina, Ludovic Margerin, Zongbo Xu, Wanbo Xiao, and Philippe Lognonné
Sabrina menina, Ludovic Margerin, Taïchi kawamura, Philippe Lognonné, Zongbo Xu, Keisuke Onodera, Marie Calvet, and Raphaël Garcia

Since its deployment at the surface of Mars in 2019, the SEIS seismometer (Seismic Experiment for Interior Structure) of the NASA-InSight mission (Interior exploration using Seismic Investigations, Geodesy and Heat Transport) has detected hundreds of small magnitude seismic events. But recently, the seismic event S1222a of magnitude 4.9 which primarily contains energy above 0.1 Hz was reported by the Mars Quake Service. It is characterized by two temporally separated arrivals with a gradual beginning, a broad maximum and a very long decay lasting up to 1 hour. In this work, we propose to use seismic coda analysis on this specific event to explore the attenuation properties of Martian crust and to provide an Earth-Mars-Moon comparison.


Fig 1: Envelope characteristics of the S1222a Broadband event. The data were collected on the VBB sensors with sampling rates of 20 samples per second. (a) Base 10 logarithm of energy envelopes normalized by the noise level as a function of time and fit of the coda decay (black dashed lines). The gray dashed line indicates signal-to-noise ratio (SNR) = 2. The different colors correspond to different frequency bands indicated in inset. The time t = 0 s corresponds to the first arrival. tP and tS refer to the arrival times of the first and second wave packets (red dashed lines). (b) Frequency dependence of Qc in the 0.4 –7.5 Hz frequency band. The color bar indicates the SNR measured on the horizontal components. The red and black dashed line shows the power-law fit Qc = 900f0.36and Qc = 655f0.94 respectively (f is the frequency.).

We employ basic multiple-scattering concepts to examine and to compare the energetic characteristics of the Earth, Martian and Lunar crust. We measured the coda quality factor Qc, that characterizes the rate of decay of the seismogram envelope, to have a first view of the attenuation characteristics [1] (see fig 1.a). We performed a fit of the frequency dependence of Qc using a parameterization of the form Qc = Q0 fn. Based on linear least-squares applied to the logarithm of Qc (see fig.1.b), we obtain two very distinct regimes of Qc as a function of the frequency f. In the [0.4 Hz - 1.5 Hz] range Q0 = 900 and n  0.36, when at higher frequency, Q0 = 655 and n  0.94 This result is specific to Mars and could be related to the stratification of the attenuation properties in the Martian crust. The quasi-linear frequency increase of Qc seen in [2.5 Hz – 7.5 Hz] range suggests that the coda decay is governed by a frequency independent attenuation time. Similar values of Q0 and n, in this frequency range, were reported by using a Very High Frequency group events [2] (see Fig 2 in orange and red).  


Fig 2:  Coda attenuation as a function of the frequency for the Earth-cratons regions (green), the Earth-tectonics regions (bleu), the Moon (gray), Mars ( α= 0 in red and α= 1 in orange) and the S1222a BB event (black).

    In the case of the regional earthquakes on Earth, the frequency dependence of Qc is radically different between tectonic areas and shields. Q0varies from 100 to 400 in tectonic areas (blue) [e.g. 3] and from 600 to 1000 in shields areas (green) [e.g. 4], where the frequency exponent n is equal to 1 and 0.1-0.5 respectively. [5] report coda Qc values of 2400 at 0.5 Hz and possibly as large as 10000 at 7 Hz on the Moon (gray). This result suggests that Q0 is a factor of 2 to 4 larger than on Mars when the frequency exponent n is more likely in the range of 0.4 – 0.6. To summarize, our Qc measurements exhibit two different regimes of frequency dependence. At higher frequency, the Q0 and n valuesare very close to the one found on Very High Frequency group events. The unique combination of high Q0 with a frequency exponent n =1 also distinguishes very clearly Mars from the Earth. By the time of the meeting, we hope to present the first results of envelope modeling based on elastic radiative transport in a stratified medium to quantify the contribution of the scattering and the absorption to the total attenuation.




 [1] Aki, K. (1969). Analysis of the seismic coda of local earthquakes as scattered waves. Journal of geophysical research74(2), 615-631.

 [2] Menina et al. (2021), Energy Envelope and Attenuation Characteristics of High‐Frequency (HF) and Very‐High‐Frequency (VF) Martian Events, Bulletin of the Seismological Society of America (2021) 111 (6): 3016–3034.

 [3] Mayeda, K., Koyanagi, S., Hoshiba, M., Aki, K., & Zeng, Y. (1992). A comparative study of scattering, intrinsic, and coda Q− 1 for Hawaii, Long Valley, and central California between 1.5 and 15.0 Hz. Journal of Geophysical Research: Solid Earth97(B5), 6643-6659.

[4] Mitchell, B, J et al. (2008). A continent-wide map of 1-Hz Lg coda Q variation across Eurasia and its relation to lithospheric evolution. Journal of Geophysical Research: Solid Earth, 113(B4).

[5] Gillet, K., et al (2017). Scattering attenuation profile of the Moon: Implications for shallow moonquakes and the structure of the megaregolith. Physics of the Earth and Planetary Interiors262, 28-40.


How to cite: menina, S., Margerin, L., kawamura, T., Lognonné, P., Xu, Z., Onodera, K., Calvet, M., and Garcia, R.: Scattering and attenuation properties of the Martian crust through the S1222a seismic event, Europlanet Science Congress 2022, Granada, Spain, 18–23 Sep 2022, EPSC2022-805,, 2022.

Ebru Bozdag, Ana-Catalina Plesa, Quancheng Huang, Daniel Peter, Caio Ciardelli, Nienke Brinkman, Martin Knapmeyer, Savas Ceylan, Mark Wieczorek, Brigitte Knapmeyer-Endrun, Attilio Rivoldini, Amir Khan, Benjamin Fernando, Tarje Nissen-Meyer, and Philippe Lognonné

The InSight (Interior exploration using Seismic Investigations, Geodesy and Heat Transport) lander provides the unique opportunity to explore the seismicity and the present-day seismic structure of Mars (Banerdt et al., 2020; Lognone et al., 2020; Giardini et al. 2020). Since the deployment of the seismometer SEIS, hundreds of seismic events have been recorded, and information about these events is available in the seismic catalogs produced by the Mars Quake Service (MQS) (Clinton et al., 2021; MQS Catalog V9).