Seismic constraints from a Mars impact experiment using InSight and Perseverance

NASA’s InSight (Interior Exploration using Seismic Investigations, Geodesy and Heat Transport) mission has operated a sophisticated suite of seismology and geophysics instruments on the surface of Mars since its arrival in 2018. On 18 February 2021, we attempted to detect the seismic and acoustic waves produced by the entry, descent and landing of the Perseverance rover using the sensors onboard the InSight lander. Similar observations have been made on Earth using data from both crewed1,2 and uncrewed3,4 spacecraft, and on the Moon during the Apollo era5, but never before on Mars or another planet. This was the only seismic event to occur on Mars since InSight began operations that had an a priori known and independently constrained timing and location. It therefore had the potential to be used as a calibration for other marsquakes recorded by InSight. Here we report that no signal from Perseverance’s entry, descent and landing is identifiable in the InSight data. Nonetheless, measurements made during the landing window enable us to place constraints on the distance–amplitude relationships used to predict the amplitude of seismic waves produced by planetary impacts and place in situ constraints on Martian impact seismic efficiency (the fraction of the impactor kinetic energy converted into seismic energy). As the Perseverance rover landed on the Martian surface, the sensors on NASA’s InSight Mars lander picked up no seismic or acoustic waves. This non-detection provides information on the crust and atmosphere of Mars.


NASA's InSight (Interior Exploration using Seismic Investigations, Geodesy and Heat Transport) mission has operated a sophisticated suite of seismology and geophysics instruments on the surface of Mars since its arrival in 2018.
On 18 February 2021, we attempted to detect the seismic and acoustic waves produced by the entry, descent and landing of the Perseverance rover using the sensors onboard the InSight lander. Similar observations have been made on Earth using data from both crewed 1,2 and uncrewed 3,4 spacecraft, and on the Moon during the Apollo era 5 , but never before on Mars or another planet. This was the only seismic event to occur on Mars since InSight began operations that had an a priori known and independently constrained timing and location. It therefore had the potential to be used as a calibration for other marsquakes recorded by InSight. Here we report that no signal from Perseverance's entry, descent and landing is identifiable in the InSight data. Nonetheless, measurements made during the landing window enable us to place constraints on the distance-amplitude relationships used to predict the amplitude of seismic waves produced by planetary impacts and place in situ constraints on Martian impact seismic efficiency (the fraction of the impactor kinetic energy converted into seismic energy).
Pre-landing predictions based on modelling and on data from Mars Science Laboratory (Curiosity)'s entry, descent and landing (EDL) in 2012 suggested that the impact of the spacecraft's two cruise balance mass devices (CBMDs, also referred to as cruise mass balance devices, or CMBDs) might produce seismic waves of high enough amplitude to be detectable at the position of InSight (Interior Exploration using Seismic Investigations, Geodesy and Heat Transport), approximately 3,450 km to the east 6,7 .
The CBMDs are solid, 77.5 kg tungsten blocks used to adjust the spacecraft's lift-to-drag ratio during EDL. Data from the spacecraft's computer indicate that the commands were generated to fire the pyrotechnic releases holding them in place at 20:28:27 utc, within 1 s of each other. At this time, the spacecraft was projected to be at an altitude of 1,253 km and travelling at a planet-relative velocity of 4,753 m s −1 .
NASA's trajectory reconstruction indicates that the CBMDs hit the surface at 20:40:33 ± 3 s utc (around five minutes before the rover's touchdown), at a speed of approximately 3,816 m s −1 and an oblique angle of around 10° from the horizontal.
The impact craters from the CBMDs were imaged by the CTX (Context Camera) 8 and HiRISE (High Resolution Imaging Science Experiment) instruments 9 on board NASA's Mars Reconnaissance Orbiter on 3 May 2021. Images identified the craters at a position around 18.9° N 76.2° E (Fig. 1).
The most promising candidate seismic phase for detection at InSight (4.5° N, 135.6° E) was expected to be a ballistic compressional (P) wave excited by the CBMD impact ( Fig. 5 in ref. 7 ). The P-wave arrival time was predicted to be approximately 420 ± 20 s after impact, or 20:47 utc (19:50 InSight Local Mars Solar Time, lmst, on sol 793).
The two variables determining whether or not this signal would be identifiable in InSight data were the seismic noise during the arrival window (constrained after landing from the seismometer recordings), and the amplitude of the impact-induced P wave. The former was particularly low during the arrival window, as it occurred during the part of InSight's day when the atmospheric noise is lowest. The latter was predicted using distance-amplitude scaling curves 7,10,11 . These relationships come with large uncertainties, as they are calibrated using only terrestrial and lunar data at closer range.

