Influence of Body Waves,Instrumentation Resonances,and Prior Assumptions on Rayleigh Wave Ellipticity Inversion for Shallow Structure at the InSight Landing Site

Based on an updated model of the regolith’s elastic properties, we simulate the ambient vibrations background wavefield recorded by InSight’s Seismic Experiment for Interior Structure (SEIS) on Mars to characterise the influence of the regolith and invert SEIS data for shallow subsurface structure. By approximately scaling the synthetics based on seismic signals of a terrestrial dust devil, we find that the high-frequency atmospheric background wavefield should be above the self-noise of SEIS’s SP sensors, even if the signals are not produced within 100–200 m of the station. We compare horizontal-to-vertical spectral ratios and Rayleigh wave ellipticity curves for a surface-wave based simulation on the one hand with synthetics explicitly considering body waves on the other hand and do not find any striking differences. Inverting the data, we find that the results are insensitive to assumptions on density. By contrast, assumptions on the velocity range in the upper-most layer have a strong influence on the results also at larger depth. Wrong assumptions can lead to results far from the true model in this case. Additional information on the general shape of the curve, i.e. single or dual peak, could help to mitigate this effect, even if it cannot directly be included into the inversion. We find that the ellipticity curves can provide stronger constraints on the minimum thickness and velocity of the second layer of the model than on the maximum values. We also consider the effect of instrumentation resonances caused by the lander flexible modes, solar panels, and the SEIS levelling system. Both the levelling system resonances and the lander flexible modes occur at significantly higher frequencies than the expected structural response, i.e. above 35 Hz and 20 Hz, respectively. While the lander and solar panel resonances might be too weak in amplitude to be recorded by SEIS, the levelling system resonances will show up clearly in horizontal spectra, the H/V and ellipticity curves. They are not removed by trying to extract only Rayleigh-wave dominated parts of the data. However, they can be distinguished from any subsurface response by their exceptionally low damping ratios of 1% or less as determined by random decrement analysis. The same applies to lander-generated signals observed in actual data from a Moon analogue experiment, so we expect this analysis will be useful in identifying instrumentation resonances in SEIS data.     

Crustal S-Wave Velocity from Apparent Incidence Angles: A Case Study in Preparation for InSight

Retrieval of crustal structure and thickness of Mars is among the main goals of InSight. Here we investigate which constraints on the crust at the landing site can be provided by apparent P-wave incidence angles derived from P-receiver functions. We consider receiver functions for six different Mars models, calculated from synthetic seismograms generated via Instaseis from the Green’s function databases of the Marsquake Service, in detail. To allow for a larger range of crustal thicknesses and structures, we additionally analyze data from five broad-band stations across Central Europe. We find that the likely usable epicentral distance range for P-wave receiver functions on Mars lies between \(35^{\circ}\) and the core shadow, and can be extended to more than \(150^{\circ}\) by also using the PP-phase. Comparison to models for the spatial distribution of Martian seismicity indicates that sufficient seismicity should occur within the P-wave distance range around InSight within the nominal mission duration to allow for the application of our method. Apparent P-wave incidence angles are derived from the amplitudes of vertical and radial receiver functions at the P-wave onset within a range of period bands, up to 120 s. The apparent incidence angles are directly related to apparent S-wave velocities, which are inverted for the subsurface S-wave velocity structure via a grid search. The veracity of the forward calculated receiver functions and apparent S-wave velocities is ensured by benchmarking various algorithms against the Instaseis synthetics. Results indicate that apparent S-wave velocity curves provide valuable constraints on crustal thickness and structure, even without any additional constraints, and considering the location uncertainty and limited data quantity of InSight. S-wave velocities in the upper half of the crust are constrained best, but if reliable measurements at long periods are available, the curves also provide constraints down to the uppermost mantle. Besides, it is demonstrated that the apparent velocity curves can differentiate between crustal velocity models that are indistinguishable by other methods.     

Planned Products of the Mars Structure Service for the InSight Mission to Mars

