Search

Abstract: On 4 May 2022 the InSight seismometer SEIS‐VBB recorded the largest marsquake ever observed, S1222a, with an initial magnitude estimate of 4.6. Understanding the depth and source properties of this event has important implications for the nature of tectonic activity on Mars. Located ∼37° to the southeast of InSight, S1222a is one of the few non‐impact marsquakes that exhibits prominent surface waves. We use waveform modeling of body waves (P and S) and surface waves (Rayleigh and Love) to constrain the focal mechanism, assuming a double‐couple source, and find that S1222a likely resulted from reverse faulting in the crust (source depth near 22 km). We estimate the scalar moment to be 2.5 × 1015–3.5 × 1015 Nm (magnitude MW 4.2–4.3). Our results suggest active compressional tectonics near the dichotomy boundary on Mars, likely due to thermal contraction from planetary cooling.

The detection of deep reflected S waves on Mars inferred a core size of 1,830 ± 40 km (ref. 1 ), requiring light-element contents that are incompatible with experimental petrological constraints. This estimate assumes a compositionally homogeneous Martian mantle, at odds with recent measurements of anomalously slow propagating P waves diffracted along the core-mantle boundary 2 . An alternative hypothesis is that Mars's mantle is heterogeneous as a consequence of an early magma ocean that solidified to form a basal layer enriched in iron and heat-producing elements. Such enrichment results in the formation of a molten silicate layer above the core, overlain by a partially molten layer 3 . Here we show that this structure is compatible with all geophysical data, notably (1) deep reflected and diffracted mantle seismic phases, (2) weak shear attenuation at seismic frequency and (3) Mars's dissipative nature at Phobos tides. The core size in this scenario is 1,650 ± 20 km, implying a density of 6.5 g cm -3 , 5-8% larger than previous seismic estimates, and can be explained by fewer, and less abundant, alloying light elements than previously required, in amounts compatible with experimental and cosmochemical constraints. Finally, the layered mantle structure requires external sources to generate the magnetic signatures recorded in Mars's crust., (© 2023. The Author(s).)

We present the first observations of seismic waves propagating through the core of Mars. These observations, made using seismic data collected by the InSight geophysical mission, have allowed us to construct the first seismically constrained models for the elastic properties of Mars' core. We observe core-transiting seismic phase SKS from two farside seismic events detected on Mars and measure the travel times of SKS relative to mantle traversing body waves. SKS travels through the core as a compressional wave, providing information about bulk modulus and density. We perform probabilistic inversions using the core-sensitive relative travel times together with gross geophysical data and travel times from other, more proximal, seismic events to seek the equation of state parameters that best describe the liquid iron-alloy core. Our inversions provide constraints on the velocities in Mars' core and are used to develop the first seismically based estimates of its composition. We show that models informed by our SKS data favor a somewhat smaller (median core radius = 1,780 to 1,810 km) and denser (core density = 6.2 to 6.3 g/cm 3 ) core compared to previous estimates, with a P-wave velocity of 4.9 to 5.0 km/s at the core-mantle boundary, with the composition and structure of the mantle as a dominant source of uncertainty. We infer from our models that Mars' core contains a median of 20 to 22 wt% light alloying elements when we consider sulfur, oxygen, carbon, and hydrogen. These data can be used to inform models of planetary accretion, composition, and evolution.

Seismic discontinuities in the mantle are indicators of its thermo-chemical state and offer clues to its dynamics. Ray-based seismic methods, though limited by the approximations made, have mapped mantle transition zone discontinuities in detail, but have yet to offer definitive conclusions on the presence and nature of mid-mantle discontinuities. Here, we show how to use a wave-equation-based imaging method, reverse-time migration of precursors to surface-reflected seismic body waves, to uncover both mantle transition zone and mid-mantle discontinuities, and interpret their physical nature. We observe a thinned mantle transition zone southeast of Hawaii, and a reduction in impedance contrast around 410 km depth in the same area, suggesting a hotter-than-average mantle in the region. Here, we furthermore reveal a 4000-5000 km-wide reflector in new images of the mid mantle below the central Pacific, at 950-1050 km depth. This deep discontinuity exhibits strong topography and generates reflections with polarity opposite to those originating at the 660 km discontinuity, implying an impedance reversal near 1000 km. We link this mid-mantle discontinuity to the upper reaches of deflected mantle plumes upwelling in the region. Reverse-time migration full-waveform imaging is a powerful approach to imaging Earth's interior, capable of broadening our understanding of its structure and dynamics and shrinking modeling uncertainties., (© 2023. The Author(s).)