NATurE ASTrONOMy
Data were recorded throughout the night (in utc) of [18][19] February 2021 on InSight's Very Broad Band (VBB) seismometer 12 , and the pressure and wind sensors of the Auxiliary Payload Sensor Suite (APSS) 13 .
Data from the arrival window are plotted in Fig. 2. The main signals observed are a large marsquake around 01:25 utc, and irregularly spaced glitches, which are artificial in origin. No other signal is observed that cannot be explained as noise excited by atmospheric phenomena, as corroborated using the wind speed and pressure measurements.
Therefore, we conclude that no signal associated with the CBMD impact is identifiable above the detection threshold (as defined in Fig. 2b), precluding a P-wave amplitude larger than 1.1 × 10 −10 m s −1 .
As the marsquake detected at 01:25 utc fits within the projected time window (4-5 h after landing) for infrasound (low-frequency sound) waves to arrive from Perseverance's EDL, we will briefly outline why we believe it is unrelated and simply coincidental.
This assessment is based on the signal's amplitude, shape and frequency content, which are entirely different from those of an infrasound wave. EDL-related infrasounds at these distances should have amplitudes below the noise floor 7 , whilst the high frequencies observed in this event (up to ~30 Hz) preclude an airborne propagation path, as the high CO 2 concentration in the Martian atmosphere would rapidly attenuate them 14 . Rather, this event ('S0794a') is a standard 'very high-frequency' marsquake 15 , in this case with an origin at a distance of ~1,100 km, much closer than Perseverance. As of the end of April 2021, 40 of these events have been recorded by InSight, often in the late evening 16 . These events are probably caused by tectonic or otherwise internal geological processes.
Having established that there is no signal from the EDL recorded in the data, we will now consider what this non-detection can be used to infer about impact processes on Mars.
There are a number of approaches to using distance-amplitude scaling relationships to predict peak P-wave amplitudes from impacts.
Depending on whether the impactor energy, total momentum or vertical component of the momentum is used to scale the amplitude, and whether terrestrial missile-impact data or lunar spacecraft-impact data are the basis, predicted amplitudes vary by up to two orders of magnitude when extrapolated to distances of 3,500 km. Five such standard scaling relationships are shown in Fig. 3-one based on impactor energy 10 , and four based on impactor momentum 11 .
Of the latter, two are based on data from artificial lunar impacts, which occurred into almost cohesionless material (black curves, ref. 5 ), and two on terrestrial missile impacts into weakly cohesive regolith soils (green curves, ref. 17 ). In each case, one curve uses the total impactor momentum as the determinant of peak P-wave amplitude (solid lines), whilst the other uses the vertical component of the impactor momentum (dashed lines).
The derived upper bound on the peak P-wave amplitude of 1.1 × 10 −10 m s −1 is shown as a horizontal grey line in Fig. 3, and at the distance in question lies below the solid green scaling curve (total impactor momentum; terrestrial missile data).
This result indicates that the distance-amplitude scaling relationship based on terrestrial missile impacts into weakly cohesive soil, and the assumption that seismic wave amplitudes scale with total impactor momentum, are not appropriate in this case. There are three possible implications.
First, for highly oblique impacts such as this, this result suggests that the vertical component of the momentum may be a more appropriate quantity to scale by.
Alternatively, as the impact data on which this model was based were collected at much closer distances and on Earth, this may indicate a stronger attenuation of seismic wave amplitudes with distance on Mars. A third possibility is that artificial impacts on the lunar surface represent a better analogue for the seismic response of small impacts on Mars than those that occur into terrestrial soils.
The upper bound on the peak P-wave amplitude may be used to place a joint constraint on the impact seismic efficiency (k s , which is site and impact specific) and the average mantle attenuation (Q μ ) along the path from Perseverance's landing site to InSight.
Because k s and Q μ are entirely independent of each other (the former being related to generation of seismic waves on a local scale, and the latter to their propagation on a global scale), jointly constraining their values in this way is valid. k s is particularly poorly constrained, due to its high sensitivity to local conditions and the lack of relevant in situ measurements 18 . All estimates for Mars thus far have therefore used modelling, simulation or material analogues to estimate its value, and no in situ measurements other than those in this paper exist.
Q μ and Q κ are the two quality factors used to describe viscoelastic (intrinsic/inelastic) attenuation properties within the solid part of the planet. The former is associated with shear properties, and the latter with bulk properties.
For typical Martian mantle rheologies, Q μ ≫ Q κ , hence Q κ has little influence on P-wave amplitudes. The P-wave attenuation Q P can then be approximated as Q P = 9 4 Qμ. Tidal observations suggest that Q μ does not vary strongly through the mantle 19 , so we assume that Q μ is a reasonable descriptor of the average attenuation along the source-receiver path (Methods).
We do not treat the crust separately to the mantle, as the P-wave propagates almost entirely within the latter and previous studies show comparably negligible crustal attenuation 15,20 . To relate our observed upper bound on the P-wave amplitude to Q μ and k s , we first estimate the seismic moment associated with the CBMD impact using the following empirically derived relationship from ref. 21 : where E k is the impactor kinetic energy (in this case, 1.1 GJ). Note that the separation between the CBMDs in space and time at impact (~1 km and around 1 s) is large enough that the impact processes In both panels, glitches in the system are recorded as sharp vertical features in the spectrogram and peaks in the r.m.s. envelope, with one exemplar highlighted. These glitches are easily distinguished from seismic signals 33 . The arrival window for the CBMD P wave is highlighted in red, as is an unrelated marsquake observed around 01:45 utc. c, The wind speed and the r.m.s. envelope of the vertical seismometer velocity in the 3.9-4.5 Hz frequency band-the latter, in this frequency range, contains a known oscillation mode of the spacecraft, which is excited by the wind and can be used as a proxy for the wind speed 34 . The absence of wind measurements around 21:00 utc occurs where the wind speed drops below the instrument threshold 13 . d, The r.m.s. envelope of the atmospheric pressure in the 0.1-4 Hz band, which is also correlated to wind speed 35 . e,f, A detail of the grey area in a and b, respectively, where the largest-amplitude event is a marked glitch just before 20:50 utc, and our defined 'detection threshold' occur independently of each other. E k is thus the kinetic energy of a single CBMD.
As the P wave reduces in amplitude during its propagation through both attenuation and geometric spreading of the wavefront, making an estimate of Q μ requires us to quantify the effects of the latter.
To do this, we undertake full-waveform seismic simulations using the AxiSEM method 22 in a purely elastic (non-attenuating) medium. Thus, the only energy loss in these simulations along the source-receiver path occurs due to the spreading of the wavefront. The effects of attenuation using an average Q μ may then be applied as a post-simulation correction (Methods). Now, the values of k s and Q μ may be independently varied (using equations (1) and (3), respectively) to determine their joint effect on the amplitude recorded by InSight. Figure 4 shows joint constraints that may be placed on k s for this impact scenario and Q μ : as increasing the value of either parameter leads to larger predicted amplitudes, these are co-constrained such that the resultant P-wave amplitude does not exceed our detection threshold.
The quality factor has previously been determined to lie in the range 300-1,000, with the lower bound derived from tidal observations 19 and the upper from observations of the spectra of marsquakes 20,23 . Thus, we constrain the maximum value of k s to be 3%, corresponding to the rightmost corner of the green zone in Fig. 4.
This upper limit on k s is compatible with experimental (terrestrial and lunar) and simulated (terrestrial, lunar and Martian) measurements, where k s ranges from 10 −6 (in porous sand/regolith) to 10 −3 (in stronger, non-porous materials) 5,11,17,24,25 . Values up to 2-3% have been observed in underground nuclear explosions in stronger target materials 26 . Figure 1 shows one of the impact craters. The surrounding morphology is representative of this region of Mars, and consists of aeolian bedforms and unconsolidated/poorly consolidated regolith. The thermal inertia of the surface at the impact location is ~250 J m −2 s −1/2 K −1 (ref. 27 ), indicative of a dominantly sandy surface. Although a quantitative relationship between k s and target properties has yet to be derived, our upper bound on k s is consistent with a surface material of this type.
The true seismic efficiency is likely to be lower than the maximum that we derived. Nonetheless, this does still demonstrate a practical method through which it may be constrained, and is notable that such a constraint has been derived on Mars. From this, we may draw the robust conclusion that the conversion of kinetic energy into seismic energy on Mars is no more efficient than the most efficient such terrestrial coupling.
We used the non-detection of the seismic waves from the impact of Perseverance's CBMDs to show that the total impactor momentum is a poor predictor of amplitude at this distance and impact angle, if a terrestrial-based scaling is assumed. In this case, the lunar-based scaling or a relationship based on impactor energy may be more appropriate. This result could also indicate that effects of attenuation on impact-generated seismic waves are stronger than previously estimated.
We also used the non-detection to constrain the impact seismic efficiency to be less than 3%, which is compatible with geological analysis of the impact site, and commensurate with previous estimates, which used modelling or proxies in place of in situ data.
The methodology presented here provides a basis for using seismic (non-) detection of artificial impacts to infer subsurface properties, and could in future be applied during geophysical missions to any of the Solar System's icy or rocky bodies.