Space Science Reviews - The InSight lander will deliver geophysical instruments to Mars in 2018, including seismometers installed directly on the surface (Seismic Experiment for Interior Structure,...     

The Marsquake Service: Securing Daily Analysis of SEIS Data and Building the Martian Seismicity Catalogue for InSight

Space Science Reviews - In recent decades, volcanic and cryovolcanic activity on moons within the Solar System has been recognised as an important source of cosmic dust. Two moons, Jupiter’s...     

Rayleigh Wave Ellipticity Modeling and Inversion for Shallow Structure at the Proposed InSight Landing Site in Elysium Planitia,Mars

The SEIS (Seismic Experiment for Interior Structure) instrument onboard the InSight mission will be the first seismometer directly deployed on the surface of Mars. From studies on the Earth and the Moon, it is well known that site amplification in low-velocity sediments on top of more competent rocks has a strong influence on seismic signals, but can also be used to constrain the subsurface structure. Here we simulate ambient vibration wavefields in a model of the shallow sub-surface at the InSight landing site in Elysium Planitia and demonstrate how the high-frequency Rayleigh wave ellipticity can be extracted from these data and inverted for shallow structure. We find that, depending on model parameters, higher mode ellipticity information can be extracted from single-station data, which significantly reduces uncertainties in inversion. Though the data are most sensitive to properties of the upper-most layer and show a strong trade-off between layer depth and velocity, it is possible to estimate the velocity and thickness of the sub-regolith layer by using reasonable constraints on regolith properties. Model parameters are best constrained if either higher mode data can be used or additional constraints on regolith properties from seismic analysis of the hammer strokes of InSight’s heat flow probe HP³ are available. In addition, the Rayleigh wave ellipticity can distinguish between models with a constant regolith velocity and models with a velocity increase in the regolith, information which is difficult to obtain otherwise.     

A Numerical Model of the SEIS Leveling System Transfer Matrix and Resonances: Application to SEIS Rotational Seismology and Dynamic Ground Interaction

Both sensors of the SEIS instrument (VBBs and SPs) are mounted on the mechanical leveling system (LVL), which has to ensure a level placement on the Martian ground under currently unknown local conditions, and provide the mechanical coupling of the seismometers to the ground. We developed a simplified analytical model of the LVL structure in order to reproduce its mechanical behavior by predicting its resonances and transfer function. This model is implemented numerically and allows to estimate the effects of the LVL on the data recorded by the VBBs and SPs on Mars. The model is validated through comparison with the horizontal resonances (between 35 and 50 Hz) observed in laboratory measurements. These modes prove to be highly dependent of the ground horizontal stiffness and torque. For this reason, an inversion study is performed and the results are compared with some experimental measurements of the LVL feet’s penetration in a martian regolith analog. This comparison shows that the analytical model can be used to estimate the elastic ground properties of the InSight landing site. Another application consists in modeling the 6 sensors on the LVL at their real positions, also considering their sensitivity axes, to study the performances of the global SEIS instrument in translation and rotation. It is found that the high frequency ground rotation can be measured by SEIS and, when compared to the ground acceleration, can provide ways to estimate the phase velocity of the seismic surface waves at shallow depths. Finally, synthetic data from the active seismic experiment made during the HP³ penetration and SEIS rotation noise are compared and used for an inversion of the Rayleigh phase velocity. This confirms the perspectives for rotational seismology with SEIS which will be developed with the SEIS data acquired during the commissioning phase after landing.     

Geology and Physical Properties Investigations by the InSight Lander

Although not the prime focus of the InSight mission, the near-surface geology and physical properties investigations provide critical information for both placing the instruments (seismometer and heat flow probe with mole) on the surface and for understanding the nature of the shallow subsurface and its effect on recorded seismic waves. Two color cameras on the lander will obtain multiple stereo images of the surface and its interaction with the spacecraft. Images will be used to identify the geologic materials and features present, quantify their areal coverage, help determine the basic geologic evolution of the area, and provide ground truth for orbital remote sensing data. A radiometer will measure the hourly temperature of the surface in two spots, which will determine the thermal inertia of the surface materials present and their particle size and/or cohesion. Continuous measurements of wind speed and direction offer a unique opportunity to correlate dust devils and high winds with eolian changes imaged at the surface and to determine the threshold friction wind stress for grain motion on Mars. During the first two weeks after landing, these investigations will support the selection of instrument placement locations that are relatively smooth, flat, free of small rocks and load bearing. Soil mechanics parameters and elastic properties of near surface materials will be determined from mole penetration and thermal conductivity measurements from the surface to 3–5 m depth, the measurement of seismic waves during mole hammering, passive monitoring of seismic waves, and experiments with the arm and scoop of the lander (indentations, scraping and trenching). These investigations will determine and test the presence and mechanical properties of the expected 3–17 m thick fragmented regolith (and underlying fractured material) built up by impact and eolian processes on top of Hesperian lava flows and determine its seismic properties for the seismic investigation of Mars’ interior.     

Analysis of Regolith Properties Using Seismic Signals Generated by InSight’s HP<Superscript>3</Superscript> Penetrator

InSight’s Seismic Experiment for Interior Structure (SEIS) provides a unique and unprecedented opportunity to conduct the first geotechnical survey of the Martian soil by taking advantage of the repeated seismic signals that will be generated by the mole of the Heat Flow and Physical Properties Package (HP³). Knowledge of the elastic properties of the Martian regolith have implications to material strength and can constrain models of water content, and provide context to geological processes and history that have acted on the landing site in western Elysium Planitia. Moreover, it will help to reduce travel-time errors introduced into the analysis of seismic data due to poor knowledge of the shallow subsurface. The challenge faced by the InSight team is to overcome the limited temporal resolution of the sharp hammer signals, which have significantly higher frequency content than the SEIS 100 Hz sampling rate. Fortunately, since the mole propagates at a rate of \(\sim1~\mbox{mm}\) per stroke down to 5 m depth, we anticipate thousands of seismic signals, which will vary very gradually as the mole travels.Using a combination of field measurements and modeling we simulate a seismic data set that mimics the InSight HP3-SEIS scenario, and the resolution of the InSight seismometer data. We demonstrate that the direct signal, and more importantly an anticipated reflected signal from the interface between the bottom of the regolith layer and an underlying lava flow, are likely to be observed both by Insight’s Very Broad Band (VBB) seismometer and Short Period (SP) seismometer. We have outlined several strategies to increase the signal temporal resolution using the multitude of hammer stroke and internal timing information to stack and interpolate multiple signals, and demonstrated that in spite of the low resolution, the key parameters—seismic velocities and regolith depth—can be retrieved with a high degree of confidence.     

A Pre-Landing Assessment of Regolith Properties at the InSight Landing Site

This article discusses relevant physical properties of the regolith at the Mars InSight landing site as understood prior to landing of the spacecraft. InSight will land in the northern lowland plains of Mars, close to the equator, where the regolith is estimated to be \(\geq3\mbox{--}5~\mbox{m}\) thick. These investigations of physical properties have relied on data collected from Mars orbital measurements, previously collected lander and rover data, results of studies of data and samples from Apollo lunar missions, laboratory measurements on regolith simulants, and theoretical studies. The investigations include changes in properties with depth and temperature. Mechanical properties investigated include density, grain-size distribution, cohesion, and angle of internal friction. Thermophysical properties include thermal inertia, surface emissivity and albedo, thermal conductivity and diffusivity, and specific heat. Regolith elastic properties not only include parameters that control seismic wave velocities in the immediate vicinity of the Insight lander but also coupling of the lander and other potential noise sources to the InSight broadband seismometer. The related properties include Poisson’s ratio, P- and S-wave velocities, Young’s modulus, and seismic attenuation. Finally, mass diffusivity was investigated to estimate gas movements in the regolith driven by atmospheric pressure changes. Physical properties presented here are all to some degree speculative. However, they form a basis for interpretation of the early data to be returned from the InSight mission.     

Flexible Mode Modelling of the InSight Lander and Consequences for the SEIS Instrument

We present an updated model for estimating the lander mechanical noise on the InSight seismometer SEIS, taking into account the flexible modes of the InSight lander. This new flexible mode model uses the Satellite Dynamics Toolbox to compute the direct and the inverse dynamic model of a satellite composed of a main body fitted with one or several dynamic appendages. Through a detailed study of the sensitivity of our results to key environment parameters we find that the frequencies of the six dominant lander resonant modes increase logarithmically with increasing ground stiffness. On the other hand, the wind strength and the incoming wind angle modify only the signal amplitude but not the frequencies of the resonances. For the baseline parameters chosen for this study, the lander mechanical noise on the SEIS instrument is not expected to exceed the instrument total noise requirements. However, in the case that the lander mechanical noise is observable in the seismic data acquired by SEIS, this may provide a complementary method for studying the ground and wind properties on Mars.