Terrestrial Single-Station Analog for Constraining the Martian Core and Deep Interior: Implications for InSight

We used a terrestrial single-station seismometer to quantify the uncertainty of InSight (INterior explorations using Seismic Investigations, Geodesy and Heat Transport) data for determining Martian core size. To mimic Martian seismicity, we formed a catalog using 917 terrestrial earthquakes, from which we randomly selected events. We stacked ScS amplitudes on modeled arrival times and searched for where ScS produced coherent seismic amplitudes. A core detection was defined by a coherent peak with small offset between predicted and user-selected arrival times. Iterating the detection algorithm with varying signal-to-noise (SNR) ranges and quantity of events determined the selection frequency of each model and quantified core depth uncertainty. Increasing the quantity of events reduced core depth uncertainty while increasing the recovery rate, while increasing event SNR had little effect. Including ScS2 multiples increased the recovery rate and reduced core depth uncertainty when we used low quantities of events. The most-frequent core depths varied by back azimuth, suggesting our method is sensitive to the presence of mantle heterogeneities. When we added 1° in source distance errors, core depth uncertainty increased by up to 11 km and recovery rates decreased by <5%. Altering epicentral distances by 25% added ~35 km of uncertainty and reduced recovery rates to <50% in some cases. From these experiments, we estimate that if InSight can detect five events with high location precision (<10 % epicentral distance errors), that there is at least an 88% chance of core depth recovery using ScS alone with uncertainty in core depth approaching 18 km and decreasing as more events are located.

Deformation of (Mg,Fe)SiO 3 Post-Perovskite and D'' Anisotropy

Polycrystalline (Mg(0.9),Fe(0.1))SiO3 post-perovskite was plastically deformed in the diamond anvil cell between 145 and 157 gigapascals. The lattice-preferred orientations obtained in the sample suggest that slip on planes near (100) and (110) dominate plastic deformation under these conditions. Assuming similar behavior at lower mantle conditions, we simulated plastic strains and the contribution of post-perovskite to anisotropy in the D'' region at the Earth core-mantle boundary using numerical convection and viscoplastic polycrystal plasticity models. We find a significant depth dependence of the anisotropy that only develops near and beyond the turning point of a downwelling slab. Our calculated anisotropies are strongly dependent on the choice of elastic moduli and remain hard to reconcile with seismic observations.

A Post-Perovskite Lens and D'' Heat Flux Beneath the Central Pacific

Temperature gradients in a low-shear-velocity province in the lowermost mantle (D'' region) beneath the central Pacific Ocean were inferred from the observation of a rapid S-wave velocity increase overlying a rapid decrease. These paired seismic discontinuities are attributed to a phase change from perovskite to post-perovskite and then back to perovskite as the temperature increases with depth. Iron enrichment could explain the occurrence of post-perovskite several hundred kilometers above the core-mantle boundary in this warm, chemically distinct province. The double phase-boundary crossing directly constrains the lowermost mantle temperature gradients. Assuming a standard but unconstrained choice of thermal conductivity, the regional core-mantle boundary heat flux (approximately 85 +/- 25 milliwatts per square meter), comparable to the average at Earth's surface, was estimated, along with a lower bound on global core-mantle boundary heat flow in the range of 13 +/- 4 terawatts. Mapped velocity-contrast variations indicate that the lens of post-perovskite minerals thins and vanishes over 1000 kilometers laterally toward the margin of the chemical distinct region as a result of a approximately 500-kelvin temperature increase.

Inferences on Flow at the Base of Earth's Mantle Based on Seismic Anisotropy

We applied global waveform tomography to model radial anisotropy in the whole mantle. We found that in the last few hundred kilometers near the core-mantle boundary, horizontally polarized S-wave velocities (VSH) are, on average, faster (by approximately 1%) than vertically polarized S-wave velocities (VSV), suggesting a large-scale predominance of horizontal shear. This confirms that the D" region at the base of the mantle is also a mechanical boundary layer for mantle convection. A notable exception to this average signature can be found at the base of the two broad low-velocity regions under the Pacific Ocean and under Africa, often referred to as "superplumes," where the anisotropic pattern indicates the onset of vertical flow.