Methods
We use an interior velocity and density model called 'TAYAK' , which combines geochemical data with geodetic constraints 28,29 . This was chosen on the basis of its

NATurE ASTrONOMy
good fit to other marsquake data 30 . The source is represented as an explosion just below the surface, as per ref. 18 . Simulations are conducted using the spectral element solver AxiSEM 22 in combination with Instaseis 31 , with a dominant source period of 1 s, and with attenuation switched off. This allows us to account for amplitude decreases due to the geometric spreading of the wavefront.
Once the amplitude in the non-attenuating case A 0 (where A 0 = A 0 (M 0 )) has been found, we account for the effects of attenuation, assuming an effective Q P,eff that averages over the propagation path S: The amplitude in a frequency band between f 1 and f 2 (here 0.2 and 0.9 Hz) is estimated through application of the following equation: where A is the amplitude at InSight. Note that Q P = 9/4Q μ , assuming a Poisson solid with a standard acoustic speed to shear speed ratio v p /v s of √ 3. t = 420 s is the predicted travel time of a P wave. Q is assumed to be frequency independent.
Combining equations (3), as derived above, and (1) allows us to determine which combinations of k s and Q P,eff (and hence k s and Q μ ) produce a permissible A at InSight's position, that is, one below the detection threshold.

Data availability
InSight APSS/TWINS/PS data can be found at https://atmos.nmsu.edu/data_ and_services/atmospheres_data/INSIGHT/insight.html. InSight SEIS data are available in the form of a seismic event catalogue and waveform data (https://doi. org/10.18715/SEIS.INSIGHT.XB_2016) that are publicly available from the IPGP Data Center and IRIS-DMC, as well as raw data available in the PDS (https:// pds-geosciences.wustl.edu/missions/insight/seis.htm). Data used here can be found in version 7 of the Mars Quake Service catalogue (https://doi.org/10.12686/a12). HiRISE data are publicly available through the Planetary Data System at https://www.uahirise.org/. Seismic modelling used the open-source Martian interior model TAYAK from refs. 29,32 , available at https://instaseis.ethz.ch, and peak amplitudes were computed using the open-source AxiSEM method of ref. 22 .