Measuring Fundamental and Higher Mode Surface Wave Dispersion on Mars From Seismic Waveforms

One of the goals of the Interior Exploration using Seismic Investigations, Geodesy and Heat Transport (InSight) mission is to constrain the interior structure of Mars. We present a hierarchical transdimensional Bayesian approach to extract phase velocity dispersion and interior shear‐wave velocity (VS) models from a single seismogram. This method was adapted to Mars from a technique recently developed for Earth (Xu & Beghein, 2019, https://doi.org/10.1093/gji/ggz133). Monte Carlo Markov Chains seek an ensemble of one dimensional (1‐D) VS models between a source and a receiver that can explain the observed waveform. The models obtained are used to calculate the phase velocities of fundamental and higher modes at selected periods, and a subsequent analysis is performed to assess which modes were reliably measured. An advantage of our approach is that it can also fit unknown data noise, which reduces the risk of overfitting the data. In addition, uncertainties in the source parameters can be propagated, yielding more accurate model parameter uncertainties. In this study, we first present our technique and discuss the challenges stemming from using a single station to characterize both structure and the source and from the absence of a Mars reference model. We then demonstrate the method feasibility using the Mars Structure Service blind test data and our own synthetic data, which included realistic noise levels based on the noise recorded by InSight.

[1]  W. Banerdt,et al.  MSS/1: Single‐Station and Single‐Event Marsquake Inversion , 2020, Earth and Space Science.

[2]  Jeroen Tromp,et al.  Initial results from the InSight mission on Mars , 2020, Nature Geoscience.

[3]  C. Russell,et al.  Constraints on the shallow elastic and anelastic structure of Mars from InSight seismic data , 2020, Nature Geoscience.

[4]  S. Kedar,et al.  The seismicity of Mars , 2020, Nature Geoscience.

[5]  Jiaxuan Li,et al.  Preparing for InSight: Evaluation of the Blind Test for Martian Seismicity , 2019, Seismological Research Letters.

[6]  C. Beghein,et al.  Measuring higher mode surface wave dispersion using a transdimensional Bayesian approach , 2019, Geophysical Journal International.

[7]  Huafeng Liu,et al.  SEIS: Insight’s Seismic Experiment for Internal Structure of Mars , 2019, Space Science Reviews.

[8]  M. Golombek,et al.  Pre-mission InSights on the Interior of Mars , 2019, Space Science Reviews.

[9]  Jeroen Tromp,et al.  Planned Products of the Mars Structure Service for the InSight Mission to Mars , 2017 .

[10]  M. van Driel,et al.  A probabilistic framework for single-station location of seismicity on Earth and Mars , 2017 .

[11]  Simon Stähler,et al.  Fully probabilistic seismic source inversion – Part 2: Modelling errors and station covariances , 2016 .

[12]  Christina Kluge,et al.  Data Reduction And Error Analysis For The Physical Sciences , 2016 .

[13]  Thorne Lay,et al.  Seismological implications of a lithospheric low seismic velocity zone in Mars , 2015 .

[14]  W. Banerdt,et al.  Verifying single-station seismic approaches using Earth-based data: Preparation for data return from the InSight mission to Mars , 2015 .

[15]  F. Nimmo,et al.  Dissipation at tidal and seismic frequencies in a melt‐free, anhydrous Mars , 2013 .

[16]  Göran Ekström,et al.  The global CMT project 2004–2010: Centroid-moment tensors for 13,017 earthquakes , 2012 .

[17]  Malcolm Sambridge,et al.  Transdimensional inversion of receiver functions and surface wave dispersion , 2012 .

[18]  G. Ekström,et al.  Automated multimode phase speed measurements for high-resolution regional-scale tomography: application to North America , 2010 .

[19]  C. Beghein.  Radial anisotropy and prior petrological constraints: A comparative study , 2010 .

[20]  A. Valentine,et al.  Reducing errors in seismic tomography: combined inversion for sources and structure , 2010 .

[21]  M. Sambridge,et al.  Seismic tomography with the reversible jump algorithm , 2009 .

[22]  A. Curtis,et al.  Global variations of temperature and water content in the mantle transition zone from higher mode surface waves , 2009 .

[23]  K. Visser.  Monte Carlo search techniques applied to the measurement of higher mode phase velocities and anisotropic surface wave tomography. Geologica Ultraiectina (285) , 2008 .

[24]  S. Lebedev,et al.  Global Love wave overtone measurements , 2007 .

[25]  B. Romanowicz,et al.  A Three-Dimensional Radially-Anisotropic Model of Shear Velocity in the Whole Mantle , 2006 .

[26]  J. Trampert,et al.  Radial anisotropy in seismic reference models of the mantle , 2006 .

[27]  Daniel E. McNamara,et al.  Ambient Noise Levels in the Continental United States , 2004 .

[28]  T. Gudkova,et al.  Mars: interior structure and excitation of free oscillations , 2004 .

[29]  É. Beucler,et al.  Surface wave higher-mode phase velocity measurements using a roller-coaster-type algorithm , 2003 .

[30]  G. Nolet,et al.  Upper mantle beneath Southeast Asia from S velocity tomography , 2003 .

[31]  A. Malinverno.  Parsimonious Bayesian Markov chain Monte Carlo inversion in a nonlinear geophysical problem , 2002 .

[32]  B. Kennett,et al.  Non-linear waveform inversion for surface waves with a neighbourhood algorithm—application to multimode dispersion measurements , 2002 .

[33]  J. Woodhouse,et al.  GLOBAL HIGH-RESOLUTION PHASE VELOCITY DISTRIBUTIONS OF OVERTONE AND FUNDAMENTAL-MODE SURFACE WAVES DETERMINED BY MODE BRANCH STRIPPING , 1999 .

[34]  Hendrik Jan van Heijst,et al.  Measuring surface-wave overtone phase velocities using a mode-branch stripping technique , 1997 .

[35]  Tilman Spohn,et al.  The interior structure of Mars: Implications from SNC meteorites , 1997 .

[36]  Barbara Romanowicz,et al.  Global mantle shear velocity model developed using nonlinear asymptotic coupling theory , 1996 .

[37]  Barbara Romanowicz,et al.  Comparison of global waveform inversions with and without considering cross-branch modal coupling , 1995 .

[38]  G. Nolet.  Partitioned waveform inversion and two‐dimensional structure under the network of autonomously recording seismographs , 1990 .

[39]  J. Montagner,et al.  Vectorial tomography—II. Application to the Indian Ocean , 1988 .

[40]  M. Cara,et al.  Waveform inversion using secondary observables , 1987 .

[41]  H. Nataf,et al.  A simple method for inverting the azimuthal anisotropy of surface waves , 1986 .

[42]  John H. Woodhouse,et al.  Determination of earthquake source parameters from waveform data for studies of global and regional seismicity , 1981 .

[43]  J. Woodhouse.  The coupling and attenuation of nearly resonant multiplets in the Earth's free oscillation spectrum , 1980 .

[44]  M. Cara.  Regional variations of higher Rayleigh-mode phase velocities: a spatial-filtering method , 1978 .

[45]  G. Nolet.  Higher Rayleigh modes in western Europe , 1975 .

[46]  Naohiro Soga,et al.  Some elastic constant data on minerals relevant to geophysics , 1968 .

[47]  F. Dahlen.  The Normal Modes of a Rotating, Elliptical Earth , 1968 .