Constraining the thermal and compositional state of the mantle is crucial for deciphering the formation and evolution of Mars. Mineral physics predicts that Mars' deep mantle is demarcated by a seismic discontinuity arising from the pressure-induced phase transformation of the mineral olivine to its higher-pressure polymorphs, making the depth of this boundary sensitive to both mantle temperature and composition. Here, we report on the seismic detection of a midmantle discontinuity using the data collected by NASA's InSight Mission to Mars that matches the expected depth and sharpness of the postolivine transition. In five teleseismic events, we observed triplicated P and S waves and constrained the depth of this discontinuity to be 1,006 [Formula: see text] 40 km by modeling the triplicated waveforms. From this depth range, we infer a mantle potential temperature of 1,605 [Formula: see text] 100 K, a result consistent with a crust that is 10 to 15 times more enriched in heat-producing elements than the underlying mantle. Our waveform fits to the data indicate a broad gradient across the boundary, implying that the Martian mantle is more enriched in iron compared to Earth. Through modeling of thermochemical evolution of Mars, we observe that only two out of the five proposed composition models are compatible with the observed boundary depth. Our geodynamic simulations suggest that the Martian mantle was relatively cold 4.5 Gyr ago (1,720 to 1,860 K) and are consistent with a present-day surface heat flow of 21 to 24 mW/m 2 .

Seismological constraints obtained from receiver function (RF) analysis provide important information about the crust and mantle structure. Here, we explore the utility of the free-surface multiple of the P-wave (PP) and the corresponding conversions in RF analysis. Using earthquake records, we demonstrate the efficacy of PPs-RFs before illustrating how they become especially useful when limited data is available in typical planetary missions. Using a transdimensional hierarchical Bayesian deconvolution approach, we compute robust P-to-S (Ps)- and PPs-RFs with InSight recordings of five marsquakes. Our Ps-RF results verify the direct Ps converted phases reported by previous RF analyses with increased coherence and reveal other phases including the primary multiple reverberating within the uppermost layer of the Martian crust. Unlike the Ps-RFs, our PPs-RFs lack an arrival at 7.2 s lag time. Whereas Ps-RFs on Mars could be equally well fit by a two- or three-layer crust, synthetic modeling shows that the disappearance of the 7.2 s phase requires a three-layer crust, and is highly sensitive to velocity and thickness of intra-crustal layers. We show that a three-layer crust is also preferred by S-to-P (Sp)-RFs. While the deepest interface of the three-layer crust represents the crust-mantle interface beneath the InSight landing site, the other two interfaces at shallower depths could represent a sharp transition between either fractured and unfractured materials or thick basaltic flows and pre-existing crustal materials. PPs-RFs can provide complementary constraints and maximize the extraction of information about crustal structure in data-constrained circumstances such as planetary missions., (© 2021. The Authors.)

