Recovery of Deep Moonquake Focal Mechanisms

Deep moonquakes are clustered not only in space but also in time: their recurrence times correspond to the durations of the anomalistic and draconic months, with some clusters preferring one of the two periods, while others are active with both periods. A key constraint for the understanding of the connection between the orbital motion of the Moon and its seismic activity is the focal mechanism: the orientation of the fault surface on which failure occurs during the quake. Due to the small aperture of the Apollo seismic network and the strong scattering of seismic waves within the lunar crust, the evaluation of P wave first motions to constrain the strike and dip of the fault planes is not feasible. Instead we evaluate the amplitude ratios of P and S waves. Seismograms are rotated into the P-SV-SH coordinate frame and amplitudes are determined as averages over short time windows after the arrival to reduce the impact of the scattering coda, which is independent of the source orientation. We allow for reversals of the fault motion, as observed for some clusters in previous studies, by taking into account the absolute amplitude only, without sign. An empirical site correction factor is applied to correct for amplitude distortions in the crust. We construct ensembles of fault plane solutions using an exhaustive grid search by accepting all orientations that reproduce the measured amplitude ratios within the observed standard deviations. Since all events of a given cluster are supposed to share the same fault plane, the combination of the individual inversion results further constrains the orientation. We evaluate 106 events from 25 different moonquake clusters. The most active cluster A001 contributes 37 events, while others contribute 1 to 9 events per cluster. Comparison of fault orientations with the variation of the tidal stress results in preferred orientations

Implications of the Global Surface Fault Distribution and of Lithospheric Cooling

This contribution presents a model that links the observed distribution of surface faults to the spatial distribution of marsquakes. The annual seismic moment budget is computed based on the as-sumption that global cooling and subsequent shrink-ing of Mars is the main source of strain today [1]. A truncated Gutenberg-Richter distribution is used to re-late the seismic moment budget to marsquake frequen-cies. We have derived a theoretical relation for the limitation of quake size by the lengths of the individual faults. This relation is used for the simulation of epi-center catalogs that may serve as input data for the development of seismological experiments

Compressive strength and elastic modulus of Comet 67P interpreted from a material science point of view

The analysis of Cometary Acoustic Surface Sounding Experiment (CASSE) data yielded values of surface compression strength and elastic modulus at the landing site Agilkia. These data are interpreted with fracture mechanical concepts from material science taking into account the high porosity of Comet 67P

Present-day Mars' seismicity predicted from 3-D thermal evolution models of interior dynamics

©2018. American Geophysical UnionThe Interior Exploration using Seismic Investigations, Geodesy and Heat Transport mission, to be launched in 2018, will perform a comprehensive geophysical investigation of Mars in situ. The Seismic Experiment for Interior Structure package aims to detect global and regional seismic events and in turn offer constraints on core size, crustal thickness, and core, mantle, and crustal composition. In this study, we estimate the present‐day amount and distribution of seismicity using 3‐D numerical thermal evolution models of Mars, taking into account contributions from convective stresses as well as from stresses associated with cooling and planetary contraction. Defining the seismogenic lithosphere by an isotherm and assuming two end‐member cases of 573 K and the 1073 K, we determine the seismogenic lithosphere thickness. Assuming a seismic efficiency between 0.025 and 1, this thickness is used to estimate the total annual seismic moment budget, and our models show values between 5.7 × 1016 and 3.9 × 1019 Nm

Re-evaluation of Apollo 17 Lunar Seismic Profiling Experiment data

We re-analyzed Apollo 17 Lunar Seismic Profiling Experiment (LSPE) data to improve our knowledge of the subsurface structure of this landing site. We use new geometrically accurate 3-D positions of the seismic equipment deployed by the astronauts, which were previously derived using high-resolution images by Lunar Reconnaissance Orbiter (LRO) in combination with Apollo astronaut photography. These include coordinates of six Explosive Packages (EPs) and four geophone stations. Re-identified P-wave arrival times are used to calculate two- and three-layer seismic velocity models. A strong increase of seismic velocity with depth can be confirmed, in particular, we suggest a more drastic increase than previously thought. For the three-layer model the P-wave velocities were calculated to 285, 580, and 1825 m/s for the uppermost, second, and third layer, respectively, with the boundaries between the layers being at 96 and 773 m depth. When compared with results obtained with previously published coordinates, we find (1) a slightly higher velocity (+4%) for the uppermost layer, and (2) lower P-wave velocities for the second and third layers, representing a decrease of 34% and 12% for second and third layer, respectively. Using P-wave arrival time readings of previous studies, we confirm that velocities increase when changing over from old to new coordinates. In the three-layer case, this means using new coordinates alone leads to thinned layers, velocities rise slightly for the uppermost layer and decrease significantly for the layers below

Modeling Approaches in Planetary Seismology

Of the many geophysical means that can be used to probe a planet's interior, seismology remains the most direct. Given that the seismic data gathered on the Moon over 40 years ago revolutionized our understanding of the Moon and are still being used today to produce new insight into the state of the lunar interior, it is no wonder that many future missions, both real and conceptual, plan to take seismometers to other planets. To best facilitate the return of high-quality data from these instruments, as well as to further our understanding of the dynamic processes that modify a planet's interior, various modeling approaches are used to quantify parameters such as the amount and distribution of seismicity, tidal deformation, and seismic structure on and of the terrestrial planets. In addition, recent advances in wavefield modeling have permitted a renewed look at seismic energy transmission and the effects of attenuation and scattering, as well as the presence and effect of a core, on recorded seismograms. In this chapter, we will review these approaches