Clues to a planet's geologic history are contained in its interior structure, particularly its core. We detected reflections of seismic waves from the core-mantle boundary of Mars using InSight seismic data and inverted these together with geodetic data to constrain the radius of the liquid metal core to 1830 ± 40 kilometers. The large core implies a martian mantle mineralogically similar to the terrestrial upper mantle and transition zone but differing from Earth by not having a bridgmanite-dominated lower mantle. We inferred a mean core density of 5.7 to 6.3 grams per cubic centimeter, which requires a substantial complement of light elements dissolved in the iron-nickel core. The seismic core shadow as seen from InSight's location covers half the surface of Mars, including the majority of potentially active regions-e.g., Tharsis-possibly limiting the number of detectable marsquakes., (Copyright © 2021 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

The Seismic Experiment for Interior Structure (SEIS) of the InSight mission to Mars, has been providing direct information on Martian interior structure and dynamics of that planet since it landed. Compared to seismic recordings on Earth, ground motion measurements acquired by SEIS on Mars are made under dramatically different ambient noise conditions, but include idiosyncratic signals that arise from coupling between different InSight sensors and spacecraft components. This work is to synthesize what is known about these signal types, illustrate how they can manifest in waveforms and noise correlations, and present pitfalls in structural interpretations based on standard seismic analysis methods. We show that glitches, a type of prominent transient signal, can produce artifacts in ambient noise correlations. Sustained signals that vary in frequency, such as lander modes which are affected by variations in temperature and wind conditions over the course of the Martian Sol, can also contaminate ambient noise results. Therefore, both types of signals have the potential to bias interpretation in terms of subsurface layering. We illustrate that signal processing in the presence of identified nonseismic signals must be informed by an understanding of the underlying physical processes in order for high fidelity waveforms of ground motion to be extracted. While the origins of most idiosyncratic signals are well understood, the 2.4 Hz resonance remains debated and the literature does not contain an explanation of its fine spectral structure. Even though the selection of idiosyncratic signal types discussed in this paper may not be exhaustive, we provide guidance on best practices for enhancing the robustness of structural interpretations., Competing Interests: Conflict of interest The authors acknowledge there are no conflicts of interested recorded.

Human activity causes vibrations that propagate into the ground as high-frequency seismic waves. Measures to mitigate the coronavirus disease 2019 (COVID-19) pandemic caused widespread changes in human activity, leading to a months-long reduction in seismic noise of up to 50%. The 2020 seismic noise quiet period is the longest and most prominent global anthropogenic seismic noise reduction on record. Although the reduction is strongest at surface seismometers in populated areas, this seismic quiescence extends for many kilometers radially and hundreds of meters in depth. This quiet period provides an opportunity to detect subtle signals from subsurface seismic sources that would have been concealed in noisier times and to benchmark sources of anthropogenic noise. A strong correlation between seismic noise and independent measurements of human mobility suggests that seismology provides an absolute, real-time estimate of human activities., (Copyright © 2020 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

The chemical and physical properties of the interiors of terrestrial planets are largely determined during their formation and differentiation. Modeling a planet's formation provides important insights into the properties of its core and mantle, and conversely, knowledge of those properties may constrain formational narratives. Here, we present a multi-stage model of Martian core formation in which we calculate core-mantle equilibration using parameterizations from high pressure-temperature metal-silicate partitioning experiments. We account for changing core-mantle boundary (CMB) conditions, composition-dependent partitioning, and partial equilibration of metal and silicate, and we evolve oxygen fugacity ( f O 2 ) self-consistently. The model successfully reproduces published meteorite-based estimates of most elemental abundances in the bulk silicate Mars, which can be used to estimate core formation conditions and core composition. This composition implies that the primordial material that formed Mars was significantly more oxidized (0.9-1.4 log units below the iron-wüstite buffer) than that of the Earth, and that core-mantle equilibration in Mars occurred at 42-60% of the evolving CMB pressure. On average, at least 84% of accreted metal and at least 40% of the mantle were equilibrated in each impact, a significantly higher degree of metal equilibration than previously reported for the Earth. In agreement with previous studies, the modeled Martian core is rich in sulfur (18-19 wt%), with less than one weight percent O and negligible Si. We have used these core and mantle compositions to produce physical models of the present-day Martian interior and evaluate the sensitivity of core radius to crustal thickness, mantle temperature, core composition, core temperature, and density of the core alloy. Trade-offs in how these properties affect observable physical parameters like planetary mass, radius, moment of inertia, and tidal Love number k 2 define a range of likely core radii: 1620-1870 km. Seismic velocity profiles for several combinations of model parameters have been used to predict seismic body-wave travel times and planetary normal mode frequencies. These results may be compared to forthcoming Martian seismic data to further constrain core formation conditions and geophysical properties.

Topography, or depth variation, of certain interfaces in the solid Earth can provide important insights into the dynamics of our planet interior. Although the intermediate- and long-range topographic variation of the 660-kilometer boundary between Earth's upper and lower mantle is well studied, small-scale measurements are far more challenging. We found a surprising amount of topography at short length scale along the 660-kilometer boundary in certain regions using scattered P'P' seismic waves. Our observations required chemical layering in regions with high short-scale roughness. By contrast, we did not see such small-scale topography along the 410-kilometer boundary in the upper mantle. Our findings support the concept of partially blocked or imperfect circulation between the upper and lower mantle., (Copyright © 2019 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

Turbulent convection of the liquid iron alloy outer core generates Earth's magnetic field and supplies heat to the mantle. The exact composition of the iron alloy is fundamentally linked to the processes powering the convection and can be constrained by its seismic properties. Discrepancies between seismic models determined using body waves and normal modes show that these properties are not yet fully agreed upon. In addition, technical challenges in experimentally measuring the equation-of-state (EoS) parameters of liquid iron alloys at high pressures and temperatures further complicate compositional inferences. We directly infer EoS parameters describing Earth's outer core from normal mode center frequency observations and present the resulting Elastic Parameters of the Outer Core (EPOC) seismic model. Unlike alternative seismic models, ours requires only three parameters and guarantees physically realistic behavior with increasing pressure for a well-mixed homogeneous material along an isentrope, consistent with the outer core's condition. We show that EPOC predicts available normal mode frequencies better than the Preliminary Reference Earth Model (PREM) while also being more consistent with body wave-derived models, eliminating a long-standing discrepancy. The velocity at the top of the outer core is lower, and increases with depth more steeply, in EPOC than in PREM, while the density in EPOC is higher than that in PREM across the outer core. The steeper profiles and higher density imply that the outer core comprises a lighter but more compressible alloy than that inferred for PREM. Furthermore, EPOC's steeper velocity gradient explains differential SmKS body wave travel times better than previous one-dimensional global models, without requiring an anomalously slow ~90- to 450-km-thick layer at the top of the outer core.