Estimation of the Seismic Moment Release Rate of Mars from InSight Seismic Data

Seismicity models for Mars usually estimate the long-term average annual seismic moment rate, and also the average annual event rate. This holds for estimations based on geological evidence (Golombek et al., 1992, Golombek, 2002, Taylor et al., 2013) as well as for models based on thermal evolution and cooling of the Martian interior (Phillips, 1991, Knapmeyer et al., 2006, Plesa et al., 2018). All studies are compatible with the conclusion based on the non-observation of any unambiguous event by Viking (Anderson et al., 1977, Goins & Lazarewicz, 1979) that Martian seismicity lies somewhere between that of the Moon and that of the Earth. We developed tools to derive reasonable estimations of the annual seismic moment rate from a number of events as small as one, provided that the observed events are beyond the global completeness threshold for observable events. Numerical tests as well as evaluation of terrestrial data shows the feasibility of the approach

SEIS: Insight's Seismic Experiment for Internal Structure of Mars.

By the end of 2018, 42 years after the landing of the two Viking seismometers on Mars, InSight will deploy onto Mars' surface the SEIS (Seismic Experiment for Internal Structure) instrument; a six-axes seismometer equipped with both a long-period three-axes Very Broad Band (VBB) instrument and a three-axes short-period (SP) instrument. These six sensors will cover a broad range of the seismic bandwidth, from 0.01 Hz to 50 Hz, with possible extension to longer periods. Data will be transmitted in the form of three continuous VBB components at 2 sample per second (sps), an estimation of the short period energy content from the SP at 1 sps and a continuous compound VBB/SP vertical axis at 10 sps. The continuous streams will be augmented by requested event data with sample rates from 20 to 100 sps. SEIS will improve upon the existing resolution of Viking's Mars seismic monitoring by a factor of ∼ 2500 at 1 Hz and ∼ 200 000 at 0.1 Hz. An additional major improvement is that, contrary to Viking, the seismometers will be deployed via a robotic arm directly onto Mars' surface and will be protected against temperature and wind by highly efficient thermal and wind shielding. Based on existing knowledge of Mars, it is reasonable to infer a moment magnitude detection threshold of M w ∼ 3 at 40 ∘ epicentral distance and a potential to detect several tens of quakes and about five impacts per year. In this paper, we first describe the science goals of the experiment and the rationale used to define its requirements. We then provide a detailed description of the hardware, from the sensors to the deployment system and associated performance, including transfer functions of the seismic sensors and temperature sensors. We conclude by describing the experiment ground segment, including data processing services, outreach and education networks and provide a description of the format to be used for future data distribution.Electronic supplementary materialThe online version of this article (10.1007/s11214-018-0574-6) contains supplementary material, which is available to authorized users

Deep Moonquake Focal Mechanisms: Recovery and Implications

A defining characteristic of deep moonquakes is their tendency to occur with tidal periodicity, prompting previous studies to infer that they are related to the buildup and release of tidal stress within the Moon. In studies of tidal forcing, a key constraint is the focal mechanism: the fault parameters describing the type of failure moonquakes represent. The quality of the lunar seismic data and the limited source/receiver geometries of the Apollo seismic network prohibit the determination of deep moonquake fault parameters using first-motion polarities, as is typically done in terrestrial seismology. Without being able to resolve tidal stress onto a known failure plane, we can examine only gross qualities of the tidal stress tensor with respect to moonquake occurrence, so we cannot fully address the role of tidal stress in moonquake generation. We will examine the extent to which shear (S) and compression (P) wave amplitude ratios can constrain moonquake fault geometry by determining whether, for a given cluster, there exists a focal mechanism that can produce a radiation pattern consistent with the amplitudes measured by the Apollo instruments. Amplitudes are read in the ray coordinate frame, directly from seismograms for which the P and S arrivals are clearly identifiable on all long-period channels of the four Apollo stations. We apply an empirical station correction to account for site effects and the differences between P- and S-wave attenuation. Instead of focusing on the best fitting solution only, we formulate the inverse problem using a falsification criterion: all source orientations that do not reproduce the observed SV/P ratios within an error margin derived from the uncertainty of amplitude readings are rejected. All others are accepted as possible solutions. The inversion is carried out using an exhaustive grid search on a regular grid with predefined step size, encompassing all possible combinations of strike, dip and slip. To assess the sensitivity of the inversion for the uncertainty of the lunar interior structure, we carry out repeated inversions with different velocity structures. Our data set consist of a total of 106 events from 25 deep moonquake clusters. The largest contribution of 37 events originates from the most active cluster, A001, while other clusters are represented by 1 to 9 events. Since the definition of a cluster implies that all events share the same source orientation, a comparison of the inversion results of all events from one cluster will reduce ambiguities of the inversion. Once we obtain a suite of fault parameters for a given source, we can attempt to further constrain the focal mechanism with refined analyses of tidal stresses and predictions based on synthetic seismograms